#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 */
double fabs(float);
int /*<<< orphan*/ printf(char *);
float betacf(float a, float b, float x) {
  int m, m2;
  float aa, c, d, del, h, qab, qam, qap;
  qab = a + b;
  qap = a + 1.0;
  qam = a - 1.0;
  c = 1.0;
  d = 1.0 - qab * x / qap;
  if (fabs(d) < 1.0e-30)
    d = 1.0e-30;
  d = 1.0 / d;
  h = d;
  for (m = 1; m <= 100; m++) {
    m2 = 2 * m;
    aa = m * (b - m) * x / ((qam + m2) * (a + m2));
    d = 1.0 + aa * d;
    if (fabs(d) < 1.0e-30)
      d = 1.0e-30;
    c = 1.0 + aa / c;
    if (fabs(c) < 1.0e-30)
      c = 1.0e-30;
    d = 1.0 / d;
    h *= d * c;
    aa = -(a + m) * (qab + m) * x / ((a + m2) * (qap + m2));
    d = 1.0 + aa * d;
    if (fabs(d) < 1.0e-30)
      d = 1.0e-30;
    c = 1.0 + aa / c;
    if (fabs(c) < 1.0e-30)
      c = 1.0e-30;
    d = 1.0 / d;
    del = d * c;
    h *= del;
    if (fabs(del - 1.0) < 3.0e-7)
      break;
  }
  if (m > 100) {
    printf("Error(betacf):a or b too big, or MAXIT too small\n");
    return (0.0);
  }
  return (h);
}


x86 assembly (clang, -O0):
<betacf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
movss  %xmm0,-0x8(%rbp)
movss  %xmm1,-0xc(%rbp)
movss  %xmm2,-0x10(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0xc(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xeb1(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe93(%rip),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  0xe8e(%rip),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x30(%rbp),%xmm0
mulss  -0x10(%rbp),%xmm0
divss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe5e(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xe34(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    11e3 <betacf+0xc3>
movss  0xe3e(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe14(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x1,-0x14(%rbp)
cmpl   $0x64,-0x14(%rbp)
jg     14ba <betacf+0x39a>
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x18(%rbp)
cvtsi2ssl -0x14(%rbp),%xmm0
movss  -0xc(%rbp),%xmm1
cvtsi2ssl -0x14(%rbp),%xmm2
subss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mulss  -0x10(%rbp),%xmm0
movss  -0x34(%rbp),%xmm1
cvtsi2ssl -0x18(%rbp),%xmm2
addss  %xmm2,%xmm1
movss  -0x8(%rbp),%xmm2
cvtsi2ssl -0x18(%rbp),%xmm3
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%rbp),%xmm0
mulss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd88(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xd5e(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    12b9 <betacf+0x199>
movss  0xd68(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x1c(%rbp),%xmm0
divss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd39(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xd0f(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1308 <betacf+0x1e8>
movss  0xd19(%rip),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xcef(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
mulss  -0x20(%rbp),%xmm0
mulss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x8(%rbp),%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
addss  %xmm1,%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  -0x30(%rbp),%xmm1
cvtsi2ssl -0x14(%rbp),%xmm2
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mulss  -0x10(%rbp),%xmm0
movss  -0x8(%rbp),%xmm1
cvtsi2ssl -0x18(%rbp),%xmm2
addss  %xmm2,%xmm1
movss  -0x38(%rbp),%xmm2
cvtsi2ssl -0x18(%rbp),%xmm3
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%rbp),%xmm0
mulss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xc5d(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xc33(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    13e4 <betacf+0x2c4>
movss  0xc3d(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x1c(%rbp),%xmm0
divss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xc0e(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xbe4(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1433 <betacf+0x313>
movss  0xbee(%rip),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xbc4(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
mulss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
mulss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xb88(%rip),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xb78(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    14a7 <betacf+0x387>
jmp    14ba <betacf+0x39a>
jmp    14ac <betacf+0x38c>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1212 <betacf+0xf2>
cmpl   $0x64,-0x14(%rbp)
jle    14dd <betacf+0x3bd>
lea    0xb55(%rip),%rdi
call   1030 <printf@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    14e7 <betacf+0x3c7>
movss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<betacf>:
push   %rbp
push   %rbx
sub    $0x68,%rsp
movaps %xmm0,%xmm3
movss  %xmm1,0x48(%rsp)
addss  %xmm1,%xmm3
movss  0xec5(%rip),%xmm4
movss  %xmm0,0xc(%rsp)
movaps %xmm0,%xmm1
addss  %xmm4,%xmm1
movss  %xmm3,0x44(%rsp)
movss  %xmm2,0x1c(%rsp)
mulss  %xmm2,%xmm3
movss  %xmm1,0x40(%rsp)
divss  %xmm1,%xmm3
movaps %xmm4,%xmm0
subss  %xmm3,%xmm0
movss  %xmm0,0x20(%rsp)
call   1040 <fabs@plt>
movss  0xe8c(%rip),%xmm1
movsd  0xe8c(%rip),%xmm2
ucomisd %xmm0,%xmm2
movss  0xe70(%rip),%xmm0
ja     119b <betacf+0x7b>
movaps %xmm0,%xmm1
divss  0x20(%rsp),%xmm1
movss  0xe61(%rip),%xmm2
addss  0xc(%rsp),%xmm2
movss  %xmm2,0x3c(%rsp)
mov    $0x1,%ebx
mov    $0x2,%ebp
movss  %xmm0,(%rsp)
movss  %xmm1,0x8(%rsp)
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xorps  %xmm0,%xmm0
cvtsi2ss %ebx,%xmm0
movss  0x48(%rsp),%xmm3
subss  %xmm0,%xmm3
xorps  %xmm2,%xmm2
cvtsi2ss %ebp,%xmm2
movaps %xmm0,0x20(%rsp)
mulss  %xmm0,%xmm3
mulss  0x1c(%rsp),%xmm3
movss  0x3c(%rsp),%xmm0
addss  %xmm2,%xmm0
movss  %xmm2,0x50(%rsp)
addss  0xc(%rsp),%xmm2
movss  %xmm2,0x4c(%rsp)
mulss  %xmm2,%xmm0
divss  %xmm0,%xmm3
movss  %xmm3,0x4(%rsp)
mulss  %xmm3,%xmm1
movss  0xdd3(%rip),%xmm0
addss  %xmm0,%xmm1
movss  %xmm1,0x10(%rsp)
movaps %xmm1,%xmm0
call   1040 <fabs@plt>
movsd  %xmm0,0x60(%rsp)
movss  0x4(%rsp),%xmm0
divss  (%rsp),%xmm0
addss  0xda8(%rip),%xmm0
movss  %xmm0,0x4(%rsp)
call   1040 <fabs@plt>
movsd  0xda5(%rip),%xmm1
ucomisd %xmm0,%xmm1
movss  0xd95(%rip),%xmm0
movss  %xmm0,(%rsp)
ja     1289 <betacf+0x169>
movss  0x4(%rsp),%xmm0
movss  %xmm0,(%rsp)
movsd  0xd7f(%rip),%xmm0
ucomisd 0x60(%rsp),%xmm0
movss  0xd69(%rip),%xmm1
ja     12af <betacf+0x18f>
movss  0xd57(%rip),%xmm1
divss  0x10(%rsp),%xmm1
movss  0x44(%rsp),%xmm0
movaps 0x20(%rsp),%xmm3
addss  %xmm3,%xmm0
addss  0xc(%rsp),%xmm3
xorps  0xd55(%rip),%xmm3
mulss  %xmm0,%xmm3
mulss  0x1c(%rsp),%xmm3
movss  0x50(%rsp),%xmm0
addss  0x40(%rsp),%xmm0
movss  0x4c(%rsp),%xmm2
mulss  %xmm0,%xmm2
divss  %xmm2,%xmm3
movaps %xmm3,0x20(%rsp)
movaps %xmm3,%xmm2
movss  %xmm1,0x4(%rsp)
mulss  %xmm1,%xmm2
movss  0xcf7(%rip),%xmm0
addss  %xmm0,%xmm2
movaps %xmm2,0x50(%rsp)
movaps %xmm2,%xmm0
call   1040 <fabs@plt>
movsd  %xmm0,0x10(%rsp)
movaps 0x20(%rsp),%xmm0
divss  (%rsp),%xmm0
addss  0xcce(%rip),%xmm0
movaps %xmm0,0x20(%rsp)
call   1040 <fabs@plt>
movsd  0xccc(%rip),%xmm1
ucomisd %xmm0,%xmm1
movss  0xcbc(%rip),%xmm2
ja     1357 <betacf+0x237>
movaps 0x20(%rsp),%xmm2
movsd  0xcb1(%rip),%xmm0
ucomisd 0x10(%rsp),%xmm0
movss  0xc9b(%rip),%xmm0
movss  0x8(%rsp),%xmm1
ja     1383 <betacf+0x263>
movss  0xc83(%rip),%xmm0
divss  0x50(%rsp),%xmm0
movss  0x4(%rsp),%xmm3
mulss  (%rsp),%xmm3
mulss  %xmm3,%xmm1
movss  %xmm0,0x10(%rsp)
movss  %xmm2,(%rsp)
mulss  %xmm2,%xmm0
mulss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
addss  0xc51(%rip),%xmm0
call   1040 <fabs@plt>
movsd  0xc58(%rip),%xmm1
ucomisd %xmm0,%xmm1
ja     13ef <betacf+0x2cf>
add    $0x1,%ebx
add    $0x2,%ebp
cmp    $0xca,%ebp
movss  0x10(%rsp),%xmm1
jne    11d0 <betacf+0xb0>
lea    0xc4b(%rip),%rdi
call   1030 <puts@plt>
xorps  %xmm0,%xmm0
jmp    13f5 <betacf+0x2d5>
movss  0x8(%rsp),%xmm0
add    $0x68,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
@.str = [50 x i8] c"Error(betacf):a or b too big, or MAXIT too small\0A\00"
define float @betacf(float %0, float %1, float %2) {
%4 = alloca float
%5 = alloca float
%6 = alloca float
%7 = alloca float
%8 = alloca i32
%9 = alloca i32
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
store float %0, float* %5
store float %1, float* %6
store float %2, float* %7
%18 = load float, float* %5
%19 = load float, float* %6
%20 = fadd float %18, %19
store float %20, float* %15
%21 = load float, float* %5
%22 = fpext float %21 to double
%23 = fadd double %22, 1.000000e+00
%24 = fptrunc double %23 to float
store float %24, float* %17
%25 = load float, float* %5
%26 = fpext float %25 to double
%27 = fsub double %26, 1.000000e+00
%28 = fptrunc double %27 to float
store float %28, float* %16
store float 1.000000e+00, float* %11
%29 = load float, float* %15
%30 = load float, float* %7
%31 = fmul float %29, %30
%32 = load float, float* %17
%33 = fdiv float %31, %32
%34 = fpext float %33 to double
%35 = fsub double 1.000000e+00, %34
%36 = fptrunc double %35 to float
store float %36, float* %12
%37 = load float, float* %12
%38 = call double @fabs(float %37)
%39 = fcmp olt double %38, 1.000000e-30
br i1 %39, label %40, label %41
40: ; preds = %3
store float 0x39B4484C00000000, float* %12
br label %41
41: ; preds = %40, %3
%42 = load float, float* %12
%43 = fpext float %42 to double
%44 = fdiv double 1.000000e+00, %43
%45 = fptrunc double %44 to float
store float %45, float* %12
%46 = load float, float* %12
store float %46, float* %14
store i32 1, i32* %8
br label %47
47: ; preds = %165, %41
%48 = load i32, i32* %8
%49 = icmp sle i32 %48, 100
br i1 %49, label %50, label %168
50: ; preds = %47
%51 = load i32, i32* %8
%52 = mul nsw i32 2, %51
store i32 %52, i32* %9
%53 = load i32, i32* %8
%54 = sitofp i32 %53 to float
%55 = load float, float* %6
%56 = load i32, i32* %8
%57 = sitofp i32 %56 to float
%58 = fsub float %55, %57
%59 = fmul float %54, %58
%60 = load float, float* %7
%61 = fmul float %59, %60
%62 = load float, float* %16
%63 = load i32, i32* %9
%64 = sitofp i32 %63 to float
%65 = fadd float %62, %64
%66 = load float, float* %5
%67 = load i32, i32* %9
%68 = sitofp i32 %67 to float
%69 = fadd float %66, %68
%70 = fmul float %65, %69
%71 = fdiv float %61, %70
store float %71, float* %10
%72 = load float, float* %10
%73 = load float, float* %12
%74 = fmul float %72, %73
%75 = fpext float %74 to double
%76 = fadd double 1.000000e+00, %75
%77 = fptrunc double %76 to float
store float %77, float* %12
%78 = load float, float* %12
%79 = call double @fabs(float %78)
%80 = fcmp olt double %79, 1.000000e-30
br i1 %80, label %81, label %82
81: ; preds = %50
store float 0x39B4484C00000000, float* %12
br label %82
82: ; preds = %81, %50
%83 = load float, float* %10
%84 = load float, float* %11
%85 = fdiv float %83, %84
%86 = fpext float %85 to double
%87 = fadd double 1.000000e+00, %86
%88 = fptrunc double %87 to float
store float %88, float* %11
%89 = load float, float* %11
%90 = call double @fabs(float %89)
%91 = fcmp olt double %90, 1.000000e-30
br i1 %91, label %92, label %93
92: ; preds = %82
store float 0x39B4484C00000000, float* %11
br label %93
93: ; preds = %92, %82
%94 = load float, float* %12
%95 = fpext float %94 to double
%96 = fdiv double 1.000000e+00, %95
%97 = fptrunc double %96 to float
store float %97, float* %12
%98 = load float, float* %12
%99 = load float, float* %11
%100 = fmul float %98, %99
%101 = load float, float* %14
%102 = fmul float %101, %100
store float %102, float* %14
%103 = load float, float* %5
%104 = load i32, i32* %8
%105 = sitofp i32 %104 to float
%106 = fadd float %103, %105
%107 = fneg float %106
%108 = load float, float* %15
%109 = load i32, i32* %8
%110 = sitofp i32 %109 to float
%111 = fadd float %108, %110
%112 = fmul float %107, %111
%113 = load float, float* %7
%114 = fmul float %112, %113
%115 = load float, float* %5
%116 = load i32, i32* %9
%117 = sitofp i32 %116 to float
%118 = fadd float %115, %117
%119 = load float, float* %17
%120 = load i32, i32* %9
%121 = sitofp i32 %120 to float
%122 = fadd float %119, %121
%123 = fmul float %118, %122
%124 = fdiv float %114, %123
store float %124, float* %10
%125 = load float, float* %10
%126 = load float, float* %12
%127 = fmul float %125, %126
%128 = fpext float %127 to double
%129 = fadd double 1.000000e+00, %128
%130 = fptrunc double %129 to float
store float %130, float* %12
%131 = load float, float* %12
%132 = call double @fabs(float %131)
%133 = fcmp olt double %132, 1.000000e-30
br i1 %133, label %134, label %135
134: ; preds = %93
store float 0x39B4484C00000000, float* %12
br label %135
135: ; preds = %134, %93
%136 = load float, float* %10
%137 = load float, float* %11
%138 = fdiv float %136, %137
%139 = fpext float %138 to double
%140 = fadd double 1.000000e+00, %139
%141 = fptrunc double %140 to float
store float %141, float* %11
%142 = load float, float* %11
%143 = call double @fabs(float %142)
%144 = fcmp olt double %143, 1.000000e-30
br i1 %144, label %145, label %146
145: ; preds = %135
store float 0x39B4484C00000000, float* %11
br label %146
146: ; preds = %145, %135
%147 = load float, float* %12
%148 = fpext float %147 to double
%149 = fdiv double 1.000000e+00, %148
%150 = fptrunc double %149 to float
store float %150, float* %12
%151 = load float, float* %12
%152 = load float, float* %11
%153 = fmul float %151, %152
store float %153, float* %13
%154 = load float, float* %13
%155 = load float, float* %14
%156 = fmul float %155, %154
store float %156, float* %14
%157 = load float, float* %13
%158 = fpext float %157 to double
%159 = fsub double %158, 1.000000e+00
%160 = fptrunc double %159 to float
%161 = call double @fabs(float %160)
%162 = fcmp olt double %161, 3.000000e-07
br i1 %162, label %163, label %164
163: ; preds = %146
br label %168
164: ; preds = %146
br label %165
165: ; preds = %164
%166 = load i32, i32* %8
%167 = add nsw i32 %166, 1
store i32 %167, i32* %8
br label %47
168: ; preds = %163, %47
%169 = load i32, i32* %8
%170 = icmp sgt i32 %169, 100
br i1 %170, label %171, label %173
171: ; preds = %168
%172 = call i32 @printf(i8* getelementptr inbounds ([50 x i8], [50 x i8]* @.str, i64 0, i64 0))
store float 0.000000e+00, float* %4
br label %175
173: ; preds = %168
%174 = load float, float* %14
store float %174, float* %4
br label %175
175: ; preds = %173, %171
%176 = load float, float* %4
ret float %176
}

LLVM IR (-O2):
@str = [49 x i8] c"Error(betacf):a or b too big, or MAXIT too small\00"
define float @betacf(float %0, float %1, float %2) {
%4 = fadd float %0, %1
%5 = fadd float %0, 1.000000e+00
%6 = fadd float %0, -1.000000e+00
%7 = fmul float %4, %2
%8 = fdiv float %7, %5
%9 = fsub float 1.000000e+00, %8
%10 = tail call double @fabs(float %9)
%11 = fcmp olt double %10, 1.000000e-30
%12 = fdiv float 1.000000e+00, %9
%13 = select i1 %11, float 0x46293E5940000000, float %12
br label %17
14: ; preds = %17
%15 = add nuw nsw i32 %18, 1
%16 = icmp eq i32 %15, 101
br i1 %16, label %69, label %17
17: ; preds = %3, %14
%18 = phi i32 [ 1, %3 ], [ %15, %14 ]
%19 = phi float [ %13, %3 ], [ %65, %14 ]
%20 = phi float [ 1.000000e+00, %3 ], [ %61, %14 ]
%21 = phi float [ %13, %3 ], [ %63, %14 ]
%22 = shl nuw nsw i32 %18, 1
%23 = sitofp i32 %18 to float
%24 = fsub float %1, %23
%25 = fmul float %24, %23
%26 = fmul float %25, %2
%27 = sitofp i32 %22 to float
%28 = fadd float %6, %27
%29 = fadd float %27, %0
%30 = fmul float %28, %29
%31 = fdiv float %26, %30
%32 = fmul float %21, %31
%33 = fadd float %32, 1.000000e+00
%34 = tail call double @fabs(float %33)
%35 = fcmp olt double %34, 1.000000e-30
%36 = fdiv float %31, %20
%37 = fadd float %36, 1.000000e+00
%38 = tail call double @fabs(float %37)
%39 = fcmp olt double %38, 1.000000e-30
%40 = select i1 %39, float 0x39B4484C00000000, float %37
%41 = fdiv float 1.000000e+00, %33
%42 = select i1 %35, float 0x46293E5940000000, float %41
%43 = fmul float %42, %40
%44 = fmul float %19, %43
%45 = fadd float %23, %0
%46 = fneg float %45
%47 = fadd float %4, %23
%48 = fmul float %47, %46
%49 = fmul float %48, %2
%50 = fadd float %5, %27
%51 = fmul float %29, %50
%52 = fdiv float %49, %51
%53 = fmul float %52, %42
%54 = fadd float %53, 1.000000e+00
%55 = tail call double @fabs(float %54)
%56 = fcmp olt double %55, 1.000000e-30
%57 = fdiv float %52, %40
%58 = fadd float %57, 1.000000e+00
%59 = tail call double @fabs(float %58)
%60 = fcmp olt double %59, 1.000000e-30
%61 = select i1 %60, float 0x39B4484C00000000, float %58
%62 = fdiv float 1.000000e+00, %54
%63 = select i1 %56, float 0x46293E5940000000, float %62
%64 = fmul float %63, %61
%65 = fmul float %44, %64
%66 = fadd float %64, -1.000000e+00
%67 = tail call double @fabs(float %66)
%68 = fcmp olt double %67, 3.000000e-07
br i1 %68, label %71, label %14
69: ; preds = %14
%70 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([49 x i8], [49 x i8]* @str, i64 0, i64 0))
br label %71
71: ; preds = %17, %69
%72 = phi float [ 0.000000e+00, %69 ], [ %65, %17 ]
ret float %72
}