#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 = (long)(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    18a1 <brcmp1+0x6d1>
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    1c19 <brcmp1+0xa49>
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    180e <brcmp1+0x63e>
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    1c19 <brcmp1+0xa49>
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    1c19 <brcmp1+0xa49>
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,%rax
mov    %eax,0x2b1f(%rip)
cmpl   $0x1,0x2b18(%rip)
jge    160f <brcmp1+0x43f>
jmp    16b4 <brcmp1+0x4e4>
movsd  0x9f9(%rip),%xmm0
movsd  %xmm0,0x2a99(%rip)
movl   $0x1,0x2aef(%rip)
mov    0x2ae9(%rip),%eax
cmp    0x2ae7(%rip),%eax
jg     1697 <brcmp1+0x4c7>
movsd  0x2a6d(%rip),%xmm0
movsd  0x9c5(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2a59(%rip)
movsd  0x2a51(%rip),%xmm0
movsd  0x2a39(%rip),%xmm1
addsd  0x2a41(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  0x2a3d(%rip),%xmm0
movsd  %xmm0,0x2a35(%rip)
mov    0x2a8f(%rip),%eax
add    $0x1,%eax
mov    %eax,0x2a86(%rip)
jmp    1629 <brcmp1+0x459>
movsd  0x2a19(%rip),%xmm0
call   10b0 <log@plt>
addsd  0x2a44(%rip),%xmm0
movsd  %xmm0,0x2a3c(%rip)
movsd  0x2a34(%rip),%xmm1
movsd  0x2a4c(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2a40(%rip)
movsd  0x29d8(%rip),%xmm0
movsd  0x930(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x29c4(%rip)
movsd  0x29ac(%rip),%xmm0
addsd  0x29b4(%rip),%xmm0
movsd  %xmm0,0x29a4(%rip)
movsd  0x299c(%rip),%xmm0
movsd  0x8fc(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1723 <brcmp1+0x553>
jmp    174b <brcmp1+0x57b>
lea    0x297e(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x8d6(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x29a2(%rip)
jmp    179a <brcmp1+0x5ca>
movsd  0x294d(%rip),%xmm0
addsd  0x2955(%rip),%xmm0
movsd  0x8ad(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2981(%rip)
lea    0x297a(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x88a(%rip),%xmm0
addsd  %xmm1,%xmm0
divsd  0x2916(%rip),%xmm0
movsd  %xmm0,0x294e(%rip)
movsd  0x28fe(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rdi
lea    0x295e(%rip),%rsi
call   1050 <esum@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
lea    0x28e1(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm2
movsd  -0x40(%rbp),%xmm0
movsd  0x82c(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  0x28f4(%rip),%xmm0
movsd  %xmm0,0x289c(%rip)
movsd  0x2894(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c19 <brcmp1+0xa49>
lea    0x288b(%rip),%rdi
call   10f0 <gamln1@plt>
movsd  %xmm0,-0x58(%rbp)
lea    0x287a(%rip),%rdi
lea    0x2883(%rip),%rsi
call   1090 <algdiv@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x28aa(%rip)
movsd  0x28c2(%rip),%xmm0
subsd  0x289a(%rip),%xmm0
movsd  %xmm0,0x28d2(%rip)
movsd  0x283a(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rdi
lea    0x28ba(%rip),%rsi
call   1050 <esum@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x2809(%rip)
movsd  0x2801(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c19 <brcmp1+0xa49>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
ucomisd (%rax),%xmm0
jbe    18bc <brcmp1+0x6ec>
jmp    1968 <brcmp1+0x798>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x27f4(%rip)
movsd  0x27ec(%rip),%xmm0
movsd  0x72c(%rip),%xmm1
addsd  0x27dc(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2808(%rip)
movsd  0x710(%rip),%xmm1
addsd  0x27c0(%rip),%xmm1
movsd  0x700(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x27ec(%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,0x276d(%rip)
jmp    1a0f <brcmp1+0x83f>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x2748(%rip)
movsd  0x688(%rip),%xmm1
addsd  0x2738(%rip),%xmm1
movsd  0x678(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x275c(%rip)
movsd  0x271c(%rip),%xmm0
movsd  0x65c(%rip),%xmm1
addsd  0x270c(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2740(%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,0x26c1(%rip)
movsd  0x26b9(%rip),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2682(%rip)
movsd  0x267a(%rip),%xmm0
movaps 0x5d3(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x5bf(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1a68 <brcmp1+0x898>
jmp    1a81 <brcmp1+0x8b1>
lea    0x2651(%rip),%rdi
call   1060 <rlog1@plt>
movsd  %xmm0,0x2674(%rip)
jmp    1ab7 <brcmp1+0x8e7>
movsd  0x2637(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x2662(%rip),%xmm0
call   10b0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2639(%rip)
movsd  0x2611(%rip),%xmm0
mov    -0x20(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x25f1(%rip)
movsd  0x25e9(%rip),%xmm0
movaps 0x542(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x52e(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1af9 <brcmp1+0x929>
jmp    1b12 <brcmp1+0x942>
lea    0x25c0(%rip),%rdi
call   1060 <rlog1@plt>
movsd  %xmm0,0x25eb(%rip)
jmp    1b48 <brcmp1+0x978>
movsd  0x25a6(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x25d9(%rip),%xmm0
call   10b0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x25b0(%rip)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0x2598(%rip),%xmm2
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  0x2590(%rip),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x25a9(%rip)
mov    -0x10(%rbp),%rdi
lea    0x259e(%rip),%rsi
call   1050 <esum@plt>
movsd  %xmm0,0x2569(%rip)
movsd  0x24d9(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  0x253c(%rip),%xmm0
call   10c0 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  0x2533(%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,0x248c(%rip)
movsd  0x2484(%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    158d <brcmp1+0x3bd>
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    1557 <brcmp1+0x387>
movsd  (%r15),%xmm0
movsd  0x0(%r13),%xmm1
call   1070 <fifdmax1@plt>
movsd  %xmm0,0x2bc0(%rip)
ucomisd 0xb10(%rip),%xmm0
jae    15a3 <brcmp1+0x3d3>
ucomisd 0xb0a(%rip),%xmm0
jbe    1608 <brcmp1+0x438>
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,%rax
test   %eax,%eax
jle    172e <brcmp1+0x55e>
movsd  0x2b57(%rip),%xmm5
cmp    $0x1,%eax
jne    168c <brcmp1+0x4bc>
movsd  0xab6(%rip),%xmm0
jmp    16e5 <brcmp1+0x515>
mov    %r15,%rdi
mov    %r13,%rsi
call   10d0 <betaln@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
lea    0x2b77(%rip),%rsi
movsd  0x2b6f(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x2b63(%rip)
mov    %r14,%rdi
call   1050 <esum@plt>
movsd  %xmm0,0x2afb(%rip)
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
lea    0x2aee(%rip),%rbx
mov    %rbx,%rdi
call   10f0 <gamln1@plt>
movsd  %xmm0,(%rsp)
lea    0x2aea(%rip),%rsi
mov    %rbx,%rdi
call   1090 <algdiv@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,0x2afd(%rip)
movsd  0x2b0d(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x2b19(%rip)
lea    0x2b12(%rip),%rsi
movsd  0x2aa2(%rip),%xmm0
movsd  %xmm0,(%rsp)
mov    %r14,%rdi
call   1050 <esum@plt>
jmp    1462 <brcmp1+0x292>
lea    0x2ad9(%rip),%rsi
mov    %r14,%rdi
call   1050 <esum@plt>
movsd  %xmm0,0x2a71(%rip)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    162f <brcmp1+0x45f>
jnp    1595 <brcmp1+0x3c5>
movsd  (%r15),%xmm0
addsd  0x0(%r13),%xmm0
movsd  %xmm0,0x2a5e(%rip)
ucomisd 0x9be(%rip),%xmm0
jbe    17fa <brcmp1+0x62a>
movsd  (%r15),%xmm0
addsd  0x0(%r13),%xmm0
addsd  0x9c5(%rip),%xmm0
movsd  %xmm0,0x2a65(%rip)
lea    0x2a5e(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x989(%rip),%xmm0
divsd  0x2a19(%rip),%xmm0
jmp    180e <brcmp1+0x63e>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movsd  0x96f(%rip),%xmm0
movsd  0x987(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
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    16b0 <brcmp1+0x4e0>
test   $0x1,%al
je     1701 <brcmp1+0x531>
addsd  0x937(%rip),%xmm1
addsd  %xmm1,%xmm5
movapd %xmm1,%xmm2
divsd  %xmm5,%xmm2
mulsd  %xmm2,%xmm0
movsd  %xmm1,0x299f(%rip)
call   10b0 <log@plt>
addsd  0x29ba(%rip),%xmm0
movsd  %xmm0,0x29b2(%rip)
movsd  0x2982(%rip),%xmm2
addsd  0x8fa(%rip),%xmm2
movsd  0x29b2(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x29a6(%rip)
movsd  %xmm2,0x295e(%rip)
addsd  0x2946(%rip),%xmm2
movsd  %xmm2,0x2946(%rip)
ucomisd 0x8a6(%rip),%xmm2
jbe    1792 <brcmp1+0x5c2>
addsd  0x8bc(%rip),%xmm2
movsd  %xmm2,0x295c(%rip)
lea    0x2955(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x880(%rip),%xmm0
divsd  0x2910(%rip),%xmm0
jmp    17a6 <brcmp1+0x5d6>
lea    0x2907(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x862(%rip),%xmm0
movsd  %xmm0,0x291a(%rip)
movsd  0x28e2(%rip),%xmm0
movsd  %xmm0,(%rsp)
lea    0x2926(%rip),%rsi
mov    %r14,%rdi
call   1050 <esum@plt>
mulsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
lea    0x28cd(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x820(%rip),%xmm0
mulsd  (%rsp),%xmm0
divsd  0x28d3(%rip),%xmm0
jmp    158d <brcmp1+0x3bd>
lea    0x289f(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x7fa(%rip),%xmm0
movsd  %xmm0,0x28d2(%rip)
mov    %r15,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x7de(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
mov    %r13,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x7c1(%rip),%xmm3
addsd  %xmm3,%xmm1
mulsd  (%rsp),%xmm1
movsd  0x2840(%rip),%xmm0
unpcklpd %xmm0,%xmm1
movsd  0x2884(%rip),%xmm2
movhpd 0x283c(%rip),%xmm2
divpd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mulsd  0x2814(%rip),%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm0
jmp    158d <brcmp1+0x3bd>

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 %218
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 %340
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 %206
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 %340
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 %340
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 i64
%144 = trunc i64 %143 to i32
store i32 %144, i32* @brcmp1.n
%145 = load i32, i32* @brcmp1.n
%146 = icmp slt i32 %145, 1
br i1 %146, label %147, label %148
147: ; preds = %139
br label %171
148: ; preds = %139
store double 1.000000e+00, double* @brcmp1.c
store i32 1, i32* @brcmp1.i
br label %149
149: ; preds = %163, %148
%150 = load i32, i32* @brcmp1.i
%151 = load i32, i32* @brcmp1.n
%152 = icmp sle i32 %150, %151
br i1 %152, label %153, label %166
153: ; preds = %149
%154 = load double, double* @brcmp1.b0
%155 = fsub double %154, 1.000000e+00
store double %155, double* @brcmp1.b0
%156 = load double, double* @brcmp1.b0
%157 = load double, double* @brcmp1.a0
%158 = load double, double* @brcmp1.b0
%159 = fadd double %157, %158
%160 = fdiv double %156, %159
%161 = load double, double* @brcmp1.c
%162 = fmul double %161, %160
store double %162, double* @brcmp1.c
br label %163
163: ; preds = %153
%164 = load i32, i32* @brcmp1.i
%165 = add nsw i32 %164, 1
store i32 %165, i32* @brcmp1.i
br label %149
166: ; preds = %149
%167 = load double, double* @brcmp1.c
%168 = call double @log(double %167)
%169 = load double, double* @brcmp1.u
%170 = fadd double %168, %169
store double %170, double* @brcmp1.u
br label %171
171: ; preds = %166, %147
%172 = load double, double* @brcmp1.u
%173 = load double, double* @brcmp1.z
%174 = fsub double %173, %172
store double %174, double* @brcmp1.z
%175 = load double, double* @brcmp1.b0
%176 = fsub double %175, 1.000000e+00
store double %176, double* @brcmp1.b0
%177 = load double, double* @brcmp1.a0
%178 = load double, double* @brcmp1.b0
%179 = fadd double %177, %178
store double %179, double* @brcmp1.apb
%180 = load double, double* @brcmp1.apb
%181 = fcmp ogt double %180, 1.000000e+00
br i1 %181, label %182, label %183
182: ; preds = %171
br label %186
183: ; preds = %171
%184 = call double @gam1(double* @brcmp1.apb)
%185 = fadd double 1.000000e+00, %184
store double %185, double* @brcmp1.t
br label %195
186: ; preds = %182
%187 = load double, double* @brcmp1.a0
%188 = load double, double* @brcmp1.b0
%189 = fadd double %187, %188
%190 = fsub double %189, 1.000000e+00
store double %190, double* @brcmp1.u
%191 = call double @gam1(double* @brcmp1.u)
%192 = fadd double 1.000000e+00, %191
%193 = load double, double* @brcmp1.apb
%194 = fdiv double %192, %193
store double %194, double* @brcmp1.t
br label %195
195: ; preds = %186, %183
%196 = load double, double* @brcmp1.a0
%197 = load i32*, i32** %7
%198 = call double @esum(i32* %197, double* @brcmp1.z)
%199 = fmul double %196, %198
%200 = call double @gam1(double* @brcmp1.b0)
%201 = fadd double 1.000000e+00, %200
%202 = fmul double %199, %201
%203 = load double, double* @brcmp1.t
%204 = fdiv double %202, %203
store double %204, double* @brcmp1.brcmp1
%205 = load double, double* @brcmp1.brcmp1
store double %205, double* %6
br label %340
206: ; preds = %83
%207 = call double @gamln1(double* @brcmp1.a0)
%208 = call double @algdiv(double* @brcmp1.a0, double* @brcmp1.b0)
%209 = fadd double %207, %208
store double %209, double* @brcmp1.u
%210 = load double, double* @brcmp1.z
%211 = load double, double* @brcmp1.u
%212 = fsub double %210, %211
store double %212, double* @brcmp1.T3
%213 = load double, double* @brcmp1.a0
%214 = load i32*, i32** %7
%215 = call double @esum(i32* %214, double* @brcmp1.T3)
%216 = fmul double %213, %215
store double %216, double* @brcmp1.brcmp1
%217 = load double, double* @brcmp1.brcmp1
store double %217, double* %6
br label %340
218: ; preds = %19
%219 = load double*, double** %8
%220 = load double, double* %219
%221 = load double*, double** %9
%222 = load double, double* %221
%223 = fcmp ogt double %220, %222
br i1 %223, label %224, label %225
224: ; preds = %218
br label %249
225: ; preds = %218
%226 = load double*, double** %8
%227 = load double, double* %226
%228 = load double*, double** %9
%229 = load double, double* %228
%230 = fdiv double %227, %229
store double %230, double* @brcmp1.h
%231 = load double, double* @brcmp1.h
%232 = load double, double* @brcmp1.h
%233 = fadd double 1.000000e+00, %232
%234 = fdiv double %231, %233
store double %234, double* @brcmp1.x0
%235 = load double, double* @brcmp1.h
%236 = fadd double 1.000000e+00, %235
%237 = fdiv double 1.000000e+00, %236
store double %237, double* @brcmp1.y0
%238 = load double*, double** %8
%239 = load double, double* %238
%240 = load double*, double** %8
%241 = load double, double* %240
%242 = load double*, double** %9
%243 = load double, double* %242
%244 = fadd double %241, %243
%245 = load double*, double** %10
%246 = load double, double* %245
%247 = fneg double %244
%248 = call double @llvm.fmuladd.f64(double %247, double %246, double %239)
store double %248, double* @brcmp1.lambda
br label %273
249: ; preds = %224
%250 = load double*, double** %9
%251 = load double, double* %250
%252 = load double*, double** %8
%253 = load double, double* %252
%254 = fdiv double %251, %253
store double %254, double* @brcmp1.h
%255 = load double, double* @brcmp1.h
%256 = fadd double 1.000000e+00, %255
%257 = fdiv double 1.000000e+00, %256
store double %257, double* @brcmp1.x0
%258 = load double, double* @brcmp1.h
%259 = load double, double* @brcmp1.h
%260 = fadd double 1.000000e+00, %259
%261 = fdiv double %258, %260
store double %261, double* @brcmp1.y0
%262 = load double*, double** %8
%263 = load double, double* %262
%264 = load double*, double** %9
%265 = load double, double* %264
%266 = fadd double %263, %265
%267 = load double*, double** %11
%268 = load double, double* %267
%269 = load double*, double** %9
%270 = load double, double* %269
%271 = fneg double %270
%272 = call double @llvm.fmuladd.f64(double %266, double %268, double %271)
store double %272, double* @brcmp1.lambda
br label %273
273: ; preds = %249, %225
%274 = load double, double* @brcmp1.lambda
%275 = load double*, double** %8
%276 = load double, double* %275
%277 = fdiv double %274, %276
%278 = fneg double %277
store double %278, double* @brcmp1.e
%279 = load double, double* @brcmp1.e
%280 = call double @llvm.fabs.f64(double %279)
%281 = fcmp ogt double %280, 6.000000e-01
br i1 %281, label %282, label %283
282: ; preds = %273
br label %285
283: ; preds = %273
%284 = call double @rlog1(double* @brcmp1.e)
store double %284, double* @brcmp1.u
br label %293
285: ; preds = %282
%286 = load double, double* @brcmp1.e
%287 = load double*, double** %10
%288 = load double, double* %287
%289 = load double, double* @brcmp1.x0
%290 = fdiv double %288, %289
%291 = call double @log(double %290)
%292 = fsub double %286, %291
store double %292, double* @brcmp1.u
br label %293
293: ; preds = %285, %283
%294 = load double, double* @brcmp1.lambda
%295 = load double*, double** %9
%296 = load double, double* %295
%297 = fdiv double %294, %296
store double %297, double* @brcmp1.e
%298 = load double, double* @brcmp1.e
%299 = call double @llvm.fabs.f64(double %298)
%300 = fcmp ogt double %299, 6.000000e-01
br i1 %300, label %301, label %302
301: ; preds = %293
br label %304
302: ; preds = %293
%303 = call double @rlog1(double* @brcmp1.e)
store double %303, double* @brcmp1.v
br label %312
304: ; preds = %301
%305 = load double, double* @brcmp1.e
%306 = load double*, double** %11
%307 = load double, double* %306
%308 = load double, double* @brcmp1.y0
%309 = fdiv double %307, %308
%310 = call double @log(double %309)
%311 = fsub double %305, %310
store double %311, double* @brcmp1.v
br label %312
312: ; preds = %304, %302
%313 = load double*, double** %8
%314 = load double, double* %313
%315 = load double, double* @brcmp1.u
%316 = load double*, double** %9
%317 = load double, double* %316
%318 = load double, double* @brcmp1.v
%319 = fmul double %317, %318
%320 = call double @llvm.fmuladd.f64(double %314, double %315, double %319)
%321 = fneg double %320
store double %321, double* @brcmp1.T4
%322 = load i32*, i32** %7
%323 = call double @esum(i32* %322, double* @brcmp1.T4)
store double %323, double* @brcmp1.z
%324 = load double, double* @brcmp1.Const
%325 = load double*, double** %9
%326 = load double, double* %325
%327 = load double, double* @brcmp1.x0
%328 = fmul double %326, %327
%329 = call double @sqrt(double %328)
%330 = fmul double %324, %329
%331 = load double, double* @brcmp1.z
%332 = fmul double %330, %331
%333 = load double*, double** %8
%334 = load double*, double** %9
%335 = call i32 @bcorr(double* %333, double* %334)
%336 = sub nsw i32 0, %335
%337 = call double @exp(i32 %336)
%338 = fmul double %332, %337
store double %338, double* @brcmp1.brcmp1
%339 = load double, double* @brcmp1.brcmp1
store double %339, double* %6
br label %340
340: ; preds = %312, %206, %195, %118, %93, %65
%341 = load double, double* %6
ret double %341
}

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 %176
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 %247
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 %167
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 %249, 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 %247
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 i64
%99 = trunc i64 %98 to i32
%100 = icmp slt i32 %99, 1
br i1 %100, label %140, label %101
101: ; preds = %94
%102 = load double, double* @brcmp1.a0
%103 = and i32 %99, 1
%104 = icmp eq i32 %99, 1
br i1 %104, label %121, label %105
105: ; preds = %101
%106 = and i32 %99, -2
br label %107
107: ; preds = %107, %105
%108 = phi double [ %96, %105 ], [ %115, %107 ]
%109 = phi double [ 1.000000e+00, %105 ], [ %118, %107 ]
%110 = phi i32 [ 0, %105 ], [ %119, %107 ]
%111 = fadd double %108, -1.000000e+00
%112 = fadd double %111, %102
%113 = fdiv double %111, %112
%114 = fmul double %109, %113
%115 = fadd double %111, -1.000000e+00
%116 = fadd double %115, %102
%117 = fdiv double %115, %116
%118 = fmul double %114, %117
%119 = add i32 %110, 2
%120 = icmp eq i32 %119, %106
br i1 %120, label %121, label %107
121: ; preds = %107, %101
%122 = phi double [ undef, %101 ], [ %115, %107 ]
%123 = phi double [ undef, %101 ], [ %118, %107 ]
%124 = phi double [ %96, %101 ], [ %115, %107 ]
%125 = phi double [ 1.000000e+00, %101 ], [ %118, %107 ]
%126 = icmp eq i32 %103, 0
br i1 %126, label %132, label %127
127: ; preds = %121
%128 = fadd double %124, -1.000000e+00
%129 = fadd double %128, %102
%130 = fdiv double %128, %129
%131 = fmul double %125, %130
br label %132
132: ; preds = %121, %127
%133 = phi double [ %122, %121 ], [ %128, %127 ]
%134 = phi double [ %123, %121 ], [ %131, %127 ]
store double %133, double* @brcmp1.b0
%135 = tail call double @log(double %134)
%136 = load double, double* @brcmp1.u
%137 = fadd double %135, %136
store double %137, double* @brcmp1.u
%138 = load double, double* @brcmp1.b0
%139 = fadd double %138, -1.000000e+00
br label %140
140: ; preds = %94, %132
%141 = phi double [ %97, %94 ], [ %139, %132 ]
%142 = phi double [ %95, %94 ], [ %137, %132 ]
%143 = load double, double* @brcmp1.z
%144 = fsub double %143, %142
store double %144, double* @brcmp1.z
store double %141, double* @brcmp1.b0
%145 = load double, double* @brcmp1.a0
%146 = fadd double %141, %145
store double %146, double* @brcmp1.apb
%147 = fcmp ogt double %146, 1.000000e+00
br i1 %147, label %151, label %148
148: ; preds = %140
%149 = tail call double @gam1(double* @brcmp1.apb)
%150 = fadd double %149, 1.000000e+00
br label %157
151: ; preds = %140
%152 = fadd double %146, -1.000000e+00
store double %152, double* @brcmp1.u
%153 = tail call double @gam1(double* @brcmp1.u)
%154 = fadd double %153, 1.000000e+00
%155 = load double, double* @brcmp1.apb
%156 = fdiv double %154, %155
br label %157
157: ; preds = %151, %148
%158 = phi double [ %150, %148 ], [ %156, %151 ]
store double %158, double* @brcmp1.t
%159 = load double, double* @brcmp1.a0
%160 = tail call double @esum(i32* %0, double* @brcmp1.z)
%161 = fmul double %159, %160
%162 = tail call double @gam1(double* @brcmp1.b0)
%163 = fadd double %162, 1.000000e+00
%164 = fmul double %161, %163
%165 = load double, double* @brcmp1.t
%166 = fdiv double %164, %165
br label %247
167: ; preds = %45
%168 = tail call double @gamln1(double* @brcmp1.a0)
%169 = tail call double @algdiv(double* @brcmp1.a0, double* @brcmp1.b0)
%170 = fadd double %168, %169
store double %170, double* @brcmp1.u
%171 = load double, double* @brcmp1.z
%172 = fsub double %171, %170
store double %172, double* @brcmp1.T3
%173 = load double, double* @brcmp1.a0
%174 = tail call double @esum(i32* %0, double* @brcmp1.T3)
%175 = fmul double %173, %174
br label %247
176: ; preds = %5
%177 = load double, double* %1
%178 = load double, double* %2
%179 = fcmp ogt double %177, %178
br i1 %179, label %187, label %180
180: ; preds = %176
%181 = fdiv double %177, %178
%182 = fadd double %181, 1.000000e+00
%183 = fdiv double %181, %182
store double %183, double* @brcmp1.x0
%184 = fdiv double 1.000000e+00, %182
store double %184, double* @brcmp1.y0
%185 = fadd double %177, %178
%186 = fneg double %185
br label %194
187: ; preds = %176
%188 = fdiv double %178, %177
%189 = fadd double %188, 1.000000e+00
%190 = fdiv double 1.000000e+00, %189
store double %190, double* @brcmp1.x0
%191 = fdiv double %188, %189
store double %191, double* @brcmp1.y0
%192 = fadd double %177, %178
%193 = fneg double %178
br label %194
194: ; preds = %187, %180
%195 = phi double [ %193, %187 ], [ %177, %180 ]
%196 = phi double* [ %4, %187 ], [ %3, %180 ]
%197 = phi double [ %192, %187 ], [ %186, %180 ]
%198 = phi double [ %190, %187 ], [ %183, %180 ]
%199 = load double, double* %196
%200 = tail call double @llvm.fmuladd.f64(double %197, double %199, double %195)
store double %200, double* @brcmp1.lambda
%201 = fneg double %200
%202 = fdiv double %201, %177
store double %202, double* @brcmp1.e
%203 = tail call double @llvm.fabs.f64(double %202)
%204 = fcmp ogt double %203, 6.000000e-01
br i1 %204, label %207, label %205
205: ; preds = %194
%206 = tail call double @rlog1(double* @brcmp1.e)
br label %212
207: ; preds = %194
%208 = load double, double* %3
%209 = fdiv double %208, %198
%210 = tail call double @log(double %209)
%211 = fsub double %202, %210
br label %212
212: ; preds = %207, %205
%213 = phi double [ %206, %205 ], [ %211, %207 ]
store double %213, double* @brcmp1.u
%214 = load double, double* @brcmp1.lambda
%215 = load double, double* %2
%216 = fdiv double %214, %215
store double %216, double* @brcmp1.e
%217 = tail call double @llvm.fabs.f64(double %216)
%218 = fcmp ogt double %217, 6.000000e-01
br i1 %218, label %221, label %219
219: ; preds = %212
%220 = tail call double @rlog1(double* @brcmp1.e)
br label %227
221: ; preds = %212
%222 = load double, double* %4
%223 = load double, double* @brcmp1.y0
%224 = fdiv double %222, %223
%225 = tail call double @log(double %224)
%226 = fsub double %216, %225
br label %227
227: ; preds = %221, %219
%228 = phi double [ %220, %219 ], [ %226, %221 ]
%229 = load double, double* %1
%230 = load double, double* @brcmp1.u
%231 = load double, double* %2
%232 = fmul double %228, %231
%233 = tail call double @llvm.fmuladd.f64(double %229, double %230, double %232)
%234 = fneg double %233
store double %234, double* @brcmp1.T4
%235 = tail call double @esum(i32* %0, double* @brcmp1.T4)
store double %235, double* @brcmp1.z
%236 = load double, double* %2
%237 = load double, double* @brcmp1.x0
%238 = fmul double %236, %237
%239 = tail call double @sqrt(double %238)
%240 = fmul double %239, 0x3FD9884533D43656
%241 = load double, double* @brcmp1.z
%242 = fmul double %240, %241
%243 = tail call i32 @bcorr(double* %1, double* %2)
%244 = sub nsw i32 0, %243
%245 = tail call double @exp(i32 %244)
%246 = fmul double %242, %245
br label %247
247: ; preds = %39, %72, %157, %167, %227
%248 = phi double [ %246, %227 ], [ %175, %167 ], [ %166, %157 ], [ %93, %72 ], [ %44, %39 ]
store double %248, double* @brcmp1.brcmp1
br label %249
249: ; preds = %247, %52
%250 = phi double [ %53, %52 ], [ %248, %247 ]
ret double %250
}