#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);
float betacf(float a, float b, float x) {
  void nrerror(char error_text[]);
  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)
    nrerror("a or b too big, or MAXIT too small in betacf");
  return h;
}


x86 assembly (clang, -O0):
<betacf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
movss  %xmm0,-0x4(%rbp)
movss  %xmm1,-0x8(%rbp)
movss  %xmm2,-0xc(%rbp)
movss  -0x4(%rbp),%xmm0
addss  -0x8(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xeb1(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe93(%rip),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  0xe8e(%rip),%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x2c(%rbp),%xmm0
mulss  -0xc(%rbp),%xmm0
divss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe5e(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%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,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe14(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x1,-0x10(%rbp)
cmpl   $0x64,-0x10(%rbp)
jg     14ba <betacf+0x39a>
mov    -0x10(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x14(%rbp)
cvtsi2ssl -0x10(%rbp),%xmm0
movss  -0x8(%rbp),%xmm1
cvtsi2ssl -0x10(%rbp),%xmm2
subss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mulss  -0xc(%rbp),%xmm0
movss  -0x30(%rbp),%xmm1
cvtsi2ssl -0x14(%rbp),%xmm2
addss  %xmm2,%xmm1
movss  -0x4(%rbp),%xmm2
cvtsi2ssl -0x14(%rbp),%xmm3
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
mulss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd88(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%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,-0x20(%rbp)
movss  -0x18(%rbp),%xmm0
divss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd39(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%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,-0x1c(%rbp)
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xcef(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
mulss  -0x1c(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x4(%rbp),%xmm0
cvtsi2ssl -0x10(%rbp),%xmm1
addss  %xmm1,%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  -0x2c(%rbp),%xmm1
cvtsi2ssl -0x10(%rbp),%xmm2
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mulss  -0xc(%rbp),%xmm0
movss  -0x4(%rbp),%xmm1
cvtsi2ssl -0x14(%rbp),%xmm2
addss  %xmm2,%xmm1
movss  -0x34(%rbp),%xmm2
cvtsi2ssl -0x14(%rbp),%xmm3
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
mulss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xc5d(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%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,-0x20(%rbp)
movss  -0x18(%rbp),%xmm0
divss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xc0e(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%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,-0x1c(%rbp)
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xbc4(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
mulss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x24(%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    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1212 <betacf+0xf2>
cmpl   $0x64,-0x10(%rbp)
jle    14d0 <betacf+0x3b0>
lea    0xb55(%rip),%rdi
call   1030 <nrerror@plt>
movss  -0x28(%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     13ea <betacf+0x2ca>
add    $0x1,%ebx
add    $0x2,%ebp
cmp    $0xca,%ebp
movss  0x10(%rsp),%xmm1
jne    11d0 <betacf+0xb0>
lea    0xc4b(%rip),%rdi
call   1030 <nrerror@plt>
movss  0x8(%rsp),%xmm0
add    $0x68,%rsp
pop    %rbx
pop    %rbp
ret

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

LLVM IR (-O2):
@.str = [45 x i8] c"a or b too big, or MAXIT too small in betacf\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 [ 1.000000e+00, %3 ], [ %61, %14 ]
%20 = phi float [ %13, %3 ], [ %65, %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, %19
%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 %20, %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 %70, label %14
69: ; preds = %14
tail call void @nrerror(i8* getelementptr inbounds ([45 x i8], [45 x i8]* @.str, i64 0, i64 0))
br label %70
70: ; preds = %17, %69
ret float %65
}