void cdflib_beta_inc(const double *a, const double *b, const double *x,
                     const double *y, double *w, double *w1, int *ierr) {
  static const int K1 = 1;
  double a0, b0, eps, lambda, t, x0, y0, z;
  int ierr1, ind, n;
  double T2, T3, T4, T5;
  eps = dpmpar(K1);
  *w = *w1 = 0.0e0;
  if (*a < 0.0e0 || *b < 0.0e0) {
    goto S270;
  }
  if (*a == 0.0e0 && *b == 0.0e0) {
    goto S280;
  }
  if (*x < 0.0e0 || *x > 1.0e0) {
    goto S290;
  }
  if (*y < 0.0e0 || *y > 1.0e0) {
    goto S300;
  }
  z = *x + *y - 0.5e0 - 0.5e0;
  if (fabs(z) > 3.0e0 * eps) {
    goto S310;
  }
  *ierr = 0;
  if (*x == 0.0e0) {
    goto S210;
  }
  if (*y == 0.0e0) {
    goto S230;
  }
  if (*a == 0.0e0) {
    goto S240;
  }
  if (*b == 0.0e0) {
    goto S220;
  }
  eps = cdflib_fifdmax1(eps, 1.e-15);
  if (cdflib_fifdmax1(*a, *b) < 1.e-3 * eps) {
    goto S260;
  }
  ind = 0;
  a0 = *a;
  b0 = *b;
  x0 = *x;
  y0 = *y;
  if (cdflib_fifdmin1(a0, b0) > 1.0e0) {
    goto S40;
  }
  if (*x <= 0.5e0) {
    goto S10;
  }
  ind = 1;
  a0 = *b;
  b0 = *a;
  x0 = *y;
  y0 = *x;
S10:
  if (b0 < cdflib_fifdmin1(eps, eps * a0)) {
    goto S90;
  }
  if (a0 < cdflib_fifdmin1(eps, eps * b0) && b0 * x0 <= 1.0e0) {
    goto S100;
  }
  if (cdflib_fifdmax1(a0, b0) > 1.0e0) {
    goto S20;
  }
  if (a0 >= cdflib_fifdmin1(0.2e0, b0)) {
    goto S110;
  }
  if (pow(x0, a0) <= 0.9e0) {
    goto S110;
  }
  if (x0 >= 0.3e0) {
    goto S120;
  }
  n = 20;
  goto S140;
S20:
  if (b0 <= 1.0e0) {
    goto S110;
  }
  if (x0 >= 0.3e0) {
    goto S120;
  }
  if (x0 >= 0.1e0) {
    goto S30;
  }
  if (pow(x0 * b0, a0) <= 0.7e0) {
    goto S110;
  }
S30:
  if (b0 > 15.0e0) {
    goto S150;
  }
  n = 20;
  goto S140;
S40:
  if (*a > *b) {
    goto S50;
  }
  lambda = *a - (*a + *b) * *x;
  goto S60;
S50:
  lambda = (*a + *b) * *y - *b;
S60:
  if (lambda >= 0.0e0) {
    goto S70;
  }
  ind = 1;
  a0 = *b;
  b0 = *a;
  x0 = *y;
  y0 = *x;
  lambda = fabs(lambda);
S70:
  if (b0 < 40.0e0 && b0 * x0 <= 0.7e0) {
    goto S110;
  }
  if (b0 < 40.0e0) {
    goto S160;
  }
  if (a0 > b0) {
    goto S80;
  }
  if (a0 <= 100.0e0) {
    goto S130;
  }
  if (lambda > 0.03e0 * a0) {
    goto S130;
  }
  goto S200;
S80:
  if (b0 <= 100.0e0) {
    goto S130;
  }
  if (lambda > 0.03e0 * b0) {
    goto S130;
  }
  goto S200;
S90:
  *w = cdflib_fpser(&a0, &b0, &x0, &eps);
  *w1 = 0.5e0 + (0.5e0 - *w);
  goto S250;
S100:
  *w1 = cdflib_apser(&a0, &b0, &x0, &eps);
  *w = 0.5e0 + (0.5e0 - *w1);
  goto S250;
S110:
  *w = cdflib_beta_pser(&a0, &b0, &x0, &eps);
  *w1 = 0.5e0 + (0.5e0 - *w);
  goto S250;
S120:
  *w1 = cdflib_beta_pser(&b0, &a0, &y0, &eps);
  *w = 0.5e0 + (0.5e0 - *w1);
  goto S250;
S130:
  T2 = 15.0e0 * eps;
  *w = cdflib_beta_frac(&a0, &b0, &x0, &y0, &lambda, &T2);
  *w1 = 0.5e0 + (0.5e0 - *w);
  goto S250;
S140:
  *w1 = cdflib_beta_up(&b0, &a0, &y0, &x0, &n, &eps);
  b0 = b0 + (double)n;
S150:
  T3 = 15.0e0 * eps;
  cdflib_beta_grat(&b0, &a0, &y0, &x0, w1, &T3, &ierr1);
  *w = 0.5e0 + (0.5e0 - *w1);
  goto S250;
S160:
  n = (int)b0;
  b0 -= (double)n;
  if (b0 != 0.0e0) {
    goto S170;
  }
  n -= 1;
  b0 = 1.0e0;
S170:
  *w = cdflib_beta_up(&b0, &a0, &y0, &x0, &n, &eps);
  if (x0 > 0.7e0) {
    goto S180;
  }
  *w = *w + cdflib_beta_pser(&a0, &b0, &x0, &eps);
  *w1 = 0.5e0 + (0.5e0 - *w);
  goto S250;
S180:
  if (a0 > 15.0e0) {
    goto S190;
  }
  n = 20;
  *w = *w + cdflib_beta_up(&a0, &b0, &x0, &y0, &n, &eps);
  a0 = a0 + (double)n;
S190:
  T4 = 15.0e0 * eps;
  cdflib_beta_grat(&a0, &b0, &x0, &y0, w, &T4, &ierr1);
  *w1 = 0.5e0 + (0.5e0 - *w);
  goto S250;
S200:
  T5 = 100.0e0 * eps;
  *w = cdflib_beta_asym(&a0, &b0, &lambda, &T5);
  *w1 = 0.5e0 + (0.5e0 - *w);
  goto S250;
S210:
  if (*a == 0.0e0) {
    goto S320;
  }
S220:
  *w = 0.0e0;
  *w1 = 1.0e0;
  return;
S230:
  if (*b == 0.0e0) {
    goto S330;
  }
S240:
  *w = 1.0e0;
  *w1 = 0.0e0;
  return;
S250:
  if (ind == 0) {
    return;
  }
  t = *w;
  *w = *w1;
  *w1 = t;
  return;
S260:
  *w = *b / (*a + *b);
  *w1 = *a / (*a + *b);
  return;
S270:
  *ierr = 1;
  return;
S280:
  *ierr = 2;
  return;
S290:
  *ierr = 3;
  return;
S300:
  *ierr = 4;
  return;
S310:
  *ierr = 5;
  return;
S320:
  *ierr = 6;
  return;
S330:
  *ierr = 7;
  return;
}


x86 assembly (clang, -O0):
<cdflib_beta_inc>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%rsp
mov    0x10(%rbp),%rax
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    $0x1,%edi
mov    $0x0,%al
call   1050 <dpmpar@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
ja     1224 <cdflib_beta_inc+0x74>
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jbe    1229 <cdflib_beta_inc+0x79>
jmp    1e42 <cdflib_beta_inc+0xc92>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1264 <cdflib_beta_inc+0xb4>
jp     1264 <cdflib_beta_inc+0xb4>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1264 <cdflib_beta_inc+0xb4>
jp     1264 <cdflib_beta_inc+0xb4>
jmp    1e51 <cdflib_beta_inc+0xca1>
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
ja     128f <cdflib_beta_inc+0xdf>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xd7b(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1294 <cdflib_beta_inc+0xe4>
jmp    1e60 <cdflib_beta_inc+0xcb0>
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
ja     12bf <cdflib_beta_inc+0x10f>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xd4b(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12c4 <cdflib_beta_inc+0x114>
jmp    1e6f <cdflib_beta_inc+0xcbf>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
addsd  (%rax),%xmm0
movsd  0xd34(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xd28(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movaps 0xd73(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0xcff(%rip),%xmm1
mulsd  -0x48(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    131d <cdflib_beta_inc+0x16d>
jmp    1e7e <cdflib_beta_inc+0xcce>
mov    0x10(%rbp),%rax
movl   $0x0,(%rax)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1347 <cdflib_beta_inc+0x197>
jp     1347 <cdflib_beta_inc+0x197>
jmp    1d2d <cdflib_beta_inc+0xb7d>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1367 <cdflib_beta_inc+0x1b7>
jp     1367 <cdflib_beta_inc+0x1b7>
jmp    1d72 <cdflib_beta_inc+0xbc2>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1387 <cdflib_beta_inc+0x1d7>
jp     1387 <cdflib_beta_inc+0x1d7>
jmp    1d97 <cdflib_beta_inc+0xbe7>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    13a7 <cdflib_beta_inc+0x1f7>
jp     13a7 <cdflib_beta_inc+0x1f7>
jmp    1d52 <cdflib_beta_inc+0xba2>
movsd  -0x48(%rbp),%xmm0
movsd  0xc6c(%rip),%xmm1
mov    $0x2,%al
call   1080 <cdflib_fifdmax1@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm1
mov    $0x2,%al
call   1080 <cdflib_fifdmax1@plt>
cvtsi2sd %eax,%xmm1
movsd  0xc31(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    13fb <cdflib_beta_inc+0x24b>
jmp    1df5 <cdflib_beta_inc+0xc45>
movl   $0x0,-0x78(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
cvtsi2sd %eax,%xmm0
movsd  0xbad(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1462 <cdflib_beta_inc+0x2b2>
jmp    16b4 <cdflib_beta_inc+0x504>
mov    -0x18(%rbp),%rax
movsd  0xba2(%rip),%xmm0
ucomisd (%rax),%xmm0
jb     147d <cdflib_beta_inc+0x2cd>
jmp    14b8 <cdflib_beta_inc+0x308>
movl   $0x1,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
movsd  -0xa8(%rbp),%xmm1
ff
cvtsi2sd %eax,%xmm0
ucomisd %xmm1,%xmm0
jbe    14f6 <cdflib_beta_inc+0x346>
jmp    18a0 <cdflib_beta_inc+0x6f0>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
movsd  -0xb0(%rbp),%xmm1
ff
cvtsi2sd %eax,%xmm0
ucomisd %xmm1,%xmm0
jbe    1550 <cdflib_beta_inc+0x3a0>
movsd  -0x40(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
movsd  0xabf(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     1550 <cdflib_beta_inc+0x3a0>
jmp    18ec <cdflib_beta_inc+0x73c>
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
mov    $0x2,%al
call   1080 <cdflib_fifdmax1@plt>
cvtsi2sd %eax,%xmm0
movsd  0xa93(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    157c <cdflib_beta_inc+0x3cc>
jmp    1609 <cdflib_beta_inc+0x459>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x40(%rbp),%xmm1
movsd  0xa92(%rip),%xmm0
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
movsd  -0xb8(%rbp),%xmm0
ff
cvtsi2sd %eax,%xmm1
ucomisd %xmm1,%xmm0
jb     15b8 <cdflib_beta_inc+0x408>
jmp    1938 <cdflib_beta_inc+0x788>
movsd  -0x60(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
call   10a0 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xa5e(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     15e1 <cdflib_beta_inc+0x431>
jmp    1938 <cdflib_beta_inc+0x788>
movsd  -0x60(%rbp),%xmm0
movsd  0xa4a(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     15fd <cdflib_beta_inc+0x44d>
jmp    1984 <cdflib_beta_inc+0x7d4>
movl   $0x14,-0x7c(%rbp)
jmp    1a3c <cdflib_beta_inc+0x88c>
movsd  0x9ef(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jb     1621 <cdflib_beta_inc+0x471>
jmp    1938 <cdflib_beta_inc+0x788>
movsd  -0x60(%rbp),%xmm0
movsd  0xa0a(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     163d <cdflib_beta_inc+0x48d>
jmp    1984 <cdflib_beta_inc+0x7d4>
movsd  -0x60(%rbp),%xmm0
movsd  0x9f6(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1659 <cdflib_beta_inc+0x4a9>
jmp    168c <cdflib_beta_inc+0x4dc>
movsd  -0x60(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
call   10a0 <pow@plt>
movaps %xmm0,%xmm1
movsd  0x9d0(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     1687 <cdflib_beta_inc+0x4d7>
jmp    1938 <cdflib_beta_inc+0x788>
jmp    168c <cdflib_beta_inc+0x4dc>
movsd  -0x40(%rbp),%xmm0
movsd  0x9b7(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    16a8 <cdflib_beta_inc+0x4f8>
jmp    1a7a <cdflib_beta_inc+0x8ca>
movl   $0x14,-0x7c(%rbp)
jmp    1a3c <cdflib_beta_inc+0x88c>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
ucomisd (%rax),%xmm0
jbe    16cf <cdflib_beta_inc+0x51f>
jmp    1718 <cdflib_beta_inc+0x568>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
addsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    175c <cdflib_beta_inc+0x5ac>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
addsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1773 <cdflib_beta_inc+0x5c3>
jmp    17c3 <cdflib_beta_inc+0x613>
movl   $0x1,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm0
movaps 0x8b6(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0x88d(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    17f7 <cdflib_beta_inc+0x647>
movsd  -0x40(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
movsd  0x860(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     17f7 <cdflib_beta_inc+0x647>
jmp    1938 <cdflib_beta_inc+0x788>
movsd  0x859(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    180f <cdflib_beta_inc+0x65f>
jmp    1ae2 <cdflib_beta_inc+0x932>
movsd  -0x38(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1824 <cdflib_beta_inc+0x674>
jmp    1862 <cdflib_beta_inc+0x6b2>
movsd  0x834(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jb     183c <cdflib_beta_inc+0x68c>
jmp    19d0 <cdflib_beta_inc+0x820>
movsd  -0x50(%rbp),%xmm0
movsd  0x81f(%rip),%xmm1
mulsd  -0x38(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    185d <cdflib_beta_inc+0x6ad>
jmp    19d0 <cdflib_beta_inc+0x820>
jmp    1cc9 <cdflib_beta_inc+0xb19>
movsd  0x7f6(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jb     187a <cdflib_beta_inc+0x6ca>
jmp    19d0 <cdflib_beta_inc+0x820>
movsd  -0x50(%rbp),%xmm0
movsd  0x7e1(%rip),%xmm1
mulsd  -0x40(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    189b <cdflib_beta_inc+0x6eb>
jmp    19d0 <cdflib_beta_inc+0x820>
jmp    1cc9 <cdflib_beta_inc+0xb19>
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x60(%rbp),%rdx
lea    -0x48(%rbp),%rcx
mov    $0x0,%al
call   1090 <cdflib_fpser@plt>
cvtsi2sd %eax,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  0x741(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x735(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x60(%rbp),%rdx
lea    -0x48(%rbp),%rcx
mov    $0x0,%al
call   10c0 <cdflib_apser@plt>
cvtsi2sd %eax,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
movsd  0x6f5(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x6e9(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x60(%rbp),%rdx
lea    -0x48(%rbp),%rcx
mov    $0x0,%al
call   10d0 <cdflib_beta_pser@plt>
cvtsi2sd %eax,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  0x6a9(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x69d(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
lea    -0x40(%rbp),%rdi
lea    -0x38(%rbp),%rsi
lea    -0x68(%rbp),%rdx
lea    -0x48(%rbp),%rcx
mov    $0x0,%al
call   10d0 <cdflib_beta_pser@plt>
cvtsi2sd %eax,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
movsd  0x65d(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x651(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
movsd  0x678(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x60(%rbp),%rdx
lea    -0x68(%rbp),%rcx
lea    -0x50(%rbp),%r8
lea    -0x88(%rbp),%r9
mov    $0x0,%al
call   1060 <cdflib_beta_frac@plt>
cvtsi2sd %eax,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  0x5f1(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x5e5(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
lea    -0x40(%rbp),%rdi
lea    -0x38(%rbp),%rsi
lea    -0x68(%rbp),%rdx
lea    -0x60(%rbp),%rcx
lea    -0x7c(%rbp),%r8
lea    -0x48(%rbp),%r9
mov    $0x0,%al
call   1040 <cdflib_beta_up@plt>
cvtsi2sd %eax,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x40(%rbp),%xmm0
cvtsi2sdl -0x7c(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0x5ce(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x30(%rbp),%r8
lea    -0x40(%rbp),%rdi
lea    -0x38(%rbp),%rsi
lea    -0x68(%rbp),%rdx
lea    -0x60(%rbp),%rcx
lea    -0x90(%rbp),%r9
lea    -0x74(%rbp),%rax
mov    %rax,(%rsp)
mov    $0x0,%al
call   1070 <cdflib_beta_grat@plt>
mov    -0x30(%rbp),%rax
movsd  0x54b(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x53f(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
cvttsd2si -0x40(%rbp),%eax
mov    %eax,-0x7c(%rbp)
cvtsi2sdl -0x7c(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1b1a <cdflib_beta_inc+0x96a>
jp     1b1a <cdflib_beta_inc+0x96a>
jmp    1b1f <cdflib_beta_inc+0x96f>
jmp    1b35 <cdflib_beta_inc+0x985>
mov    -0x7c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x7c(%rbp)
movsd  0x4d0(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
lea    -0x40(%rbp),%rdi
lea    -0x38(%rbp),%rsi
lea    -0x68(%rbp),%rdx
lea    -0x60(%rbp),%rcx
lea    -0x7c(%rbp),%r8
lea    -0x48(%rbp),%r9
mov    $0x0,%al
call   1040 <cdflib_beta_up@plt>
cvtsi2sd %eax,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x60(%rbp),%xmm0
movsd  0x4db(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1b7c <cdflib_beta_inc+0x9cc>
jmp    1be4 <cdflib_beta_inc+0xa34>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x60(%rbp),%rdx
lea    -0x48(%rbp),%rcx
mov    $0x0,%al
call   10d0 <cdflib_beta_pser@plt>
movsd  -0xc0(%rbp),%xmm0
ff
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  0x449(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x43d(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
movsd  -0x38(%rbp),%xmm0
movsd  0x45f(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1c00 <cdflib_beta_inc+0xa50>
jmp    1c61 <cdflib_beta_inc+0xab1>
movl   $0x14,-0x7c(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x60(%rbp),%rdx
lea    -0x68(%rbp),%rcx
lea    -0x7c(%rbp),%r8
lea    -0x48(%rbp),%r9
mov    $0x0,%al
call   1040 <cdflib_beta_up@plt>
movsd  -0xc8(%rbp),%xmm0
ff
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x38(%rbp),%xmm0
cvtsi2sdl -0x7c(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0x3e7(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x28(%rbp),%r8
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x60(%rbp),%rdx
lea    -0x68(%rbp),%rcx
lea    -0x98(%rbp),%r9
lea    -0x74(%rbp),%rax
mov    %rax,(%rsp)
mov    $0x0,%al
call   1070 <cdflib_beta_grat@plt>
mov    -0x28(%rbp),%rax
movsd  0x364(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x358(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
movsd  0x38f(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
lea    -0x38(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x50(%rbp),%rdx
lea    -0xa0(%rbp),%rcx
mov    $0x0,%al
call   10b0 <cdflib_beta_asym@plt>
cvtsi2sd %eax,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  0x300(%rip),%xmm1
subsd  (%rax),%xmm1
movsd  0x2f4(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1db7 <cdflib_beta_inc+0xc07>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1d4d <cdflib_beta_inc+0xb9d>
jp     1d4d <cdflib_beta_inc+0xb9d>
jmp    1e8d <cdflib_beta_inc+0xcdd>
jmp    1d52 <cdflib_beta_inc+0xba2>
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
movsd  0x297(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1d92 <cdflib_beta_inc+0xbe2>
jp     1d92 <cdflib_beta_inc+0xbe2>
jmp    1e9c <cdflib_beta_inc+0xcec>
jmp    1d97 <cdflib_beta_inc+0xbe7>
mov    -0x28(%rbp),%rax
movsd  0x25d(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
cmpl   $0x0,-0x78(%rbp)
jne    1dc6 <cdflib_beta_inc+0xc16>
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x58(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x10(%rbp),%rax
addsd  (%rax),%xmm1
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x10(%rbp),%rax
addsd  (%rax),%xmm1
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    0x10(%rbp),%rax
movl   $0x1,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    0x10(%rbp),%rax
movl   $0x2,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    0x10(%rbp),%rax
movl   $0x3,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    0x10(%rbp),%rax
movl   $0x4,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    0x10(%rbp),%rax
movl   $0x5,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    0x10(%rbp),%rax
movl   $0x6,(%rax)
jmp    1ea6 <cdflib_beta_inc+0xcf6>
mov    0x10(%rbp),%rax
movl   $0x7,(%rax)
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cdflib_beta_inc>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9,%r15
mov    %r8,%r12
mov    %rcx,0x28(%rsp)
mov    %rdx,%r13
mov    %rsi,%rbx
mov    %rdi,%rbp
mov    0xb0(%rsp),%r14
mov    $0x1,%edi
xor    %eax,%eax
call   1050 <dpmpar@plt>
movq   $0x0,(%r15)
movq   $0x0,(%r12)
movsd  0x0(%rbp),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
ja     120e <cdflib_beta_inc+0x5e>
movsd  (%rbx),%xmm2
ucomisd %xmm2,%xmm0
jbe    1224 <cdflib_beta_inc+0x74>
movl   $0x1,(%r14)
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    1241 <cdflib_beta_inc+0x91>
jp     1241 <cdflib_beta_inc+0x91>
ucomisd %xmm0,%xmm2
jne    1241 <cdflib_beta_inc+0x91>
jp     1241 <cdflib_beta_inc+0x91>
movl   $0x2,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
movsd  0x0(%r13),%xmm3
ucomisd %xmm3,%xmm0
ja     12b9 <cdflib_beta_inc+0x109>
ucomisd 0xdab(%rip),%xmm3
ja     12b9 <cdflib_beta_inc+0x109>
mov    0x28(%rsp),%rcx
movsd  (%rcx),%xmm4
xorpd  %xmm0,%xmm0
ucomisd %xmm4,%xmm0
ja     12c5 <cdflib_beta_inc+0x115>
ucomisd 0xd8e(%rip),%xmm4
ja     12c5 <cdflib_beta_inc+0x115>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movapd %xmm3,%xmm5
addsd  %xmm4,%xmm5
movsd  0xd7d(%rip),%xmm6
addsd  %xmm6,%xmm5
addsd  %xmm6,%xmm5
andpd  0xde5(%rip),%xmm5
movsd  0xd6d(%rip),%xmm6
mulsd  %xmm0,%xmm6
ucomisd %xmm6,%xmm5
jbe    12d1 <cdflib_beta_inc+0x121>
movl   $0x5,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
movl   $0x3,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
movl   $0x4,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
movl   $0x0,(%r14)
xorpd  %xmm5,%xmm5
ucomisd %xmm5,%xmm3
jne    12f8 <cdflib_beta_inc+0x148>
jp     12f8 <cdflib_beta_inc+0x148>
ucomisd %xmm5,%xmm1
jne    1346 <cdflib_beta_inc+0x196>
jp     1346 <cdflib_beta_inc+0x196>
movl   $0x6,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
ucomisd %xmm5,%xmm4
jne    1318 <cdflib_beta_inc+0x168>
jp     1318 <cdflib_beta_inc+0x168>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm2
jne    1324 <cdflib_beta_inc+0x174>
jp     1324 <cdflib_beta_inc+0x174>
movl   $0x7,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm1
jne    133e <cdflib_beta_inc+0x18e>
jp     133e <cdflib_beta_inc+0x18e>
movabs $0x3ff0000000000000,%rax
mov    %rax,(%r12)
movq   $0x0,(%r15)
jmp    1215 <cdflib_beta_inc+0x65>
ucomisd %xmm3,%xmm2
jne    1360 <cdflib_beta_inc+0x1b0>
jp     1360 <cdflib_beta_inc+0x1b0>
movq   $0x0,(%r12)
movabs $0x3ff0000000000000,%rax
mov    %rax,(%r15)
jmp    1215 <cdflib_beta_inc+0x65>
mov    %r15,%r14
movsd  0xcad(%rip),%xmm1
mov    $0x2,%al
call   1080 <cdflib_fifdmax1@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x58(%rsp)
movsd  %xmm0,0x30(%rsp)
movsd  0x0(%rbp),%xmm0
movsd  (%rbx),%xmm1
mov    $0x2,%al
call   1080 <cdflib_fifdmax1@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0xc7c(%rip),%xmm1
mulsd  0x58(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jbe    13e2 <cdflib_beta_inc+0x232>
movsd  (%rbx),%xmm0
movsd  0x0(%rbp),%xmm1
addsd  %xmm0,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r12)
movsd  0x0(%rbp),%xmm0
movsd  (%rbx),%xmm1
addsd  %xmm0,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
movsd  0x0(%rbp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  (%rbx),%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  0x0(%r13),%xmm2
movsd  %xmm2,0x48(%rsp)
movsd  %xmm2,0x20(%rsp)
mov    0x28(%rsp),%rax
movsd  (%rax),%xmm2
movsd  %xmm2,0x40(%rsp)
movsd  %xmm0,0x50(%rsp)
movsd  %xmm1,0x38(%rsp)
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
cmp    $0x1,%eax
jle    1454 <cdflib_beta_inc+0x2a4>
movsd  0x0(%rbp),%xmm0
movsd  (%rbx),%xmm1
ucomisd %xmm1,%xmm0
movapd 0xc4b(%rip),%xmm4
ja     14f7 <cdflib_beta_inc+0x347>
movapd %xmm0,%xmm3
jmp    14ff <cdflib_beta_inc+0x34f>
movsd  0x0(%r13),%xmm0
xor    %r15d,%r15d
movsd  0xbc3(%rip),%xmm1
ucomisd %xmm0,%xmm1
jae    14ad <cdflib_beta_inc+0x2fd>
movsd  (%rbx),%xmm1
movsd  %xmm1,0x50(%rsp)
movsd  %xmm1,0x10(%rsp)
movsd  0x0(%rbp),%xmm1
movsd  %xmm1,0x38(%rsp)
movsd  %xmm1,0x8(%rsp)
mov    0x28(%rsp),%rax
movsd  (%rax),%xmm1
movsd  %xmm1,0x48(%rsp)
movsd  %xmm1,0x20(%rsp)
movsd  %xmm0,0x40(%rsp)
mov    $0x1,%r15d
movsd  0x50(%rsp),%xmm1
movsd  0x58(%rsp),%xmm0
mulsd  %xmm0,%xmm1
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
ucomisd 0x38(%rsp),%xmm0
jbe    157b <cdflib_beta_inc+0x3cb>
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x20(%rsp),%rdx
lea    0x30(%rsp),%rcx
xor    %eax,%eax
call   1090 <cdflib_fpser@plt>
jmp    1a17 <cdflib_beta_inc+0x867>
movapd %xmm1,%xmm3
xorpd  %xmm4,%xmm3
movapd %xmm0,%xmm2
addsd  %xmm1,%xmm2
mov    %r13,%rax
cmova  0x28(%rsp),%rax
ja     1516 <cdflib_beta_inc+0x366>
xorpd  %xmm4,%xmm2
mulsd  (%rax),%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,0x68(%rsp)
xor    %r15d,%r15d
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm2
jae    15fa <cdflib_beta_inc+0x44a>
movsd  %xmm1,0x10(%rsp)
movsd  %xmm0,0x8(%rsp)
mov    0x28(%rsp),%rax
movsd  (%rax),%xmm3
movsd  %xmm3,0x48(%rsp)
movsd  %xmm3,0x20(%rsp)
movsd  0x0(%r13),%xmm3
movsd  %xmm3,0x40(%rsp)
andpd  0xb16(%rip),%xmm2
movlpd %xmm2,0x68(%rsp)
mov    $0x1,%r15d
jmp    1606 <cdflib_beta_inc+0x456>
movsd  0x38(%rsp),%xmm1
movsd  0x58(%rsp),%xmm0
mulsd  %xmm0,%xmm1
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
ucomisd 0x50(%rsp),%xmm0
movsd  0x38(%rsp),%xmm0
mulsd  0x48(%rsp),%xmm0
movsd  %xmm0,0x28(%rsp)
jbe    17a3 <cdflib_beta_inc+0x5f3>
movsd  0xa41(%rip),%xmm0
ucomisd 0x28(%rsp),%xmm0
jb     17a3 <cdflib_beta_inc+0x5f3>
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x20(%rsp),%rdx
lea    0x30(%rsp),%rcx
xor    %eax,%eax
call   10c0 <cdflib_apser@plt>
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  %xmm2,(%r14)
mov    %r12,%rax
jmp    1a27 <cdflib_beta_inc+0x877>
movsd  0x50(%rsp),%xmm1
movsd  0x38(%rsp),%xmm0
movsd  0xa52(%rip),%xmm3
ucomisd %xmm0,%xmm3
jbe    1770 <cdflib_beta_inc+0x5c0>
movsd  0x48(%rsp),%xmm2
mulsd  %xmm0,%xmm2
movsd  0xa26(%rip),%xmm1
ucomisd %xmm2,%xmm1
jae    19fc <cdflib_beta_inc+0x84c>
cvttsd2si %xmm0,%eax
cvttpd2dq %xmm0,%xmm1
mov    %eax,0x1c(%rsp)
cvtdq2pd %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rsp)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1670 <cdflib_beta_inc+0x4c0>
jp     1670 <cdflib_beta_inc+0x4c0>
add    $0xffffffff,%eax
mov    %eax,0x1c(%rsp)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x8(%rsp)
lea    0x8(%rsp),%rdi
lea    0x10(%rsp),%rsi
lea    0x40(%rsp),%rdx
lea    0x20(%rsp),%rcx
lea    0x1c(%rsp),%r8
lea    0x30(%rsp),%r9
xor    %eax,%eax
call   1040 <cdflib_beta_up@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  %xmm0,(%r12)
movsd  0x20(%rsp),%xmm0
ucomisd 0x99a(%rip),%xmm0
jbe    18c1 <cdflib_beta_inc+0x711>
movsd  0x10(%rsp),%xmm0
ucomisd 0x98e(%rip),%xmm0
ja     1725 <cdflib_beta_inc+0x575>
movl   $0x14,0x1c(%rsp)
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x20(%rsp),%rdx
lea    0x40(%rsp),%rcx
lea    0x1c(%rsp),%r8
lea    0x30(%rsp),%r9
xor    %eax,%eax
call   1040 <cdflib_beta_up@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x28(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r12)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x1c(%rsp),%xmm0
addsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x30(%rsp),%xmm0
mulsd  0x925(%rip),%xmm0
movsd  %xmm0,0x60(%rsp)
lea    0x74(%rsp),%rax
mov    %rax,(%rsp)
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x20(%rsp),%rdx
lea    0x40(%rsp),%rcx
lea    0x60(%rsp),%r9
mov    %r12,%r8
xor    %eax,%eax
call   1070 <cdflib_beta_grat@plt>
movsd  (%r12),%xmm2
jmp    1a24 <cdflib_beta_inc+0x874>
ucomisd %xmm0,%xmm1
jbe    1833 <cdflib_beta_inc+0x683>
movsd  0x8e6(%rip),%xmm1
ucomisd %xmm0,%xmm1
jae    1883 <cdflib_beta_inc+0x6d3>
mulsd  0x8dc(%rip),%xmm0
ucomisd %xmm0,%xmm2
jbe    184f <cdflib_beta_inc+0x69f>
jmp    1883 <cdflib_beta_inc+0x6d3>
movsd  0x50(%rsp),%xmm0
movsd  0x38(%rsp),%xmm1
mov    $0x2,%al
call   1080 <cdflib_fifdmax1@plt>
cmp    $0x1,%eax
jle    18f2 <cdflib_beta_inc+0x742>
movsd  0x839(%rip),%xmm0
ucomisd 0x38(%rsp),%xmm0
jae    19fc <cdflib_beta_inc+0x84c>
movsd  0x48(%rsp),%xmm0
ucomisd 0x85f(%rip),%xmm0
jae    1a63 <cdflib_beta_inc+0x8b3>
movsd  0x48(%rsp),%xmm0
ucomisd 0x853(%rip),%xmm0
jae    181a <cdflib_beta_inc+0x66a>
movsd  0x28(%rsp),%xmm0
movsd  0x50(%rsp),%xmm1
call   10a0 <pow@plt>
movsd  0x840(%rip),%xmm1
ucomisd %xmm0,%xmm1
jae    19fc <cdflib_beta_inc+0x84c>
movsd  0x38(%rsp),%xmm0
ucomisd 0x830(%rip),%xmm0
jbe    1955 <cdflib_beta_inc+0x7a5>
jmp    19af <cdflib_beta_inc+0x7ff>
movsd  0x82d(%rip),%xmm0
ucomisd %xmm1,%xmm0
jae    1883 <cdflib_beta_inc+0x6d3>
mulsd  0x827(%rip),%xmm1
ucomisd %xmm1,%xmm2
ja     1883 <cdflib_beta_inc+0x6d3>
movsd  0x58(%rsp),%xmm0
mulsd  0x80b(%rip),%xmm0
movsd  %xmm0,0x60(%rsp)
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x68(%rsp),%rdx
lea    0x60(%rsp),%rcx
xor    %eax,%eax
call   10b0 <cdflib_beta_asym@plt>
jmp    1a17 <cdflib_beta_inc+0x867>
movsd  0x58(%rsp),%xmm0
mulsd  0x7c7(%rip),%xmm0
movsd  %xmm0,0x60(%rsp)
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x20(%rsp),%rdx
lea    0x40(%rsp),%rcx
lea    0x68(%rsp),%r8
lea    0x60(%rsp),%r9
xor    %eax,%eax
call   1060 <cdflib_beta_frac@plt>
jmp    1a17 <cdflib_beta_inc+0x867>
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x20(%rsp),%rdx
lea    0x30(%rsp),%rcx
xor    %eax,%eax
call   10d0 <cdflib_beta_pser@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x28(%rsp),%xmm2
addsd  %xmm0,%xmm2
jmp    1a1e <cdflib_beta_inc+0x86e>
movsd  0x736(%rip),%xmm0
movsd  0x38(%rsp),%xmm1
mov    $0x2,%al
call   1030 <cdflib_fifdmin1@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x50(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jae    19fc <cdflib_beta_inc+0x84c>
movsd  0x48(%rsp),%xmm0
movsd  0x50(%rsp),%xmm1
call   10a0 <pow@plt>
movsd  0x701(%rip),%xmm1
ucomisd %xmm0,%xmm1
jae    19fc <cdflib_beta_inc+0x84c>
movsd  0x48(%rsp),%xmm0
ucomisd 0x6f1(%rip),%xmm0
jae    1a63 <cdflib_beta_inc+0x8b3>
movl   $0x14,0x1c(%rsp)
lea    0x8(%rsp),%rdi
lea    0x10(%rsp),%rsi
lea    0x40(%rsp),%rdx
lea    0x20(%rsp),%rcx
lea    0x1c(%rsp),%r8
lea    0x30(%rsp),%r9
xor    %eax,%eax
call   1040 <cdflib_beta_up@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r14)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x1c(%rsp),%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x58(%rsp)
movsd  0x58(%rsp),%xmm0
mulsd  0x69b(%rip),%xmm0
movsd  %xmm0,0x68(%rsp)
lea    0x74(%rsp),%rax
mov    %rax,(%rsp)
lea    0x8(%rsp),%rdi
lea    0x10(%rsp),%rsi
lea    0x40(%rsp),%rdx
lea    0x20(%rsp),%rcx
lea    0x68(%rsp),%r9
mov    %r14,%rbx
mov    %r14,%r8
xor    %eax,%eax
call   1070 <cdflib_beta_grat@plt>
movsd  (%r14),%xmm2
jmp    15f2 <cdflib_beta_inc+0x442>
lea    0x10(%rsp),%rdi
lea    0x8(%rsp),%rsi
lea    0x20(%rsp),%rdx
lea    0x30(%rsp),%rcx
xor    %eax,%eax
call   10d0 <cdflib_beta_pser@plt>
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  %xmm2,(%r12)
mov    %r14,%rax
movsd  0x5f9(%rip),%xmm0
movapd %xmm0,%xmm1
subsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rax)
test   %r15d,%r15d
je     1215 <cdflib_beta_inc+0x65>
movsd  (%r12),%xmm0
movsd  (%r14),%xmm1
movsd  %xmm1,(%r12)
movsd  %xmm0,(%r14)
jmp    1215 <cdflib_beta_inc+0x65>
lea    0x8(%rsp),%rdi
lea    0x10(%rsp),%rsi
lea    0x40(%rsp),%rdx
lea    0x30(%rsp),%rcx
xor    %eax,%eax
call   10d0 <cdflib_beta_pser@plt>
jmp    15e6 <cdflib_beta_inc+0x436>

LLVM IR (-O0):
@cdflib_beta_inc.K1 = internal constant i32 1
define void @cdflib_beta_inc(double* %0, double* %1, double* %2, double* %3, double* %4, double* %5, i32* %6) {
%8 = alloca double*
%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
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
store double* %0, double** %8
store double* %1, double** %9
store double* %2, double** %10
store double* %3, double** %11
store double* %4, double** %12
store double* %5, double** %13
store i32* %6, i32** %14
%30 = call i32 (i32, ...) bitcast (i32 (...)* @dpmpar to i32 (i32, ...)*)(i32 1)
%31 = sitofp i32 %30 to double
store double %31, double* %17
%32 = load double*, double** %13
store double 0.000000e+00, double* %32
%33 = load double*, double** %12
store double 0.000000e+00, double* %33
%34 = load double*, double** %8
%35 = load double, double* %34
%36 = fcmp olt double %35, 0.000000e+00
br i1 %36, label %41, label %37
37: ; preds = %7
%38 = load double*, double** %9
%39 = load double, double* %38
%40 = fcmp olt double %39, 0.000000e+00
br i1 %40, label %41, label %42
41: ; preds = %37, %7
br label %491
42: ; preds = %37
%43 = load double*, double** %8
%44 = load double, double* %43
%45 = fcmp oeq double %44, 0.000000e+00
br i1 %45, label %46, label %51
46: ; preds = %42
%47 = load double*, double** %9
%48 = load double, double* %47
%49 = fcmp oeq double %48, 0.000000e+00
br i1 %49, label %50, label %51
50: ; preds = %46
br label %493
51: ; preds = %46, %42
%52 = load double*, double** %10
%53 = load double, double* %52
%54 = fcmp olt double %53, 0.000000e+00
br i1 %54, label %59, label %55
55: ; preds = %51
%56 = load double*, double** %10
%57 = load double, double* %56
%58 = fcmp ogt double %57, 1.000000e+00
br i1 %58, label %59, label %60
59: ; preds = %55, %51
br label %495
60: ; preds = %55
%61 = load double*, double** %11
%62 = load double, double* %61
%63 = fcmp olt double %62, 0.000000e+00
br i1 %63, label %68, label %64
64: ; preds = %60
%65 = load double*, double** %11
%66 = load double, double* %65
%67 = fcmp ogt double %66, 1.000000e+00
br i1 %67, label %68, label %69
68: ; preds = %64, %60
br label %497
69: ; preds = %64
%70 = load double*, double** %10
%71 = load double, double* %70
%72 = load double*, double** %11
%73 = load double, double* %72
%74 = fadd double %71, %73
%75 = fsub double %74, 5.000000e-01
%76 = fsub double %75, 5.000000e-01
store double %76, double* %22
%77 = load double, double* %22
%78 = call double @llvm.fabs.f64(double %77)
%79 = load double, double* %17
%80 = fmul double 3.000000e+00, %79
%81 = fcmp ogt double %78, %80
br i1 %81, label %82, label %83
82: ; preds = %69
br label %499
83: ; preds = %69
%84 = load i32*, i32** %14
store i32 0, i32* %84
%85 = load double*, double** %10
%86 = load double, double* %85
%87 = fcmp oeq double %86, 0.000000e+00
br i1 %87, label %88, label %89
88: ; preds = %83
br label %442
89: ; preds = %83
%90 = load double*, double** %11
%91 = load double, double* %90
%92 = fcmp oeq double %91, 0.000000e+00
br i1 %92, label %93, label %94
93: ; preds = %89
br label %451
94: ; preds = %89
%95 = load double*, double** %8
%96 = load double, double* %95
%97 = fcmp oeq double %96, 0.000000e+00
br i1 %97, label %98, label %99
98: ; preds = %94
br label %457
99: ; preds = %94
%100 = load double*, double** %9
%101 = load double, double* %100
%102 = fcmp oeq double %101, 0.000000e+00
br i1 %102, label %103, label %104
103: ; preds = %99
br label %448
104: ; preds = %99
%105 = load double, double* %17
%106 = call i32 (double, double, ...) bitcast (i32 (...)* @cdflib_fifdmax1 to i32 (double, double, ...)*)(double %105, double 1.000000e-15)
%107 = sitofp i32 %106 to double
store double %107, double* %17
%108 = load double*, double** %8
%109 = load double, double* %108
%110 = load double*, double** %9
%111 = load double, double* %110
%112 = call i32 (double, double, ...) bitcast (i32 (...)* @cdflib_fifdmax1 to i32 (double, double, ...)*)(double %109, double %111)
%113 = sitofp i32 %112 to double
%114 = load double, double* %17
%115 = fmul double 1.000000e-03, %114
%116 = fcmp olt double %113, %115
br i1 %116, label %117, label %118
117: ; preds = %104
br label %472
118: ; preds = %104
store i32 0, i32* %24
%119 = load double*, double** %8
%120 = load double, double* %119
store double %120, double* %15
%121 = load double*, double** %9
%122 = load double, double* %121
store double %122, double* %16
%123 = load double*, double** %10
%124 = load double, double* %123
store double %124, double* %20
%125 = load double*, double** %11
%126 = load double, double* %125
store double %126, double* %21
%127 = load double, double* %15
%128 = load double, double* %16
%129 = call i32 (double, double, ...) bitcast (i32 (...)* @cdflib_fifdmin1 to i32 (double, double, ...)*)(double %127, double %128)
%130 = sitofp i32 %129 to double
%131 = fcmp ogt double %130, 1.000000e+00
br i1 %131, label %132, label %133
132: ; preds = %118
br label %223
133: ; preds = %118
%134 = load double*, double** %10
%135 = load double, double* %134
%136 = fcmp ole double %135, 5.000000e-01
br i1 %136, label %137, label %138
137: ; preds = %133
br label %147
138: ; preds = %133
store i32 1, i32* %24
%139 = load double*, double** %9
%140 = load double, double* %139
store double %140, double* %15
%141 = load double*, double** %8
%142 = load double, double* %141
store double %142, double* %16
%143 = load double*, double** %11
%144 = load double, double* %143
store double %144, double* %20
%145 = load double*, double** %10
%146 = load double, double* %145
store double %146, double* %21
br label %147
147: ; preds = %138, %137
%148 = load double, double* %16
%149 = load double, double* %17
%150 = load double, double* %17
%151 = load double, double* %15
%152 = fmul double %150, %151
%153 = call i32 (double, double, ...) bitcast (i32 (...)* @cdflib_fifdmin1 to i32 (double, double, ...)*)(double %149, double %152)
%154 = sitofp i32 %153 to double
%155 = fcmp olt double %148, %154
br i1 %155, label %156, label %157
156: ; preds = %147
br label %309
157: ; preds = %147
%158 = load double, double* %15
%159 = load double, double* %17
%160 = load double, double* %17
%161 = load double, double* %16
%162 = fmul double %160, %161
%163 = call i32 (double, double, ...) bitcast (i32 (...)* @cdflib_fifdmin1 to i32 (double, double, ...)*)(double %159, double %162)
%164 = sitofp i32 %163 to double
%165 = fcmp olt double %158, %164
br i1 %165, label %166, label %172
166: ; preds = %157
%167 = load double, double* %16
%168 = load double, double* %20
%169 = fmul double %167, %168
%170 = fcmp ole double %169, 1.000000e+00
br i1 %170, label %171, label %172
171: ; preds = %166
br label %318
172: ; preds = %166, %157
%173 = load double, double* %15
%174 = load double, double* %16
%175 = call i32 (double, double, ...) bitcast (i32 (...)* @cdflib_fifdmax1 to i32 (double, double, ...)*)(double %173, double %174)
%176 = sitofp i32 %175 to double
%177 = fcmp ogt double %176, 1.000000e+00
br i1 %177, label %178, label %179
178: ; preds = %172
br label %197
179: ; preds = %172
%180 = load double, double* %15
%181 = load double, double* %16
%182 = call i32 (double, double, ...) bitcast (i32 (...)* @cdflib_fifdmin1 to i32 (double, double, ...)*)(double 2.000000e-01, double %181)
%183 = sitofp i32 %182 to double
%184 = fcmp oge double %180, %183
br i1 %184, label %185, label %186
185: ; preds = %179
br label %327
186: ; preds = %179
%187 = load double, double* %20
%188 = load double, double* %15
%189 = call double @pow(double %187, double %188)
%190 = fcmp ole double %189, 9.000000e-01
br i1 %190, label %191, label %192
191: ; preds = %186
br label %327
192: ; preds = %186
%193 = load double, double* %20
%194 = fcmp oge double %193, 3.000000e-01
br i1 %194, label %195, label %196
195: ; preds = %192
br label %336
196: ; preds = %192
store i32 20, i32* %25
br label %356
197: ; preds = %178
%198 = load double, double* %16
%199 = fcmp ole double %198, 1.000000e+00
br i1 %199, label %200, label %201
200: ; preds = %197
br label %327
201: ; preds = %197
%202 = load double, double* %20
%203 = fcmp oge double %202, 3.000000e-01
br i1 %203, label %204, label %205
204: ; preds = %201
br label %336
205: ; preds = %201
%206 = load double, double* %20
%207 = fcmp oge double %206, 1.000000e-01
br i1 %207, label %208, label %209
208: ; preds = %205
br label %218
209: ; preds = %205
%210 = load double, double* %20
%211 = load double, double* %16
%212 = fmul double %210, %211
%213 = load double, double* %15
%214 = call double @pow(double %212, double %213)
%215 = fcmp ole double %214, 0x3FE6666666666666
br i1 %215, label %216, label %217
216: ; preds = %209
br label %327
217: ; preds = %209
br label %218
218: ; preds = %217, %208
%219 = load double, double* %16
%220 = fcmp ogt double %219, 1.500000e+01
br i1 %220, label %221, label %222
221: ; preds = %218
br label %364
222: ; preds = %218
store i32 20, i32* %25
br label %356
223: ; preds = %132
%224 = load double*, double** %8
%225 = load double, double* %224
%226 = load double*, double** %9
%227 = load double, double* %226
%228 = fcmp ogt double %225, %227
br i1 %228, label %229, label %230
229: ; preds = %223
br label %242
230: ; preds = %223
%231 = load double*, double** %8
%232 = load double, double* %231
%233 = load double*, double** %8
%234 = load double, double* %233
%235 = load double*, double** %9
%236 = load double, double* %235
%237 = fadd double %234, %236
%238 = load double*, double** %10
%239 = load double, double* %238
%240 = fneg double %237
%241 = call double @llvm.fmuladd.f64(double %240, double %239, double %232)
store double %241, double* %18
br label %254
242: ; preds = %229
%243 = load double*, double** %8
%244 = load double, double* %243
%245 = load double*, double** %9
%246 = load double, double* %245
%247 = fadd double %244, %246
%248 = load double*, double** %11
%249 = load double, double* %248
%250 = load double*, double** %9
%251 = load double, double* %250
%252 = fneg double %251
%253 = call double @llvm.fmuladd.f64(double %247, double %249, double %252)
store double %253, double* %18
br label %254
254: ; preds = %242, %230
%255 = load double, double* %18
%256 = fcmp oge double %255, 0.000000e+00
br i1 %256, label %257, label %258
257: ; preds = %254
br label %269
258: ; preds = %254
store i32 1, i32* %24
%259 = load double*, double** %9
%260 = load double, double* %259
store double %260, double* %15
%261 = load double*, double** %8
%262 = load double, double* %261
store double %262, double* %16
%263 = load double*, double** %11
%264 = load double, double* %263
store double %264, double* %20
%265 = load double*, double** %10
%266 = load double, double* %265
store double %266, double* %21
%267 = load double, double* %18
%268 = call double @llvm.fabs.f64(double %267)
store double %268, double* %18
br label %269
269: ; preds = %258, %257
%270 = load double, double* %16
%271 = fcmp olt double %270, 4.000000e+01
br i1 %271, label %272, label %278
272: ; preds = %269
%273 = load double, double* %16
%274 = load double, double* %20
%275 = fmul double %273, %274
%276 = fcmp ole double %275, 0x3FE6666666666666
br i1 %276, label %277, label %278
277: ; preds = %272
br label %327
278: ; preds = %272, %269
%279 = load double, double* %16
%280 = fcmp olt double %279, 4.000000e+01
br i1 %280, label %281, label %282
281: ; preds = %278
br label %374
282: ; preds = %278
%283 = load double, double* %15
%284 = load double, double* %16
%285 = fcmp ogt double %283, %284
br i1 %285, label %286, label %287
286: ; preds = %282
br label %298
287: ; preds = %282
%288 = load double, double* %15
%289 = fcmp ole double %288, 1.000000e+02
br i1 %289, label %290, label %291
290: ; preds = %287
br label %345
291: ; preds = %287
%292 = load double, double* %18
%293 = load double, double* %15
%294 = fmul double 3.000000e-02, %293
%295 = fcmp ogt double %292, %294
br i1 %295, label %296, label %297
296: ; preds = %291
br label %345
297: ; preds = %291
br label %431
298: ; preds = %286
%299 = load double, double* %16
%300 = fcmp ole double %299, 1.000000e+02
br i1 %300, label %301, label %302
301: ; preds = %298
br label %345
302: ; preds = %298
%303 = load double, double* %18
%304 = load double, double* %16
%305 = fmul double 3.000000e-02, %304
%306 = fcmp ogt double %303, %305
br i1 %306, label %307, label %308
307: ; preds = %302
br label %345
308: ; preds = %302
br label %431
309: ; preds = %156
%310 = call i32 (double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_fpser to i32 (double*, double*, double*, double*, ...)*)(double* %15, double* %16, double* %20, double* %17)
%311 = sitofp i32 %310 to double
%312 = load double*, double** %12
store double %311, double* %312
%313 = load double*, double** %12
%314 = load double, double* %313
%315 = fsub double 5.000000e-01, %314
%316 = fadd double 5.000000e-01, %315
%317 = load double*, double** %13
store double %316, double* %317
br label %460
318: ; preds = %171
%319 = call i32 (double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_apser to i32 (double*, double*, double*, double*, ...)*)(double* %15, double* %16, double* %20, double* %17)
%320 = sitofp i32 %319 to double
%321 = load double*, double** %13
store double %320, double* %321
%322 = load double*, double** %13
%323 = load double, double* %322
%324 = fsub double 5.000000e-01, %323
%325 = fadd double 5.000000e-01, %324
%326 = load double*, double** %12
store double %325, double* %326
br label %460
327: ; preds = %277, %216, %200, %191, %185
%328 = call i32 (double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_beta_pser to i32 (double*, double*, double*, double*, ...)*)(double* %15, double* %16, double* %20, double* %17)
%329 = sitofp i32 %328 to double
%330 = load double*, double** %12
store double %329, double* %330
%331 = load double*, double** %12
%332 = load double, double* %331
%333 = fsub double 5.000000e-01, %332
%334 = fadd double 5.000000e-01, %333
%335 = load double*, double** %13
store double %334, double* %335
br label %460
336: ; preds = %204, %195
%337 = call i32 (double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_beta_pser to i32 (double*, double*, double*, double*, ...)*)(double* %16, double* %15, double* %21, double* %17)
%338 = sitofp i32 %337 to double
%339 = load double*, double** %13
store double %338, double* %339
%340 = load double*, double** %13
%341 = load double, double* %340
%342 = fsub double 5.000000e-01, %341
%343 = fadd double 5.000000e-01, %342
%344 = load double*, double** %12
store double %343, double* %344
br label %460
345: ; preds = %307, %301, %296, %290
%346 = load double, double* %17
%347 = fmul double 1.500000e+01, %346
store double %347, double* %26
%348 = call i32 (double*, double*, double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_beta_frac to i32 (double*, double*, double*, double*, double*, double*, ...)*)(double* %15, double* %16, double* %20, double* %21, double* %18, double* %26)
%349 = sitofp i32 %348 to double
%350 = load double*, double** %12
store double %349, double* %350
%351 = load double*, double** %12
%352 = load double, double* %351
%353 = fsub double 5.000000e-01, %352
%354 = fadd double 5.000000e-01, %353
%355 = load double*, double** %13
store double %354, double* %355
br label %460
356: ; preds = %222, %196
%357 = call i32 (double*, double*, double*, double*, i32*, double*, ...) bitcast (i32 (...)* @cdflib_beta_up to i32 (double*, double*, double*, double*, i32*, double*, ...)*)(double* %16, double* %15, double* %21, double* %20, i32* %25, double* %17)
%358 = sitofp i32 %357 to double
%359 = load double*, double** %13
store double %358, double* %359
%360 = load double, double* %16
%361 = load i32, i32* %25
%362 = sitofp i32 %361 to double
%363 = fadd double %360, %362
store double %363, double* %16
br label %364
364: ; preds = %356, %221
%365 = load double, double* %17
%366 = fmul double 1.500000e+01, %365
store double %366, double* %27
%367 = load double*, double** %13
%368 = call i32 (double*, double*, double*, double*, double*, double*, i32*, ...) bitcast (i32 (...)* @cdflib_beta_grat to i32 (double*, double*, double*, double*, double*, double*, i32*, ...)*)(double* %16, double* %15, double* %21, double* %20, double* %367, double* %27, i32* %23)
%369 = load double*, double** %13
%370 = load double, double* %369
%371 = fsub double 5.000000e-01, %370
%372 = fadd double 5.000000e-01, %371
%373 = load double*, double** %12
store double %372, double* %373
br label %460
374: ; preds = %281
%375 = load double, double* %16
%376 = fptosi double %375 to i32
store i32 %376, i32* %25
%377 = load i32, i32* %25
%378 = sitofp i32 %377 to double
%379 = load double, double* %16
%380 = fsub double %379, %378
store double %380, double* %16
%381 = load double, double* %16
%382 = fcmp une double %381, 0.000000e+00
br i1 %382, label %383, label %384
383: ; preds = %374
br label %387
384: ; preds = %374
%385 = load i32, i32* %25
%386 = sub nsw i32 %385, 1
store i32 %386, i32* %25
store double 1.000000e+00, double* %16
br label %387
387: ; preds = %384, %383
%388 = call i32 (double*, double*, double*, double*, i32*, double*, ...) bitcast (i32 (...)* @cdflib_beta_up to i32 (double*, double*, double*, double*, i32*, double*, ...)*)(double* %16, double* %15, double* %21, double* %20, i32* %25, double* %17)
%389 = sitofp i32 %388 to double
%390 = load double*, double** %12
store double %389, double* %390
%391 = load double, double* %20
%392 = fcmp ogt double %391, 0x3FE6666666666666
br i1 %392, label %393, label %394
393: ; preds = %387
br label %406
394: ; preds = %387
%395 = load double*, double** %12
%396 = load double, double* %395
%397 = call i32 (double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_beta_pser to i32 (double*, double*, double*, double*, ...)*)(double* %15, double* %16, double* %20, double* %17)
%398 = sitofp i32 %397 to double
%399 = fadd double %396, %398
%400 = load double*, double** %12
store double %399, double* %400
%401 = load double*, double** %12
%402 = load double, double* %401
%403 = fsub double 5.000000e-01, %402
%404 = fadd double 5.000000e-01, %403
%405 = load double*, double** %13
store double %404, double* %405
br label %460
406: ; preds = %393
%407 = load double, double* %15
%408 = fcmp ogt double %407, 1.500000e+01
br i1 %408, label %409, label %410
409: ; preds = %406
br label %421
410: ; preds = %406
store i32 20, i32* %25
%411 = load double*, double** %12
%412 = load double, double* %411
%413 = call i32 (double*, double*, double*, double*, i32*, double*, ...) bitcast (i32 (...)* @cdflib_beta_up to i32 (double*, double*, double*, double*, i32*, double*, ...)*)(double* %15, double* %16, double* %20, double* %21, i32* %25, double* %17)
%414 = sitofp i32 %413 to double
%415 = fadd double %412, %414
%416 = load double*, double** %12
store double %415, double* %416
%417 = load double, double* %15
%418 = load i32, i32* %25
%419 = sitofp i32 %418 to double
%420 = fadd double %417, %419
store double %420, double* %15
br label %421
421: ; preds = %410, %409
%422 = load double, double* %17
%423 = fmul double 1.500000e+01, %422
store double %423, double* %28
%424 = load double*, double** %12
%425 = call i32 (double*, double*, double*, double*, double*, double*, i32*, ...) bitcast (i32 (...)* @cdflib_beta_grat to i32 (double*, double*, double*, double*, double*, double*, i32*, ...)*)(double* %15, double* %16, double* %20, double* %21, double* %424, double* %28, i32* %23)
%426 = load double*, double** %12
%427 = load double, double* %426
%428 = fsub double 5.000000e-01, %427
%429 = fadd double 5.000000e-01, %428
%430 = load double*, double** %13
store double %429, double* %430
br label %460
431: ; preds = %308, %297
%432 = load double, double* %17
%433 = fmul double 1.000000e+02, %432
store double %433, double* %29
%434 = call i32 (double*, double*, double*, double*, ...) bitcast (i32 (...)* @cdflib_beta_asym to i32 (double*, double*, double*, double*, ...)*)(double* %15, double* %16, double* %18, double* %29)
%435 = sitofp i32 %434 to double
%436 = load double*, double** %12
store double %435, double* %436
%437 = load double*, double** %12
%438 = load double, double* %437
%439 = fsub double 5.000000e-01, %438
%440 = fadd double 5.000000e-01, %439
%441 = load double*, double** %13
store double %440, double* %441
br label %460
442: ; preds = %88
%443 = load double*, double** %8
%444 = load double, double* %443
%445 = fcmp oeq double %444, 0.000000e+00
br i1 %445, label %446, label %447
446: ; preds = %442
br label %501
447: ; preds = %442
br label %448
448: ; preds = %447, %103
%449 = load double*, double** %12
store double 0.000000e+00, double* %449
%450 = load double*, double** %13
store double 1.000000e+00, double* %450
br label %505
451: ; preds = %93
%452 = load double*, double** %9
%453 = load double, double* %452
%454 = fcmp oeq double %453, 0.000000e+00
br i1 %454, label %455, label %456
455: ; preds = %451
br label %503
456: ; preds = %451
br label %457
457: ; preds = %456, %98
%458 = load double*, double** %12
store double 1.000000e+00, double* %458
%459 = load double*, double** %13
store double 0.000000e+00, double* %459
br label %505
460: ; preds = %431, %421, %394, %364, %345, %336, %327, %318, %309
%461 = load i32, i32* %24
%462 = icmp eq i32 %461, 0
br i1 %462, label %463, label %464
463: ; preds = %460
br label %505
464: ; preds = %460
%465 = load double*, double** %12
%466 = load double, double* %465
store double %466, double* %19
%467 = load double*, double** %13
%468 = load double, double* %467
%469 = load double*, double** %12
store double %468, double* %469
%470 = load double, double* %19
%471 = load double*, double** %13
store double %470, double* %471
br label %505
472: ; preds = %117
%473 = load double*, double** %9
%474 = load double, double* %473
%475 = load double*, double** %8
%476 = load double, double* %475
%477 = load double*, double** %9
%478 = load double, double* %477
%479 = fadd double %476, %478
%480 = fdiv double %474, %479
%481 = load double*, double** %12
store double %480, double* %481
%482 = load double*, double** %8
%483 = load double, double* %482
%484 = load double*, double** %8
%485 = load double, double* %484
%486 = load double*, double** %9
%487 = load double, double* %486
%488 = fadd double %485, %487
%489 = fdiv double %483, %488
%490 = load double*, double** %13
store double %489, double* %490
br label %505
491: ; preds = %41
%492 = load i32*, i32** %14
store i32 1, i32* %492
br label %505
493: ; preds = %50
%494 = load i32*, i32** %14
store i32 2, i32* %494
br label %505
495: ; preds = %59
%496 = load i32*, i32** %14
store i32 3, i32* %496
br label %505
497: ; preds = %68
%498 = load i32*, i32** %14
store i32 4, i32* %498
br label %505
499: ; preds = %82
%500 = load i32*, i32** %14
store i32 5, i32* %500
br label %505
501: ; preds = %446
%502 = load i32*, i32** %14
store i32 6, i32* %502
br label %505
503: ; preds = %455
%504 = load i32*, i32** %14
store i32 7, i32* %504
br label %505
505: ; preds = %503, %501, %499, %497, %495, %493, %491, %472, %464, %463, %457, %448
ret void
}

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