#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 *);
double bcorr(double *, double *);
scalar_t__ betaln(double *, double *);
double exp(double);
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 brcomp(double *a, double *b, double *x, double *y) {
  static double Const = .398942280401433e0;
  static double brcomp, a0, apb, b0, c, e, h, lambda, lnx, lny, t, u, v, x0, y0,
      z;
  static int i, n;
  static double T1, T2;
  brcomp = 0.0e0;
  if (*x == 0.0e0 || *y == 0.0e0)
    return brcomp;
  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);
  brcomp = exp(z);
  return brcomp;
S40:
  b0 = fifdmax1(*a, *b);
  if (b0 >= 8.0e0)
    goto S120;
  if (b0 > 1.0e0)
    goto S70;
  brcomp = exp(z);
  if (brcomp == 0.0e0)
    return brcomp;
  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;
  brcomp = brcomp * (a0 * c) / (1.0e0 + a0 / b0);
  return brcomp;
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:
  brcomp = a0 * exp(z) * (1.0e0 + gam1(&b0)) / t;
  return brcomp;
S120:
  u = gamln1(&a0) + algdiv(&a0, &b0);
  brcomp = a0 * exp(z - u);
  return brcomp;
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:
  z = exp(-(*a * u + *b * v));
  brcomp = Const * sqrt(*b * x0) * z * exp(-bcorr(a, b));
  return brcomp;
}


x86 assembly (clang, -O0):
<brcomp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,0x2ead(%rip)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1203 <brcomp+0x43>
jp     1203 <brcomp+0x43>
jmp    121e <brcomp+0x5e>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1230 <brcomp+0x70>
jp     1230 <brcomp+0x70>
movsd  0x2e6a(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3d <brcomp+0xa7d>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
call   1070 <fifdmin1@plt>
movsd  %xmm0,0x2e4b(%rip)
movsd  0x2e43(%rip),%xmm0
movsd  0xda3(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     126c <brcomp+0xac>
jmp    18c5 <brcomp+0x705>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xd8c(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    128b <brcomp+0xcb>
jmp    12e0 <brcomp+0x120>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
call   10a0 <log@plt>
movsd  %xmm0,0x2e30(%rip)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2e51(%rip)
lea    0x2e4a(%rip),%rdi
call   1040 <alnrel@plt>
movsd  %xmm0,0x2dfd(%rip)
jmp    137e <brcomp+0x1be>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xd18(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12ff <brcomp+0x13f>
jmp    1354 <brcomp+0x194>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2dfa(%rip)
lea    0x2df3(%rip),%rdi
call   1040 <alnrel@plt>
movsd  %xmm0,0x2d96(%rip)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
call   10a0 <log@plt>
movsd  %xmm0,0x2d89(%rip)
jmp    137e <brcomp+0x1be>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
call   10a0 <log@plt>
movsd  %xmm0,0x2d67(%rip)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
call   10a0 <log@plt>
movsd  %xmm0,0x2d5a(%rip)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0x2d42(%rip),%xmm2
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  0x2d3a(%rip),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x2d5a(%rip)
movsd  0xc5a(%rip),%xmm0
ucomisd 0x2cda(%rip),%xmm0
jbe    13c9 <brcomp+0x209>
jmp    1416 <brcomp+0x256>
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   10c0 <betaln@plt>
cvtsi2sd %rax,%xmm1
movsd  0x2d25(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2d19(%rip)
movsd  0x2d11(%rip),%xmm0
call   1090 <exp@plt>
movsd  %xmm0,0x2c8c(%rip)
movsd  0x2c84(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3d <brcomp+0xa7d>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
call   1060 <fifdmax1@plt>
movsd  %xmm0,0x2c75(%rip)
movsd  0x2c6d(%rip),%xmm0
movsd  0xbbd(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1452 <brcomp+0x292>
jmp    1845 <brcomp+0x685>
movsd  0x2c4e(%rip),%xmm0
movsd  0xbae(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1471 <brcomp+0x2b1>
jmp    1605 <brcomp+0x445>
movsd  0x2c8f(%rip),%xmm0
call   1090 <exp@plt>
movsd  %xmm0,0x2c0a(%rip)
movsd  0x2c02(%rip),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    14b3 <brcomp+0x2f3>
jp     14b3 <brcomp+0x2f3>
movsd  0x2be7(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3d <brcomp+0xa7d>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,0x2bd5(%rip)
movsd  0x2bcd(%rip),%xmm0
movsd  0xb35(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    14ea <brcomp+0x32a>
jmp    1512 <brcomp+0x352>
lea    0x2baf(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0xb0f(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x2bfb(%rip)
jmp    1561 <brcomp+0x3a1>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
addsd  (%rax),%xmm0
movsd  0xae6(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2bb2(%rip)
lea    0x2bab(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0xac3(%rip),%xmm0
addsd  %xmm1,%xmm0
divsd  0x2b47(%rip),%xmm0
movsd  %xmm0,0x2ba7(%rip)
mov    -0x10(%rbp),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0xa9b(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm2
movsd  -0x30(%rbp),%xmm0
movsd  0xa79(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  0x2b61(%rip),%xmm0
movsd  %xmm0,0x2b01(%rip)
movsd  0x2ad9(%rip),%xmm0
movsd  0x2ad9(%rip),%xmm1
mulsd  0x2ae9(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0x2ac5(%rip),%xmm2
divsd  0x2acd(%rip),%xmm2
movsd  0xa2d(%rip),%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2a9d(%rip)
movsd  0x2a95(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3d <brcomp+0xa7d>
lea    0x2a8c(%rip),%rdi
call   10e0 <gamln1@plt>
movsd  %xmm0,0x2acf(%rip)
movsd  0x2a87(%rip),%xmm0
movsd  0x9e7(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,0x2add(%rip)
cmpl   $0x1,0x2ad6(%rip)
jge    1649 <brcomp+0x489>
jmp    16ee <brcomp+0x52e>
movsd  0x9bf(%rip),%xmm0
movsd  %xmm0,0x2a57(%rip)
movl   $0x1,0x2aad(%rip)
mov    0x2aa7(%rip),%eax
cmp    0x2aa5(%rip),%eax
jg     16d1 <brcomp+0x511>
movsd  0x2a2b(%rip),%xmm0
movsd  0x98b(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2a17(%rip)
movsd  0x2a0f(%rip),%xmm0
movsd  0x29f7(%rip),%xmm1
addsd  0x29ff(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  0x29fb(%rip),%xmm0
movsd  %xmm0,0x29f3(%rip)
mov    0x2a4d(%rip),%eax
add    $0x1,%eax
mov    %eax,0x2a44(%rip)
jmp    1663 <brcomp+0x4a3>
movsd  0x29d7(%rip),%xmm0
call   10a0 <log@plt>
addsd  0x2a02(%rip),%xmm0
movsd  %xmm0,0x29fa(%rip)
movsd  0x29f2(%rip),%xmm1
movsd  0x2a0a(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x29fe(%rip)
movsd  0x2996(%rip),%xmm0
movsd  0x8f6(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2982(%rip)
movsd  0x296a(%rip),%xmm0
addsd  0x2972(%rip),%xmm0
movsd  %xmm0,0x2962(%rip)
movsd  0x295a(%rip),%xmm0
movsd  0x8c2(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    175d <brcomp+0x59d>
jmp    1785 <brcomp+0x5c5>
lea    0x293c(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x89c(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x2960(%rip)
jmp    17d4 <brcomp+0x614>
movsd  0x290b(%rip),%xmm0
addsd  0x2913(%rip),%xmm0
movsd  0x873(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x293f(%rip)
lea    0x2938(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x850(%rip),%xmm0
addsd  %xmm1,%xmm0
divsd  0x28d4(%rip),%xmm0
movsd  %xmm0,0x290c(%rip)
movsd  0x28bc(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0x291f(%rip),%xmm0
call   1090 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x40(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
lea    0x28a2(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm2
movsd  -0x38(%rbp),%xmm0
movsd  0x7f5(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  0x28b5(%rip),%xmm0
movsd  %xmm0,0x285d(%rip)
movsd  0x2855(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3d <brcomp+0xa7d>
lea    0x284c(%rip),%rdi
call   10e0 <gamln1@plt>
movsd  %xmm0,-0x50(%rbp)
lea    0x283b(%rip),%rdi
lea    0x2844(%rip),%rsi
call   1080 <algdiv@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x286b(%rip)
movsd  0x2813(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0x2876(%rip),%xmm0
subsd  0x284e(%rip),%xmm0
call   1090 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x27dd(%rip)
movsd  0x27d5(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3d <brcomp+0xa7d>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
ucomisd (%rax),%xmm0
jbe    18e0 <brcomp+0x720>
jmp    198c <brcomp+0x7cc>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x27c8(%rip)
movsd  0x27c0(%rip),%xmm0
movsd  0x708(%rip),%xmm1
addsd  0x27b0(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x27dc(%rip)
movsd  0x6ec(%rip),%xmm1
addsd  0x2794(%rip),%xmm1
movsd  0x6dc(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x27c0(%rip)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
addsd  (%rax),%xmm0
mov    -0x20(%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,0x2741(%rip)
jmp    1a33 <brcomp+0x873>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x271c(%rip)
movsd  0x664(%rip),%xmm1
addsd  0x270c(%rip),%xmm1
movsd  0x654(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2730(%rip)
movsd  0x26f0(%rip),%xmm0
movsd  0x638(%rip),%xmm1
addsd  0x26e0(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2714(%rip)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
addsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x18(%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,0x2695(%rip)
movsd  0x268d(%rip),%xmm0
mov    -0x10(%rbp),%rax
divsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2656(%rip)
movsd  0x264e(%rip),%xmm0
movaps 0x5af(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x59b(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1a8c <brcomp+0x8cc>
jmp    1aa5 <brcomp+0x8e5>
lea    0x2625(%rip),%rdi
call   1050 <rlog1@plt>
movsd  %xmm0,0x2648(%rip)
jmp    1adb <brcomp+0x91b>
movsd  0x260b(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x2636(%rip),%xmm0
call   10a0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x260d(%rip)
movsd  0x25e5(%rip),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x25c5(%rip)
movsd  0x25bd(%rip),%xmm0
movaps 0x51e(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x50a(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1b1d <brcomp+0x95d>
jmp    1b36 <brcomp+0x976>
lea    0x2594(%rip),%rdi
call   1050 <rlog1@plt>
movsd  %xmm0,0x25bf(%rip)
jmp    1b6c <brcomp+0x9ac>
movsd  0x257a(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x25ad(%rip),%xmm0
call   10a0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2584(%rip)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0x256c(%rip),%xmm2
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  0x2564(%rip),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1090 <exp@plt>
movsd  %xmm0,0x2550(%rip)
movsd  0x24c0(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  0x2523(%rip),%xmm0
call   10b0 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  0x251a(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   10d0 <bcorr@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1090 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x2460(%rip)
movsd  0x2458(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<brcomp>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x18,%rsp
movq   $0x0,0x2eb2(%rip)
movsd  (%rdx),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    11ea <brcomp+0x2a>
jnp    1592 <brcomp+0x3d2>
mov    %rcx,%r12
movsd  (%rcx),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    1201 <brcomp+0x41>
jnp    1592 <brcomp+0x3d2>
mov    %rdx,%rbx
mov    %rsi,%r15
mov    %rdi,%r14
movsd  (%rdi),%xmm0
movsd  (%rsi),%xmm1
call   1070 <fifdmin1@plt>
movsd  %xmm0,0x2e71(%rip)
ucomisd 0xdd9(%rip),%xmm0
jae    1255 <brcomp+0x95>
movsd  (%rbx),%xmm0
ucomisd 0xdeb(%rip),%xmm0
jbe    12b2 <brcomp+0xf2>
movsd  (%r12),%xmm1
ucomisd 0xddb(%rip),%xmm1
jbe    1478 <brcomp+0x2b8>
call   10a0 <log@plt>
jmp    1494 <brcomp+0x2d4>
movsd  (%r14),%xmm0
movsd  (%r15),%xmm4
ucomisd %xmm4,%xmm0
jbe    12e4 <brcomp+0x124>
movapd %xmm4,%xmm2
divsd  %xmm0,%xmm2
movsd  0xd93(%rip),%xmm1
movapd %xmm2,%xmm3
addsd  %xmm1,%xmm3
divsd  %xmm3,%xmm1
movsd  %xmm1,0x2e47(%rip)
divsd  %xmm3,%xmm2
movsd  %xmm2,0x2e43(%rip)
movapd %xmm0,%xmm3
addsd  %xmm4,%xmm3
xorpd  0xd8b(%rip),%xmm4
movapd %xmm4,%xmm2
mov    %r12,%rax
movapd %xmm3,%xmm4
jmp    1327 <brcomp+0x167>
call   10a0 <log@plt>
movsd  %xmm0,0x2df9(%rip)
movsd  (%rbx),%xmm0
xorpd  0xd65(%rip),%xmm0
movlpd %xmm0,0x2e15(%rip)
lea    0x2e0e(%rip),%rdi
call   1040 <alnrel@plt>
jmp    14a7 <brcomp+0x2e7>
movapd %xmm0,%xmm1
divsd  %xmm4,%xmm1
movsd  0xd14(%rip),%xmm2
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
divsd  %xmm3,%xmm1
movsd  %xmm1,0x2dc8(%rip)
divsd  %xmm3,%xmm2
movsd  %xmm2,0x2dc4(%rip)
addsd  %xmm0,%xmm4
xorpd  0xd10(%rip),%xmm4
movapd %xmm0,%xmm2
mov    %rbx,%rax
mulsd  (%rax),%xmm4
addsd  %xmm2,%xmm4
movsd  %xmm4,0x2d79(%rip)
xorpd  0xcf1(%rip),%xmm4
divsd  %xmm0,%xmm4
movsd  %xmm4,0x2d5d(%rip)
movapd 0xced(%rip),%xmm0
andpd  %xmm4,%xmm0
ucomisd 0xcb1(%rip),%xmm0
jbe    137e <brcomp+0x1be>
movsd  (%rbx),%xmm0
divsd  %xmm1,%xmm0
movapd %xmm4,(%rsp)
call   10a0 <log@plt>
movapd (%rsp),%xmm1
subsd  %xmm0,%xmm1
jmp    138e <brcomp+0x1ce>
lea    0x2d23(%rip),%rdi
call   1050 <rlog1@plt>
movapd %xmm0,%xmm1
movsd  %xmm1,0x2d32(%rip)
movsd  0x2d12(%rip),%xmm1
divsd  (%r15),%xmm1
movsd  %xmm1,0x2cfd(%rip)
movapd 0xc8d(%rip),%xmm0
andpd  %xmm1,%xmm0
ucomisd 0xc51(%rip),%xmm0
jbe    13e4 <brcomp+0x224>
movsd  (%r12),%xmm0
divsd  0x2d09(%rip),%xmm0
movapd %xmm1,(%rsp)
call   10a0 <log@plt>
movapd (%rsp),%xmm1
subsd  %xmm0,%xmm1
jmp    13f4 <brcomp+0x234>
lea    0x2cbd(%rip),%rdi
call   1050 <rlog1@plt>
movapd %xmm0,%xmm1
movsd  (%r14),%xmm0
mulsd  (%r15),%xmm1
mulsd  0x2cc2(%rip),%xmm0
addsd  %xmm1,%xmm0
xorpd  0xc1e(%rip),%xmm0
call   1090 <exp@plt>
movsd  %xmm0,0x2cc1(%rip)
movsd  (%r15),%xmm0
mulsd  0x2ca4(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     143c <brcomp+0x27c>
sqrtsd %xmm0,%xmm0
jmp    1441 <brcomp+0x281>
call   10b0 <sqrt@plt>
mulsd  0xbcf(%rip),%xmm0
mulsd  0x2c8f(%rip),%xmm0
movsd  %xmm0,(%rsp)
mov    %r14,%rdi
mov    %r15,%rsi
call   10d0 <bcorr@plt>
xorpd  0xbc7(%rip),%xmm0
call   1090 <exp@plt>
mulsd  (%rsp),%xmm0
jmp    158a <brcomp+0x3ca>
xorpd  0xbb0(%rip),%xmm1
movlpd %xmm1,0x2c68(%rip)
lea    0x2c61(%rip),%rdi
call   1040 <alnrel@plt>
movsd  %xmm0,0x2c1c(%rip)
movsd  (%r12),%xmm0
call   10a0 <log@plt>
movsd  (%r14),%xmm2
movsd  (%r15),%xmm1
mulsd  %xmm1,%xmm0
movsd  0x2bfb(%rip),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,0x2c13(%rip)
movsd  0xb33(%rip),%xmm0
ucomisd 0x2bb3(%rip),%xmm0
jbe    155e <brcomp+0x39e>
movapd %xmm2,%xmm0
call   1060 <fifdmax1@plt>
movsd  %xmm0,0x2bb0(%rip)
ucomisd 0xb08(%rip),%xmm0
jae    159e <brcomp+0x3de>
ucomisd 0xb02(%rip),%xmm0
jbe    15f0 <brcomp+0x430>
lea    0x2b7d(%rip),%rdi
call   10e0 <gamln1@plt>
movsd  %xmm0,0x2ba8(%rip)
movsd  0x2b78(%rip),%xmm1
movsd  0xaf8(%rip),%xmm2
addsd  %xmm1,%xmm2
cvttsd2si %xmm2,%eax
test   %eax,%eax
jle    16ff <brcomp+0x53f>
movsd  0x2b48(%rip),%xmm5
cmp    $0x1,%eax
jne    166c <brcomp+0x4ac>
movsd  0xaaf(%rip),%xmm0
jmp    16b6 <brcomp+0x4f6>
mov    %r14,%rdi
mov    %r15,%rsi
call   10c0 <betaln@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %rax,%xmm1
movsd  0x2b67(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2b5b(%rip)
call   1090 <exp@plt>
movsd  %xmm0,0x2af6(%rip)
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
lea    0x2aeb(%rip),%rbx
mov    %rbx,%rdi
call   10e0 <gamln1@plt>
movsd  %xmm0,(%rsp)
lea    0x2ae7(%rip),%rsi
mov    %rbx,%rdi
call   1080 <algdiv@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,0x2afa(%rip)
movsd  0x2aba(%rip),%xmm1
movsd  %xmm1,(%rsp)
movsd  0x2afd(%rip),%xmm1
subsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
jmp    1469 <brcomp+0x2a9>
movsd  0x2ae8(%rip),%xmm0
call   1090 <exp@plt>
movsd  %xmm0,0x2a83(%rip)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1611 <brcomp+0x451>
jnp    1592 <brcomp+0x3d2>
movsd  (%r14),%xmm0
addsd  (%r15),%xmm0
movsd  %xmm0,0x2a75(%rip)
ucomisd 0x9dd(%rip),%xmm0
jbe    17c9 <brcomp+0x609>
movsd  (%r14),%xmm0
addsd  (%r15),%xmm0
addsd  0x9e5(%rip),%xmm0
movsd  %xmm0,0x2a7d(%rip)
lea    0x2a76(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x9a9(%rip),%xmm0
divsd  0x2a31(%rip),%xmm0
jmp    17dd <brcomp+0x61d>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movsd  0x98f(%rip),%xmm0
movsd  0x9a7(%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    1681 <brcomp+0x4c1>
test   $0x1,%al
je     16d2 <brcomp+0x512>
addsd  0x966(%rip),%xmm1
addsd  %xmm1,%xmm5
movapd %xmm1,%xmm2
divsd  %xmm5,%xmm2
mulsd  %xmm2,%xmm0
movsd  %xmm1,0x29c6(%rip)
call   10a0 <log@plt>
addsd  0x29e1(%rip),%xmm0
movsd  %xmm0,0x29d9(%rip)
movsd  0x29a9(%rip),%xmm2
addsd  0x929(%rip),%xmm2
movsd  0x29d9(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x29cd(%rip)
movsd  %xmm2,0x2985(%rip)
addsd  0x296d(%rip),%xmm2
movsd  %xmm2,0x296d(%rip)
ucomisd 0x8d5(%rip),%xmm2
jbe    1763 <brcomp+0x5a3>
addsd  0x8eb(%rip),%xmm2
movsd  %xmm2,0x2983(%rip)
lea    0x297c(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x8af(%rip),%xmm0
divsd  0x2937(%rip),%xmm0
jmp    1777 <brcomp+0x5b7>
lea    0x292e(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x891(%rip),%xmm0
movsd  %xmm0,0x2941(%rip)
movsd  0x2909(%rip),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x294c(%rip),%xmm0
call   1090 <exp@plt>
mulsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
lea    0x28f6(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x851(%rip),%xmm0
mulsd  (%rsp),%xmm0
divsd  0x28fc(%rip),%xmm0
jmp    158a <brcomp+0x3ca>
lea    0x28c8(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x82b(%rip),%xmm0
movsd  %xmm0,0x28fb(%rip)
mov    %r14,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x80f(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
mov    %r15,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x7f2(%rip),%xmm3
addsd  %xmm3,%xmm1
mulsd  (%rsp),%xmm1
movsd  0x2869(%rip),%xmm0
unpcklpd %xmm0,%xmm1
movsd  0x28ad(%rip),%xmm2
movhpd 0x2865(%rip),%xmm2
divpd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mulsd  0x283d(%rip),%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm0
jmp    158a <brcomp+0x3ca>

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

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