#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 = (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:
  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    1c3e <brcomp+0xa7e>
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    18c6 <brcomp+0x706>
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    1c3e <brcomp+0xa7e>
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    1846 <brcomp+0x686>
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    1c3e <brcomp+0xa7e>
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    1c3e <brcomp+0xa7e>
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,%rax
mov    %eax,0x2adc(%rip)
cmpl   $0x1,0x2ad5(%rip)
jge    164a <brcomp+0x48a>
jmp    16ef <brcomp+0x52f>
movsd  0x9be(%rip),%xmm0
movsd  %xmm0,0x2a56(%rip)
movl   $0x1,0x2aac(%rip)
mov    0x2aa6(%rip),%eax
cmp    0x2aa4(%rip),%eax
jg     16d2 <brcomp+0x512>
movsd  0x2a2a(%rip),%xmm0
movsd  0x98a(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2a16(%rip)
movsd  0x2a0e(%rip),%xmm0
movsd  0x29f6(%rip),%xmm1
addsd  0x29fe(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  0x29fa(%rip),%xmm0
movsd  %xmm0,0x29f2(%rip)
mov    0x2a4c(%rip),%eax
add    $0x1,%eax
mov    %eax,0x2a43(%rip)
jmp    1664 <brcomp+0x4a4>
movsd  0x29d6(%rip),%xmm0
call   10a0 <log@plt>
addsd  0x2a01(%rip),%xmm0
movsd  %xmm0,0x29f9(%rip)
movsd  0x29f1(%rip),%xmm1
movsd  0x2a09(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x29fd(%rip)
movsd  0x2995(%rip),%xmm0
movsd  0x8f5(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2981(%rip)
movsd  0x2969(%rip),%xmm0
addsd  0x2971(%rip),%xmm0
movsd  %xmm0,0x2961(%rip)
movsd  0x2959(%rip),%xmm0
movsd  0x8c1(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    175e <brcomp+0x59e>
jmp    1786 <brcomp+0x5c6>
lea    0x293b(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x89b(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x295f(%rip)
jmp    17d5 <brcomp+0x615>
movsd  0x290a(%rip),%xmm0
addsd  0x2912(%rip),%xmm0
movsd  0x872(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x293e(%rip)
lea    0x2937(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm1
movsd  0x84f(%rip),%xmm0
addsd  %xmm1,%xmm0
divsd  0x28d3(%rip),%xmm0
movsd  %xmm0,0x290b(%rip)
movsd  0x28bb(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0x291e(%rip),%xmm0
call   1090 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x40(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
lea    0x28a1(%rip),%rdi
call   1030 <gam1@plt>
movaps %xmm0,%xmm2
movsd  -0x38(%rbp),%xmm0
movsd  0x7f4(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  0x28b4(%rip),%xmm0
movsd  %xmm0,0x285c(%rip)
movsd  0x2854(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3e <brcomp+0xa7e>
lea    0x284b(%rip),%rdi
call   10e0 <gamln1@plt>
movsd  %xmm0,-0x50(%rbp)
lea    0x283a(%rip),%rdi
lea    0x2843(%rip),%rsi
call   1080 <algdiv@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x286a(%rip)
movsd  0x2812(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0x2875(%rip),%xmm0
subsd  0x284d(%rip),%xmm0
call   1090 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x27dc(%rip)
movsd  0x27d4(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1c3e <brcomp+0xa7e>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
ucomisd (%rax),%xmm0
jbe    18e1 <brcomp+0x721>
jmp    198d <brcomp+0x7cd>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x27c7(%rip)
movsd  0x27bf(%rip),%xmm0
movsd  0x707(%rip),%xmm1
addsd  0x27af(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x27db(%rip)
movsd  0x6eb(%rip),%xmm1
addsd  0x2793(%rip),%xmm1
movsd  0x6db(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x27bf(%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,0x2740(%rip)
jmp    1a34 <brcomp+0x874>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x271b(%rip)
movsd  0x663(%rip),%xmm1
addsd  0x270b(%rip),%xmm1
movsd  0x653(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x272f(%rip)
movsd  0x26ef(%rip),%xmm0
movsd  0x637(%rip),%xmm1
addsd  0x26df(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2713(%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,0x2694(%rip)
movsd  0x268c(%rip),%xmm0
mov    -0x10(%rbp),%rax
divsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,0x2655(%rip)
movsd  0x264d(%rip),%xmm0
movaps 0x5ae(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x59a(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1a8d <brcomp+0x8cd>
jmp    1aa6 <brcomp+0x8e6>
lea    0x2624(%rip),%rdi
call   1050 <rlog1@plt>
movsd  %xmm0,0x2647(%rip)
jmp    1adc <brcomp+0x91c>
movsd  0x260a(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x2635(%rip),%xmm0
call   10a0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x260c(%rip)
movsd  0x25e4(%rip),%xmm0
mov    -0x18(%rbp),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,0x25c4(%rip)
movsd  0x25bc(%rip),%xmm0
movaps 0x51d(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x509(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1b1e <brcomp+0x95e>
jmp    1b37 <brcomp+0x977>
lea    0x2593(%rip),%rdi
call   1050 <rlog1@plt>
movsd  %xmm0,0x25be(%rip)
jmp    1b6d <brcomp+0x9ad>
movsd  0x2579(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
divsd  0x25ac(%rip),%xmm0
call   10a0 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2583(%rip)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0x256b(%rip),%xmm2
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  0x2563(%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,0x254f(%rip)
movsd  0x24bf(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  0x2522(%rip),%xmm0
call   10b0 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  0x2519(%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,0x245f(%rip)
movsd  0x2457(%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    1597 <brcomp+0x3d7>
mov    %rcx,%r12
movsd  (%rcx),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    1201 <brcomp+0x41>
jnp    1597 <brcomp+0x3d7>
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    158f <brcomp+0x3cf>
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    1563 <brcomp+0x3a3>
movapd %xmm2,%xmm0
call   1060 <fifdmax1@plt>
movsd  %xmm0,0x2bac(%rip)
ucomisd 0xb04(%rip),%xmm0
jae    15a3 <brcomp+0x3e3>
ucomisd 0xafe(%rip),%xmm0
jbe    15f5 <brcomp+0x435>
lea    0x2b79(%rip),%rdi
call   10e0 <gamln1@plt>
movsd  %xmm0,0x2ba4(%rip)
movsd  0x2b74(%rip),%xmm1
movsd  0xaf4(%rip),%xmm2
addsd  %xmm1,%xmm2
cvttsd2si %xmm2,%rax
test   %eax,%eax
jle    1704 <brcomp+0x544>
movsd  0x2b43(%rip),%xmm5
cmp    $0x1,%eax
jne    1671 <brcomp+0x4b1>
movsd  0xaaa(%rip),%xmm0
jmp    16bb <brcomp+0x4fb>
mov    %r14,%rdi
mov    %r15,%rsi
call   10c0 <betaln@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %rax,%xmm1
movsd  0x2b62(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x2b56(%rip)
call   1090 <exp@plt>
movsd  %xmm0,0x2af1(%rip)
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
lea    0x2ae6(%rip),%rbx
mov    %rbx,%rdi
call   10e0 <gamln1@plt>
movsd  %xmm0,(%rsp)
lea    0x2ae2(%rip),%rsi
mov    %rbx,%rdi
call   1080 <algdiv@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,0x2af5(%rip)
movsd  0x2ab5(%rip),%xmm1
movsd  %xmm1,(%rsp)
movsd  0x2af8(%rip),%xmm1
subsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
jmp    1469 <brcomp+0x2a9>
movsd  0x2ae3(%rip),%xmm0
call   1090 <exp@plt>
movsd  %xmm0,0x2a7e(%rip)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1616 <brcomp+0x456>
jnp    1597 <brcomp+0x3d7>
movsd  (%r14),%xmm0
addsd  (%r15),%xmm0
movsd  %xmm0,0x2a70(%rip)
ucomisd 0x9d8(%rip),%xmm0
jbe    17ce <brcomp+0x60e>
movsd  (%r14),%xmm0
addsd  (%r15),%xmm0
addsd  0x9e0(%rip),%xmm0
movsd  %xmm0,0x2a78(%rip)
lea    0x2a71(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x9a4(%rip),%xmm0
divsd  0x2a2c(%rip),%xmm0
jmp    17e2 <brcomp+0x622>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movsd  0x98a(%rip),%xmm0
movsd  0x9a2(%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    1686 <brcomp+0x4c6>
test   $0x1,%al
je     16d7 <brcomp+0x517>
addsd  0x961(%rip),%xmm1
addsd  %xmm1,%xmm5
movapd %xmm1,%xmm2
divsd  %xmm5,%xmm2
mulsd  %xmm2,%xmm0
movsd  %xmm1,0x29c1(%rip)
call   10a0 <log@plt>
addsd  0x29dc(%rip),%xmm0
movsd  %xmm0,0x29d4(%rip)
movsd  0x29a4(%rip),%xmm2
addsd  0x924(%rip),%xmm2
movsd  0x29d4(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x29c8(%rip)
movsd  %xmm2,0x2980(%rip)
addsd  0x2968(%rip),%xmm2
movsd  %xmm2,0x2968(%rip)
ucomisd 0x8d0(%rip),%xmm2
jbe    1768 <brcomp+0x5a8>
addsd  0x8e6(%rip),%xmm2
movsd  %xmm2,0x297e(%rip)
lea    0x2977(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x8aa(%rip),%xmm0
divsd  0x2932(%rip),%xmm0
jmp    177c <brcomp+0x5bc>
lea    0x2929(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x88c(%rip),%xmm0
movsd  %xmm0,0x293c(%rip)
movsd  0x2904(%rip),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x2947(%rip),%xmm0
call   1090 <exp@plt>
mulsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
lea    0x28f1(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x84c(%rip),%xmm0
mulsd  (%rsp),%xmm0
divsd  0x28f7(%rip),%xmm0
jmp    158f <brcomp+0x3cf>
lea    0x28c3(%rip),%rdi
call   1030 <gam1@plt>
addsd  0x826(%rip),%xmm0
movsd  %xmm0,0x28f6(%rip)
mov    %r14,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x80a(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
mov    %r15,%rdi
call   1030 <gam1@plt>
movapd %xmm0,%xmm1
movsd  0x7ed(%rip),%xmm3
addsd  %xmm3,%xmm1
mulsd  (%rsp),%xmm1
movsd  0x2864(%rip),%xmm0
unpcklpd %xmm0,%xmm1
movsd  0x28a8(%rip),%xmm2
movhpd 0x2860(%rip),%xmm2
divpd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mulsd  0x2838(%rip),%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm0
jmp    158f <brcomp+0x3cf>

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

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