#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double algdiv(double *, double *);
double alnrel(double *);
int /*<<< orphan*/ bcorr(double *, double *);
scalar_t__ betaln(double *, double *);
double esum(int *, double *);
double exp(int /*<<< orphan*/);
double fabs(double);
double fifdmax1(double, double);
double fifdmin1(double, double);
double gam1(double *);
double gamln1(double *);
double log(double);
double rlog1(double *);
double sqrt(double);
double brcmp1(int *mu, double *a, double *b, double *x, double *y) {
  static double Const = .398942280401433e0;
  static double brcmp1, a0, apb, b0, c, e, h, lambda, lnx, lny, t, u, v, x0, y0,
      z;
  static int i, n;
  static double T1, T2, T3, T4;
  a0 = fifdmin1(*a, *b);
  if (a0 >= 8.0e0)
    goto S130;
  if (*x > 0.375e0)
    goto S10;
  lnx = log(*x);
  T1 = -*x;
  lny = alnrel(&T1);
  goto S30;
S10:
  if (*y > 0.375e0)
    goto S20;
  T2 = -*y;
  lnx = alnrel(&T2);
  lny = log(*y);
  goto S30;
S20:
  lnx = log(*x);
  lny = log(*y);
S30:
  z = *a * lnx + *b * lny;
  if (a0 < 1.0e0)
    goto S40;
  z -= betaln(a, b);
  brcmp1 = esum(mu, &z);
  return brcmp1;
S40:
  b0 = fifdmax1(*a, *b);
  if (b0 >= 8.0e0)
    goto S120;
  if (b0 > 1.0e0)
    goto S70;
  brcmp1 = esum(mu, &z);
  if (brcmp1 == 0.0e0)
    return brcmp1;
  apb = *a + *b;
  if (apb > 1.0e0)
    goto S50;
  z = 1.0e0 + gam1(&apb);
  goto S60;
S50:
  u = *a + *b - 1.e0;
  z = (1.0e0 + gam1(&u)) / apb;
S60:
  c = (1.0e0 + gam1(a)) * (1.0e0 + gam1(b)) / z;
  brcmp1 = brcmp1 * (a0 * c) / (1.0e0 + a0 / b0);
  return brcmp1;
S70:
  u = gamln1(&a0);
  n = b0 - 1.0e0;
  if (n < 1)
    goto S90;
  c = 1.0e0;
  for (i = 1; i <= n; i++) {
    b0 -= 1.0e0;
    c *= (b0 / (a0 + b0));
  }
  u = log(c) + u;
S90:
  z -= u;
  b0 -= 1.0e0;
  apb = a0 + b0;
  if (apb > 1.0e0)
    goto S100;
  t = 1.0e0 + gam1(&apb);
  goto S110;
S100:
  u = a0 + b0 - 1.e0;
  t = (1.0e0 + gam1(&u)) / apb;
S110:
  brcmp1 = a0 * esum(mu, &z) * (1.0e0 + gam1(&b0)) / t;
  return brcmp1;
S120:
  u = gamln1(&a0) + algdiv(&a0, &b0);
  T3 = z - u;
  brcmp1 = a0 * esum(mu, &T3);
  return brcmp1;
S130:
  if (*a > *b)
    goto S140;
  h = *a / *b;
  x0 = h / (1.0e0 + h);
  y0 = 1.0e0 / (1.0e0 + h);
  lambda = *a - (*a + *b) * *x;
  goto S150;
S140:
  h = *b / *a;
  x0 = 1.0e0 / (1.0e0 + h);
  y0 = h / (1.0e0 + h);
  lambda = (*a + *b) * *y - *b;
S150:
  e = -(lambda / *a);
  if (fabs(e) > 0.6e0)
    goto S160;
  u = rlog1(&e);
  goto S170;
S160:
  u = e - log(*x / x0);
S170:
  e = lambda / *b;
  if (fabs(e) > 0.6e0)
    goto S180;
  v = rlog1(&e);
  goto S190;
S180:
  v = e - log(*y / y0);
S190:
  T4 = -(*a * u + *b * v);
  z = esum(mu, &T4);
  brcmp1 = Const * sqrt(*b * x0) * z * exp(-bcorr(a, b));
  return brcmp1;
}


x86 assembly (clang, -O0):
<brcmp1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
call   1080 <fifdmin1@plt>
movsd  %xmm0,0x2e94(%rip)
movsd  0x2e8c(%rip),%xmm0
movsd  0xde4(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     122b <brcmp1+0x5b>
jmp    18a0 <brcmp1+0x6d0>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xdcd(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    124a <brcmp1+0x7a>
jmp    129f <brcmp1+0xcf>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
call   10b0 <log@plt>
movsd  %xmm0,0x2e79(%rip)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2e9a(%rip)
lea    0x2e93(%rip),%rdi
call   1040 <alnrel@plt>
movsd  %xmm0,0x2e46(%rip)
jmp    133d <brcmp1+0x16d>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xd59(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12be <brcmp1+0xee>
jmp    1313 <brcmp1+0x143>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2e43(%rip)
lea    0x2e3c(%rip),%rdi
call   1040 <alnrel@plt>
movsd  %xmm0,0x2ddf(%rip)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
call   10b0 <log@plt>
movsd  %xmm0,0x2dd2(%rip)
jmp    133d <brcmp1+0x16d>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
call   10b0 <log@plt>
movsd  %xmm0,0x2db0(%rip)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
call   10b0 <log@plt>
movsd  %xmm0,0x2da3(%rip)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0x2d8b(%rip),%xmm2
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  0x2d83(%rip),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x2da3(%rip)
movsd  0xc9b(%rip),%xmm0
ucomisd 0x2d23(%rip),%xmm0
jbe    1388 <brcmp1+0x1b8>
jmp    13d8 <brcmp1+0x208>
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
call   10d0 <betaln@plt>
cvtsi2sd %rax,%xmm1
movsd  0x2d6e(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2d62(%rip)
mov    -0x10(%rbp),%rdi
lea    0x2d57(%rip),%rsi
call   1050 <esum@plt>
movsd  %xmm0,0x2cd2(%rip)
movsd  0x2cca(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c18 <brcmp1+0xa48>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
call   1070 <fifdmax1@plt>
movsd  %xmm0,0x2cbb(%rip)
movsd  0x2cb3(%rip),%xmm0
movsd  0xbfb(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1414 <brcmp1+0x244>
jmp    180d <brcmp1+0x63d>
movsd  0x2c94(%rip),%xmm0
movsd  0xbec(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1433 <brcmp1+0x263>
jmp    15ca <brcmp1+0x3fa>
mov    -0x10(%rbp),%rdi
lea    0x2cd2(%rip),%rsi
call   1050 <esum@plt>
movsd  %xmm0,0x2c4d(%rip)
movsd  0x2c45(%rip),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1478 <brcmp1+0x2a8>
jp     1478 <brcmp1+0x2a8>
movsd  0x2c2a(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c18 <brcmp1+0xa48>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,0x2c18(%rip)
movsd  0x2c10(%rip),%xmm0
movsd  0xb70(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    14af <brcmp1+0x2df>
jmp    14d7 <brcmp1+0x307>
lea    0x2bf2(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0xb4a(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x2c3e(%rip)
jmp    1526 <brcmp1+0x356>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
addsd  (%rax),%xmm0
movsd  0xb21(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2bf5(%rip)
lea    0x2bee(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0xafe(%rip),%xmm0
addsd  %xmm1,%xmm0
divsd  0x2b8a(%rip),%xmm0
movsd  %xmm0,0x2bea(%rip)
mov    -0x18(%rbp),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0xad6(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm2
movsd  -0x38(%rbp),%xmm0
movsd  0xab4(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  0x2ba4(%rip),%xmm0
movsd  %xmm0,0x2b44(%rip)
movsd  0x2b1c(%rip),%xmm0
movsd  0x2b1c(%rip),%xmm1
mulsd  0x2b2c(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0x2b08(%rip),%xmm2
divsd  0x2b10(%rip),%xmm2
movsd  0xa68(%rip),%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2ae0(%rip)
movsd  0x2ad8(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c18 <brcmp1+0xa48>
lea    0x2acf(%rip),%rdi
call   10f0 <gamln1@plt>
movsd  %xmm0,0x2b12(%rip)
movsd  0x2aca(%rip),%xmm0
movsd  0xa22(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,0x2b20(%rip)
cmpl   $0x1,0x2b19(%rip)
jge    160e <brcmp1+0x43e>
jmp    16b3 <brcmp1+0x4e3>
movsd  0x9fa(%rip),%xmm0
movsd  %xmm0,0x2a9a(%rip)
movl   $0x1,0x2af0(%rip)
mov    0x2aea(%rip),%eax
cmp    0x2ae8(%rip),%eax
jg     1696 <brcmp1+0x4c6>
movsd  0x2a6e(%rip),%xmm0
movsd  0x9c6(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2a5a(%rip)
movsd  0x2a52(%rip),%xmm0
movsd  0x2a3a(%rip),%xmm1
addsd  0x2a42(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  0x2a3e(%rip),%xmm0
movsd  %xmm0,0x2a36(%rip)
mov    0x2a90(%rip),%eax
add    $0x1,%eax
mov    %eax,0x2a87(%rip)
jmp    1628 <brcmp1+0x458>
movsd  0x2a1a(%rip),%xmm0
call   10b0 <log@plt>
addsd  0x2a45(%rip),%xmm0
movsd  %xmm0,0x2a3d(%rip)
movsd  0x2a35(%rip),%xmm1
movsd  0x2a4d(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2a41(%rip)
movsd  0x29d9(%rip),%xmm0
movsd  0x931(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x29c5(%rip)
movsd  0x29ad(%rip),%xmm0
addsd  0x29b5(%rip),%xmm0
movsd  %xmm0,0x29a5(%rip)
movsd  0x299d(%rip),%xmm0
movsd  0x8fd(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1722 <brcmp1+0x552>
jmp    174a <brcmp1+0x57a>
lea    0x297f(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x8d7(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x29a3(%rip)
jmp    1799 <brcmp1+0x5c9>
movsd  0x294e(%rip),%xmm0
addsd  0x2956(%rip),%xmm0
movsd  0x8ae(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2982(%rip)
lea    0x297b(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x88b(%rip),%xmm0
addsd  %xmm1,%xmm0
divsd  0x2917(%rip),%xmm0
movsd  %xmm0,0x294f(%rip)
movsd  0x28ff(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rdi
lea    0x295f(%rip),%rsi
call   1050 <esum@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
lea    0x28e2(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm2
movsd  -0x40(%rbp),%xmm0
movsd  0x82d(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  0x28f5(%rip),%xmm0
movsd  %xmm0,0x289d(%rip)
movsd  0x2895(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c18 <brcmp1+0xa48>
lea    0x288c(%rip),%rdi
call   10f0 <gamln1@plt>
movsd  %xmm0,-0x58(%rbp)
lea    0x287b(%rip),%rdi
lea    0x2884(%rip),%rsi
call   1090 <algdiv@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x28ab(%rip)
movsd  0x28c3(%rip),%xmm0
subsd  0x289b(%rip),%xmm0
movsd  %xmm0,0x28d3(%rip)
movsd  0x283b(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rdi
lea    0x28bb(%rip),%rsi
call   1050 <esum@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x280a(%rip)
movsd  0x2802(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c18 <brcmp1+0xa48>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
ucomisd (%rax),%xmm0
jbe    18bb <brcmp1+0x6eb>
jmp    1967 <brcmp1+0x797>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x27f5(%rip)
movsd  0x27ed(%rip),%xmm0
movsd  0x72d(%rip),%xmm1
addsd  0x27dd(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2809(%rip)
movsd  0x711(%rip),%xmm1
addsd  0x27c1(%rip),%xmm1
movsd  0x701(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x27ed(%rip)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
addsd  (%rax),%xmm0
mov    -0x28(%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,0x276e(%rip)
jmp    1a0e <brcmp1+0x83e>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x2749(%rip)
movsd  0x689(%rip),%xmm1
addsd  0x2739(%rip),%xmm1
movsd  0x679(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x275d(%rip)
movsd  0x271d(%rip),%xmm0
movsd  0x65d(%rip),%xmm1
addsd  0x270d(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2741(%rip)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
addsd  (%rax),%xmm0
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x20(%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,0x26c2(%rip)
movsd  0x26ba(%rip),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2683(%rip)
movsd  0x267b(%rip),%xmm0
movaps 0x5d4(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x5c0(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1a67 <brcmp1+0x897>
jmp    1a80 <brcmp1+0x8b0>
lea    0x2652(%rip),%rdi
call   1060 <rlog1@plt>
movsd  %xmm0,0x2675(%rip)
jmp    1ab6 <brcmp1+0x8e6>
movsd  0x2638(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x2663(%rip),%xmm0
call   10b0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x263a(%rip)
movsd  0x2612(%rip),%xmm0
mov    -0x20(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x25f2(%rip)
movsd  0x25ea(%rip),%xmm0
movaps 0x543(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x52f(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1af8 <brcmp1+0x928>
jmp    1b11 <brcmp1+0x941>
lea    0x25c1(%rip),%rdi
call   1060 <rlog1@plt>
movsd  %xmm0,0x25ec(%rip)
jmp    1b47 <brcmp1+0x977>
movsd  0x25a7(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x25da(%rip),%xmm0
call   10b0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x25b1(%rip)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0x2599(%rip),%xmm2
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  0x2591(%rip),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x25aa(%rip)
mov    -0x10(%rbp),%rdi
lea    0x259f(%rip),%rsi
call   1050 <esum@plt>
movsd  %xmm0,0x256a(%rip)
movsd  0x24da(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  0x253d(%rip),%xmm0
call   10c0 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  0x2534(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
call   10e0 <bcorr@plt>
xor    %edi,%edi
sub    %eax,%edi
call   10a0 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x248d(%rip)
movsd  0x2485(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<brcmp1>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %r8,%r12
mov    %rcx,%rbx
mov    %rdx,%r13
mov    %rsi,%r15
mov    %rdi,%r14
movsd  (%rsi),%xmm0
movsd  (%rdx),%xmm1
call   1080 <fifdmin1@plt>
movsd  %xmm0,0x2e97(%rip)
ucomisd 0xdf7(%rip),%xmm0
jae    1237 <brcmp1+0x67>
movsd  (%rbx),%xmm0
ucomisd 0xe09(%rip),%xmm0
jbe    1295 <brcmp1+0xc5>
movsd  (%r12),%xmm1
ucomisd 0xdf9(%rip),%xmm1
jbe    146c <brcmp1+0x29c>
call   10b0 <log@plt>
jmp    1488 <brcmp1+0x2b8>
movsd  (%r15),%xmm0
movsd  0x0(%r13),%xmm4
ucomisd %xmm4,%xmm0
jbe    12c7 <brcmp1+0xf7>
movapd %xmm4,%xmm2
divsd  %xmm0,%xmm2
movsd  0xdb0(%rip),%xmm1
movapd %xmm2,%xmm3
addsd  %xmm1,%xmm3
divsd  %xmm3,%xmm1
movsd  %xmm1,0x2e6c(%rip)
divsd  %xmm3,%xmm2
movsd  %xmm2,0x2e68(%rip)
movapd %xmm0,%xmm3
addsd  %xmm4,%xmm3
xorpd  0xda8(%rip),%xmm4
movapd %xmm4,%xmm2
mov    %r12,%rax
movapd %xmm3,%xmm4
jmp    130a <brcmp1+0x13a>
call   10b0 <log@plt>
movsd  %xmm0,0x2e1e(%rip)
movsd  (%rbx),%xmm0
xorpd  0xd82(%rip),%xmm0
movlpd %xmm0,0x2e3a(%rip)
lea    0x2e33(%rip),%rdi
call   1040 <alnrel@plt>
jmp    149b <brcmp1+0x2cb>
movapd %xmm0,%xmm1
divsd  %xmm4,%xmm1
movsd  0xd31(%rip),%xmm2
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
divsd  %xmm3,%xmm1
movsd  %xmm1,0x2ded(%rip)
divsd  %xmm3,%xmm2
movsd  %xmm2,0x2de9(%rip)
addsd  %xmm0,%xmm4
xorpd  0xd2d(%rip),%xmm4
movapd %xmm0,%xmm2
mov    %rbx,%rax
mulsd  (%rax),%xmm4
addsd  %xmm2,%xmm4
movsd  %xmm4,0x2d9e(%rip)
xorpd  0xd0e(%rip),%xmm4
divsd  %xmm0,%xmm4
movsd  %xmm4,0x2d82(%rip)
movapd 0xd0a(%rip),%xmm0
andpd  %xmm4,%xmm0
ucomisd 0xcce(%rip),%xmm0
jbe    1361 <brcmp1+0x191>
movsd  (%rbx),%xmm0
divsd  %xmm1,%xmm0
movapd %xmm4,(%rsp)
call   10b0 <log@plt>
movapd (%rsp),%xmm1
subsd  %xmm0,%xmm1
jmp    1371 <brcmp1+0x1a1>
lea    0x2d48(%rip),%rdi
call   1060 <rlog1@plt>
movapd %xmm0,%xmm1
movsd  %xmm1,0x2d57(%rip)
movsd  0x2d37(%rip),%xmm1
divsd  0x0(%r13),%xmm1
movsd  %xmm1,0x2d21(%rip)
movapd 0xca9(%rip),%xmm0
andpd  %xmm1,%xmm0
ucomisd 0xc6d(%rip),%xmm0
jbe    13c8 <brcmp1+0x1f8>
movsd  (%r12),%xmm0
divsd  0x2d2d(%rip),%xmm0
movapd %xmm1,(%rsp)
call   10b0 <log@plt>
movapd (%rsp),%xmm1
subsd  %xmm0,%xmm1
jmp    13d8 <brcmp1+0x208>
lea    0x2ce1(%rip),%rdi
call   1060 <rlog1@plt>
movapd %xmm0,%xmm1
mulsd  0x0(%r13),%xmm1
movsd  (%r15),%xmm0
mulsd  0x2ce5(%rip),%xmm0
addsd  %xmm1,%xmm0
xorpd  0xc39(%rip),%xmm0
movlpd %xmm0,0x2d09(%rip)
lea    0x2d02(%rip),%rsi
mov    %r14,%rdi
call   1050 <esum@plt>
movsd  %xmm0,0x2cd2(%rip)
movsd  0x0(%r13),%xmm0
mulsd  0x2cb4(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1434 <brcmp1+0x264>
sqrtsd %xmm0,%xmm0
jmp    1439 <brcmp1+0x269>
call   10c0 <sqrt@plt>
mulsd  0xbd7(%rip),%xmm0
mulsd  0x2c9f(%rip),%xmm0
movsd  %xmm0,(%rsp)
mov    %r15,%rdi
mov    %r13,%rsi
call   10e0 <bcorr@plt>
neg    %eax
mov    %eax,%edi
call   10a0 <exp@plt>
mulsd  (%rsp),%xmm0
jmp    158c <brcmp1+0x3bc>
xorpd  0xbbc(%rip),%xmm1
movlpd %xmm1,0x2c7c(%rip)
lea    0x2c75(%rip),%rdi
call   1040 <alnrel@plt>
movsd  %xmm0,0x2c30(%rip)
movsd  (%r12),%xmm0
call   10b0 <log@plt>
movsd  (%r15),%xmm1
mulsd  0x0(%r13),%xmm0
mulsd  0x2c12(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x2c2e(%rip)
movsd  0xb46(%rip),%xmm0
ucomisd 0x2bce(%rip),%xmm0
jbe    1556 <brcmp1+0x386>
movsd  (%r15),%xmm0
movsd  0x0(%r13),%xmm1
call   1070 <fifdmax1@plt>
movsd  %xmm0,0x2bc0(%rip)
ucomisd 0xb10(%rip),%xmm0
jae    15a2 <brcmp1+0x3d2>
ucomisd 0xb0a(%rip),%xmm0
jbe    1607 <brcmp1+0x437>
lea    0x2b8d(%rip),%rdi
call   10f0 <gamln1@plt>
movsd  %xmm0,0x2bb8(%rip)
movsd  0x2b88(%rip),%xmm1
movsd  0xb00(%rip),%xmm2
addsd  %xmm1,%xmm2
cvttsd2si %xmm2,%eax
test   %eax,%eax
jle    171e <brcmp1+0x54e>
movsd  0x2b58(%rip),%xmm5
cmp    $0x1,%eax
jne    168b <brcmp1+0x4bb>
movsd  0xab7(%rip),%xmm0
jmp    16d5 <brcmp1+0x505>
mov    %r15,%rdi
mov    %r13,%rsi
call   10d0 <betaln@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
lea    0x2b78(%rip),%rsi
movsd  0x2b70(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x2b64(%rip)
mov    %r14,%rdi
call   1050 <esum@plt>
movsd  %xmm0,0x2afc(%rip)
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
lea    0x2aef(%rip),%rbx
mov    %rbx,%rdi
call   10f0 <gamln1@plt>
movsd  %xmm0,(%rsp)
lea    0x2aeb(%rip),%rsi
mov    %rbx,%rdi
call   1090 <algdiv@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,0x2afe(%rip)
movsd  0x2b0e(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x2b1a(%rip)
lea    0x2b13(%rip),%rsi
movsd  0x2aa3(%rip),%xmm0
movsd  %xmm0,(%rsp)
mov    %r14,%rdi
call   1050 <esum@plt>
jmp    1462 <brcmp1+0x292>
lea    0x2ada(%rip),%rsi
mov    %r14,%rdi
call   1050 <esum@plt>
movsd  %xmm0,0x2a72(%rip)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    162e <brcmp1+0x45e>
jnp    1594 <brcmp1+0x3c4>
movsd  (%r15),%xmm0
addsd  0x0(%r13),%xmm0
movsd  %xmm0,0x2a5f(%rip)
ucomisd 0x9bf(%rip),%xmm0
jbe    17ea <brcmp1+0x61a>
movsd  (%r15),%xmm0
addsd  0x0(%r13),%xmm0
addsd  0x9c6(%rip),%xmm0
movsd  %xmm0,0x2a66(%rip)
lea    0x2a5f(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x98a(%rip),%xmm0
divsd  0x2a1a(%rip),%xmm0
jmp    17fe <brcmp1+0x62e>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movsd  0x970(%rip),%xmm0
movsd  0x988(%rip),%xmm2
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm3
addsd  %xmm5,%xmm3
movapd %xmm1,%xmm4
divsd  %xmm3,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm3
addsd  %xmm5,%xmm3
movapd %xmm1,%xmm0
divsd  %xmm3,%xmm0
mulsd  %xmm4,%xmm0
add    $0xfffffffe,%ecx
jne    16a0 <brcmp1+0x4d0>
test   $0x1,%al
je     16f1 <brcmp1+0x521>
addsd  0x947(%rip),%xmm1
addsd  %xmm1,%xmm5
movapd %xmm1,%xmm2
divsd  %xmm5,%xmm2
mulsd  %xmm2,%xmm0
movsd  %xmm1,0x29af(%rip)
call   10b0 <log@plt>
addsd  0x29ca(%rip),%xmm0
movsd  %xmm0,0x29c2(%rip)
movsd  0x2992(%rip),%xmm2
addsd  0x90a(%rip),%xmm2
movsd  0x29c2(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x29b6(%rip)
movsd  %xmm2,0x296e(%rip)
addsd  0x2956(%rip),%xmm2
movsd  %xmm2,0x2956(%rip)
ucomisd 0x8b6(%rip),%xmm2
jbe    1782 <brcmp1+0x5b2>
addsd  0x8cc(%rip),%xmm2
movsd  %xmm2,0x296c(%rip)
lea    0x2965(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x890(%rip),%xmm0
divsd  0x2920(%rip),%xmm0
jmp    1796 <brcmp1+0x5c6>
lea    0x2917(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x872(%rip),%xmm0
movsd  %xmm0,0x292a(%rip)
movsd  0x28f2(%rip),%xmm0
movsd  %xmm0,(%rsp)
lea    0x2936(%rip),%rsi
mov    %r14,%rdi
call   1050 <esum@plt>
mulsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
lea    0x28dd(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x830(%rip),%xmm0
mulsd  (%rsp),%xmm0
divsd  0x28e3(%rip),%xmm0
jmp    158c <brcmp1+0x3bc>
lea    0x28af(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x80a(%rip),%xmm0
movsd  %xmm0,0x28e2(%rip)
mov    %r15,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x7ee(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
mov    %r13,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x7d1(%rip),%xmm3
addsd  %xmm3,%xmm1
mulsd  (%rsp),%xmm1
movsd  0x2850(%rip),%xmm0
unpcklpd %xmm0,%xmm1
movsd  0x2894(%rip),%xmm2
movhpd 0x284c(%rip),%xmm2
divpd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mulsd  0x2824(%rip),%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm0
jmp    158c <brcmp1+0x3bc>

LLVM IR (-O0):
@brcmp1.Const = internal global double 0x3FD9884533D43656
@brcmp1.brcmp1 = internal global double 0.000000e+00
@brcmp1.a0 = internal global double 0.000000e+00
@brcmp1.apb = internal global double 0.000000e+00
@brcmp1.b0 = internal global double 0.000000e+00
@brcmp1.c = internal global double 0.000000e+00
@brcmp1.e = internal global double 0.000000e+00
@brcmp1.h = internal global double 0.000000e+00
@brcmp1.lambda = internal global double 0.000000e+00
@brcmp1.lnx = internal global double 0.000000e+00
@brcmp1.lny = internal global double 0.000000e+00
@brcmp1.t = internal global double 0.000000e+00
@brcmp1.u = internal global double 0.000000e+00
@brcmp1.v = internal global double 0.000000e+00
@brcmp1.x0 = internal global double 0.000000e+00
@brcmp1.y0 = internal global double 0.000000e+00
@brcmp1.z = internal global double 0.000000e+00
@brcmp1.i = internal global i32 0
@brcmp1.n = internal global i32 0
@brcmp1.T1 = internal global double 0.000000e+00
@brcmp1.T2 = internal global double 0.000000e+00
@brcmp1.T3 = internal global double 0.000000e+00
@brcmp1.T4 = internal global double 0.000000e+00
define double @brcmp1(i32* %0, double* %1, double* %2, double* %3, double* %4) {
%6 = alloca double
%7 = alloca i32*
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
store i32* %0, i32** %7
store double* %1, double** %8
store double* %2, double** %9
store double* %3, double** %10
store double* %4, double** %11
%12 = load double*, double** %8
%13 = load double, double* %12
%14 = load double*, double** %9
%15 = load double, double* %14
%16 = call double @fifdmin1(double %13, double %15)
store double %16, double* @brcmp1.a0
%17 = load double, double* @brcmp1.a0
%18 = fcmp oge double %17, 8.000000e+00
br i1 %18, label %19, label %20
19: ; preds = %5
br label %217
20: ; preds = %5
%21 = load double*, double** %10
%22 = load double, double* %21
%23 = fcmp ogt double %22, 3.750000e-01
br i1 %23, label %24, label %25
24: ; preds = %20
br label %33
25: ; preds = %20
%26 = load double*, double** %10
%27 = load double, double* %26
%28 = call double @log(double %27)
store double %28, double* @brcmp1.lnx
%29 = load double*, double** %10
%30 = load double, double* %29
%31 = fneg double %30
store double %31, double* @brcmp1.T1
%32 = call double @alnrel(double* @brcmp1.T1)
store double %32, double* @brcmp1.lny
br label %53
33: ; preds = %24
%34 = load double*, double** %11
%35 = load double, double* %34
%36 = fcmp ogt double %35, 3.750000e-01
br i1 %36, label %37, label %38
37: ; preds = %33
br label %46
38: ; preds = %33
%39 = load double*, double** %11
%40 = load double, double* %39
%41 = fneg double %40
store double %41, double* @brcmp1.T2
%42 = call double @alnrel(double* @brcmp1.T2)
store double %42, double* @brcmp1.lnx
%43 = load double*, double** %11
%44 = load double, double* %43
%45 = call double @log(double %44)
store double %45, double* @brcmp1.lny
br label %53
46: ; preds = %37
%47 = load double*, double** %10
%48 = load double, double* %47
%49 = call double @log(double %48)
store double %49, double* @brcmp1.lnx
%50 = load double*, double** %11
%51 = load double, double* %50
%52 = call double @log(double %51)
store double %52, double* @brcmp1.lny
br label %53
53: ; preds = %46, %38, %25
%54 = load double*, double** %8
%55 = load double, double* %54
%56 = load double, double* @brcmp1.lnx
%57 = load double*, double** %9
%58 = load double, double* %57
%59 = load double, double* @brcmp1.lny
%60 = fmul double %58, %59
%61 = call double @llvm.fmuladd.f64(double %55, double %56, double %60)
store double %61, double* @brcmp1.z
%62 = load double, double* @brcmp1.a0
%63 = fcmp olt double %62, 1.000000e+00
br i1 %63, label %64, label %65
64: ; preds = %53
br label %75
65: ; preds = %53
%66 = load double*, double** %8
%67 = load double*, double** %9
%68 = call i64 @betaln(double* %66, double* %67)
%69 = sitofp i64 %68 to double
%70 = load double, double* @brcmp1.z
%71 = fsub double %70, %69
store double %71, double* @brcmp1.z
%72 = load i32*, i32** %7
%73 = call double @esum(i32* %72, double* @brcmp1.z)
store double %73, double* @brcmp1.brcmp1
%74 = load double, double* @brcmp1.brcmp1
store double %74, double* %6
br label %339
75: ; preds = %64
%76 = load double*, double** %8
%77 = load double, double* %76
%78 = load double*, double** %9
%79 = load double, double* %78
%80 = call double @fifdmax1(double %77, double %79)
store double %80, double* @brcmp1.b0
%81 = load double, double* @brcmp1.b0
%82 = fcmp oge double %81, 8.000000e+00
br i1 %82, label %83, label %84
83: ; preds = %75
br label %205
84: ; preds = %75
%85 = load double, double* @brcmp1.b0
%86 = fcmp ogt double %85, 1.000000e+00
br i1 %86, label %87, label %88
87: ; preds = %84
br label %139
88: ; preds = %84
%89 = load i32*, i32** %7
%90 = call double @esum(i32* %89, double* @brcmp1.z)
store double %90, double* @brcmp1.brcmp1
%91 = load double, double* @brcmp1.brcmp1
%92 = fcmp oeq double %91, 0.000000e+00
br i1 %92, label %93, label %95
93: ; preds = %88
%94 = load double, double* @brcmp1.brcmp1
store double %94, double* %6
br label %339
95: ; preds = %88
%96 = load double*, double** %8
%97 = load double, double* %96
%98 = load double*, double** %9
%99 = load double, double* %98
%100 = fadd double %97, %99
store double %100, double* @brcmp1.apb
%101 = load double, double* @brcmp1.apb
%102 = fcmp ogt double %101, 1.000000e+00
br i1 %102, label %103, label %104
103: ; preds = %95
br label %107
104: ; preds = %95
%105 = call double @gam1(double* @brcmp1.apb)
%106 = fadd double 1.000000e+00, %105
store double %106, double* @brcmp1.z
br label %118
107: ; preds = %103
%108 = load double*, double** %8
%109 = load double, double* %108
%110 = load double*, double** %9
%111 = load double, double* %110
%112 = fadd double %109, %111
%113 = fsub double %112, 1.000000e+00
store double %113, double* @brcmp1.u
%114 = call double @gam1(double* @brcmp1.u)
%115 = fadd double 1.000000e+00, %114
%116 = load double, double* @brcmp1.apb
%117 = fdiv double %115, %116
store double %117, double* @brcmp1.z
br label %118
118: ; preds = %107, %104
%119 = load double*, double** %8
%120 = call double @gam1(double* %119)
%121 = fadd double 1.000000e+00, %120
%122 = load double*, double** %9
%123 = call double @gam1(double* %122)
%124 = fadd double 1.000000e+00, %123
%125 = fmul double %121, %124
%126 = load double, double* @brcmp1.z
%127 = fdiv double %125, %126
store double %127, double* @brcmp1.c
%128 = load double, double* @brcmp1.brcmp1
%129 = load double, double* @brcmp1.a0
%130 = load double, double* @brcmp1.c
%131 = fmul double %129, %130
%132 = fmul double %128, %131
%133 = load double, double* @brcmp1.a0
%134 = load double, double* @brcmp1.b0
%135 = fdiv double %133, %134
%136 = fadd double 1.000000e+00, %135
%137 = fdiv double %132, %136
store double %137, double* @brcmp1.brcmp1
%138 = load double, double* @brcmp1.brcmp1
store double %138, double* %6
br label %339
139: ; preds = %87
%140 = call double @gamln1(double* @brcmp1.a0)
store double %140, double* @brcmp1.u
%141 = load double, double* @brcmp1.b0
%142 = fsub double %141, 1.000000e+00
%143 = fptosi double %142 to i32
store i32 %143, i32* @brcmp1.n
%144 = load i32, i32* @brcmp1.n
%145 = icmp slt i32 %144, 1
br i1 %145, label %146, label %147
146: ; preds = %139
br label %170
147: ; preds = %139
store double 1.000000e+00, double* @brcmp1.c
store i32 1, i32* @brcmp1.i
br label %148
148: ; preds = %162, %147
%149 = load i32, i32* @brcmp1.i
%150 = load i32, i32* @brcmp1.n
%151 = icmp sle i32 %149, %150
br i1 %151, label %152, label %165
152: ; preds = %148
%153 = load double, double* @brcmp1.b0
%154 = fsub double %153, 1.000000e+00
store double %154, double* @brcmp1.b0
%155 = load double, double* @brcmp1.b0
%156 = load double, double* @brcmp1.a0
%157 = load double, double* @brcmp1.b0
%158 = fadd double %156, %157
%159 = fdiv double %155, %158
%160 = load double, double* @brcmp1.c
%161 = fmul double %160, %159
store double %161, double* @brcmp1.c
br label %162
162: ; preds = %152
%163 = load i32, i32* @brcmp1.i
%164 = add nsw i32 %163, 1
store i32 %164, i32* @brcmp1.i
br label %148
165: ; preds = %148
%166 = load double, double* @brcmp1.c
%167 = call double @log(double %166)
%168 = load double, double* @brcmp1.u
%169 = fadd double %167, %168
store double %169, double* @brcmp1.u
br label %170
170: ; preds = %165, %146
%171 = load double, double* @brcmp1.u
%172 = load double, double* @brcmp1.z
%173 = fsub double %172, %171
store double %173, double* @brcmp1.z
%174 = load double, double* @brcmp1.b0
%175 = fsub double %174, 1.000000e+00
store double %175, double* @brcmp1.b0
%176 = load double, double* @brcmp1.a0
%177 = load double, double* @brcmp1.b0
%178 = fadd double %176, %177
store double %178, double* @brcmp1.apb
%179 = load double, double* @brcmp1.apb
%180 = fcmp ogt double %179, 1.000000e+00
br i1 %180, label %181, label %182
181: ; preds = %170
br label %185
182: ; preds = %170
%183 = call double @gam1(double* @brcmp1.apb)
%184 = fadd double 1.000000e+00, %183
store double %184, double* @brcmp1.t
br label %194
185: ; preds = %181
%186 = load double, double* @brcmp1.a0
%187 = load double, double* @brcmp1.b0
%188 = fadd double %186, %187
%189 = fsub double %188, 1.000000e+00
store double %189, double* @brcmp1.u
%190 = call double @gam1(double* @brcmp1.u)
%191 = fadd double 1.000000e+00, %190
%192 = load double, double* @brcmp1.apb
%193 = fdiv double %191, %192
store double %193, double* @brcmp1.t
br label %194
194: ; preds = %185, %182
%195 = load double, double* @brcmp1.a0
%196 = load i32*, i32** %7
%197 = call double @esum(i32* %196, double* @brcmp1.z)
%198 = fmul double %195, %197
%199 = call double @gam1(double* @brcmp1.b0)
%200 = fadd double 1.000000e+00, %199
%201 = fmul double %198, %200
%202 = load double, double* @brcmp1.t
%203 = fdiv double %201, %202
store double %203, double* @brcmp1.brcmp1
%204 = load double, double* @brcmp1.brcmp1
store double %204, double* %6
br label %339
205: ; preds = %83
%206 = call double @gamln1(double* @brcmp1.a0)
%207 = call double @algdiv(double* @brcmp1.a0, double* @brcmp1.b0)
%208 = fadd double %206, %207
store double %208, double* @brcmp1.u
%209 = load double, double* @brcmp1.z
%210 = load double, double* @brcmp1.u
%211 = fsub double %209, %210
store double %211, double* @brcmp1.T3
%212 = load double, double* @brcmp1.a0
%213 = load i32*, i32** %7
%214 = call double @esum(i32* %213, double* @brcmp1.T3)
%215 = fmul double %212, %214
store double %215, double* @brcmp1.brcmp1
%216 = load double, double* @brcmp1.brcmp1
store double %216, double* %6
br label %339
217: ; preds = %19
%218 = load double*, double** %8
%219 = load double, double* %218
%220 = load double*, double** %9
%221 = load double, double* %220
%222 = fcmp ogt double %219, %221
br i1 %222, label %223, label %224
223: ; preds = %217
br label %248
224: ; preds = %217
%225 = load double*, double** %8
%226 = load double, double* %225
%227 = load double*, double** %9
%228 = load double, double* %227
%229 = fdiv double %226, %228
store double %229, double* @brcmp1.h
%230 = load double, double* @brcmp1.h
%231 = load double, double* @brcmp1.h
%232 = fadd double 1.000000e+00, %231
%233 = fdiv double %230, %232
store double %233, double* @brcmp1.x0
%234 = load double, double* @brcmp1.h
%235 = fadd double 1.000000e+00, %234
%236 = fdiv double 1.000000e+00, %235
store double %236, double* @brcmp1.y0
%237 = load double*, double** %8
%238 = load double, double* %237
%239 = load double*, double** %8
%240 = load double, double* %239
%241 = load double*, double** %9
%242 = load double, double* %241
%243 = fadd double %240, %242
%244 = load double*, double** %10
%245 = load double, double* %244
%246 = fneg double %243
%247 = call double @llvm.fmuladd.f64(double %246, double %245, double %238)
store double %247, double* @brcmp1.lambda
br label %272
248: ; preds = %223
%249 = load double*, double** %9
%250 = load double, double* %249
%251 = load double*, double** %8
%252 = load double, double* %251
%253 = fdiv double %250, %252
store double %253, double* @brcmp1.h
%254 = load double, double* @brcmp1.h
%255 = fadd double 1.000000e+00, %254
%256 = fdiv double 1.000000e+00, %255
store double %256, double* @brcmp1.x0
%257 = load double, double* @brcmp1.h
%258 = load double, double* @brcmp1.h
%259 = fadd double 1.000000e+00, %258
%260 = fdiv double %257, %259
store double %260, double* @brcmp1.y0
%261 = load double*, double** %8
%262 = load double, double* %261
%263 = load double*, double** %9
%264 = load double, double* %263
%265 = fadd double %262, %264
%266 = load double*, double** %11
%267 = load double, double* %266
%268 = load double*, double** %9
%269 = load double, double* %268
%270 = fneg double %269
%271 = call double @llvm.fmuladd.f64(double %265, double %267, double %270)
store double %271, double* @brcmp1.lambda
br label %272
272: ; preds = %248, %224
%273 = load double, double* @brcmp1.lambda
%274 = load double*, double** %8
%275 = load double, double* %274
%276 = fdiv double %273, %275
%277 = fneg double %276
store double %277, double* @brcmp1.e
%278 = load double, double* @brcmp1.e
%279 = call double @llvm.fabs.f64(double %278)
%280 = fcmp ogt double %279, 6.000000e-01
br i1 %280, label %281, label %282
281: ; preds = %272
br label %284
282: ; preds = %272
%283 = call double @rlog1(double* @brcmp1.e)
store double %283, double* @brcmp1.u
br label %292
284: ; preds = %281
%285 = load double, double* @brcmp1.e
%286 = load double*, double** %10
%287 = load double, double* %286
%288 = load double, double* @brcmp1.x0
%289 = fdiv double %287, %288
%290 = call double @log(double %289)
%291 = fsub double %285, %290
store double %291, double* @brcmp1.u
br label %292
292: ; preds = %284, %282
%293 = load double, double* @brcmp1.lambda
%294 = load double*, double** %9
%295 = load double, double* %294
%296 = fdiv double %293, %295
store double %296, double* @brcmp1.e
%297 = load double, double* @brcmp1.e
%298 = call double @llvm.fabs.f64(double %297)
%299 = fcmp ogt double %298, 6.000000e-01
br i1 %299, label %300, label %301
300: ; preds = %292
br label %303
301: ; preds = %292
%302 = call double @rlog1(double* @brcmp1.e)
store double %302, double* @brcmp1.v
br label %311
303: ; preds = %300
%304 = load double, double* @brcmp1.e
%305 = load double*, double** %11
%306 = load double, double* %305
%307 = load double, double* @brcmp1.y0
%308 = fdiv double %306, %307
%309 = call double @log(double %308)
%310 = fsub double %304, %309
store double %310, double* @brcmp1.v
br label %311
311: ; preds = %303, %301
%312 = load double*, double** %8
%313 = load double, double* %312
%314 = load double, double* @brcmp1.u
%315 = load double*, double** %9
%316 = load double, double* %315
%317 = load double, double* @brcmp1.v
%318 = fmul double %316, %317
%319 = call double @llvm.fmuladd.f64(double %313, double %314, double %318)
%320 = fneg double %319
store double %320, double* @brcmp1.T4
%321 = load i32*, i32** %7
%322 = call double @esum(i32* %321, double* @brcmp1.T4)
store double %322, double* @brcmp1.z
%323 = load double, double* @brcmp1.Const
%324 = load double*, double** %9
%325 = load double, double* %324
%326 = load double, double* @brcmp1.x0
%327 = fmul double %325, %326
%328 = call double @sqrt(double %327)
%329 = fmul double %323, %328
%330 = load double, double* @brcmp1.z
%331 = fmul double %329, %330
%332 = load double*, double** %8
%333 = load double*, double** %9
%334 = call i32 @bcorr(double* %332, double* %333)
%335 = sub nsw i32 0, %334
%336 = call double @exp(i32 %335)
%337 = fmul double %331, %336
store double %337, double* @brcmp1.brcmp1
%338 = load double, double* @brcmp1.brcmp1
store double %338, double* %6
br label %339
339: ; preds = %311, %205, %194, %118, %93, %65
%340 = load double, double* %6
ret double %340
}

LLVM IR (-O2):
@brcmp1.brcmp1 = internal unnamed_addr global double 0.000000e+00
@brcmp1.a0 = internal global double 0.000000e+00
@brcmp1.apb = internal global double 0.000000e+00
@brcmp1.b0 = internal global double 0.000000e+00
@brcmp1.e = internal global double 0.000000e+00
@brcmp1.lambda = internal unnamed_addr global double 0.000000e+00
@brcmp1.lnx = internal unnamed_addr global double 0.000000e+00
@brcmp1.t = internal unnamed_addr global double 0.000000e+00
@brcmp1.u = internal global double 0.000000e+00
@brcmp1.x0 = internal unnamed_addr global double 0.000000e+00
@brcmp1.y0 = internal unnamed_addr global double 0.000000e+00
@brcmp1.z = internal global double 0.000000e+00
@brcmp1.T1 = internal global double 0.000000e+00
@brcmp1.T2 = internal global double 0.000000e+00
@brcmp1.T3 = internal global double 0.000000e+00
@brcmp1.T4 = internal global double 0.000000e+00
define double @brcmp1(i32* %0, double* %1, double* %2, double* nocapture %3, double* nocapture %4) {
%6 = load double, double* %1
%7 = load double, double* %2
%8 = tail call double @fifdmin1(double %6, double %7)
store double %8, double* @brcmp1.a0
%9 = fcmp ult double %8, 8.000000e+00
br i1 %9, label %10, label %175
10: ; preds = %5
%11 = load double, double* %3
%12 = fcmp ogt double %11, 3.750000e-01
br i1 %12, label %18, label %13
13: ; preds = %10
%14 = tail call double @log(double %11)
store double %14, double* @brcmp1.lnx
%15 = load double, double* %3
%16 = fneg double %15
store double %16, double* @brcmp1.T1
%17 = tail call double @alnrel(double* @brcmp1.T1)
br label %30
18: ; preds = %10
%19 = load double, double* %4
%20 = fcmp ogt double %19, 3.750000e-01
br i1 %20, label %26, label %21
21: ; preds = %18
%22 = fneg double %19
store double %22, double* @brcmp1.T2
%23 = tail call double @alnrel(double* @brcmp1.T2)
store double %23, double* @brcmp1.lnx
%24 = load double, double* %4
%25 = tail call double @log(double %24)
br label %30
26: ; preds = %18
%27 = tail call double @log(double %11)
store double %27, double* @brcmp1.lnx
%28 = load double, double* %4
%29 = tail call double @log(double %28)
br label %30
30: ; preds = %26, %21, %13
%31 = phi double [ %29, %26 ], [ %25, %21 ], [ %17, %13 ]
%32 = load double, double* %1
%33 = load double, double* @brcmp1.lnx
%34 = load double, double* %2
%35 = fmul double %34, %31
%36 = tail call double @llvm.fmuladd.f64(double %32, double %33, double %35)
store double %36, double* @brcmp1.z
%37 = load double, double* @brcmp1.a0
%38 = fcmp olt double %37, 1.000000e+00
br i1 %38, label %45, label %39
39: ; preds = %30
%40 = tail call i64 @betaln(double* %1, double* %2)
%41 = sitofp i64 %40 to double
%42 = load double, double* @brcmp1.z
%43 = fsub double %42, %41
store double %43, double* @brcmp1.z
%44 = tail call double @esum(i32* %0, double* @brcmp1.z)
br label %246
45: ; preds = %30
%46 = load double, double* %1
%47 = load double, double* %2
%48 = tail call double @fifdmax1(double %46, double %47)
store double %48, double* @brcmp1.b0
%49 = fcmp ult double %48, 8.000000e+00
br i1 %49, label %50, label %166
50: ; preds = %45
%51 = fcmp ogt double %48, 1.000000e+00
br i1 %51, label %94, label %52
52: ; preds = %50
%53 = tail call double @esum(i32* %0, double* @brcmp1.z)
store double %53, double* @brcmp1.brcmp1
%54 = fcmp oeq double %53, 0.000000e+00
br i1 %54, label %248, label %55
55: ; preds = %52
%56 = load double, double* %1
%57 = load double, double* %2
%58 = fadd double %56, %57
store double %58, double* @brcmp1.apb
%59 = fcmp ogt double %58, 1.000000e+00
br i1 %59, label %63, label %60
60: ; preds = %55
%61 = tail call double @gam1(double* @brcmp1.apb)
%62 = fadd double %61, 1.000000e+00
br label %72
63: ; preds = %55
%64 = load double, double* %1
%65 = load double, double* %2
%66 = fadd double %64, %65
%67 = fadd double %66, -1.000000e+00
store double %67, double* @brcmp1.u
%68 = tail call double @gam1(double* @brcmp1.u)
%69 = fadd double %68, 1.000000e+00
%70 = load double, double* @brcmp1.apb
%71 = fdiv double %69, %70
br label %72
72: ; preds = %63, %60
%73 = phi double [ %62, %60 ], [ %71, %63 ]
store double %73, double* @brcmp1.z
%74 = tail call double @gam1(double* %1)
%75 = fadd double %74, 1.000000e+00
%76 = tail call double @gam1(double* %2)
%77 = fadd double %76, 1.000000e+00
%78 = fmul double %75, %77
%79 = load double, double* @brcmp1.z
%80 = load double, double* @brcmp1.brcmp1
%81 = load double, double* @brcmp1.a0
%82 = load double, double* @brcmp1.b0
%83 = insertelement <2 x double> poison, double %78, i64 0
%84 = insertelement <2 x double> %83, double %81, i64 1
%85 = insertelement <2 x double> poison, double %79, i64 0
%86 = insertelement <2 x double> %85, double %82, i64 1
%87 = fdiv <2 x double> %84, %86
%88 = extractelement <2 x double> %87, i64 0
%89 = fmul double %88, %81
%90 = fmul double %80, %89
%91 = extractelement <2 x double> %87, i64 1
%92 = fadd double %91, 1.000000e+00
%93 = fdiv double %90, %92
br label %246
94: ; preds = %50
%95 = tail call double @gamln1(double* @brcmp1.a0)
store double %95, double* @brcmp1.u
%96 = load double, double* @brcmp1.b0
%97 = fadd double %96, -1.000000e+00
%98 = fptosi double %97 to i32
%99 = icmp slt i32 %98, 1
br i1 %99, label %139, label %100
100: ; preds = %94
%101 = load double, double* @brcmp1.a0
%102 = and i32 %98, 1
%103 = icmp eq i32 %98, 1
br i1 %103, label %120, label %104
104: ; preds = %100
%105 = and i32 %98, -2
br label %106
106: ; preds = %106, %104
%107 = phi double [ %96, %104 ], [ %114, %106 ]
%108 = phi double [ 1.000000e+00, %104 ], [ %117, %106 ]
%109 = phi i32 [ 0, %104 ], [ %118, %106 ]
%110 = fadd double %107, -1.000000e+00
%111 = fadd double %110, %101
%112 = fdiv double %110, %111
%113 = fmul double %108, %112
%114 = fadd double %110, -1.000000e+00
%115 = fadd double %114, %101
%116 = fdiv double %114, %115
%117 = fmul double %113, %116
%118 = add i32 %109, 2
%119 = icmp eq i32 %118, %105
br i1 %119, label %120, label %106
120: ; preds = %106, %100
%121 = phi double [ undef, %100 ], [ %114, %106 ]
%122 = phi double [ undef, %100 ], [ %117, %106 ]
%123 = phi double [ %96, %100 ], [ %114, %106 ]
%124 = phi double [ 1.000000e+00, %100 ], [ %117, %106 ]
%125 = icmp eq i32 %102, 0
br i1 %125, label %131, label %126
126: ; preds = %120
%127 = fadd double %123, -1.000000e+00
%128 = fadd double %127, %101
%129 = fdiv double %127, %128
%130 = fmul double %124, %129
br label %131
131: ; preds = %120, %126
%132 = phi double [ %121, %120 ], [ %127, %126 ]
%133 = phi double [ %122, %120 ], [ %130, %126 ]
store double %132, double* @brcmp1.b0
%134 = tail call double @log(double %133)
%135 = load double, double* @brcmp1.u
%136 = fadd double %134, %135
store double %136, double* @brcmp1.u
%137 = load double, double* @brcmp1.b0
%138 = fadd double %137, -1.000000e+00
br label %139
139: ; preds = %94, %131
%140 = phi double [ %97, %94 ], [ %138, %131 ]
%141 = phi double [ %95, %94 ], [ %136, %131 ]
%142 = load double, double* @brcmp1.z
%143 = fsub double %142, %141
store double %143, double* @brcmp1.z
store double %140, double* @brcmp1.b0
%144 = load double, double* @brcmp1.a0
%145 = fadd double %140, %144
store double %145, double* @brcmp1.apb
%146 = fcmp ogt double %145, 1.000000e+00
br i1 %146, label %150, label %147
147: ; preds = %139
%148 = tail call double @gam1(double* @brcmp1.apb)
%149 = fadd double %148, 1.000000e+00
br label %156
150: ; preds = %139
%151 = fadd double %145, -1.000000e+00
store double %151, double* @brcmp1.u
%152 = tail call double @gam1(double* @brcmp1.u)
%153 = fadd double %152, 1.000000e+00
%154 = load double, double* @brcmp1.apb
%155 = fdiv double %153, %154
br label %156
156: ; preds = %150, %147
%157 = phi double [ %149, %147 ], [ %155, %150 ]
store double %157, double* @brcmp1.t
%158 = load double, double* @brcmp1.a0
%159 = tail call double @esum(i32* %0, double* @brcmp1.z)
%160 = fmul double %158, %159
%161 = tail call double @gam1(double* @brcmp1.b0)
%162 = fadd double %161, 1.000000e+00
%163 = fmul double %160, %162
%164 = load double, double* @brcmp1.t
%165 = fdiv double %163, %164
br label %246
166: ; preds = %45
%167 = tail call double @gamln1(double* @brcmp1.a0)
%168 = tail call double @algdiv(double* @brcmp1.a0, double* @brcmp1.b0)
%169 = fadd double %167, %168
store double %169, double* @brcmp1.u
%170 = load double, double* @brcmp1.z
%171 = fsub double %170, %169
store double %171, double* @brcmp1.T3
%172 = load double, double* @brcmp1.a0
%173 = tail call double @esum(i32* %0, double* @brcmp1.T3)
%174 = fmul double %172, %173
br label %246
175: ; preds = %5
%176 = load double, double* %1
%177 = load double, double* %2
%178 = fcmp ogt double %176, %177
br i1 %178, label %186, label %179
179: ; preds = %175
%180 = fdiv double %176, %177
%181 = fadd double %180, 1.000000e+00
%182 = fdiv double %180, %181
store double %182, double* @brcmp1.x0
%183 = fdiv double 1.000000e+00, %181
store double %183, double* @brcmp1.y0
%184 = fadd double %176, %177
%185 = fneg double %184
br label %193
186: ; preds = %175
%187 = fdiv double %177, %176
%188 = fadd double %187, 1.000000e+00
%189 = fdiv double 1.000000e+00, %188
store double %189, double* @brcmp1.x0
%190 = fdiv double %187, %188
store double %190, double* @brcmp1.y0
%191 = fadd double %176, %177
%192 = fneg double %177
br label %193
193: ; preds = %186, %179
%194 = phi double [ %192, %186 ], [ %176, %179 ]
%195 = phi double* [ %4, %186 ], [ %3, %179 ]
%196 = phi double [ %191, %186 ], [ %185, %179 ]
%197 = phi double [ %189, %186 ], [ %182, %179 ]
%198 = load double, double* %195
%199 = tail call double @llvm.fmuladd.f64(double %196, double %198, double %194)
store double %199, double* @brcmp1.lambda
%200 = fneg double %199
%201 = fdiv double %200, %176
store double %201, double* @brcmp1.e
%202 = tail call double @llvm.fabs.f64(double %201)
%203 = fcmp ogt double %202, 6.000000e-01
br i1 %203, label %206, label %204
204: ; preds = %193
%205 = tail call double @rlog1(double* @brcmp1.e)
br label %211
206: ; preds = %193
%207 = load double, double* %3
%208 = fdiv double %207, %197
%209 = tail call double @log(double %208)
%210 = fsub double %201, %209
br label %211
211: ; preds = %206, %204
%212 = phi double [ %205, %204 ], [ %210, %206 ]
store double %212, double* @brcmp1.u
%213 = load double, double* @brcmp1.lambda
%214 = load double, double* %2
%215 = fdiv double %213, %214
store double %215, double* @brcmp1.e
%216 = tail call double @llvm.fabs.f64(double %215)
%217 = fcmp ogt double %216, 6.000000e-01
br i1 %217, label %220, label %218
218: ; preds = %211
%219 = tail call double @rlog1(double* @brcmp1.e)
br label %226
220: ; preds = %211
%221 = load double, double* %4
%222 = load double, double* @brcmp1.y0
%223 = fdiv double %221, %222
%224 = tail call double @log(double %223)
%225 = fsub double %215, %224
br label %226
226: ; preds = %220, %218
%227 = phi double [ %219, %218 ], [ %225, %220 ]
%228 = load double, double* %1
%229 = load double, double* @brcmp1.u
%230 = load double, double* %2
%231 = fmul double %227, %230
%232 = tail call double @llvm.fmuladd.f64(double %228, double %229, double %231)
%233 = fneg double %232
store double %233, double* @brcmp1.T4
%234 = tail call double @esum(i32* %0, double* @brcmp1.T4)
store double %234, double* @brcmp1.z
%235 = load double, double* %2
%236 = load double, double* @brcmp1.x0
%237 = fmul double %235, %236
%238 = tail call double @sqrt(double %237)
%239 = fmul double %238, 0x3FD9884533D43656
%240 = load double, double* @brcmp1.z
%241 = fmul double %239, %240
%242 = tail call i32 @bcorr(double* %1, double* %2)
%243 = sub nsw i32 0, %242
%244 = tail call double @exp(i32 %243)
%245 = fmul double %241, %244
br label %246
246: ; preds = %39, %72, %156, %166, %226
%247 = phi double [ %245, %226 ], [ %174, %166 ], [ %165, %156 ], [ %93, %72 ], [ %44, %39 ]
store double %247, double* @brcmp1.brcmp1
br label %248
248: ; preds = %246, %52
%249 = phi double [ %53, %52 ], [ %247, %246 ]
ret double %249
}