#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 */
int /*<<< orphan*/ abort();
double const fabs(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ stderr;
double dstrb_betacfd(double a, double b, double x) {
  int m, m2;
  double aa, c, d, del, h;
  const int max_iterations = 100;
  const double eps = 3.0e-14;
  const double fp_min = 1.0e-44;
  const double qab = a + b;
  const double qap = a + 1.0;
  const double qam = a - 1.0;
  c = 1.0;
  d = 1.0 - qab * x / qap;
  if (fabs(d) < fp_min)
    d = fp_min;
  d = 1.0 / d;
  h = d;
  for (m = 1; m <= max_iterations; m++) {
    m2 = 2 * m;
    aa = m * (b - m) * x / ((qam + m2) * (a + m2));
    d = 1.0 + aa * d;
    if (fabs(d) < fp_min)
      d = fp_min;
    c = 1.0 + aa / c;
    if (fabs(c) < fp_min)
      c = fp_min;
    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) < fp_min)
      d = fp_min;
    c = 1.0 + aa / c;
    if (fabs(c) < fp_min)
      c = fp_min;
    d = 1.0 / d;
    del = d * c;
    h *= del;
    if (fabs(del - 1.0) < eps)
      break;
  }
  if (m > max_iterations) {
    fprintf(stderr, "a or b too big, or max_iterations too small in betacf\n");
    abort();
  }
  return h;
}


x86 assembly (clang, -O0):
<dstrb_betacfd>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movl   $0x64,-0x4c(%rbp)
movsd  0xeb7(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xe9a(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x8(%rbp),%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  0xe86(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  0xe6f(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xe5e(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x68(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
divsd  -0x70(%rbp),%xmm1
movsd  0xe42(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xe1c(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    11fb <dstrb_betacfd+0xcb>
movsd  0xe0a(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe05(%rip),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x1c(%rbp)
cmpl   $0x64,-0x1c(%rbp)
jg     1498 <dstrb_betacfd+0x368>
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x20(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x1c(%rbp),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm1
cvtsi2sdl -0x20(%rbp),%xmm2
addsd  %xmm2,%xmm1
movsd  -0x8(%rbp),%xmm2
cvtsi2sdl -0x20(%rbp),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd7c(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xd56(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    12c1 <dstrb_betacfd+0x191>
movsd  0xd44(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm1
divsd  -0x30(%rbp),%xmm1
movsd  0xd35(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xd0f(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1308 <dstrb_betacfd+0x1d8>
movsd  0xcfd(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xcf8(%rip),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x8(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x68(%rbp),%xmm1
cvtsi2sdl -0x1c(%rbp),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
cvtsi2sdl -0x20(%rbp),%xmm2
addsd  %xmm2,%xmm1
movsd  -0x70(%rbp),%xmm2
cvtsi2sdl -0x20(%rbp),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xc5f(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xc39(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    13de <dstrb_betacfd+0x2ae>
movsd  0xc27(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm1
divsd  -0x30(%rbp),%xmm1
movsd  0xc18(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xbf2(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1425 <dstrb_betacfd+0x2f5>
movsd  0xbe0(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xbdb(%rip),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  0xba6(%rip),%xmm1
subsd  %xmm1,%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xb9a(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1485 <dstrb_betacfd+0x355>
jmp    1498 <dstrb_betacfd+0x368>
jmp    148a <dstrb_betacfd+0x35a>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    121e <dstrb_betacfd+0xee>
cmpl   $0x64,-0x1c(%rbp)
jle    14be <dstrb_betacfd+0x38e>
mov    0x2b37(%rip),%rax
mov    (%rax),%edi
lea    0xb66(%rip),%rsi
call   1050 <fprintf@plt>
mov    $0x0,%al
call   1030 <abort@plt>
movsd  -0x48(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dstrb_betacfd>:
sub    $0x18,%rsp
movapd %xmm0,%xmm3
movsd  %xmm1,0x10(%rsp)
addsd  %xmm1,%xmm3
movsd  0xec5(%rip),%xmm11
movapd %xmm0,%xmm10
addsd  %xmm11,%xmm10
movsd  %xmm3,0x8(%rsp)
mulsd  %xmm2,%xmm3
divsd  %xmm10,%xmm3
movapd %xmm11,%xmm1
subsd  %xmm3,%xmm1
movapd 0xecb(%rip),%xmm3
andpd  %xmm1,%xmm3
movsd  0xe9f(%rip),%xmm5
movsd  0xe9f(%rip),%xmm7
ucomisd %xmm3,%xmm7
ja     1188 <dstrb_betacfd+0x68>
movapd %xmm11,%xmm5
divsd  %xmm1,%xmm5
movsd  0xe77(%rip),%xmm15
addsd  %xmm0,%xmm15
mov    $0x1,%eax
mov    $0x2,%ecx
movapd 0xe87(%rip),%xmm9
movapd %xmm11,%xmm4
movapd %xmm5,%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorps  %xmm12,%xmm12
cvtsi2sd %eax,%xmm12
movsd  0x10(%rsp),%xmm1
subsd  %xmm12,%xmm1
mulsd  %xmm12,%xmm1
mulsd  %xmm2,%xmm1
xorps  %xmm14,%xmm14
cvtsi2sd %ecx,%xmm14
movapd %xmm15,%xmm6
addsd  %xmm14,%xmm6
movapd %xmm14,%xmm13
addsd  %xmm0,%xmm13
mulsd  %xmm13,%xmm6
divsd  %xmm6,%xmm1
mulsd  %xmm1,%xmm5
addsd  %xmm11,%xmm5
movapd %xmm5,%xmm6
andpd  %xmm9,%xmm6
ucomisd %xmm6,%xmm7
movsd  0xdee(%rip),%xmm8
ja     122e <dstrb_betacfd+0x10e>
movapd %xmm11,%xmm8
divsd  %xmm5,%xmm8
movsd  0x8(%rsp),%xmm5
addsd  %xmm12,%xmm5
addsd  %xmm0,%xmm12
xorpd  0xdf9(%rip),%xmm12
mulsd  %xmm5,%xmm12
mulsd  %xmm2,%xmm12
addsd  %xmm10,%xmm14
mulsd  %xmm14,%xmm13
divsd  %xmm13,%xmm12
movapd %xmm12,%xmm6
mulsd  %xmm8,%xmm6
addsd  %xmm11,%xmm6
movapd %xmm6,%xmm5
andpd  %xmm9,%xmm5
ucomisd %xmm5,%xmm7
movsd  0xd8c(%rip),%xmm5
ja     128f <dstrb_betacfd+0x16f>
movapd %xmm11,%xmm5
divsd  %xmm6,%xmm5
divsd  %xmm4,%xmm1
addsd  %xmm11,%xmm1
movapd %xmm1,%xmm4
andpd  %xmm9,%xmm4
cmpltsd %xmm7,%xmm4
movapd %xmm4,%xmm6
andnpd %xmm1,%xmm6
andpd  %xmm7,%xmm4
orpd   %xmm6,%xmm4
mulsd  %xmm4,%xmm8
mulsd  %xmm8,%xmm3
divsd  %xmm4,%xmm12
addsd  %xmm11,%xmm12
movapd %xmm12,%xmm4
andpd  %xmm9,%xmm4
cmpltsd %xmm7,%xmm4
movapd %xmm4,%xmm1
andnpd %xmm12,%xmm1
andpd  %xmm7,%xmm4
orpd   %xmm1,%xmm4
movapd %xmm5,%xmm1
mulsd  %xmm4,%xmm1
mulsd  %xmm1,%xmm3
addsd  0xd0a(%rip),%xmm1
andpd  0xd2a(%rip),%xmm1
movsd  0xd12(%rip),%xmm6
ucomisd %xmm1,%xmm6
ja     1342 <dstrb_betacfd+0x222>
add    $0x1,%eax
add    $0x2,%ecx
cmp    $0xca,%ecx
jne    11c0 <dstrb_betacfd+0xa0>
mov    0x2cb3(%rip),%rax
mov    (%rax),%edi
lea    0xd1a(%rip),%rsi
call   1040 <fprintf@plt>
xor    %eax,%eax
call   1030 <abort@plt>
movapd %xmm3,%xmm0
add    $0x18,%rsp
ret

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

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