#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 bcorr(double, double);
double erfc1(int, double);
double exp(double);
double fabs(double);
double rlog1(double);
double sqrt(double);
double basym(double a, double b, double lambda, double eps) {
  static int num = 20;
  static double e0 = 1.12837916709551;
  static double e1 = .353553390593274;
  int i__1, i__2, i__3, i__4;
  double ret_val, r__1, r__2;
  double c__[21], d__[21], f, h__;
  int i__, j, m, n;
  double r__, s, t, u, w, z__, a0[21], b0[21], j0, j1, h2, r0, r1, t0, t1, w0,
      z0, z2, hn, zn;
  int im1, mm1, np1, imj, mmj;
  double sum, znm1, bsum, dsum;
  ret_val = 0.0;
  if (a >= b) {
    goto L10;
  }
  h__ = a / b;
  r0 = 1.0 / (h__ + 1.0);
  r1 = (b - a) / b;
  w0 = 1.0 / sqrt(a * (h__ + 1.0));
  goto L20;
L10:
  h__ = b / a;
  r0 = 1.0 / (h__ + 1.0);
  r1 = (b - a) / a;
  w0 = 1.0 / sqrt(b * (h__ + 1.0));
L20:
  r__1 = -lambda / a;
  r__2 = lambda / b;
  f = a * rlog1(r__1) + b * rlog1(r__2);
  t = exp(-f);
  if (t == 0.0) {
    return ret_val;
  }
  z0 = sqrt(f);
  z__ = z0 / e1 * 0.5;
  z2 = f + f;
  a0[0] = r1 * .66666666666666663;
  c__[0] = a0[0] * -0.5;
  d__[0] = -c__[0];
  j0 = 0.5 / e0 * erfc1(1, z0);
  j1 = e1;
  sum = j0 + d__[0] * w0 * j1;
  s = 1.0;
  h2 = h__ * h__;
  hn = 1.0;
  w = w0;
  znm1 = z__;
  zn = z2;
  i__1 = num;
  for (n = 2; n <= i__1; n += 2) {
    hn = h2 * hn;
    a0[n - 1] = r0 * 2.0 * (h__ * hn + 1.0) / (n + 2.0);
    np1 = n + 1;
    s += hn;
    a0[np1 - 1] = r1 * 2.0 * s / (n + 3.0);
    i__2 = np1;
    for (i__ = n; i__ <= i__2; ++i__) {
      r__ = (i__ + 1.0) * -0.5;
      b0[0] = r__ * a0[0];
      i__3 = i__;
      for (m = 2; m <= i__3; ++m) {
        bsum = 0.0;
        mm1 = m - 1;
        i__4 = mm1;
        for (j = 1; j <= i__4; ++j) {
          mmj = m - j;
          bsum += (j * r__ - mmj) * a0[j - 1] * b0[mmj - 1];
        }
        b0[m - 1] = r__ * a0[m - 1] + bsum / m;
      }
      c__[i__ - 1] = b0[i__ - 1] / (i__ + 1.0);
      dsum = 0.0;
      im1 = i__ - 1;
      i__3 = im1;
      for (j = 1; j <= i__3; ++j) {
        imj = i__ - j;
        dsum += d__[imj - 1] * c__[j - 1];
      }
      d__[i__ - 1] = -(dsum + c__[i__ - 1]);
    }
    j0 = e1 * znm1 + (n - 1.0) * j0;
    j1 = e1 * zn + n * j1;
    znm1 = z2 * znm1;
    zn = z2 * zn;
    w = w0 * w;
    t0 = d__[n - 1] * w * j0;
    w = w0 * w;
    t1 = d__[np1 - 1] * w * j1;
    sum += t0 + t1;
    if (fabs(t0) + fabs(t1) <= eps * sum) {
      goto L60;
    }
  }
L60:
  u = exp(-bcorr(a, b));
  ret_val = e0 * t * u * sum;
  return ret_val;
}


x86 assembly (clang, -O0):
<basym>:
push   %rbp
mov    %rsp,%rbp
sub    $0x420,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  %xmm3,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x10(%rbp),%xmm0
ucomisd -0x18(%rbp),%xmm0
jb     118c <basym+0x3c>
jmp    1213 <basym+0xc3>
movsd  -0x10(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  0xe5a(%rip),%xmm1
addsd  -0x1c0(%rbp),%xmm1
ff
movsd  0xe4a(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x380(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x388(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
movsd  0xe1a(%rip),%xmm1
addsd  -0x1c0(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
call   1060 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xdfe(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x3a0(%rbp)
ff
jmp    1295 <basym+0x145>
movsd  -0x18(%rbp),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  0xdd3(%rip),%xmm1
addsd  -0x1c0(%rbp),%xmm1
ff
movsd  0xdc3(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x380(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x388(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  0xd93(%rip),%xmm1
addsd  -0x1c0(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
call   1060 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xd77(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x3a0(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x20(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x400(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
call   1030 <rlog1@plt>
movsd  %xmm0,-0x408(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x410(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
call   1030 <rlog1@plt>
movsd  -0x410(%rbp),%xmm1
ff
movsd  -0x408(%rbp),%xmm2
ff
movaps %xmm0,%xmm3
movsd  -0x400(%rbp),%xmm0
ff
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movsd  -0x1b8(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,-0x1e8(%rbp)
ff
movsd  -0x1e8(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1385 <basym+0x235>
jp     1385 <basym+0x235>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1b21 <basym+0x9d1>
movsd  -0x1b8(%rbp),%xmm0
ff
call   1060 <sqrt@plt>
movsd  %xmm0,-0x3a8(%rbp)
ff
movsd  -0x3a8(%rbp),%xmm0
ff
divsd  0x2cae(%rip),%xmm0
movsd  0xc56(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x200(%rbp)
ff
movsd  -0x1b8(%rbp),%xmm0
ff
addsd  -0x1b8(%rbp),%xmm0
ff
movsd  %xmm0,-0x3b0(%rbp)
ff
movsd  0xc3a(%rip),%xmm0
mulsd  -0x388(%rbp),%xmm0
ff
movsd  %xmm0,-0x2b0(%rbp)
ff
movsd  0xc1a(%rip),%xmm0
mulsd  -0x2b0(%rbp),%xmm0
ff
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x100(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x1b0(%rbp)
ff
movsd  0xbd3(%rip),%xmm0
divsd  0x2c13(%rip),%xmm0
movsd  %xmm0,-0x418(%rbp)
ff
movsd  -0x3a8(%rbp),%xmm0
ff
mov    $0x1,%edi
call   1050 <erfc1@plt>
movaps %xmm0,%xmm1
movsd  -0x418(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x368(%rbp)
ff
movsd  0x2be2(%rip),%xmm0
movsd  %xmm0,-0x370(%rbp)
ff
movsd  -0x368(%rbp),%xmm1
ff
movsd  -0x1b0(%rbp),%xmm0
ff
mulsd  -0x3a0(%rbp),%xmm0
ff
movsd  -0x370(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x3e0(%rbp)
ff
movsd  0xb4a(%rip),%xmm0
movsd  %xmm0,-0x1e0(%rbp)
ff
movsd  -0x1c0(%rbp),%xmm0
ff
mulsd  -0x1c0(%rbp),%xmm0
ff
movsd  %xmm0,-0x378(%rbp)
ff
movsd  0xb22(%rip),%xmm0
movsd  %xmm0,-0x3b8(%rbp)
ff
movsd  -0x3a0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1f8(%rbp)
ff
movsd  -0x200(%rbp),%xmm0
ff
movsd  %xmm0,-0x3e8(%rbp)
ff
movsd  -0x3b0(%rbp),%xmm0
ff
movsd  %xmm0,-0x3c0(%rbp)
ff
mov    0x2b2c(%rip),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x2,-0x1d0(%rbp)
mov    -0x1d0(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1aba <basym+0x96a>
movsd  -0x378(%rbp),%xmm0
ff
mulsd  -0x3b8(%rbp),%xmm0
ff
movsd  %xmm0,-0x3b8(%rbp)
ff
movsd  0xad0(%rip),%xmm0
mulsd  -0x380(%rbp),%xmm0
ff
movsd  -0x1c0(%rbp),%xmm1
ff
movsd  -0x3b8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm1
movsd  0xa84(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x1d0(%rbp),%xmm1
ff
movsd  0xa94(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x1d0(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  %xmm0,-0x2b0(%rbp,%rax,8)
ff ff
mov    -0x1d0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3cc(%rbp)
movsd  -0x3b8(%rbp),%xmm0
ff
addsd  -0x1e0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1e0(%rbp)
ff
movsd  0xa49(%rip),%xmm0
mulsd  -0x388(%rbp),%xmm0
ff
mulsd  -0x1e0(%rbp),%xmm0
ff
cvtsi2sdl -0x1d0(%rbp),%xmm1
ff
movsd  0xa21(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x3cc(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  %xmm0,-0x2b0(%rbp,%rax,8)
ff ff
mov    -0x3cc(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x1d0(%rbp),%eax
mov    %eax,-0x1c4(%rbp)
mov    -0x1c4(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     1917 <basym+0x7c7>
cvtsi2sdl -0x1c4(%rbp),%xmm0
ff
movsd  0x9b1(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0x9b5(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  -0x1d8(%rbp),%xmm0
ff
mulsd  -0x2b0(%rbp),%xmm0
ff
movsd  %xmm0,-0x360(%rbp)
ff
mov    -0x1c4(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x2,-0x1cc(%rbp)
mov    -0x1cc(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     17da <basym+0x68a>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x3f0(%rbp)
ff
mov    -0x1cc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x3c8(%rbp)
mov    -0x3c8(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x1,-0x1c8(%rbp)
mov    -0x1c8(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     177a <basym+0x62a>
mov    -0x1cc(%rbp),%eax
sub    -0x1c8(%rbp),%eax
mov    %eax,-0x3d4(%rbp)
cvtsi2sdl -0x1c8(%rbp),%xmm0
ff
movsd  -0x1d8(%rbp),%xmm2
ff
cvtsi2sdl -0x3d4(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x1c8(%rbp),%eax
sub    $0x1,%eax
cltq
mulsd  -0x2b0(%rbp,%rax,8),%xmm0
ff ff
mov    -0x3d4(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x360(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x3f0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x3f0(%rbp)
ff
mov    -0x1c8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c8(%rbp)
jmp    16ce <basym+0x57e>
movsd  -0x1d8(%rbp),%xmm0
ff
mov    -0x1cc(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x2b0(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x3f0(%rbp),%xmm1
ff
cvtsi2sdl -0x1cc(%rbp),%xmm3
ff
divsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x1cc(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  %xmm0,-0x360(%rbp,%rax,8)
ff ff
mov    -0x1cc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1cc(%rbp)
jmp    1692 <basym+0x542>
mov    -0x1c4(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x360(%rbp,%rax,8),%xmm0
ff ff
cvtsi2sdl -0x1c4(%rbp),%xmm1
ff
movsd  0x802(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x1c4(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  %xmm0,-0x100(%rbp,%rax,8)
ff ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x3f8(%rbp)
ff
mov    -0x1c4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x3c4(%rbp)
mov    -0x3c4(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x1,-0x1c8(%rbp)
mov    -0x1c8(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     18bc <basym+0x76c>
mov    -0x1c4(%rbp),%eax
sub    -0x1c8(%rbp),%eax
mov    %eax,-0x3d0(%rbp)
mov    -0x3d0(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x1b0(%rbp,%rax,8),%xmm0
ff ff
mov    -0x1c8(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x100(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x3f8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x3f8(%rbp)
ff
mov    -0x1c8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c8(%rbp)
jmp    1847 <basym+0x6f7>
movsd  -0x3f8(%rbp),%xmm0
ff
mov    -0x1c4(%rbp),%eax
sub    $0x1,%eax
cltq
addsd  -0x100(%rbp,%rax,8),%xmm0
ff ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x1c4(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  %xmm0,-0x1b0(%rbp,%rax,8)
ff ff
mov    -0x1c4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c4(%rbp)
jmp    1630 <basym+0x4e0>
movsd  0x2739(%rip),%xmm0
movsd  -0x3e8(%rbp),%xmm2
ff
cvtsi2sdl -0x1d0(%rbp),%xmm1
ff
movsd  0x6c9(%rip),%xmm3
subsd  %xmm3,%xmm1
mulsd  -0x368(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x368(%rbp)
ff
movsd  0x26fd(%rip),%xmm0
movsd  -0x3c0(%rbp),%xmm2
ff
cvtsi2sdl -0x1d0(%rbp),%xmm1
ff
mulsd  -0x370(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x370(%rbp)
ff
movsd  -0x3b0(%rbp),%xmm0
ff
mulsd  -0x3e8(%rbp),%xmm0
ff
movsd  %xmm0,-0x3e8(%rbp)
ff
movsd  -0x3b0(%rbp),%xmm0
ff
mulsd  -0x3c0(%rbp),%xmm0
ff
movsd  %xmm0,-0x3c0(%rbp)
ff
movsd  -0x3a0(%rbp),%xmm0
ff
mulsd  -0x1f8(%rbp),%xmm0
ff
movsd  %xmm0,-0x1f8(%rbp)
ff
mov    -0x1d0(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x1b0(%rbp,%rax,8),%xmm0
ff ff
mulsd  -0x1f8(%rbp),%xmm0
ff
mulsd  -0x368(%rbp),%xmm0
ff
movsd  %xmm0,-0x390(%rbp)
ff
movsd  -0x3a0(%rbp),%xmm0
ff
mulsd  -0x1f8(%rbp),%xmm0
ff
movsd  %xmm0,-0x1f8(%rbp)
ff
mov    -0x3cc(%rbp),%eax
sub    $0x1,%eax
cltq
movsd  -0x1b0(%rbp,%rax,8),%xmm0
ff ff
mulsd  -0x1f8(%rbp),%xmm0
ff
mulsd  -0x370(%rbp),%xmm0
ff
movsd  %xmm0,-0x398(%rbp)
ff
movsd  -0x390(%rbp),%xmm0
ff
addsd  -0x398(%rbp),%xmm0
ff
addsd  -0x3e0(%rbp),%xmm0
ff
movsd  %xmm0,-0x3e0(%rbp)
ff
movsd  -0x390(%rbp),%xmm1
ff
movaps 0x5c6(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  -0x398(%rbp),%xmm0
ff
movaps 0x5b3(%rip),%xmm2
pand   %xmm2,%xmm0
addsd  %xmm0,%xmm1
movsd  -0x28(%rbp),%xmm0
mulsd  -0x3e0(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
jb     1aa1 <basym+0x951>
jmp    1abf <basym+0x96f>
jmp    1aa6 <basym+0x956>
mov    -0x1d0(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1d0(%rbp)
jmp    1529 <basym+0x3d9>
jmp    1abf <basym+0x96f>
movsd  -0x10(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
call   1070 <bcorr@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1040 <exp@plt>
movsd  %xmm0,-0x1f0(%rbp)
ff
movsd  0x2556(%rip),%xmm0
mulsd  -0x1e8(%rbp),%xmm0
ff
mulsd  -0x1f0(%rbp),%xmm0
ff
mulsd  -0x3e0(%rbp),%xmm0
ff
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x420,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<basym>:
push   %rbx
sub    $0x360,%rsp
movsd  %xmm3,0x78(%rsp)
movaps %xmm2,0x20(%rsp)
movapd %xmm0,%xmm2
ucomisd %xmm1,%xmm0
jae    11a9 <basym+0x59>
movapd %xmm2,%xmm0
divsd  %xmm1,%xmm0
movsd  0xe83(%rip),%xmm3
movsd  %xmm0,0x48(%rsp)
addsd  %xmm0,%xmm3
movapd %xmm1,%xmm0
subsd  %xmm2,%xmm0
divsd  %xmm1,%xmm0
movapd %xmm0,0x60(%rsp)
movsd  %xmm3,0x50(%rsp)
movapd %xmm3,%xmm0
mulsd  %xmm2,%xmm0
jmp    11e3 <basym+0x93>
movapd %xmm1,%xmm0
divsd  %xmm2,%xmm0
movsd  0xe47(%rip),%xmm3
movsd  %xmm0,0x48(%rsp)
addsd  %xmm0,%xmm3
movapd %xmm1,%xmm0
subsd  %xmm2,%xmm0
divsd  %xmm2,%xmm0
movapd %xmm0,0x60(%rsp)
movsd  %xmm3,0x50(%rsp)
movapd %xmm3,%xmm0
mulsd  %xmm1,%xmm0
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
movsd  %xmm1,0x40(%rsp)
movsd  %xmm2,0x38(%rsp)
jb     11ff <basym+0xaf>
sqrtsd %xmm0,%xmm0
jmp    1210 <basym+0xc0>
call   1060 <sqrt@plt>
movsd  0x38(%rsp),%xmm2
movsd  0x40(%rsp),%xmm1
movsd  %xmm0,0x8(%rsp)
movapd 0x20(%rsp),%xmm3
movapd %xmm3,%xmm0
xorpd  0xe18(%rip),%xmm0
divsd  %xmm2,%xmm0
divsd  %xmm1,%xmm3
movapd %xmm3,0x20(%rsp)
call   1030 <rlog1@plt>
movapd %xmm0,0x10(%rsp)
movapd 0x20(%rsp),%xmm0
call   1030 <rlog1@plt>
mulsd  0x40(%rsp),%xmm0
movapd 0x10(%rsp),%xmm1
mulsd  0x38(%rsp),%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
movapd 0xdd0(%rip),%xmm0
xorpd  %xmm1,%xmm0
call   1040 <exp@plt>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
jne    1289 <basym+0x139>
jnp    17fd <basym+0x6ad>
movsd  %xmm0,0x70(%rsp)
movsd  0xd69(%rip),%xmm0
movapd %xmm0,%xmm1
divsd  0x8(%rsp),%xmm1
movsd  %xmm1,0x8(%rsp)
divsd  0x50(%rsp),%xmm0
movapd %xmm0,0x50(%rsp)
xorpd  %xmm0,%xmm0
movapd 0x10(%rsp),%xmm3
ucomisd %xmm0,%xmm3
jb     12cc <basym+0x17c>
xorps  %xmm0,%xmm0
sqrtsd %xmm3,%xmm0
jmp    12db <basym+0x18b>
movapd %xmm3,%xmm0
call   1060 <sqrt@plt>
movapd 0x10(%rsp),%xmm3
movapd %xmm0,%xmm4
divsd  0xd21(%rip),%xmm4
movapd %xmm4,%xmm1
unpcklpd %xmm3,%xmm1
movsd  0xd19(%rip),%xmm2
mulsd  %xmm2,%xmm4
unpcklpd %xmm3,%xmm2
addpd  %xmm1,%xmm2
movapd %xmm2,0x10(%rsp)
shufpd $0x2,%xmm2,%xmm4
movapd %xmm4,0x90(%rsp)
movsd  0xcf9(%rip),%xmm1
mulsd  0x60(%rsp),%xmm1
movsd  %xmm1,0x2b0(%rsp)
movsd  0xcea(%rip),%xmm2
movsd  %xmm1,0x20(%rsp)
mulsd  %xmm1,%xmm2
movapd %xmm2,0x80(%rsp)
movsd  %xmm2,0x150(%rsp)
movapd 0xce6(%rip),%xmm1
xorpd  %xmm2,%xmm1
movlpd %xmm1,0x200(%rsp)
mov    $0x1,%edi
call   1050 <erfc1@plt>
movapd 0x90(%rsp),%xmm12
mulsd  0xca5(%rip),%xmm0
movsd  0x8(%rsp),%xmm7
movapd 0x80(%rsp),%xmm2
mulsd  %xmm7,%xmm2
mulsd  0xc6a(%rip),%xmm2
movapd %xmm0,%xmm4
subsd  %xmm2,%xmm4
movsd  0x48(%rsp),%xmm14
mulsd  %xmm14,%xmm14
movapd 0x50(%rsp),%xmm11
unpcklpd 0x60(%rsp),%xmm11
addpd  %xmm11,%xmm11
shufpd $0x2,0xc82(%rip),%xmm0
movapd 0x10(%rsp),%xmm2
unpckhpd %xmm2,%xmm2
movapd %xmm2,0x10(%rsp)
lea    0x208(%rsp),%r9
mov    $0x3,%ecx
mov    $0x2,%r11d
xor    %r8d,%r8d
movapd 0xc73(%rip),%xmm13
xor    %r10d,%r10d
movsd  0xbf7(%rip),%xmm15
movapd %xmm15,%xmm9
xchg   %ax,%ax
mulsd  %xmm14,%xmm9
movsd  0x48(%rsp),%xmm1
mulsd  %xmm9,%xmm1
movsd  0xbd8(%rip),%xmm2
addsd  %xmm2,%xmm1
cvtsi2sd %r11d,%xmm10
addsd  %xmm9,%xmm15
unpcklpd %xmm15,%xmm1
mulpd  %xmm11,%xmm1
movapd %xmm10,%xmm3
unpcklpd %xmm10,%xmm3
addpd  0xc0e(%rip),%xmm3
divpd  %xmm3,%xmm1
movupd %xmm1,0x2a8(%rsp,%r11,8)
movapd %xmm10,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm3
mulsd  0xbab(%rip),%xmm3
movapd %xmm3,%xmm1
mulsd  0x20(%rsp),%xmm1
movsd  %xmm1,0xa0(%rsp)
mov    $0x1,%edi
mov    $0x2,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorpd  %xmm1,%xmm1
mov    %rdi,%rsi
xor    %edx,%edx
nopl   0x0(%rax)
lea    0x1(%rdx),%rbx
xorps  %xmm5,%xmm5
cvtsi2sd %ebx,%xmm5
xorps  %xmm8,%xmm8
cvtsi2sd %esi,%xmm8
mulsd  %xmm3,%xmm5
subsd  %xmm8,%xmm5
mulsd  0x2b0(%rsp,%rdx,8),%xmm5
mulsd  0x98(%rsp,%rsi,8),%xmm5
addsd  %xmm5,%xmm1
mov    %rbx,%rdx
add    $0xffffffffffffffff,%rsi
jne    14b0 <basym+0x360>
xorps  %xmm5,%xmm5
cvtsi2sd %eax,%xmm5
divsd  %xmm5,%xmm1
movsd  0x2a8(%rsp,%rax,8),%xmm5
mulsd  %xmm3,%xmm5
addsd  %xmm1,%xmm5
movsd  %xmm5,0x98(%rsp,%rax,8)
add    $0x1,%rax
add    $0x1,%rdi
cmp    %rcx,%rax
jne    14a0 <basym+0x350>
movsd  0x98(%rsp,%r11,8),%xmm3
divsd  %xmm6,%xmm3
movsd  %xmm3,0x148(%rsp,%r11,8)
mov    $0x1,%esi
movapd %xmm4,%xmm8
test   %r10,%r10
je     15a0 <basym+0x450>
xorpd  %xmm5,%xmm5
mov    $0x1,%eax
mov    %r8,%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x1f8(%rsp,%rdi,8),%xmm1
movsd  0x200(%rsp,%rdi,8),%xmm6
mulsd  0x148(%rsp,%rax,8),%xmm6
mulsd  0x150(%rsp,%rax,8),%xmm1
addsd  %xmm5,%xmm6
addsd  %xmm6,%xmm1
add    $0x2,%rax
movapd %xmm1,%xmm5
add    $0xfffffffffffffffe,%rdi
jne    1560 <basym+0x410>
jmp    15a9 <basym+0x459>
nopl   0x0(%rax)
xorpd  %xmm1,%xmm1
mov    $0x1,%eax
mov    %rax,%rdx
not    %rdx
add    %r11,%rdx
movsd  0x200(%rsp,%rdx,8),%xmm5
mulsd  0x148(%rsp,%rax,8),%xmm5
addsd  %xmm1,%xmm5
addsd  %xmm5,%xmm3
xorpd  0xa6c(%rip),%xmm3
movlpd %xmm3,0x1f8(%rsp,%r11,8)
mov    %r11d,%eax
or     $0x1,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  0xa0d(%rip),%xmm1
movapd %xmm1,%xmm6
mulsd  0xa21(%rip),%xmm6
movapd %xmm6,%xmm5
mulsd  0x20(%rsp),%xmm5
movsd  %xmm5,0xa0(%rsp)
mov    $0x2,%eax
nopw   0x0(%rax,%rax,1)
xorpd  %xmm5,%xmm5
xor    %edx,%edx
mov    %rsi,%rdi
nopl   0x0(%rax)
lea    0x1(%rdx),%rbx
xorps  %xmm2,%xmm2
cvtsi2sd %ebx,%xmm2
xorps  %xmm4,%xmm4
cvtsi2sd %edi,%xmm4
mulsd  %xmm6,%xmm2
subsd  %xmm4,%xmm2
mulsd  0x2b0(%rsp,%rdx,8),%xmm2
mulsd  0x98(%rsp,%rdi,8),%xmm2
addsd  %xmm2,%xmm5
mov    %rbx,%rdx
add    $0xffffffffffffffff,%rdi
jne    1630 <basym+0x4e0>
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
divsd  %xmm2,%xmm5
movsd  0x2a8(%rsp,%rax,8),%xmm2
mulsd  %xmm6,%xmm2
addsd  %xmm5,%xmm2
movsd  %xmm2,0x98(%rsp,%rax,8)
lea    0x1(%rax),%rdx
add    $0x1,%rsi
cmp    %rcx,%rax
mov    %rdx,%rax
jne    1620 <basym+0x4d0>
movsd  0xa0(%rsp,%r11,8),%xmm6
divsd  %xmm1,%xmm6
movsd  %xmm6,0x150(%rsp,%r11,8)
xorpd  %xmm1,%xmm1
mov    %r9,%rax
xor    %esi,%esi
nop
movsd  (%rax),%xmm2
mulsd  0x150(%rsp,%rsi,8),%xmm2
addsd  %xmm1,%xmm2
movsd  -0x8(%rax),%xmm1
mulsd  0x158(%rsp,%rsi,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x2,%rsi
add    $0xfffffffffffffff0,%rax
cmp    %rsi,%r11
jne    16c0 <basym+0x570>
addsd  %xmm1,%xmm6
xorpd  0x944(%rip),%xmm6
movlpd %xmm6,0x200(%rsp,%r11,8)
movapd %xmm10,%xmm1
addsd  0x91d(%rip),%xmm1
unpcklpd %xmm10,%xmm1
mulpd  %xmm1,%xmm0
movapd %xmm12,%xmm1
mulpd  %xmm13,%xmm1
addpd  %xmm0,%xmm1
movapd %xmm7,%xmm2
movsd  0x8(%rsp),%xmm7
mulsd  %xmm7,%xmm2
mulsd  %xmm2,%xmm7
unpcklpd %xmm7,%xmm2
unpcklpd %xmm6,%xmm3
mulpd  %xmm2,%xmm3
mulpd  %xmm1,%xmm3
movapd %xmm3,%xmm1
unpckhpd %xmm3,%xmm1
addsd  %xmm3,%xmm1
movapd %xmm8,%xmm4
addsd  %xmm1,%xmm4
andpd  0x917(%rip),%xmm3
pshufd $0xee,%xmm3,%xmm1
addsd  %xmm3,%xmm1
movapd %xmm4,%xmm2
mulsd  0x78(%rsp),%xmm2
ucomisd %xmm1,%xmm2
jae    17bd <basym+0x66d>
movapd 0x10(%rsp),%xmm1
mulpd  %xmm12,%xmm1
mulpd  %xmm13,%xmm12
addpd  %xmm12,%xmm0
lea    0x2(%r11),%rax
add    $0x2,%rcx
add    $0x1,%r10
add    $0x2,%r8
add    $0x10,%r9
movapd %xmm1,%xmm12
cmp    $0x13,%r11
mov    %rax,%r11
jb     1410 <basym+0x2c0>
movsd  0x38(%rsp),%xmm0
movsd  0x40(%rsp),%xmm1
movapd %xmm4,0x20(%rsp)
call   1070 <bcorr@plt>
xorpd  0x864(%rip),%xmm0
call   1040 <exp@plt>
movapd 0x20(%rsp),%xmm2
movsd  0x70(%rsp),%xmm1
mulsd  0x843(%rip),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm1,%xmm2
movapd %xmm2,%xmm0
add    $0x360,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@basym.num = internal global i32 20
@basym.e0 = internal global double 0x3FF20DD750429B61
@basym.e1 = internal global double 0x3FD6A09E667F3BD1
define double @basym(double %0, double %1, double %2, double %3) {
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca [21 x double]
%18 = alloca [21 x double]
%19 = alloca double
%20 = alloca double
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca [21 x double]
%32 = alloca [21 x double]
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca double
%38 = alloca double
%39 = alloca double
%40 = alloca double
%41 = alloca double
%42 = alloca double
%43 = alloca double
%44 = alloca double
%45 = alloca i32
%46 = alloca i32
%47 = alloca i32
%48 = alloca i32
%49 = alloca i32
%50 = alloca double
%51 = alloca double
%52 = alloca double
%53 = 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* %14
%54 = load double, double* %6
%55 = load double, double* %7
%56 = fcmp oge double %54, %55
br i1 %56, label %57, label %58
57: ; preds = %4
br label %76
58: ; preds = %4
%59 = load double, double* %6
%60 = load double, double* %7
%61 = fdiv double %59, %60
store double %61, double* %20
%62 = load double, double* %20
%63 = fadd double %62, 1.000000e+00
%64 = fdiv double 1.000000e+00, %63
store double %64, double* %36
%65 = load double, double* %7
%66 = load double, double* %6
%67 = fsub double %65, %66
%68 = load double, double* %7
%69 = fdiv double %67, %68
store double %69, double* %37
%70 = load double, double* %6
%71 = load double, double* %20
%72 = fadd double %71, 1.000000e+00
%73 = fmul double %70, %72
%74 = call double @sqrt(double %73)
%75 = fdiv double 1.000000e+00, %74
store double %75, double* %40
br label %94
76: ; preds = %57
%77 = load double, double* %7
%78 = load double, double* %6
%79 = fdiv double %77, %78
store double %79, double* %20
%80 = load double, double* %20
%81 = fadd double %80, 1.000000e+00
%82 = fdiv double 1.000000e+00, %81
store double %82, double* %36
%83 = load double, double* %7
%84 = load double, double* %6
%85 = fsub double %83, %84
%86 = load double, double* %6
%87 = fdiv double %85, %86
store double %87, double* %37
%88 = load double, double* %7
%89 = load double, double* %20
%90 = fadd double %89, 1.000000e+00
%91 = fmul double %88, %90
%92 = call double @sqrt(double %91)
%93 = fdiv double 1.000000e+00, %92
store double %93, double* %40
br label %94
94: ; preds = %76, %58
%95 = load double, double* %8
%96 = fneg double %95
%97 = load double, double* %6
%98 = fdiv double %96, %97
store double %98, double* %15
%99 = load double, double* %8
%100 = load double, double* %7
%101 = fdiv double %99, %100
store double %101, double* %16
%102 = load double, double* %6
%103 = load double, double* %15
%104 = call double @rlog1(double %103)
%105 = load double, double* %7
%106 = load double, double* %16
%107 = call double @rlog1(double %106)
%108 = fmul double %105, %107
%109 = call double @llvm.fmuladd.f64(double %102, double %104, double %108)
store double %109, double* %19
%110 = load double, double* %19
%111 = fneg double %110
%112 = call double @exp(double %111)
store double %112, double* %27
%113 = load double, double* %27
%114 = fcmp oeq double %113, 0.000000e+00
br i1 %114, label %115, label %117
115: ; preds = %94
%116 = load double, double* %14
store double %116, double* %5
br label %408
117: ; preds = %94
%118 = load double, double* %19
%119 = call double @sqrt(double %118)
store double %119, double* %41
%120 = load double, double* %41
%121 = load double, double* @basym.e1
%122 = fdiv double %120, %121
%123 = fmul double %122, 5.000000e-01
store double %123, double* %30
%124 = load double, double* %19
%125 = load double, double* %19
%126 = fadd double %124, %125
store double %126, double* %42
%127 = load double, double* %37
%128 = fmul double %127, 0x3FE5555555555555
%129 = getelementptr inbounds [21 x double], [21 x double]* %31, i64 0, i64 0
store double %128, double* %129
%130 = getelementptr inbounds [21 x double], [21 x double]* %31, i64 0, i64 0
%131 = load double, double* %130
%132 = fmul double %131, -5.000000e-01
%133 = getelementptr inbounds [21 x double], [21 x double]* %17, i64 0, i64 0
store double %132, double* %133
%134 = getelementptr inbounds [21 x double], [21 x double]* %17, i64 0, i64 0
%135 = load double, double* %134
%136 = fneg double %135
%137 = getelementptr inbounds [21 x double], [21 x double]* %18, i64 0, i64 0
store double %136, double* %137
%138 = load double, double* @basym.e0
%139 = fdiv double 5.000000e-01, %138
%140 = load double, double* %41
%141 = call double @erfc1(i32 1, double %140)
%142 = fmul double %139, %141
store double %142, double* %33
%143 = load double, double* @basym.e1
store double %143, double* %34
%144 = load double, double* %33
%145 = getelementptr inbounds [21 x double], [21 x double]* %18, i64 0, i64 0
%146 = load double, double* %145
%147 = load double, double* %40
%148 = fmul double %146, %147
%149 = load double, double* %34
%150 = call double @llvm.fmuladd.f64(double %148, double %149, double %144)
store double %150, double* %50
store double 1.000000e+00, double* %26
%151 = load double, double* %20
%152 = load double, double* %20
%153 = fmul double %151, %152
store double %153, double* %35
store double 1.000000e+00, double* %43
%154 = load double, double* %40
store double %154, double* %29
%155 = load double, double* %30
store double %155, double* %51
%156 = load double, double* %42
store double %156, double* %44
%157 = load i32, i32* @basym.num
store i32 %157, i32* %10
store i32 2, i32* %24
br label %158
158: ; preds = %390, %117
%159 = load i32, i32* %24
%160 = load i32, i32* %10
%161 = icmp sle i32 %159, %160
br i1 %161, label %162, label %393
162: ; preds = %158
%163 = load double, double* %35
%164 = load double, double* %43
%165 = fmul double %163, %164
store double %165, double* %43
%166 = load double, double* %36
%167 = fmul double %166, 2.000000e+00
%168 = load double, double* %20
%169 = load double, double* %43
%170 = call double @llvm.fmuladd.f64(double %168, double %169, double 1.000000e+00)
%171 = fmul double %167, %170
%172 = load i32, i32* %24
%173 = sitofp i32 %172 to double
%174 = fadd double %173, 2.000000e+00
%175 = fdiv double %171, %174
%176 = load i32, i32* %24
%177 = sub nsw i32 %176, 1
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds [21 x double], [21 x double]* %31, i64 0, i64 %178
store double %175, double* %179
%180 = load i32, i32* %24
%181 = add nsw i32 %180, 1
store i32 %181, i32* %47
%182 = load double, double* %43
%183 = load double, double* %26
%184 = fadd double %183, %182
store double %184, double* %26
%185 = load double, double* %37
%186 = fmul double %185, 2.000000e+00
%187 = load double, double* %26
%188 = fmul double %186, %187
%189 = load i32, i32* %24
%190 = sitofp i32 %189 to double
%191 = fadd double %190, 3.000000e+00
%192 = fdiv double %188, %191
%193 = load i32, i32* %47
%194 = sub nsw i32 %193, 1
%195 = sext i32 %194 to i64
%196 = getelementptr inbounds [21 x double], [21 x double]* %31, i64 0, i64 %195
store double %192, double* %196
%197 = load i32, i32* %47
store i32 %197, i32* %11
%198 = load i32, i32* %24
store i32 %198, i32* %21
br label %199
199: ; preds = %325, %162
%200 = load i32, i32* %21
%201 = load i32, i32* %11
%202 = icmp sle i32 %200, %201
br i1 %202, label %203, label %328
203: ; preds = %199
%204 = load i32, i32* %21
%205 = sitofp i32 %204 to double
%206 = fadd double %205, 1.000000e+00
%207 = fmul double %206, -5.000000e-01
store double %207, double* %25
%208 = load double, double* %25
%209 = getelementptr inbounds [21 x double], [21 x double]* %31, i64 0, i64 0
%210 = load double, double* %209
%211 = fmul double %208, %210
%212 = getelementptr inbounds [21 x double], [21 x double]* %32, i64 0, i64 0
store double %211, double* %212
%213 = load i32, i32* %21
store i32 %213, i32* %12
store i32 2, i32* %23
br label %214
214: ; preds = %269, %203
%215 = load i32, i32* %23
%216 = load i32, i32* %12
%217 = icmp sle i32 %215, %216
br i1 %217, label %218, label %272
218: ; preds = %214
store double 0.000000e+00, double* %52
%219 = load i32, i32* %23
%220 = sub nsw i32 %219, 1
store i32 %220, i32* %46
%221 = load i32, i32* %46
store i32 %221, i32* %13
store i32 1, i32* %22
br label %222
222: ; preds = %250, %218
%223 = load i32, i32* %22
%224 = load i32, i32* %13
%225 = icmp sle i32 %223, %224
br i1 %225, label %226, label %253
226: ; preds = %222
%227 = load i32, i32* %23
%228 = load i32, i32* %22
%229 = sub nsw i32 %227, %228
store i32 %229, i32* %49
%230 = load i32, i32* %22
%231 = sitofp i32 %230 to double
%232 = load double, double* %25
%233 = load i32, i32* %49
%234 = sitofp i32 %233 to double
%235 = fneg double %234
%236 = call double @llvm.fmuladd.f64(double %231, double %232, double %235)
%237 = load i32, i32* %22
%238 = sub nsw i32 %237, 1
%239 = sext i32 %238 to i64
%240 = getelementptr inbounds [21 x double], [21 x double]* %31, i64 0, i64 %239
%241 = load double, double* %240
%242 = fmul double %236, %241
%243 = load i32, i32* %49
%244 = sub nsw i32 %243, 1
%245 = sext i32 %244 to i64
%246 = getelementptr inbounds [21 x double], [21 x double]* %32, i64 0, i64 %245
%247 = load double, double* %246
%248 = load double, double* %52
%249 = call double @llvm.fmuladd.f64(double %242, double %247, double %248)
store double %249, double* %52
br label %250
250: ; preds = %226
%251 = load i32, i32* %22
%252 = add nsw i32 %251, 1
store i32 %252, i32* %22
br label %222
253: ; preds = %222
%254 = load double, double* %25
%255 = load i32, i32* %23
%256 = sub nsw i32 %255, 1
%257 = sext i32 %256 to i64
%258 = getelementptr inbounds [21 x double], [21 x double]* %31, i64 0, i64 %257
%259 = load double, double* %258
%260 = load double, double* %52
%261 = load i32, i32* %23
%262 = sitofp i32 %261 to double
%263 = fdiv double %260, %262
%264 = call double @llvm.fmuladd.f64(double %254, double %259, double %263)
%265 = load i32, i32* %23
%266 = sub nsw i32 %265, 1
%267 = sext i32 %266 to i64
%268 = getelementptr inbounds [21 x double], [21 x double]* %32, i64 0, i64 %267
store double %264, double* %268
br label %269
269: ; preds = %253
%270 = load i32, i32* %23
%271 = add nsw i32 %270, 1
store i32 %271, i32* %23
br label %214
272: ; preds = %214
%273 = load i32, i32* %21
%274 = sub nsw i32 %273, 1
%275 = sext i32 %274 to i64
%276 = getelementptr inbounds [21 x double], [21 x double]* %32, i64 0, i64 %275
%277 = load double, double* %276
%278 = load i32, i32* %21
%279 = sitofp i32 %278 to double
%280 = fadd double %279, 1.000000e+00
%281 = fdiv double %277, %280
%282 = load i32, i32* %21
%283 = sub nsw i32 %282, 1
%284 = sext i32 %283 to i64
%285 = getelementptr inbounds [21 x double], [21 x double]* %17, i64 0, i64 %284
store double %281, double* %285
store double 0.000000e+00, double* %53
%286 = load i32, i32* %21
%287 = sub nsw i32 %286, 1
store i32 %287, i32* %45
%288 = load i32, i32* %45
store i32 %288, i32* %12
store i32 1, i32* %22
br label %289
289: ; preds = %309, %272
%290 = load i32, i32* %22
%291 = load i32, i32* %12
%292 = icmp sle i32 %290, %291
br i1 %292, label %293, label %312
293: ; preds = %289
%294 = load i32, i32* %21
%295 = load i32, i32* %22
%296 = sub nsw i32 %294, %295
store i32 %296, i32* %48
%297 = load i32, i32* %48
%298 = sub nsw i32 %297, 1
%299 = sext i32 %298 to i64
%300 = getelementptr inbounds [21 x double], [21 x double]* %18, i64 0, i64 %299
%301 = load double, double* %300
%302 = load i32, i32* %22
%303 = sub nsw i32 %302, 1
%304 = sext i32 %303 to i64
%305 = getelementptr inbounds [21 x double], [21 x double]* %17, i64 0, i64 %304
%306 = load double, double* %305
%307 = load double, double* %53
%308 = call double @llvm.fmuladd.f64(double %301, double %306, double %307)
store double %308, double* %53
br label %309
309: ; preds = %293
%310 = load i32, i32* %22
%311 = add nsw i32 %310, 1
store i32 %311, i32* %22
br label %289
312: ; preds = %289
%313 = load double, double* %53
%314 = load i32, i32* %21
%315 = sub nsw i32 %314, 1
%316 = sext i32 %315 to i64
%317 = getelementptr inbounds [21 x double], [21 x double]* %17, i64 0, i64 %316
%318 = load double, double* %317
%319 = fadd double %313, %318
%320 = fneg double %319
%321 = load i32, i32* %21
%322 = sub nsw i32 %321, 1
%323 = sext i32 %322 to i64
%324 = getelementptr inbounds [21 x double], [21 x double]* %18, i64 0, i64 %323
store double %320, double* %324
br label %325
325: ; preds = %312
%326 = load i32, i32* %21
%327 = add nsw i32 %326, 1
store i32 %327, i32* %21
br label %199
328: ; preds = %199
%329 = load double, double* @basym.e1
%330 = load double, double* %51
%331 = load i32, i32* %24
%332 = sitofp i32 %331 to double
%333 = fsub double %332, 1.000000e+00
%334 = load double, double* %33
%335 = fmul double %333, %334
%336 = call double @llvm.fmuladd.f64(double %329, double %330, double %335)
store double %336, double* %33
%337 = load double, double* @basym.e1
%338 = load double, double* %44
%339 = load i32, i32* %24
%340 = sitofp i32 %339 to double
%341 = load double, double* %34
%342 = fmul double %340, %341
%343 = call double @llvm.fmuladd.f64(double %337, double %338, double %342)
store double %343, double* %34
%344 = load double, double* %42
%345 = load double, double* %51
%346 = fmul double %344, %345
store double %346, double* %51
%347 = load double, double* %42
%348 = load double, double* %44
%349 = fmul double %347, %348
store double %349, double* %44
%350 = load double, double* %40
%351 = load double, double* %29
%352 = fmul double %350, %351
store double %352, double* %29
%353 = load i32, i32* %24
%354 = sub nsw i32 %353, 1
%355 = sext i32 %354 to i64
%356 = getelementptr inbounds [21 x double], [21 x double]* %18, i64 0, i64 %355
%357 = load double, double* %356
%358 = load double, double* %29
%359 = fmul double %357, %358
%360 = load double, double* %33
%361 = fmul double %359, %360
store double %361, double* %38
%362 = load double, double* %40
%363 = load double, double* %29
%364 = fmul double %362, %363
store double %364, double* %29
%365 = load i32, i32* %47
%366 = sub nsw i32 %365, 1
%367 = sext i32 %366 to i64
%368 = getelementptr inbounds [21 x double], [21 x double]* %18, i64 0, i64 %367
%369 = load double, double* %368
%370 = load double, double* %29
%371 = fmul double %369, %370
%372 = load double, double* %34
%373 = fmul double %371, %372
store double %373, double* %39
%374 = load double, double* %38
%375 = load double, double* %39
%376 = fadd double %374, %375
%377 = load double, double* %50
%378 = fadd double %377, %376
store double %378, double* %50
%379 = load double, double* %38
%380 = call double @llvm.fabs.f64(double %379)
%381 = load double, double* %39
%382 = call double @llvm.fabs.f64(double %381)
%383 = fadd double %380, %382
%384 = load double, double* %9
%385 = load double, double* %50
%386 = fmul double %384, %385
%387 = fcmp ole double %383, %386
br i1 %387, label %388, label %389
388: ; preds = %328
br label %394
389: ; preds = %328
br label %390
390: ; preds = %389
%391 = load i32, i32* %24
%392 = add nsw i32 %391, 2
store i32 %392, i32* %24
br label %158
393: ; preds = %158
br label %394
394: ; preds = %393, %388
%395 = load double, double* %6
%396 = load double, double* %7
%397 = call double @bcorr(double %395, double %396)
%398 = fneg double %397
%399 = call double @exp(double %398)
store double %399, double* %28
%400 = load double, double* @basym.e0
%401 = load double, double* %27
%402 = fmul double %400, %401
%403 = load double, double* %28
%404 = fmul double %402, %403
%405 = load double, double* %50
%406 = fmul double %404, %405
store double %406, double* %14
%407 = load double, double* %14
store double %407, double* %5
br label %408
408: ; preds = %394, %115
%409 = load double, double* %5
ret double %409
}

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