#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
unsigned long OS_AsmMult(unsigned long, unsigned long, unsigned long *);
float __mulsf3(float a_fp, float b_fp) {
  unsigned long a, b, c;
  unsigned long as, af, bs, bf, cs, cf;
  long ae, be, ce;
  unsigned long hi, lo;
  a = (*(unsigned long *)&(a_fp));
  b = (*(unsigned long *)&(b_fp));
  as = a >> 31;
  ae = (a >> 23) & 0xff;
  af = 0x00800000 | (a & 0x007fffff);
  bs = b >> 31;
  be = (b >> 23) & 0xff;
  bf = 0x00800000 | (b & 0x007fffff);
  cs = as ^ bs;
  lo = OS_AsmMult(af, bf, &hi);
  cf = (hi << 9) | (lo >> 23);
  ce = ae + be - 0x80 + 1;
  if (cf == 0)
    return (*(float *)&(cf));
  while (cf & 0xff000000) {
    ++ce;
    cf >>= 1;
  }
  c = (cs << 31) | (ce << 23) | (cf & 0x007fffff);
  if (ce < 1)
    c = 0;
  return (*(float *)&(c));
}


x86 assembly (clang, -O0):
<__mulsf3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
movss  %xmm0,-0x8(%rbp)
movss  %xmm1,-0xc(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x18(%rbp)
mov    -0xc(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x18(%rbp),%rax
shr    $0x1f,%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
shr    $0x17,%rax
and    $0xff,%rax
mov    %rax,-0x60(%rbp)
mov    -0x18(%rbp),%rax
and    $0x7fffff,%rax
or     $0x800000,%rax
mov    %rax,-0x38(%rbp)
mov    -0x20(%rbp),%rax
shr    $0x1f,%rax
mov    %rax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
shr    $0x17,%rax
and    $0xff,%rax
mov    %rax,-0x68(%rbp)
mov    -0x20(%rbp),%rax
and    $0x7fffff,%rax
or     $0x800000,%rax
mov    %rax,-0x48(%rbp)
mov    -0x30(%rbp),%rax
xor    -0x40(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x38(%rbp),%rdi
mov    -0x48(%rbp),%rsi
lea    -0x78(%rbp),%rdx
call   1030 <OS_AsmMult@plt>
mov    %rax,-0x80(%rbp)
mov    -0x78(%rbp),%rax
shl    $0x9,%rax
mov    -0x80(%rbp),%rcx
shr    $0x17,%rcx
or     %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x60(%rbp),%rax
add    -0x68(%rbp),%rax
sub    $0x80,%rax
add    $0x1,%rax
mov    %rax,-0x70(%rbp)
cmpq   $0x0,-0x58(%rbp)
jne    1201 <__mulsf3+0xf1>
movss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    127d <__mulsf3+0x16d>
jmp    1206 <__mulsf3+0xf6>
movabs $0xff000000,%rax
and    -0x58(%rbp),%rax
cmp    $0x0,%rax
je     123b <__mulsf3+0x12b>
mov    -0x70(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x70(%rbp)
mov    -0x58(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    1206 <__mulsf3+0xf6>
mov    -0x50(%rbp),%rax
shl    $0x1f,%rax
mov    -0x70(%rbp),%rcx
shl    $0x17,%rcx
or     %rcx,%rax
mov    -0x58(%rbp),%rcx
and    $0x7fffff,%rcx
or     %rcx,%rax
mov    %rax,-0x28(%rbp)
cmpq   $0x1,-0x70(%rbp)
jge    1273 <__mulsf3+0x163>
movq   $0x0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<__mulsf3>:
push   %r14
push   %rbx
sub    $0x18,%rsp
movd   %xmm0,0x8(%rsp)
movss  %xmm1,(%rsp)
mov    0x8(%rsp),%r14
mov    (%rsp),%rbx
mov    %r14d,%edi
and    $0x7fffff,%edi
or     $0x800000,%rdi
mov    %ebx,%esi
and    $0x7fffff,%esi
or     $0x800000,%rsi
lea    0x10(%rsp),%rdx
call   1030 <OS_AsmMult@plt>
mov    0x10(%rsp),%rcx
shl    $0x9,%rcx
shr    $0x17,%rax
or     %rcx,%rax
je     11b4 <__mulsf3+0xa4>
mov    %r14d,%ecx
shr    $0x17,%ecx
movzbl %cl,%ecx
mov    %ebx,%edx
shr    $0x17,%edx
movzbl %dl,%edx
add    %rdx,%rcx
add    $0xffffffffffffff81,%rcx
test   $0xff000000,%eax
je     11ba <__mulsf3+0xaa>
movabs $0x1fe000000,%rsi
mov    %rax,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rcx
shr    %rdx
test   %rsi,%rax
mov    %rdx,%rax
jne    11a0 <__mulsf3+0x90>
test   %rcx,%rcx
jg     11c2 <__mulsf3+0xb2>
pxor   %xmm0,%xmm0
jmp    11dc <__mulsf3+0xcc>
mov    %rax,%rdx
test   %rcx,%rcx
jle    11b4 <__mulsf3+0xa4>
xor    %r14d,%ebx
and    $0x80000000,%ebx
shl    $0x17,%ecx
and    $0x7fffff,%edx
or     %edx,%ebx
or     %ecx,%ebx
movd   %ebx,%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define float @__mulsf3(float %0, float %1) {
%3 = alloca float
%4 = alloca float
%5 = alloca float
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
%15 = alloca i64
%16 = alloca i64
%17 = alloca i64
%18 = alloca i64
%19 = alloca i64
store float %0, float* %4
store float %1, float* %5
%20 = bitcast float* %4 to i64*
%21 = load i64, i64* %20
store i64 %21, i64* %6
%22 = bitcast float* %5 to i64*
%23 = load i64, i64* %22
store i64 %23, i64* %7
%24 = load i64, i64* %6
%25 = lshr i64 %24, 31
store i64 %25, i64* %9
%26 = load i64, i64* %6
%27 = lshr i64 %26, 23
%28 = and i64 %27, 255
store i64 %28, i64* %15
%29 = load i64, i64* %6
%30 = and i64 %29, 8388607
%31 = or i64 8388608, %30
store i64 %31, i64* %10
%32 = load i64, i64* %7
%33 = lshr i64 %32, 31
store i64 %33, i64* %11
%34 = load i64, i64* %7
%35 = lshr i64 %34, 23
%36 = and i64 %35, 255
store i64 %36, i64* %16
%37 = load i64, i64* %7
%38 = and i64 %37, 8388607
%39 = or i64 8388608, %38
store i64 %39, i64* %12
%40 = load i64, i64* %9
%41 = load i64, i64* %11
%42 = xor i64 %40, %41
store i64 %42, i64* %13
%43 = load i64, i64* %10
%44 = load i64, i64* %12
%45 = call i64 @OS_AsmMult(i64 %43, i64 %44, i64* %18)
store i64 %45, i64* %19
%46 = load i64, i64* %18
%47 = shl i64 %46, 9
%48 = load i64, i64* %19
%49 = lshr i64 %48, 23
%50 = or i64 %47, %49
store i64 %50, i64* %14
%51 = load i64, i64* %15
%52 = load i64, i64* %16
%53 = add nsw i64 %51, %52
%54 = sub nsw i64 %53, 128
%55 = add nsw i64 %54, 1
store i64 %55, i64* %17
%56 = load i64, i64* %14
%57 = icmp eq i64 %56, 0
br i1 %57, label %58, label %61
58: ; preds = %2
%59 = bitcast i64* %14 to float*
%60 = load float, float* %59
store float %60, float* %3
br label %86
61: ; preds = %2
br label %62
62: ; preds = %66, %61
%63 = load i64, i64* %14
%64 = and i64 %63, 4278190080
%65 = icmp ne i64 %64, 0
br i1 %65, label %66, label %71
66: ; preds = %62
%67 = load i64, i64* %17
%68 = add nsw i64 %67, 1
store i64 %68, i64* %17
%69 = load i64, i64* %14
%70 = lshr i64 %69, 1
store i64 %70, i64* %14
br label %62
71: ; preds = %62
%72 = load i64, i64* %13
%73 = shl i64 %72, 31
%74 = load i64, i64* %17
%75 = shl i64 %74, 23
%76 = or i64 %73, %75
%77 = load i64, i64* %14
%78 = and i64 %77, 8388607
%79 = or i64 %76, %78
store i64 %79, i64* %8
%80 = load i64, i64* %17
%81 = icmp slt i64 %80, 1
br i1 %81, label %82, label %83
82: ; preds = %71
store i64 0, i64* %8
br label %83
83: ; preds = %82, %71
%84 = bitcast i64* %8 to float*
%85 = load float, float* %84
store float %85, float* %3
br label %86
86: ; preds = %83, %58
%87 = load float, float* %3
ret float %87
}

LLVM IR (-O2):
define float @__mulsf3(float %0, float %1) {
%3 = alloca float
%4 = alloca float
%5 = alloca i64
store float %0, float* %3
store float %1, float* %4
%6 = bitcast i64* %5 to i8*
%7 = bitcast float* %3 to i64*
%8 = load i64, i64* %7
%9 = bitcast float* %4 to i64*
%10 = load i64, i64* %9
%11 = lshr i64 %8, 23
%12 = and i64 %11, 255
%13 = and i64 %8, 8388607
%14 = or i64 %13, 8388608
%15 = lshr i64 %10, 23
%16 = and i64 %15, 255
%17 = and i64 %10, 8388607
%18 = or i64 %17, 8388608
%19 = xor i64 %10, %8
%20 = and i64 %19, 2147483648
%21 = call i64 @OS_AsmMult(i64 %14, i64 %18, i64* %5)
%22 = load i64, i64* %5
%23 = shl i64 %22, 9
%24 = lshr i64 %21, 23
%25 = or i64 %23, %24
%26 = add nsw i64 %12, -127
%27 = add nsw i64 %26, %16
%28 = icmp eq i64 %25, 0
br i1 %28, label %50, label %29
29: ; preds = %2
%30 = and i64 %25, 4278190080
%31 = icmp eq i64 %30, 0
br i1 %31, label %39, label %32
32: ; preds = %29, %32
%33 = phi i64 [ %35, %32 ], [ %27, %29 ]
%34 = phi i64 [ %36, %32 ], [ %25, %29 ]
%35 = add nsw i64 %33, 1
%36 = lshr i64 %34, 1
%37 = and i64 %34, 8556380160
%38 = icmp eq i64 %37, 0
br i1 %38, label %39, label %32
39: ; preds = %32, %29
%40 = phi i64 [ %25, %29 ], [ %36, %32 ]
%41 = phi i64 [ %27, %29 ], [ %35, %32 ]
%42 = shl i64 %41, 23
%43 = and i64 %40, 8388607
%44 = or i64 %43, %20
%45 = or i64 %44, %42
%46 = icmp sgt i64 %41, 0
%47 = trunc i64 %45 to i32
%48 = bitcast i32 %47 to float
%49 = select i1 %46, float %48, float 0.000000e+00
br label %50
50: ; preds = %2, %39
%51 = phi float [ %49, %39 ], [ 0.000000e+00, %2 ]
ret float %51
}