#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(double);
double betacf(int a, int b, double x) {
  int m, m2;
  double 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;
  }
  return h;
}


x86 assembly (clang, -O0):
<betacf>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    -0x4(%rbp),%eax
add    -0x8(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm0
movsd  0xedd(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm0
movsd  0xec7(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xeb6(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x48(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
divsd  -0x58(%rbp),%xmm1
movsd  0xe9a(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm1
movaps 0xe9d(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xe71(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    11a6 <betacf+0xa6>
movsd  0xe5f(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xe5a(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x14(%rbp)
cmpl   $0x64,-0x14(%rbp)
jg     142c <betacf+0x32c>
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    -0x14(%rbp),%ecx
imul   %ecx,%eax
cvtsi2sd %eax,%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
cvtsi2sdl -0x18(%rbp),%xmm2
addsd  %xmm2,%xmm1
mov    -0x4(%rbp),%eax
add    -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xddc(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm1
movaps 0xddf(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xdb3(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1264 <betacf+0x164>
movsd  0xda1(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  0xd92(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
movaps 0xd95(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xd69(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    12ae <betacf+0x1ae>
movsd  0xd57(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xd52(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4(%rbp),%ecx
add    -0x14(%rbp),%ecx
xor    %eax,%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  -0x48(%rbp),%xmm1
cvtsi2sdl -0x14(%rbp),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x10(%rbp),%xmm0
mov    -0x4(%rbp),%eax
add    -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  -0x58(%rbp),%xmm2
cvtsi2sdl -0x18(%rbp),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xcd4(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm1
movaps 0xcd7(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xcab(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    136c <betacf+0x26c>
movsd  0xc99(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm1
divsd  -0x28(%rbp),%xmm1
movsd  0xc8a(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
movaps 0xc8d(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xc61(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    13b6 <betacf+0x2b6>
movsd  0xc4f(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xc4a(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  0xc15(%rip),%xmm0
subsd  %xmm0,%xmm1
movaps 0xc22(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xc06(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1419 <betacf+0x319>
jmp    142c <betacf+0x32c>
jmp    141e <betacf+0x31e>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11c9 <betacf+0xc9>
movsd  -0x40(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<betacf>:
lea    (%rsi,%rdi,1),%eax
cvtsi2sd %eax,%xmm1
cvtsi2sd %edi,%xmm13
movsd  0xeeb(%rip),%xmm8
movapd %xmm13,%xmm14
addsd  %xmm8,%xmm14
movsd  %xmm1,-0x8(%rsp)
mulsd  %xmm0,%xmm1
divsd  %xmm14,%xmm1
movapd %xmm8,%xmm2
subsd  %xmm1,%xmm2
movapd 0xef1(%rip),%xmm1
andpd  %xmm2,%xmm1
movsd  0xec4(%rip),%xmm9
movsd  0xec4(%rip),%xmm7
ucomisd %xmm1,%xmm7
movapd %xmm9,%xmm1
ja     1168 <betacf+0x68>
movapd %xmm8,%xmm1
divsd  %xmm2,%xmm1
addsd  0xe97(%rip),%xmm13
mov    %edi,%r8d
not    %r8d
add    $0xffffffff,%esi
mov    $0x1,%ecx
mov    $0x2,%edx
movapd 0xea3(%rip),%xmm15
movapd %xmm8,%xmm4
movapd %xmm1,%xmm6
cs nopw 0x0(%rax,%rax,1)
mov    %esi,%eax
imul   %ecx,%eax
xorps  %xmm5,%xmm5
cvtsi2sd %eax,%xmm5
xorps  %xmm10,%xmm10
cvtsi2sd %edx,%xmm10
mulsd  %xmm0,%xmm5
movapd %xmm13,%xmm2
lea    (%rdi,%rdx,1),%eax
cvtsi2sd %eax,%xmm11
addsd  %xmm10,%xmm2
mulsd  %xmm11,%xmm2
divsd  %xmm2,%xmm5
mulsd  %xmm5,%xmm6
addsd  %xmm8,%xmm6
movapd %xmm6,%xmm2
andpd  %xmm15,%xmm2
ucomisd %xmm2,%xmm7
movapd %xmm9,%xmm12
ja     11fb <betacf+0xfb>
movapd %xmm8,%xmm12
divsd  %xmm6,%xmm12
xorps  %xmm3,%xmm3
cvtsi2sd %r8d,%xmm3
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
addsd  -0x8(%rsp),%xmm2
mulsd  %xmm3,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm14,%xmm10
mulsd  %xmm11,%xmm10
divsd  %xmm10,%xmm2
movapd %xmm2,%xmm3
mulsd  %xmm12,%xmm3
addsd  %xmm8,%xmm3
movapd %xmm3,%xmm6
andpd  %xmm15,%xmm6
ucomisd %xmm6,%xmm7
movapd %xmm9,%xmm6
ja     1252 <betacf+0x152>
movapd %xmm8,%xmm6
divsd  %xmm3,%xmm6
divsd  %xmm4,%xmm5
addsd  %xmm8,%xmm5
movapd %xmm5,%xmm3
andpd  %xmm15,%xmm3
cmpltsd %xmm7,%xmm3
movapd %xmm3,%xmm4
andnpd %xmm5,%xmm4
andpd  %xmm7,%xmm3
orpd   %xmm4,%xmm3
mulsd  %xmm3,%xmm12
mulsd  %xmm12,%xmm1
divsd  %xmm3,%xmm2
addsd  %xmm8,%xmm2
movapd %xmm2,%xmm4
andpd  %xmm15,%xmm4
cmpltsd %xmm7,%xmm4
movapd %xmm4,%xmm3
andnpd %xmm2,%xmm3
andpd  %xmm7,%xmm4
orpd   %xmm3,%xmm4
movapd %xmm6,%xmm2
mulsd  %xmm4,%xmm2
mulsd  %xmm2,%xmm1
addsd  0xd4a(%rip),%xmm2
andpd  0xd6a(%rip),%xmm2
movsd  0xd52(%rip),%xmm3
ucomisd %xmm2,%xmm3
ja     12ef <betacf+0x1ef>
add    $0x1,%ecx
lea    0x2(%rdx),%eax
add    $0xffffffff,%r8d
add    $0xffffffff,%esi
cmp    $0xc8,%edx
mov    %eax,%edx
jne    11a0 <betacf+0xa0>
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @betacf(i32 %0, i32 %1, double %2) {
%4 = alloca i32
%5 = alloca i32
%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 double
%15 = alloca double
%16 = alloca double
store i32 %0, i32* %4
store i32 %1, i32* %5
store double %2, double* %6
%17 = load i32, i32* %4
%18 = load i32, i32* %5
%19 = add nsw i32 %17, %18
%20 = sitofp i32 %19 to double
store double %20, double* %14
%21 = load i32, i32* %4
%22 = sitofp i32 %21 to double
%23 = fadd double %22, 1.000000e+00
store double %23, double* %16
%24 = load i32, i32* %4
%25 = sitofp i32 %24 to double
%26 = fsub double %25, 1.000000e+00
store double %26, double* %15
store double 1.000000e+00, double* %10
%27 = load double, double* %14
%28 = load double, double* %6
%29 = fmul double %27, %28
%30 = load double, double* %16
%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 @llvm.fabs.f64(double %33)
%35 = fcmp olt double %34, 1.000000e-30
br i1 %35, label %36, label %37
36: ; preds = %3
store double 1.000000e-30, 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 = %142, %37
%42 = load i32, i32* %7
%43 = icmp sle i32 %42, 100
br i1 %43, label %44, label %145
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 = load i32, i32* %5
%49 = load i32, i32* %7
%50 = sub nsw i32 %48, %49
%51 = mul nsw i32 %47, %50
%52 = sitofp i32 %51 to double
%53 = load double, double* %6
%54 = fmul double %52, %53
%55 = load double, double* %15
%56 = load i32, i32* %8
%57 = sitofp i32 %56 to double
%58 = fadd double %55, %57
%59 = load i32, i32* %4
%60 = load i32, i32* %8
%61 = add nsw i32 %59, %60
%62 = sitofp i32 %61 to double
%63 = fmul double %58, %62
%64 = fdiv double %54, %63
store double %64, double* %9
%65 = load double, double* %9
%66 = load double, double* %11
%67 = call double @llvm.fmuladd.f64(double %65, double %66, double 1.000000e+00)
store double %67, double* %11
%68 = load double, double* %11
%69 = call double @llvm.fabs.f64(double %68)
%70 = fcmp olt double %69, 1.000000e-30
br i1 %70, label %71, label %72
71: ; preds = %44
store double 1.000000e-30, double* %11
br label %72
72: ; preds = %71, %44
%73 = load double, double* %9
%74 = load double, double* %10
%75 = fdiv double %73, %74
%76 = fadd double 1.000000e+00, %75
store double %76, double* %10
%77 = load double, double* %10
%78 = call double @llvm.fabs.f64(double %77)
%79 = fcmp olt double %78, 1.000000e-30
br i1 %79, label %80, label %81
80: ; preds = %72
store double 1.000000e-30, double* %10
br label %81
81: ; preds = %80, %72
%82 = load double, double* %11
%83 = fdiv double 1.000000e+00, %82
store double %83, double* %11
%84 = load double, double* %11
%85 = load double, double* %10
%86 = fmul double %84, %85
%87 = load double, double* %13
%88 = fmul double %87, %86
store double %88, double* %13
%89 = load i32, i32* %4
%90 = load i32, i32* %7
%91 = add nsw i32 %89, %90
%92 = sub nsw i32 0, %91
%93 = sitofp i32 %92 to double
%94 = load double, double* %14
%95 = load i32, i32* %7
%96 = sitofp i32 %95 to double
%97 = fadd double %94, %96
%98 = fmul double %93, %97
%99 = load double, double* %6
%100 = fmul double %98, %99
%101 = load i32, i32* %4
%102 = load i32, i32* %8
%103 = add nsw i32 %101, %102
%104 = sitofp i32 %103 to double
%105 = load double, double* %16
%106 = load i32, i32* %8
%107 = sitofp i32 %106 to double
%108 = fadd double %105, %107
%109 = fmul double %104, %108
%110 = fdiv double %100, %109
store double %110, double* %9
%111 = load double, double* %9
%112 = load double, double* %11
%113 = call double @llvm.fmuladd.f64(double %111, double %112, double 1.000000e+00)
store double %113, double* %11
%114 = load double, double* %11
%115 = call double @llvm.fabs.f64(double %114)
%116 = fcmp olt double %115, 1.000000e-30
br i1 %116, label %117, label %118
117: ; preds = %81
store double 1.000000e-30, double* %11
br label %118
118: ; preds = %117, %81
%119 = load double, double* %9
%120 = load double, double* %10
%121 = fdiv double %119, %120
%122 = fadd double 1.000000e+00, %121
store double %122, double* %10
%123 = load double, double* %10
%124 = call double @llvm.fabs.f64(double %123)
%125 = fcmp olt double %124, 1.000000e-30
br i1 %125, label %126, label %127
126: ; preds = %118
store double 1.000000e-30, double* %10
br label %127
127: ; preds = %126, %118
%128 = load double, double* %11
%129 = fdiv double 1.000000e+00, %128
store double %129, double* %11
%130 = load double, double* %11
%131 = load double, double* %10
%132 = fmul double %130, %131
store double %132, double* %12
%133 = load double, double* %12
%134 = load double, double* %13
%135 = fmul double %134, %133
store double %135, double* %13
%136 = load double, double* %12
%137 = fsub double %136, 1.000000e+00
%138 = call double @llvm.fabs.f64(double %137)
%139 = fcmp olt double %138, 3.000000e-07
br i1 %139, label %140, label %141
140: ; preds = %127
br label %145
141: ; preds = %127
br label %142
142: ; preds = %141
%143 = load i32, i32* %7
%144 = add nsw i32 %143, 1
store i32 %144, i32* %7
br label %41
145: ; preds = %140, %41
%146 = load double, double* %13
ret double %146
}

LLVM IR (-O2):
define double @betacf(i32 %0, i32 %1, double %2) {
%4 = add nsw i32 %1, %0
%5 = sitofp i32 %4 to double
%6 = sitofp i32 %0 to double
%7 = fadd double %6, 1.000000e+00
%8 = fadd double %6, -1.000000e+00
%9 = fmul double %5, %2
%10 = fdiv double %9, %7
%11 = fsub double 1.000000e+00, %10
%12 = tail call double @llvm.fabs.f64(double %11)
%13 = fcmp olt double %12, 1.000000e-30
%14 = fdiv double 1.000000e+00, %11
%15 = select i1 %13, double 0x46293E5939A08CE9, double %14
br label %16
16: ; preds = %16, %3
%17 = phi i32 [ 1, %3 ], [ %69, %16 ]
%18 = phi double [ 1.000000e+00, %3 ], [ %61, %16 ]
%19 = phi double [ %15, %3 ], [ %65, %16 ]
%20 = phi double [ %15, %3 ], [ %63, %16 ]
%21 = shl nuw nsw i32 %17, 1
%22 = sub nsw i32 %1, %17
%23 = mul nsw i32 %22, %17
%24 = sitofp i32 %23 to double
%25 = fmul double %24, %2
%26 = sitofp i32 %21 to double
%27 = fadd double %8, %26
%28 = add nsw i32 %21, %0
%29 = sitofp i32 %28 to double
%30 = fmul double %27, %29
%31 = fdiv double %25, %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, 1.000000e-30
%35 = fdiv double %31, %18
%36 = fadd double %35, 1.000000e+00
%37 = tail call double @llvm.fabs.f64(double %36)
%38 = fcmp olt double %37, 1.000000e-30
%39 = select i1 %38, double 1.000000e-30, double %36
%40 = fdiv double 1.000000e+00, %32
%41 = select i1 %34, double 0x46293E5939A08CE9, double %40
%42 = fmul double %41, %39
%43 = fmul double %19, %42
%44 = add nsw i32 %17, %0
%45 = sub nsw i32 0, %44
%46 = sitofp i32 %45 to double
%47 = sitofp i32 %17 to double
%48 = fadd double %5, %47
%49 = fmul double %48, %46
%50 = fmul double %49, %2
%51 = fadd double %7, %26
%52 = fmul double %51, %29
%53 = fdiv double %50, %52
%54 = tail call double @llvm.fmuladd.f64(double %53, double %41, double 1.000000e+00)
%55 = tail call double @llvm.fabs.f64(double %54)
%56 = fcmp olt double %55, 1.000000e-30
%57 = fdiv double %53, %39
%58 = fadd double %57, 1.000000e+00
%59 = tail call double @llvm.fabs.f64(double %58)
%60 = fcmp olt double %59, 1.000000e-30
%61 = select i1 %60, double 1.000000e-30, double %58
%62 = fdiv double 1.000000e+00, %54
%63 = select i1 %56, double 0x46293E5939A08CE9, double %62
%64 = fmul double %63, %61
%65 = fmul double %43, %64
%66 = fadd double %64, -1.000000e+00
%67 = tail call double @llvm.fabs.f64(double %66)
%68 = fcmp olt double %67, 3.000000e-07
%69 = add nuw nsw i32 %17, 1
%70 = icmp eq i32 %69, 101
%71 = select i1 %68, i1 true, i1 %70
br i1 %71, label %72, label %16
72: ; preds = %16
ret double %65
}