void cdflib_cumfnc(const double *f, const double *dfn, const double *dfd,
                   const double *pnonc, double *cum, double *ccum) {
  static const double half = 0.5e0;
  static const double done = 1.0e0;
  static const double eps = 1.0e-4;
  double dsum, dummy, prod, xx, yy, adn, aup, b, betdn, betup, centwt, dnterm,
      sum, upterm, xmult, xnonc;
  int i, icent, ierr;
  double T1, T2, T3, T4, T5, T6;
  if (!(*f <= 0.0e0)) {
    goto S10;
  }
  *cum = 0.0e0;
  *ccum = 1.0e0;
  return;
S10:
  if (!(*pnonc < 1.0e-10)) {
    goto S20;
  }
  cdflib_cumf(f, dfn, dfd, cum, ccum);
  return;
S20:
  xnonc = *pnonc / 2.0e0;
  icent = (int)xnonc;
  if (icent == 0) {
    icent = 1;
  }
  T1 = (double)(icent + 1);
  centwt = exp(-xnonc + (double)icent * log(xnonc) - cdflib_gamma_log(&T1));
  prod = *dfn * *f;
  dsum = *dfd + prod;
  yy = *dfd / dsum;
  if (yy > half) {
    xx = prod / dsum;
    yy = done - xx;
  } else {
    xx = done - yy;
  }
  T2 = *dfn * half + (double)icent;
  T3 = *dfd * half;
  cdflib_beta_inc(&T2, &T3, &xx, &yy, &betdn, &dummy, &ierr);
  adn = *dfn / 2.0e0 + (double)icent;
  aup = adn;
  b = *dfd / 2.0e0;
  betup = betdn;
  sum = centwt * betdn;
  xmult = centwt;
  i = icent;
  T4 = adn + b;
  T5 = adn + 1.0e0;
  dnterm = exp(cdflib_gamma_log(&T4) - cdflib_gamma_log(&T5) -
               cdflib_gamma_log(&b) + adn * log(xx) + b * log(yy));
S30:
  if ((int)(sum < 1.0e-20 || (xmult * betdn) < eps * sum) || i <= 0) {
    goto S40;
  }
  xmult *= ((double)i / xnonc);
  i -= 1;
  adn -= 1.0;
  dnterm = (adn + 1.0) / ((adn + b) * xx) * dnterm;
  betdn += dnterm;
  sum += (xmult * betdn);
  goto S30;
S40:
  i = icent + 1;
  xmult = centwt;
  if (aup - 1.0 + b == 0)
    upterm = exp(-cdflib_gamma_log(&aup) - cdflib_gamma_log(&b) +
                 (aup - 1.0) * log(xx) + b * log(yy));
  else {
    T6 = aup - 1.0 + b;
    upterm = exp(cdflib_gamma_log(&T6) - cdflib_gamma_log(&aup) -
                 cdflib_gamma_log(&b) + (aup - 1.0) * log(xx) + b * log(yy));
  }
  goto S60;
S50:
  if ((int)(sum < 1.0e-20 || (xmult * betup) < eps * sum)) {
    goto S70;
  }
S60:
  xmult *= (xnonc / (double)i);
  i += 1;
  aup += 1.0;
  upterm = (aup + b - 2.0e0) * xx / (aup - 1.0) * upterm;
  betup -= upterm;
  sum += (xmult * betup);
  goto S50;
S70:
  *cum = sum;
  *ccum = 0.5e0 + (0.5e0 - *cum);
  return;
}


x86 assembly (clang, -O0):
<cdflib_cumfnc>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1a0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jae    1189 <cdflib_cumfnc+0x39>
jmp    11a9 <cdflib_cumfnc+0x59>
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
movsd  0xe78(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    19b3 <cdflib_cumfnc+0x863>
mov    -0x20(%rbp),%rax
movsd  0xe4b(%rip),%xmm0
ucomisd (%rax),%xmm0
ja     11c4 <cdflib_cumfnc+0x74>
jmp    11e4 <cdflib_cumfnc+0x94>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x28(%rbp),%rcx
mov    -0x30(%rbp),%r8
mov    $0x0,%al
call   1040 <cdflib_cumf@plt>
jmp    19b3 <cdflib_cumfnc+0x863>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xe14(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
cvttsd2si -0xb0(%rbp),%eax
ff
mov    %eax,-0xb8(%rbp)
cmpl   $0x0,-0xb8(%rbp)
jne    1225 <cdflib_cumfnc+0xd5>
movl   $0x1,-0xb8(%rbp)
mov    -0xb8(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
cvtsi2sdl -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
call   1060 <log@plt>
movsd  -0x108(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x100(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
lea    -0xc8(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
movsd  -0xf8(%rbp),%xmm0
ff
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  0xcfe(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1342 <cdflib_cumfnc+0x1f2>
movsd  -0x48(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xce5(%rip),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    1354 <cdflib_cumfnc+0x204>
movsd  0xcce(%rip),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
cvtsi2sdl -0xb8(%rbp),%xmm1
ff
movsd  0xca4(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0xc88(%rip),%xmm0
mulsd  (%rax),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
lea    -0xd0(%rbp),%rdi
lea    -0xd8(%rbp),%rsi
lea    -0x50(%rbp),%rdx
lea    -0x58(%rbp),%rcx
lea    -0x78(%rbp),%r8
lea    -0x40(%rbp),%r9
lea    -0xbc(%rbp),%rax
mov    %rax,(%rsp)
mov    $0x0,%al
call   1030 <cdflib_beta_inc@plt>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xc34(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsi2sdl -0xb8(%rbp),%xmm1
ff
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xc05(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0xb8(%rbp),%eax
mov    %eax,-0xb4(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  0xbb7(%rip),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
lea    -0xe0(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,-0x130(%rbp)
lea    -0xe8(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ecx
mov    -0x130(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x12c(%rbp)
lea    -0x70(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ecx
mov    -0x12c(%rbp),%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
call   1060 <log@plt>
movsd  -0x128(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x120(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
call   1060 <log@plt>
movsd  -0x118(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x110(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  %xmm0,-0x90(%rbp)
ff
mov    $0x1,%al
movsd  0xadc(%rip),%xmm0
ucomisd -0x98(%rbp),%xmm0
ff
mov    %al,-0x131(%rbp)
ja     1582 <cdflib_cumfnc+0x432>
movsd  -0xa8(%rbp),%xmm1
ff
mulsd  -0x78(%rbp),%xmm1
movsd  0xabb(%rip),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x131(%rbp)
mov    -0x131(%rbp),%al
test   $0x1,%al
jne    159d <cdflib_cumfnc+0x44d>
cmpl   $0x0,-0xb4(%rbp)
jg     15a2 <cdflib_cumfnc+0x452>
jmp    1653 <cdflib_cumfnc+0x503>
cvtsi2sdl -0xb4(%rbp),%xmm0
ff
divsd  -0xb0(%rbp),%xmm0
ff
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0xb4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0xb4(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  0xa3a(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xa29(%rip),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm1
addsd  -0x70(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0x78(%rbp),%xmm2
movsd  -0x98(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
jmp    153a <cdflib_cumfnc+0x3ea>
mov    -0xb8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  0x999(%rip),%xmm1
subsd  %xmm1,%xmm0
addsd  -0x70(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    175d <cdflib_cumfnc+0x60d>
jp     175d <cdflib_cumfnc+0x60d>
lea    -0x68(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x15c(%rbp)
lea    -0x70(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ecx
mov    -0x15c(%rbp),%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x158(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  0x938(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x150(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
call   1060 <log@plt>
movsd  -0x158(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x150(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x148(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x140(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
call   1060 <log@plt>
movsd  -0x148(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x140(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  %xmm0,-0xa0(%rbp)
ff
jmp    1850 <cdflib_cumfnc+0x700>
movsd  -0x68(%rbp),%xmm0
movsd  0x8ae(%rip),%xmm1
subsd  %xmm1,%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
lea    -0xf0(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,-0x188(%rbp)
lea    -0x68(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ecx
mov    -0x188(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x184(%rbp)
lea    -0x70(%rbp),%rdi
mov    $0x0,%al
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ecx
mov    -0x184(%rbp),%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x180(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  0x840(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x178(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
call   1060 <log@plt>
movsd  -0x180(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x178(%rbp),%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x170(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x168(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
call   1060 <log@plt>
movsd  -0x170(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x168(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  %xmm0,-0xa0(%rbp)
ff
jmp    18ba <cdflib_cumfnc+0x76a>
mov    $0x1,%al
movsd  0x7c1(%rip),%xmm0
ucomisd -0x98(%rbp),%xmm0
ff
mov    %al,-0x189(%rbp)
ja     189d <cdflib_cumfnc+0x74d>
movsd  -0xa8(%rbp),%xmm1
ff
mulsd  -0x80(%rbp),%xmm1
movsd  0x7a0(%rip),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x189(%rbp)
mov    -0x189(%rbp),%al
test   $0x1,%al
jne    18b0 <cdflib_cumfnc+0x760>
jmp    18b5 <cdflib_cumfnc+0x765>
jmp    197f <cdflib_cumfnc+0x82f>
jmp    18ba <cdflib_cumfnc+0x76a>
movsd  -0xb0(%rbp),%xmm0
ff
cvtsi2sdl -0xb4(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
movsd  0x723(%rip),%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  0x6f7(%rip),%xmm1
subsd  %xmm1,%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm1
movsd  0x6f1(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm1
ff
movsd  -0x80(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm2
movsd  -0x98(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
jmp    1855 <cdflib_cumfnc+0x705>
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  0x675(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x669(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x1a0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cdflib_cumfnc>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xb8,%rsp
xorpd  %xmm0,%xmm0
ucomisd (%rdi),%xmm0
jae    11ad <cdflib_cumfnc+0x5d>
mov    %rdx,%r12
mov    %rsi,%r13
mov    %rdi,%rbx
movsd  (%rcx),%xmm0
movsd  0xe80(%rip),%xmm1
ucomisd %xmm0,%xmm1
jbe    11c6 <cdflib_cumfnc+0x76>
mov    %rbx,%rdi
mov    %r13,%rsi
mov    %r12,%rdx
mov    %r8,%rcx
mov    %r9,%r8
xor    %eax,%eax
add    $0xb8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <cdflib_cumf@plt>
movq   $0x0,(%r8)
movabs $0x3ff0000000000000,%rax
mov    %rax,(%r9)
jmp    16d7 <cdflib_cumfnc+0x587>
mov    %r8,0x70(%rsp)
mov    %r9,0x48(%rsp)
mulsd  0xe30(%rip),%xmm0
movsd  %xmm0,0x50(%rsp)
cvttsd2si %xmm0,%eax
test   %eax,%eax
mov    $0x1,%r15d
cmovne %eax,%r15d
lea    0x1(%r15),%r14d
xorps  %xmm1,%xmm1
cvtsi2sd %r14d,%xmm1
movsd  %xmm1,0xa0(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sd %r15d,%xmm1
movsd  %xmm1,0x20(%rsp)
call   1060 <log@plt>
mulsd  0x20(%rsp),%xmm0
subsd  0x50(%rsp),%xmm0
movsd  %xmm0,0x38(%rsp)
lea    0xa0(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x38(%rsp),%xmm0
subsd  %xmm1,%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x38(%rsp)
movsd  0x0(%r13),%xmm0
movsd  (%rbx),%xmm2
mulsd  %xmm0,%xmm2
movsd  (%r12),%xmm1
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
movapd %xmm1,%xmm4
divsd  %xmm3,%xmm4
movsd  %xmm4,0x40(%rsp)
ucomisd 0xd83(%rip),%xmm4
jbe    12a5 <cdflib_cumfnc+0x155>
divsd  %xmm3,%xmm2
movsd  %xmm2,0x30(%rsp)
movsd  0xd77(%rip),%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,0x40(%rsp)
jmp    12b7 <cdflib_cumfnc+0x167>
movsd  0xd63(%rip),%xmm2
subsd  %xmm4,%xmm2
movsd  %xmm2,0x30(%rsp)
movsd  0xd49(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  0x20(%rsp),%xmm0
movsd  %xmm0,0x98(%rsp)
mulsd  %xmm2,%xmm1
movsd  %xmm1,0x90(%rsp)
lea    0xac(%rsp),%rax
mov    %rax,(%rsp)
lea    0x98(%rsp),%rdi
lea    0x90(%rsp),%rsi
lea    0x30(%rsp),%rdx
lea    0x40(%rsp),%rcx
lea    0x58(%rsp),%r8
lea    0xb0(%rsp),%r9
xor    %eax,%eax
call   1030 <cdflib_beta_inc@plt>
movsd  0x0(%r13),%xmm3
movsd  0xce1(%rip),%xmm1
mulsd  %xmm1,%xmm3
addsd  0x20(%rsp),%xmm3
movsd  %xmm3,0x10(%rsp)
movsd  %xmm3,0x28(%rsp)
movsd  (%r12),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0x58(%rsp),%xmm1
movsd  0x38(%rsp),%xmm2
movsd  %xmm1,0x68(%rsp)
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x20(%rsp)
addsd  %xmm3,%xmm0
movsd  %xmm0,0x88(%rsp)
movsd  0xc92(%rip),%xmm0
addsd  %xmm3,%xmm0
movsd  %xmm0,0x80(%rsp)
lea    0x88(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ebx
lea    0x80(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ebp
lea    0x18(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
add    %ebp,%eax
sub    %eax,%ebx
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x30(%rsp),%xmm0
call   1060 <log@plt>
mulsd  0x10(%rsp),%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x60(%rsp)
movsd  0x40(%rsp),%xmm0
call   1060 <log@plt>
mulsd  0x60(%rsp),%xmm0
addsd  0x8(%rsp),%xmm0
call   1050 <exp@plt>
movsd  0x20(%rsp),%xmm3
movsd  0xbfa(%rip),%xmm12
ucomisd %xmm3,%xmm12
jbe    1431 <cdflib_cumfnc+0x2e1>
movsd  0x18(%rsp),%xmm11
jmp    14f1 <cdflib_cumfnc+0x3a1>
movsd  0x58(%rsp),%xmm2
movsd  0x18(%rsp),%xmm11
movsd  0x30(%rsp),%xmm13
movsd  0xbd2(%rip),%xmm10
movsd  0xbd1(%rip),%xmm8
movsd  0xbb0(%rip),%xmm9
movsd  0x38(%rsp),%xmm5
movsd  0x50(%rsp),%xmm4
movsd  0x10(%rsp),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm5,%xmm7
mulsd  %xmm2,%xmm7
movapd %xmm3,%xmm6
mulsd  %xmm10,%xmm6
ucomisd %xmm7,%xmm6
ja     14f1 <cdflib_cumfnc+0x3a1>
test   %r15d,%r15d
jle    14f1 <cdflib_cumfnc+0x3a1>
xorps  %xmm6,%xmm6
cvtsi2sd %r15d,%xmm6
divsd  %xmm4,%xmm6
mulsd  %xmm6,%xmm5
add    $0xffffffff,%r15d
addsd  %xmm8,%xmm1
movapd %xmm1,%xmm6
addsd  %xmm9,%xmm6
movapd %xmm1,%xmm7
addsd  %xmm11,%xmm7
mulsd  %xmm13,%xmm7
divsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm0
addsd  %xmm0,%xmm2
movsd  %xmm2,0x58(%rsp)
movapd %xmm5,%xmm6
mulsd  %xmm2,%xmm6
addsd  %xmm6,%xmm3
ucomisd %xmm3,%xmm12
jbe    1480 <cdflib_cumfnc+0x330>
mov    0x48(%rsp),%r15
movsd  0x28(%rsp),%xmm0
addsd  0xb24(%rip),%xmm0
addsd  %xmm11,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
movsd  %xmm3,0x20(%rsp)
jne    1542 <cdflib_cumfnc+0x3f2>
jp     1542 <cdflib_cumfnc+0x3f2>
lea    0x28(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ebx
lea    0x18(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
add    %ebx,%eax
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
jmp    157b <cdflib_cumfnc+0x42b>
movsd  %xmm0,0x78(%rsp)
lea    0x78(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ebx
lea    0x28(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
mov    %eax,%ebp
lea    0x18(%rsp),%rdi
xor    %eax,%eax
call   1070 <cdflib_gamma_log@plt>
add    %ebp,%eax
sub    %eax,%ebx
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x28(%rsp),%xmm0
addsd  0xa99(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x30(%rsp),%xmm0
call   1060 <log@plt>
mulsd  0x8(%rsp),%xmm0
addsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x40(%rsp),%xmm0
call   1060 <log@plt>
mulsd  0x8(%rsp),%xmm0
addsd  0x10(%rsp),%xmm0
call   1050 <exp@plt>
movsd  0x18(%rsp),%xmm12
movsd  0x30(%rsp),%xmm13
movsd  0x28(%rsp),%xmm2
movsd  0xa19(%rip),%xmm8
movsd  0xa30(%rip),%xmm10
movsd  0xa1f(%rip),%xmm11
movsd  0xa0e(%rip),%xmm9
movsd  0x50(%rsp),%xmm14
movsd  0x38(%rsp),%xmm5
movsd  0x20(%rsp),%xmm3
movsd  0x68(%rsp),%xmm1
movsd  0x9e5(%rip),%xmm6
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd %xmm5,%xmm4
movapd %xmm0,%xmm7
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
movapd %xmm14,%xmm5
divsd  %xmm0,%xmm5
mulsd  %xmm4,%xmm5
addsd  %xmm8,%xmm2
movapd %xmm2,%xmm0
addsd  %xmm12,%xmm0
addsd  %xmm10,%xmm0
mulsd  %xmm13,%xmm0
movapd %xmm2,%xmm4
addsd  %xmm11,%xmm4
divsd  %xmm4,%xmm0
mulsd  %xmm7,%xmm0
subsd  %xmm0,%xmm1
movapd %xmm5,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm4,%xmm3
ucomisd %xmm3,%xmm6
ja     16af <cdflib_cumfnc+0x55f>
add    $0x1,%r14d
movapd %xmm3,%xmm7
mulsd  %xmm9,%xmm7
ucomisd %xmm4,%xmm7
jbe    1640 <cdflib_cumfnc+0x4f0>
movsd  %xmm2,0x28(%rsp)
mov    0x70(%rsp),%rax
movsd  %xmm3,(%rax)
movsd  0x942(%rip),%xmm1
movapd %xmm1,%xmm0
subsd  %xmm3,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r15)
add    $0xb8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@cdflib_cumfnc.half = internal constant double 5.000000e-01
@cdflib_cumfnc.done = internal constant double 1.000000e+00
@cdflib_cumfnc.eps = internal constant double 1.000000e-04
define void @cdflib_cumfnc(double* %0, double* %1, double* %2, double* %3, double* %4, double* %5) {
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
%12 = alloca double*
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca double
store double* %0, double** %7
store double* %1, double** %8
store double* %2, double** %9
store double* %3, double** %10
store double* %4, double** %11
store double* %5, double** %12
%38 = load double*, double** %7
%39 = load double, double* %38
%40 = fcmp ole double %39, 0.000000e+00
br i1 %40, label %42, label %41
41: ; preds = %6
br label %45
42: ; preds = %6
%43 = load double*, double** %11
store double 0.000000e+00, double* %43
%44 = load double*, double** %12
store double 1.000000e+00, double* %44
br label %293
45: ; preds = %41
%46 = load double*, double** %10
%47 = load double, double* %46
%48 = fcmp olt double %47, 1.000000e-10
br i1 %48, label %50, label %49
49: ; preds = %45
br label %57
50: ; preds = %45
%51 = load double*, double** %7
%52 = load double*, double** %8
%53 = load double*, double** %9
%54 = load double*, double** %11
%55 = load double*, double** %12
%56 = call i32 (double*, double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_cumf to i32 (double*, double*, double*, double*, double*, ...)*)(double* %51, double* %52, double* %53, double* %54, double* %55)
br label %293
57: ; preds = %49
%58 = load double*, double** %10
%59 = load double, double* %58
%60 = fdiv double %59, 2.000000e+00
store double %60, double* %28
%61 = load double, double* %28
%62 = fptosi double %61 to i32
store i32 %62, i32* %30
%63 = load i32, i32* %30
%64 = icmp eq i32 %63, 0
br i1 %64, label %65, label %66
65: ; preds = %57
store i32 1, i32* %30
br label %66
66: ; preds = %65, %57
%67 = load i32, i32* %30
%68 = add nsw i32 %67, 1
%69 = sitofp i32 %68 to double
store double %69, double* %32
%70 = load double, double* %28
%71 = fneg double %70
%72 = load i32, i32* %30
%73 = sitofp i32 %72 to double
%74 = load double, double* %28
%75 = call double @log(double %74)
%76 = call double @llvm.fmuladd.f64(double %73, double %75, double %71)
%77 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %32)
%78 = sitofp i32 %77 to double
%79 = fsub double %76, %78
%80 = call double @exp(double %79)
store double %80, double* %23
%81 = load double*, double** %8
%82 = load double, double* %81
%83 = load double*, double** %7
%84 = load double, double* %83
%85 = fmul double %82, %84
store double %85, double* %15
%86 = load double*, double** %9
%87 = load double, double* %86
%88 = load double, double* %15
%89 = fadd double %87, %88
store double %89, double* %13
%90 = load double*, double** %9
%91 = load double, double* %90
%92 = load double, double* %13
%93 = fdiv double %91, %92
store double %93, double* %17
%94 = load double, double* %17
%95 = fcmp ogt double %94, 5.000000e-01
br i1 %95, label %96, label %102
96: ; preds = %66
%97 = load double, double* %15
%98 = load double, double* %13
%99 = fdiv double %97, %98
store double %99, double* %16
%100 = load double, double* %16
%101 = fsub double 1.000000e+00, %100
store double %101, double* %17
br label %105
102: ; preds = %66
%103 = load double, double* %17
%104 = fsub double 1.000000e+00, %103
store double %104, double* %16
br label %105
105: ; preds = %102, %96
%106 = load double*, double** %8
%107 = load double, double* %106
%108 = load i32, i32* %30
%109 = sitofp i32 %108 to double
%110 = call double @llvm.fmuladd.f64(double %107, double 5.000000e-01, double %109)
store double %110, double* %33
%111 = load double*, double** %9
%112 = load double, double* %111
%113 = fmul double %112, 5.000000e-01
store double %113, double* %34
%114 = call i32 (double*, double*, double*, double*, double*, double*, i32*, ...) bitcast (i32 (...)* @cdflib_beta_inc to i32 (double*, double*, double*, double*, double*, double*, i32*, ...)*)(double* %33, double* %34, double* %16, double* %17, double* %21, double* %14, i32* %31)
%115 = load double*, double** %8
%116 = load double, double* %115
%117 = fdiv double %116, 2.000000e+00
%118 = load i32, i32* %30
%119 = sitofp i32 %118 to double
%120 = fadd double %117, %119
store double %120, double* %18
%121 = load double, double* %18
store double %121, double* %19
%122 = load double*, double** %9
%123 = load double, double* %122
%124 = fdiv double %123, 2.000000e+00
store double %124, double* %20
%125 = load double, double* %21
store double %125, double* %22
%126 = load double, double* %23
%127 = load double, double* %21
%128 = fmul double %126, %127
store double %128, double* %25
%129 = load double, double* %23
store double %129, double* %27
%130 = load i32, i32* %30
store i32 %130, i32* %29
%131 = load double, double* %18
%132 = load double, double* %20
%133 = fadd double %131, %132
store double %133, double* %35
%134 = load double, double* %18
%135 = fadd double %134, 1.000000e+00
store double %135, double* %36
%136 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %35)
%137 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %36)
%138 = sub nsw i32 %136, %137
%139 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %20)
%140 = sub nsw i32 %138, %139
%141 = sitofp i32 %140 to double
%142 = load double, double* %18
%143 = load double, double* %16
%144 = call double @log(double %143)
%145 = call double @llvm.fmuladd.f64(double %142, double %144, double %141)
%146 = load double, double* %20
%147 = load double, double* %17
%148 = call double @log(double %147)
%149 = call double @llvm.fmuladd.f64(double %146, double %148, double %145)
%150 = call double @exp(double %149)
store double %150, double* %24
br label %151
151: ; preds = %167, %105
%152 = load double, double* %25
%153 = fcmp olt double %152, 0x3BC79CA10C924223
br i1 %153, label %161, label %154
154: ; preds = %151
%155 = load double, double* %27
%156 = load double, double* %21
%157 = fmul double %155, %156
%158 = load double, double* %25
%159 = fmul double 1.000000e-04, %158
%160 = fcmp olt double %157, %159
br label %161
161: ; preds = %154, %151
%162 = phi i1 [ true, %151 ], [ %160, %154 ]
br i1 %162, label %166, label %163
163: ; preds = %161
%164 = load i32, i32* %29
%165 = icmp sle i32 %164, 0
br i1 %165, label %166, label %167
166: ; preds = %163, %161
br label %195
167: ; preds = %163
%168 = load i32, i32* %29
%169 = sitofp i32 %168 to double
%170 = load double, double* %28
%171 = fdiv double %169, %170
%172 = load double, double* %27
%173 = fmul double %172, %171
store double %173, double* %27
%174 = load i32, i32* %29
%175 = sub nsw i32 %174, 1
store i32 %175, i32* %29
%176 = load double, double* %18
%177 = fsub double %176, 1.000000e+00
store double %177, double* %18
%178 = load double, double* %18
%179 = fadd double %178, 1.000000e+00
%180 = load double, double* %18
%181 = load double, double* %20
%182 = fadd double %180, %181
%183 = load double, double* %16
%184 = fmul double %182, %183
%185 = fdiv double %179, %184
%186 = load double, double* %24
%187 = fmul double %185, %186
store double %187, double* %24
%188 = load double, double* %24
%189 = load double, double* %21
%190 = fadd double %189, %188
store double %190, double* %21
%191 = load double, double* %27
%192 = load double, double* %21
%193 = load double, double* %25
%194 = call double @llvm.fmuladd.f64(double %191, double %192, double %193)
store double %194, double* %25
br label %151
195: ; preds = %166
%196 = load i32, i32* %30
%197 = add nsw i32 %196, 1
store i32 %197, i32* %29
%198 = load double, double* %23
store double %198, double* %27
%199 = load double, double* %19
%200 = fsub double %199, 1.000000e+00
%201 = load double, double* %20
%202 = fadd double %200, %201
%203 = fcmp oeq double %202, 0.000000e+00
br i1 %203, label %204, label %220
204: ; preds = %195
%205 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %19)
%206 = sub nsw i32 0, %205
%207 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %20)
%208 = sub nsw i32 %206, %207
%209 = sitofp i32 %208 to double
%210 = load double, double* %19
%211 = fsub double %210, 1.000000e+00
%212 = load double, double* %16
%213 = call double @log(double %212)
%214 = call double @llvm.fmuladd.f64(double %211, double %213, double %209)
%215 = load double, double* %20
%216 = load double, double* %17
%217 = call double @log(double %216)
%218 = call double @llvm.fmuladd.f64(double %215, double %217, double %214)
%219 = call double @exp(double %218)
store double %219, double* %26
br label %241
220: ; preds = %195
%221 = load double, double* %19
%222 = fsub double %221, 1.000000e+00
%223 = load double, double* %20
%224 = fadd double %222, %223
store double %224, double* %37
%225 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %37)
%226 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %19)
%227 = sub nsw i32 %225, %226
%228 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %20)
%229 = sub nsw i32 %227, %228
%230 = sitofp i32 %229 to double
%231 = load double, double* %19
%232 = fsub double %231, 1.000000e+00
%233 = load double, double* %16
%234 = call double @log(double %233)
%235 = call double @llvm.fmuladd.f64(double %232, double %234, double %230)
%236 = load double, double* %20
%237 = load double, double* %17
%238 = call double @log(double %237)
%239 = call double @llvm.fmuladd.f64(double %236, double %238, double %235)
%240 = call double @exp(double %239)
store double %240, double* %26
br label %241
241: ; preds = %220, %204
br label %256
242: ; preds = %256
%243 = load double, double* %25
%244 = fcmp olt double %243, 0x3BC79CA10C924223
br i1 %244, label %252, label %245
245: ; preds = %242
%246 = load double, double* %27
%247 = load double, double* %22
%248 = fmul double %246, %247
%249 = load double, double* %25
%250 = fmul double 1.000000e-04, %249
%251 = fcmp olt double %248, %250
br label %252
252: ; preds = %245, %242
%253 = phi i1 [ true, %242 ], [ %251, %245 ]
br i1 %253, label %254, label %255
254: ; preds = %252
br label %285
255: ; preds = %252
br label %256
256: ; preds = %255, %241
%257 = load double, double* %28
%258 = load i32, i32* %29
%259 = sitofp i32 %258 to double
%260 = fdiv double %257, %259
%261 = load double, double* %27
%262 = fmul double %261, %260
store double %262, double* %27
%263 = load i32, i32* %29
%264 = add nsw i32 %263, 1
store i32 %264, i32* %29
%265 = load double, double* %19
%266 = fadd double %265, 1.000000e+00
store double %266, double* %19
%267 = load double, double* %19
%268 = load double, double* %20
%269 = fadd double %267, %268
%270 = fsub double %269, 2.000000e+00
%271 = load double, double* %16
%272 = fmul double %270, %271
%273 = load double, double* %19
%274 = fsub double %273, 1.000000e+00
%275 = fdiv double %272, %274
%276 = load double, double* %26
%277 = fmul double %275, %276
store double %277, double* %26
%278 = load double, double* %26
%279 = load double, double* %22
%280 = fsub double %279, %278
store double %280, double* %22
%281 = load double, double* %27
%282 = load double, double* %22
%283 = load double, double* %25
%284 = call double @llvm.fmuladd.f64(double %281, double %282, double %283)
store double %284, double* %25
br label %242
285: ; preds = %254
%286 = load double, double* %25
%287 = load double*, double** %11
store double %286, double* %287
%288 = load double*, double** %11
%289 = load double, double* %288
%290 = fsub double 5.000000e-01, %289
%291 = fadd double 5.000000e-01, %290
%292 = load double*, double** %12
store double %291, double* %292
br label %293
293: ; preds = %285, %50, %42
ret void
}

LLVM IR (-O2):
define void @cdflib_cumfnc(double* %0, double* %1, double* %2, double* nocapture %3, double* %4, double* %5) {
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = bitcast double* %7 to i8*
%21 = bitcast double* %8 to i8*
%22 = bitcast double* %9 to i8*
%23 = bitcast double* %10 to i8*
%24 = bitcast double* %11 to i8*
%25 = bitcast double* %12 to i8*
%26 = bitcast i32* %13 to i8*
%27 = bitcast double* %14 to i8*
%28 = bitcast double* %15 to i8*
%29 = bitcast double* %16 to i8*
%30 = bitcast double* %17 to i8*
%31 = bitcast double* %18 to i8*
%32 = bitcast double* %19 to i8*
%33 = load double, double* %0
%34 = fcmp ugt double %33, 0.000000e+00
br i1 %34, label %36, label %35
35: ; preds = %6
store double 0.000000e+00, double* %4
store double 1.000000e+00, double* %5
br label %200
36: ; preds = %6
%37 = load double, double* %3
%38 = fcmp olt double %37, 1.000000e-10
br i1 %38, label %39, label %41
39: ; preds = %36
%40 = tail call i32 (double*, double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_cumf to i32 (double*, double*, double*, double*, double*, ...)*)(double* %0, double* %1, double* %2, double* %4, double* %5)
br label %200
41: ; preds = %36
%42 = fmul double %37, 5.000000e-01
%43 = fptosi double %42 to i32
%44 = icmp eq i32 %43, 0
%45 = select i1 %44, i32 1, i32 %43
%46 = add nsw i32 %45, 1
%47 = sitofp i32 %46 to double
store double %47, double* %14
%48 = fneg double %42
%49 = sitofp i32 %45 to double
%50 = tail call double @log(double %42)
%51 = tail call double @llvm.fmuladd.f64(double %49, double %50, double %48)
%52 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %14)
%53 = sitofp i32 %52 to double
%54 = fsub double %51, %53
%55 = call double @exp(double %54)
%56 = load double, double* %1
%57 = load double, double* %0
%58 = fmul double %56, %57
%59 = load double, double* %2
%60 = fadd double %59, %58
%61 = fdiv double %59, %60
store double %61, double* %9
%62 = fcmp ogt double %61, 5.000000e-01
br i1 %62, label %63, label %66
63: ; preds = %41
%64 = fdiv double %58, %60
store double %64, double* %8
%65 = fsub double 1.000000e+00, %64
store double %65, double* %9
br label %68
66: ; preds = %41
%67 = fsub double 1.000000e+00, %61
store double %67, double* %8
br label %68
68: ; preds = %66, %63
%69 = call double @llvm.fmuladd.f64(double %56, double 5.000000e-01, double %49)
store double %69, double* %15
%70 = fmul double %59, 5.000000e-01
store double %70, double* %16
%71 = call i32 (double*, double*, double*, double*, double*, double*, i32*, ...) bitcast (i32 (...)* @cdflib_beta_inc to i32 (double*, double*, double*, double*, double*, double*, i32*, ...)*)(double* %15, double* %16, double* %8, double* %9, double* %12, double* %7, i32* %13)
%72 = load double, double* %1
%73 = fmul double %72, 5.000000e-01
%74 = fadd double %73, %49
store double %74, double* %10
%75 = load double, double* %2
%76 = fmul double %75, 5.000000e-01
store double %76, double* %11
%77 = load double, double* %12
%78 = fmul double %55, %77
%79 = fadd double %74, %76
store double %79, double* %17
%80 = fadd double %74, 1.000000e+00
store double %80, double* %18
%81 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %17)
%82 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %18)
%83 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %11)
%84 = add i32 %82, %83
%85 = sub i32 %81, %84
%86 = sitofp i32 %85 to double
%87 = load double, double* %8
%88 = call double @log(double %87)
%89 = call double @llvm.fmuladd.f64(double %74, double %88, double %86)
%90 = load double, double* %11
%91 = load double, double* %9
%92 = call double @log(double %91)
%93 = call double @llvm.fmuladd.f64(double %90, double %92, double %89)
%94 = call double @exp(double %93)
%95 = fcmp olt double %78, 0x3BC79CA10C924223
br i1 %95, label %96, label %98
96: ; preds = %68
%97 = load double, double* %11
br label %128
98: ; preds = %68
%99 = load double, double* %12
%100 = load double, double* %11
%101 = load double, double* %8
br label %102
102: ; preds = %98, %114
%103 = phi i32 [ %45, %98 ], [ %118, %114 ]
%104 = phi double [ %74, %98 ], [ %119, %114 ]
%105 = phi double [ %55, %98 ], [ %117, %114 ]
%106 = phi double [ %78, %98 ], [ %126, %114 ]
%107 = phi double [ %94, %98 ], [ %124, %114 ]
%108 = phi double [ %99, %98 ], [ %125, %114 ]
%109 = fmul double %105, %108
%110 = fmul double %106, 1.000000e-04
%111 = fcmp olt double %109, %110
%112 = icmp slt i32 %103, 1
%113 = select i1 %111, i1 true, i1 %112
br i1 %113, label %128, label %114
114: ; preds = %102
%115 = sitofp i32 %103 to double
%116 = fdiv double %115, %42
%117 = fmul double %105, %116
%118 = add nsw i32 %103, -1
%119 = fadd double %104, -1.000000e+00
%120 = fadd double %119, 1.000000e+00
%121 = fadd double %119, %100
%122 = fmul double %121, %101
%123 = fdiv double %120, %122
%124 = fmul double %107, %123
%125 = fadd double %108, %124
store double %125, double* %12
%126 = call double @llvm.fmuladd.f64(double %117, double %125, double %106)
%127 = fcmp olt double %126, 0x3BC79CA10C924223
br i1 %127, label %128, label %102
128: ; preds = %114, %102, %96
%129 = phi double [ %97, %96 ], [ %100, %102 ], [ %100, %114 ]
%130 = phi double [ %78, %96 ], [ %126, %114 ], [ %106, %102 ]
%131 = load double, double* %10
%132 = fadd double %131, -1.000000e+00
%133 = fadd double %132, %129
%134 = fcmp oeq double %133, 0.000000e+00
br i1 %134, label %135, label %150
135: ; preds = %128
%136 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %10)
%137 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %11)
%138 = add i32 %136, %137
%139 = sub i32 0, %138
%140 = sitofp i32 %139 to double
%141 = load double, double* %10
%142 = fadd double %141, -1.000000e+00
%143 = load double, double* %8
%144 = call double @log(double %143)
%145 = call double @llvm.fmuladd.f64(double %142, double %144, double %140)
%146 = load double, double* %11
%147 = load double, double* %9
%148 = call double @log(double %147)
%149 = call double @llvm.fmuladd.f64(double %146, double %148, double %145)
br label %166
150: ; preds = %128
store double %133, double* %19
%151 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %19)
%152 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %10)
%153 = call i32 (double*, ...) bitcast (i32 (...)* @cdflib_gamma_log to i32 (double*, ...)*)(double* %11)
%154 = add i32 %152, %153
%155 = sub i32 %151, %154
%156 = sitofp i32 %155 to double
%157 = load double, double* %10
%158 = fadd double %157, -1.000000e+00
%159 = load double, double* %8
%160 = call double @log(double %159)
%161 = call double @llvm.fmuladd.f64(double %158, double %160, double %156)
%162 = load double, double* %11
%163 = load double, double* %9
%164 = call double @log(double %163)
%165 = call double @llvm.fmuladd.f64(double %162, double %164, double %161)
br label %166
166: ; preds = %150, %135
%167 = phi double [ %165, %150 ], [ %149, %135 ]
%168 = call double @exp(double %167)
%169 = load double, double* %11
%170 = load double, double* %8
%171 = load double, double* %10
br label %177
172: ; preds = %177
%173 = add nsw i32 %183, 1
%174 = fmul double %186, %194
%175 = fmul double %195, 1.000000e-04
%176 = fcmp olt double %174, %175
br i1 %176, label %197, label %177
177: ; preds = %172, %166
%178 = phi double [ %171, %166 ], [ %187, %172 ]
%179 = phi double [ %77, %166 ], [ %194, %172 ]
%180 = phi double [ %130, %166 ], [ %195, %172 ]
%181 = phi double [ %168, %166 ], [ %193, %172 ]
%182 = phi double [ %55, %166 ], [ %186, %172 ]
%183 = phi i32 [ %46, %166 ], [ %173, %172 ]
%184 = sitofp i32 %183 to double
%185 = fdiv double %42, %184
%186 = fmul double %182, %185
%187 = fadd double %178, 1.000000e+00
%188 = fadd double %187, %169
%189 = fadd double %188, -2.000000e+00
%190 = fmul double %170, %189
%191 = fadd double %187, -1.000000e+00
%192 = fdiv double %190, %191
%193 = fmul double %181, %192
%194 = fsub double %179, %193
%195 = call double @llvm.fmuladd.f64(double %186, double %194, double %180)
%196 = fcmp olt double %195, 0x3BC79CA10C924223
br i1 %196, label %197, label %172
197: ; preds = %177, %172
store double %187, double* %10
store double %195, double* %4
%198 = fsub double 5.000000e-01, %195
%199 = fadd double %198, 5.000000e-01
store double %199, double* %5
br label %200
200: ; preds = %197, %39, %35
ret void
}