int odrrj(int *n, double *d__, double *e2, int *ifirst, int *ilast,
          double *rtol, int *offset, double *w, double *werr, double *work,
          int *iwork, double *pivmin, double *spdiam, int *info) {
  int i__1, i__2;
  double d__1, d__2;
  int i__, j, k, p;
  double s;
  int i1, i2, ii;
  double fac, mid;
  int cnt;
  double tmp, left;
  int iter, nint, prev, next, savi1;
  double right, width, dplus;
  int olnint, maxitr;
  --iwork;
  --work;
  --werr;
  --w;
  --e2;
  --d__;
  *info = 0;
  maxitr = (int)((log(*spdiam + *pivmin) - log(*pivmin)) / log(2.)) + 2;
  i1 = *ifirst;
  i2 = *ilast;
  nint = 0;
  prev = 0;
  i__1 = i2;
  for (i__ = i1; i__ <= i__1; ++i__) {
    k = i__ << 1;
    ii = i__ - *offset;
    left = w[ii] - werr[ii];
    mid = w[ii];
    right = w[ii] + werr[ii];
    width = right - mid;
    d__1 = fabs(left), d__2 = fabs(right);
    tmp = fmax(d__1, d__2);
    if (width < *rtol * tmp) {
      iwork[k - 1] = -1;
      if (i__ == i1 && i__ < i2) {
        i1 = i__ + 1;
      }
      if (prev >= i1 && i__ <= i2) {
        iwork[(prev << 1) - 1] = i__ + 1;
      }
    } else {
      prev = i__;
      fac = 1.;
    L20:
      cnt = 0;
      s = left;
      dplus = d__[1] - s;
      if (dplus < 0.) {
        ++cnt;
      }
      i__2 = *n;
      for (j = 2; j <= i__2; ++j) {
        dplus = d__[j] - s - e2[j - 1] / dplus;
        if (dplus < 0.) {
          ++cnt;
        }
      }
      if (cnt > i__ - 1) {
        left -= werr[ii] * fac;
        fac *= 2.;
        goto L20;
      }
      fac = 1.;
    L50:
      cnt = 0;
      s = right;
      dplus = d__[1] - s;
      if (dplus < 0.) {
        ++cnt;
      }
      i__2 = *n;
      for (j = 2; j <= i__2; ++j) {
        dplus = d__[j] - s - e2[j - 1] / dplus;
        if (dplus < 0.) {
          ++cnt;
        }
      }
      if (cnt < i__) {
        right += werr[ii] * fac;
        fac *= 2.;
        goto L50;
      }
      ++nint;
      iwork[k - 1] = i__ + 1;
      iwork[k] = cnt;
    }
    work[k - 1] = left;
    work[k] = right;
  }
  savi1 = i1;
  iter = 0;
L80:
  prev = i1 - 1;
  i__ = i1;
  olnint = nint;
  i__1 = olnint;
  for (p = 1; p <= i__1; ++p) {
    k = i__ << 1;
    ii = i__ - *offset;
    next = iwork[k - 1];
    left = work[k - 1];
    right = work[k];
    mid = (left + right) * .5;
    width = right - mid;
    d__1 = fabs(left), d__2 = fabs(right);
    tmp = fmax(d__1, d__2);
    if (width < *rtol * tmp || iter == maxitr) {
      --nint;
      iwork[k - 1] = 0;
      if (i1 == i__) {
        i1 = next;
      } else {
        if (prev >= i1) {
          iwork[(prev << 1) - 1] = next;
        }
      }
      i__ = next;
      goto L100;
    }
    prev = i__;
    cnt = 0;
    s = mid;
    dplus = d__[1] - s;
    if (dplus < 0.) {
      ++cnt;
    }
    i__2 = *n;
    for (j = 2; j <= i__2; ++j) {
      dplus = d__[j] - s - e2[j - 1] / dplus;
      if (dplus < 0.) {
        ++cnt;
      }
    }
    if (cnt <= i__ - 1) {
      work[k - 1] = mid;
    } else {
      work[k] = mid;
    }
    i__ = next;
  L100:;
  }
  ++iter;
  if (nint > 0 && iter <= maxitr) {
    goto L80;
  }
  i__1 = *ilast;
  for (i__ = savi1; i__ <= i__1; ++i__) {
    k = i__ << 1;
    ii = i__ - *offset;
    if (iwork[k - 1] == 0) {
      w[ii] = (work[k - 1] + work[k]) * .5;
      werr[ii] = work[k] - w[ii];
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<odrrj>:
push   %rbp
mov    %rsp,%rbp
sub    $0xe0,%rsp
mov    0x48(%rbp),%rax
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    0x30(%rbp),%rax
add    $0xfffffffffffffffc,%rax
mov    %rax,0x30(%rbp)
mov    0x28(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,0x28(%rbp)
mov    0x20(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,0x20(%rbp)
mov    0x18(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,0x18(%rbp)
mov    -0x18(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x10(%rbp)
mov    0x48(%rbp),%rax
movl   $0x0,(%rax)
mov    0x40(%rbp),%rax
movsd  (%rax),%xmm0
mov    0x38(%rbp),%rax
addsd  (%rax),%xmm0
call   1030 <log@plt>
movsd  %xmm0,-0xe0(%rbp)
ff
mov    0x38(%rbp),%rax
movsd  (%rax),%xmm0
call   1030 <log@plt>
movaps %xmm0,%xmm1
movsd  -0xe0(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  0xe12(%rip),%xmm0
call   1030 <log@plt>
movaps %xmm0,%xmm1
movsd  -0xd8(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
add    $0x2,%eax
mov    %eax,-0xd0(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x64(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x68(%rbp)
movl   $0x0,-0xa0(%rbp)
movl   $0x0,-0xa4(%rbp)
mov    -0x68(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     16ba <odrrj+0x5aa>
mov    -0x4c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x4c(%rbp),%eax
mov    0x10(%rbp),%rcx
sub    (%rcx),%eax
mov    %eax,-0x6c(%rbp)
mov    0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movaps 0xd44(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
movaps 0xd2c(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm3
movsd  -0x48(%rbp),%xmm2
movaps %xmm2,%xmm4
movaps %xmm3,%xmm0
cmpunordsd %xmm0,%xmm0
movaps %xmm0,%xmm1
pand   %xmm4,%xmm1
maxsd  %xmm3,%xmm2
pandn  %xmm2,%xmm0
por    %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0xc0(%rbp),%xmm1
ff
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
jbe    13c0 <odrrj+0x2b0>
mov    0x30(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movl   $0xffffffff,(%rax,%rcx,4)
mov    -0x4c(%rbp),%eax
cmp    -0x64(%rbp),%eax
jne    1384 <odrrj+0x274>
mov    -0x4c(%rbp),%eax
cmp    -0x68(%rbp),%eax
jge    1384 <odrrj+0x274>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
mov    -0xa4(%rbp),%eax
cmp    -0x64(%rbp),%eax
jl     13bb <odrrj+0x2ab>
mov    -0x4c(%rbp),%eax
cmp    -0x68(%rbp),%eax
jg     13bb <odrrj+0x2ab>
mov    -0x4c(%rbp),%edx
add    $0x1,%edx
mov    0x30(%rbp),%rax
mov    -0xa4(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
jmp    167d <odrrj+0x56d>
mov    -0x4c(%rbp),%eax
mov    %eax,-0xa4(%rbp)
movsd  0xc3f(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x84(%rbp)
movsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xc8(%rbp),%xmm0
ff
jbe    1423 <odrrj+0x313>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x38(%rbp)
movl   $0x2,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     14aa <odrrj+0x39a>
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
divsd  -0xc8(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xc8(%rbp),%xmm0
ff
jbe    1497 <odrrj+0x387>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    149c <odrrj+0x38c>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1433 <odrrj+0x323>
mov    -0x84(%rbp),%eax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    1516 <odrrj+0x406>
mov    0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x78(%rbp),%xmm2
movsd  -0x98(%rbp),%xmm1
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xaf9(%rip),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    13d6 <odrrj+0x2c6>
movsd  0xaf2(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x84(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xc8(%rbp),%xmm0
ff
jbe    1570 <odrrj+0x460>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x38(%rbp)
movl   $0x2,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     15f7 <odrrj+0x4e7>
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
divsd  -0xc8(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xc8(%rbp),%xmm0
ff
jbe    15e4 <odrrj+0x4d4>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    15e9 <odrrj+0x4d9>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1580 <odrrj+0x470>
mov    -0x84(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jge    1647 <odrrj+0x537>
mov    0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x78(%rbp),%xmm2
movsd  -0xb8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0x9c8(%rip),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    1523 <odrrj+0x413>
mov    -0xa0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa0(%rbp)
mov    -0x4c(%rbp),%edx
add    $0x1,%edx
mov    0x30(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x84(%rbp),%edx
mov    0x30(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
movsd  -0x98(%rbp),%xmm0
ff
mov    0x28(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb8(%rbp),%xmm0
ff
mov    0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1241 <odrrj+0x131>
mov    -0x64(%rbp),%eax
mov    %eax,-0xac(%rbp)
movl   $0x0,-0x9c(%rbp)
mov    -0x64(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0xa4(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0xa0(%rbp),%eax
mov    %eax,-0xcc(%rbp)
mov    -0xcc(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     19de <odrrj+0x8ce>
mov    -0x4c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x4c(%rbp),%eax
mov    0x10(%rbp),%rcx
sub    (%rcx),%eax
mov    %eax,-0x6c(%rbp)
mov    0x30(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xa8(%rbp)
mov    0x28(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
addsd  -0xb8(%rbp),%xmm0
ff
movsd  0x88f(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movaps 0x87a(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
movaps 0x862(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm3
movsd  -0x48(%rbp),%xmm2
movaps %xmm2,%xmm4
movaps %xmm3,%xmm0
cmpunordsd %xmm0,%xmm0
movaps %xmm0,%xmm1
pand   %xmm4,%xmm1
maxsd  %xmm3,%xmm2
pandn  %xmm2,%xmm0
por    %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0xc0(%rbp),%xmm1
ff
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
ja     182b <odrrj+0x71b>
mov    -0x9c(%rbp),%eax
cmp    -0xd0(%rbp),%eax
jne    18a6 <odrrj+0x796>
mov    -0xa0(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0xa0(%rbp)
mov    0x30(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x64(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jne    1868 <odrrj+0x758>
mov    -0xa8(%rbp),%eax
mov    %eax,-0x64(%rbp)
jmp    1898 <odrrj+0x788>
mov    -0xa4(%rbp),%eax
cmp    -0x64(%rbp),%eax
jl     1893 <odrrj+0x783>
mov    -0xa8(%rbp),%edx
mov    0x30(%rbp),%rax
mov    -0xa4(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
jmp    1898 <odrrj+0x788>
mov    -0xa8(%rbp),%eax
mov    %eax,-0x4c(%rbp)
jmp    19cb <odrrj+0x8bb>
mov    -0x4c(%rbp),%eax
mov    %eax,-0xa4(%rbp)
movl   $0x0,-0x84(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xc8(%rbp),%xmm0
ff
jbe    18f9 <odrrj+0x7e9>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x38(%rbp)
movl   $0x2,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     1980 <odrrj+0x870>
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
divsd  -0xc8(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xc8(%rbp),%xmm0
ff
jbe    196d <odrrj+0x85d>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    1972 <odrrj+0x862>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1909 <odrrj+0x7f9>
mov    -0x84(%rbp),%eax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     19b0 <odrrj+0x8a0>
movsd  -0x80(%rbp),%xmm0
mov    0x28(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    19c2 <odrrj+0x8b2>
movsd  -0x80(%rbp),%xmm0
mov    0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xa8(%rbp),%eax
mov    %eax,-0x4c(%rbp)
jmp    19d0 <odrrj+0x8c0>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    16fb <odrrj+0x5eb>
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
cmpl   $0x0,-0xa0(%rbp)
jle    1a11 <odrrj+0x901>
mov    -0x9c(%rbp),%eax
cmp    -0xd0(%rbp),%eax
jg     1a11 <odrrj+0x901>
jmp    16cd <odrrj+0x5bd>
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x34(%rbp)
mov    -0xac(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     1acd <odrrj+0x9bd>
mov    -0x4c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x4c(%rbp),%eax
mov    0x10(%rbp),%rcx
sub    (%rcx),%eax
mov    %eax,-0x6c(%rbp)
mov    0x30(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
jne    1aba <odrrj+0x9aa>
mov    0x28(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0x586(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1abf <odrrj+0x9af>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1a23 <odrrj+0x913>
xor    %eax,%eax
add    $0xe0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<odrrj>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9,0x20(%rsp)
mov    %r8,%rbp
mov    %rcx,(%rsp)
mov    %rdx,%r12
mov    %rsi,%rbx
mov    %rdi,0x38(%rsp)
mov    0xa0(%rsp),%r15
mov    0x98(%rsp),%r13
mov    0xa8(%rsp),%r14
mov    0xb0(%rsp),%rax
mov    0xb8(%rsp),%rcx
movl   $0x0,(%rcx)
movsd  (%rax),%xmm0
addsd  (%r14),%xmm0
call   1030 <log@plt>
movsd  %xmm0,0x18(%rsp)
movsd  (%r14),%xmm0
call   1030 <log@plt>
movsd  0x18(%rsp),%xmm1
subsd  %xmm0,%xmm1
divsd  0xe6d(%rip),%xmm1
cvttsd2si %xmm1,%eax
mov    %eax,0xc(%rsp)
mov    (%rsp),%rax
mov    (%rax),%edx
mov    %rbp,0x40(%rsp)
movslq 0x0(%rbp),%r11
xor    %eax,%eax
mov    %rax,0x18(%rsp)
cmp    %r11d,%edx
jle    144b <odrrj+0x33b>
mov    %edx,%r8d
movl   $0x0,0x8(%rsp)
mov    0x20(%rsp),%r14
addl   $0x2,0xc(%rsp)
movsd  0xe38(%rip),%xmm8
movapd 0xe40(%rip),%xmm1
xorpd  %xmm2,%xmm2
mov    %r8,0x10(%rsp)
mov    %r8d,%r10d
mov    0x8(%rsp),%r9d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
test   %r9d,%r9d
jle    141e <odrrj+0x30e>
lea    -0x1(%r10),%r8d
mov    $0x1,%esi
mov    %r9d,%edx
mov    %r10d,%eax
jmp    1240 <odrrj+0x130>
nopw   0x0(%rax,%rax,1)
movsd  %xmm3,-0x8(%r13,%rdi,8)
mov    %ecx,%r8d
mov    0x20(%rsp),%r14
mov    (%rsp),%edx
lea    0x1(%rsi),%ecx
cmp    %r9d,%esi
mov    %ecx,%esi
je     1400 <odrrj+0x2f0>
mov    %eax,%ecx
lea    (%rcx,%rcx,1),%eax
movslq %eax,%rdi
lea    -0x1(%rdi),%r11
mov    -0x8(%r15,%rdi,4),%eax
movsd  -0x10(%r13,%rdi,8),%xmm4
movsd  -0x8(%r13,%rdi,8),%xmm5
movapd %xmm4,%xmm3
addsd  %xmm5,%xmm3
mulsd  %xmm8,%xmm3
movapd %xmm5,%xmm6
subsd  %xmm3,%xmm6
andpd  %xmm1,%xmm4
andpd  %xmm1,%xmm5
movapd %xmm4,%xmm7
cmpunordsd %xmm4,%xmm7
movapd %xmm7,%xmm0
andpd  %xmm5,%xmm0
maxsd  %xmm4,%xmm5
andnpd %xmm5,%xmm7
orpd   %xmm0,%xmm7
mulsd  (%r14),%xmm7
ucomisd %xmm6,%xmm7
ja     12f0 <odrrj+0x1e0>
mov    0x18(%rsp),%rbp
cmp    0xc(%rsp),%ebp
je     12f0 <odrrj+0x1e0>
mov    %edx,(%rsp)
movsd  (%rbx),%xmm5
subsd  %xmm3,%xmm5
xor    %edx,%edx
ucomisd %xmm5,%xmm2
seta   %dl
mov    0x38(%rsp),%rbp
mov    0x0(%rbp),%r8d
cmp    $0x2,%r8d
jl     13df <odrrj+0x2cf>
add    $0x1,%r8
mov    $0x2,%r14d
cmp    $0x3,%r8d
jne    132a <odrrj+0x21a>
movapd %xmm5,%xmm4
jmp    13a3 <odrrj+0x293>
nop
add    $0xffffffff,%edx
movl   $0x0,-0x4(%r15,%r11,4)
cmp    %ecx,%r10d
jne    1310 <odrrj+0x200>
mov    %eax,%r10d
jmp    1232 <odrrj+0x122>
nopl   0x0(%rax)
cmp    %r10d,%r8d
jl     1232 <odrrj+0x122>
lea    (%r8,%r8,1),%ecx
movslq %ecx,%rcx
mov    %eax,-0x8(%r15,%rcx,4)
jmp    1232 <odrrj+0x122>
lea    -0x2(%r8),%r15
and    $0xfffffffffffffffe,%r15
neg    %r15
mov    $0x2,%r14d
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rbx,%r14,8),%xmm0
movsd  (%rbx,%r14,8),%xmm4
subsd  %xmm3,%xmm0
movsd  -0x10(%r12,%r14,8),%xmm6
divsd  %xmm5,%xmm6
subsd  %xmm6,%xmm0
xor    %r13d,%r13d
ucomisd %xmm0,%xmm2
movsd  -0x8(%r12,%r14,8),%xmm5
seta   %r13b
add    %edx,%r13d
subsd  %xmm3,%xmm4
divsd  %xmm0,%xmm5
subsd  %xmm5,%xmm4
xor    %edx,%edx
ucomisd %xmm4,%xmm2
seta   %dl
add    %r13d,%edx
lea    (%r15,%r14,1),%rbp
add    $0x2,%rbp
add    $0x2,%r14
movapd %xmm4,%xmm5
cmp    $0x2,%rbp
jne    1340 <odrrj+0x230>
test   $0x1,%r8b
mov    0xa0(%rsp),%r15
mov    0x98(%rsp),%r13
je     13df <odrrj+0x2cf>
movsd  -0x8(%rbx,%r14,8),%xmm0
subsd  %xmm3,%xmm0
movsd  -0x10(%r12,%r14,8),%xmm5
divsd  %xmm4,%xmm5
subsd  %xmm5,%xmm0
xor    %ebp,%ebp
ucomisd %xmm0,%xmm2
seta   %bpl
add    %ebp,%edx
cmp    %ecx,%edx
jge    1220 <odrrj+0x110>
movsd  %xmm3,-0x8(%r13,%r11,8)
jmp    1227 <odrrj+0x117>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
test   %edx,%edx
jle    141e <odrrj+0x30e>
mov    0x18(%rsp),%rcx
lea    0x1(%rcx),%eax
mov    %edx,%r9d
cmp    0xc(%rsp),%ecx
mov    %rax,0x18(%rsp)
jl     1200 <odrrj+0xf0>
mov    0x40(%rsp),%rax
mov    (%rax),%eax
mov    0x10(%rsp),%rcx
cmp    %eax,%ecx
mov    0x80(%rsp),%r8
jle    17e8 <odrrj+0x6d8>
xor    %eax,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
movslq %edx,%rdi
lea    0x1(%r11),%r10d
xor    %eax,%eax
mov    %rax,0x30(%rsp)
movapd 0xbbe(%rip),%xmm10
movapd 0xbc5(%rip),%xmm8
movsd  0xb94(%rip),%xmm9
xorpd  %xmm3,%xmm3
movl   $0x0,0x8(%rsp)
mov    0x80(%rsp),%r9
mov    0x20(%rsp),%r14
mov    %r10d,0x2c(%rsp)
jmp    14fe <odrrj+0x3ee>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addl   $0x1,0x8(%rsp)
mov    %rdi,%rbp
add    $0x1,%rbp
mov    (%rsp),%rcx
lea    -0x1(%rcx),%rsi
mov    %ebp,-0x8(%r15,%rcx,4)
mov    %rbp,%rcx
mov    %edx,%r8d
mov    %rdi,0x30(%rsp)
mov    0x80(%rsp),%r9
mov    0x2c(%rsp),%r10d
movslq 0x10(%rsp),%rdx
mov    %eax,-0x4(%r15,%rdx,4)
mov    %ebp,%eax
mov    %rcx,%rdi
mov    (%rsp),%rbp
movsd  %xmm5,-0x8(%r13,%rsi,8)
movsd  %xmm4,-0x8(%r13,%rbp,8)
mov    %r8d,%edx
cmp    %r10d,%eax
je     11ca <odrrj+0xba>
movslq (%r9),%rax
mov    %rdi,%rcx
sub    %rax,%rcx
mov    0x88(%rsp),%rax
movsd  -0x8(%rax,%rcx,8),%xmm5
mov    0x90(%rsp),%rax
movsd  -0x8(%rax,%rcx,8),%xmm11
movapd %xmm5,%xmm4
addsd  %xmm11,%xmm4
movapd %xmm4,%xmm0
subsd  %xmm5,%xmm0
subsd  %xmm11,%xmm5
movapd %xmm5,%xmm1
andpd  %xmm10,%xmm1
movapd %xmm4,%xmm2
andpd  %xmm10,%xmm2
movapd %xmm1,%xmm6
cmpunordsd %xmm1,%xmm6
movapd %xmm6,%xmm7
andpd  %xmm2,%xmm7
maxsd  %xmm1,%xmm2
andnpd %xmm2,%xmm6
orpd   %xmm7,%xmm6
mulsd  (%r14),%xmm6
lea    (%rdi,%rdi,1),%eax
ucomisd %xmm0,%xmm6
movslq %eax,%rbp
jbe    15c0 <odrrj+0x4b0>
lea    -0x1(%rbp),%rsi
movl   $0xffffffff,-0x8(%r15,%rbp,4)
ff ff
lea    0x1(%rdi),%rax
cmp    %r11,%rdi
mov    %edx,%r8d
cmovl  %eax,%r8d
cmp    %edi,%edx
cmovne %edx,%r8d
mov    0x30(%rsp),%rdx
cmp    %r8d,%edx
jge    17d0 <odrrj+0x6c0>
mov    %rax,%rdi
jmp    14e4 <odrrj+0x3d4>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %eax,0x10(%rsp)
mov    %rbp,(%rsp)
movsd  (%rbx),%xmm13
mov    0x38(%rsp),%rax
mov    (%rax),%esi
movapd %xmm11,%xmm12
xorpd  %xmm8,%xmm12
lea    0x1(%rsi),%r9d
lea    -0x2(%r9),%r8
lea    -0x3(%r9),%r10
and    $0xfffffffffffffffe,%r8
neg    %r8
movapd %xmm9,%xmm2
cs nopw 0x0(%rax,%rax,1)
movapd %xmm13,%xmm7
subsd  %xmm5,%xmm7
xor    %eax,%eax
ucomisd %xmm7,%xmm3
seta   %al
cmp    $0x2,%esi
jl     16b8 <odrrj+0x5a8>
mov    $0x2,%ecx
test   %r10,%r10
je     168e <odrrj+0x57e>
movapd %xmm7,%xmm1
nopl   0x0(%rax)
movsd  -0x8(%rbx,%rcx,8),%xmm6
movsd  (%rbx,%rcx,8),%xmm7
subsd  %xmm5,%xmm6
movsd  -0x10(%r12,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
subsd  %xmm0,%xmm6
xor    %ebp,%ebp
ucomisd %xmm6,%xmm3
movsd  -0x8(%r12,%rcx,8),%xmm0
seta   %bpl
add    %eax,%ebp
subsd  %xmm5,%xmm7
divsd  %xmm6,%xmm0
subsd  %xmm0,%xmm7
xor    %eax,%eax
ucomisd %xmm7,%xmm3
seta   %al
add    %ebp,%eax
lea    (%r8,%rcx,1),%rbp
add    $0x2,%rbp
add    $0x2,%rcx
movapd %xmm7,%xmm1
cmp    $0x2,%rbp
jne    1630 <odrrj+0x520>
test   $0x1,%r9b
je     16b8 <odrrj+0x5a8>
movsd  -0x8(%rbx,%rcx,8),%xmm0
subsd  %xmm5,%xmm0
movsd  -0x10(%r12,%rcx,8),%xmm1
divsd  %xmm7,%xmm1
subsd  %xmm1,%xmm0
xor    %ecx,%ecx
ucomisd %xmm0,%xmm3
seta   %cl
add    %ecx,%eax
cltq
cmp    %rax,%rdi
jg     16e0 <odrrj+0x5d0>
movapd %xmm12,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm0,%xmm5
addsd  %xmm2,%xmm2
jmp    1600 <odrrj+0x4f0>
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm9,%xmm2
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm13,%xmm6
subsd  %xmm4,%xmm6
xor    %eax,%eax
ucomisd %xmm6,%xmm3
seta   %al
cmp    $0x2,%esi
jl     17a8 <odrrj+0x698>
mov    $0x2,%ecx
test   %r10,%r10
je     177e <odrrj+0x66e>
movapd %xmm6,%xmm1
nopl   0x0(%rax)
movsd  -0x8(%rbx,%rcx,8),%xmm0
movsd  (%rbx,%rcx,8),%xmm6
subsd  %xmm4,%xmm0
movsd  -0x10(%r12,%rcx,8),%xmm7
divsd  %xmm1,%xmm7
subsd  %xmm7,%xmm0
xor    %ebp,%ebp
ucomisd %xmm0,%xmm3
movsd  -0x8(%r12,%rcx,8),%xmm1
seta   %bpl
add    %eax,%ebp
subsd  %xmm4,%xmm6
divsd  %xmm0,%xmm1
subsd  %xmm1,%xmm6
xor    %eax,%eax
ucomisd %xmm6,%xmm3
seta   %al
add    %ebp,%eax
lea    (%r8,%rcx,1),%rbp
add    $0x2,%rbp
add    $0x2,%rcx
movapd %xmm6,%xmm1
cmp    $0x2,%rbp
jne    1720 <odrrj+0x610>
test   $0x1,%r9b
je     17a8 <odrrj+0x698>
movsd  -0x8(%rbx,%rcx,8),%xmm0
subsd  %xmm4,%xmm0
movsd  -0x10(%r12,%rcx,8),%xmm1
divsd  %xmm6,%xmm1
subsd  %xmm1,%xmm0
xor    %ecx,%ecx
ucomisd %xmm0,%xmm3
seta   %cl
add    %ecx,%eax
movslq %eax,%rcx
cmp    %rcx,%rdi
jle    14a0 <odrrj+0x390>
movapd %xmm11,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm0,%xmm4
addsd  %xmm2,%xmm2
jmp    16f0 <odrrj+0x5e0>
nopw   0x0(%rax,%rax,1)
mov    %rbp,(%rsp)
lea    (%rdx,%rdx,1),%ecx
add    $0xffffffff,%ecx
mov    %ecx,0x10(%rsp)
mov    %eax,%ebp
mov    %rax,%rcx
jmp    14d1 <odrrj+0x3c1>
movslq %ecx,%rdx
mov    0x88(%rsp),%rcx
lea    (%rcx,%rdx,8),%rbp
add    $0xfffffffffffffff8,%rbp
mov    0x90(%rsp),%rcx
lea    (%rcx,%rdx,8),%rdx
add    $0xfffffffffffffff8,%rdx
mov    0x10(%rsp),%rcx
sub    %ecx,%eax
add    %ecx,%ecx
add    $0x1,%eax
jmp    1834 <odrrj+0x724>
nopl   0x0(%rax)
add    $0x8,%rbp
add    $0x8,%rdx
add    $0x2,%ecx
add    $0xffffffff,%eax
je     143a <odrrj+0x32a>
lea    -0x1(%rcx),%esi
movslq %esi,%rsi
cmpl   $0x0,-0x4(%r15,%rsi,4)
jne    1820 <odrrj+0x710>
movslq (%r8),%rdi
movsd  -0x8(%r13,%rsi,8),%xmm0
movslq %ecx,%rsi
addsd  -0x8(%r13,%rsi,8),%xmm0
mulsd  %xmm8,%xmm0
shl    $0x3,%rdi
mov    %rcx,%rbx
mov    %rbp,%rcx
sub    %rdi,%rcx
movsd  %xmm0,(%rcx)
movsd  -0x8(%r13,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
mov    %rdx,%rcx
sub    %rdi,%rcx
movsd  %xmm1,(%rcx)
mov    %rbx,%rcx
jmp    1820 <odrrj+0x710>

LLVM IR (-O0):
define i32 @odrrj(i32* %0, double* %1, double* %2, i32* %3, i32* %4, double* %5, i32* %6, double* %7, double* %8, double* %9, i32* %10, double* %11, double* %12, i32* %13) {
%15 = alloca i32*
%16 = alloca double*
%17 = alloca double*
%18 = alloca i32*
%19 = alloca i32*
%20 = alloca double*
%21 = alloca i32*
%22 = alloca double*
%23 = alloca double*
%24 = alloca double*
%25 = alloca i32*
%26 = alloca double*
%27 = alloca double*
%28 = alloca i32*
%29 = alloca i32
%30 = alloca i32
%31 = alloca double
%32 = alloca double
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
%36 = alloca i32
%37 = alloca double
%38 = alloca i32
%39 = alloca i32
%40 = alloca i32
%41 = alloca double
%42 = alloca double
%43 = alloca i32
%44 = alloca double
%45 = alloca double
%46 = alloca i32
%47 = alloca i32
%48 = alloca i32
%49 = alloca i32
%50 = alloca i32
%51 = alloca double
%52 = alloca double
%53 = alloca double
%54 = alloca i32
%55 = alloca i32
store i32* %0, i32** %15
store double* %1, double** %16
store double* %2, double** %17
store i32* %3, i32** %18
store i32* %4, i32** %19
store double* %5, double** %20
store i32* %6, i32** %21
store double* %7, double** %22
store double* %8, double** %23
store double* %9, double** %24
store i32* %10, i32** %25
store double* %11, double** %26
store double* %12, double** %27
store i32* %13, i32** %28
%56 = load i32*, i32** %25
%57 = getelementptr inbounds i32, i32* %56, i32 -1
store i32* %57, i32** %25
%58 = load double*, double** %24
%59 = getelementptr inbounds double, double* %58, i32 -1
store double* %59, double** %24
%60 = load double*, double** %23
%61 = getelementptr inbounds double, double* %60, i32 -1
store double* %61, double** %23
%62 = load double*, double** %22
%63 = getelementptr inbounds double, double* %62, i32 -1
store double* %63, double** %22
%64 = load double*, double** %17
%65 = getelementptr inbounds double, double* %64, i32 -1
store double* %65, double** %17
%66 = load double*, double** %16
%67 = getelementptr inbounds double, double* %66, i32 -1
store double* %67, double** %16
%68 = load i32*, i32** %28
store i32 0, i32* %68
%69 = load double*, double** %27
%70 = load double, double* %69
%71 = load double*, double** %26
%72 = load double, double* %71
%73 = fadd double %70, %72
%74 = call double @log(double %73)
%75 = load double*, double** %26
%76 = load double, double* %75
%77 = call double @log(double %76)
%78 = fsub double %74, %77
%79 = call double @log(double 2.000000e+00)
%80 = fdiv double %78, %79
%81 = fptosi double %80 to i32
%82 = add nsw i32 %81, 2
store i32 %82, i32* %55
%83 = load i32*, i32** %18
%84 = load i32, i32* %83
store i32 %84, i32* %38
%85 = load i32*, i32** %19
%86 = load i32, i32* %85
store i32 %86, i32* %39
store i32 0, i32* %47
store i32 0, i32* %48
%87 = load i32, i32* %39
store i32 %87, i32* %29
%88 = load i32, i32* %38
store i32 %88, i32* %33
br label %89
89: ; preds = %329, %14
%90 = load i32, i32* %33
%91 = load i32, i32* %29
%92 = icmp sle i32 %90, %91
br i1 %92, label %93, label %332
93: ; preds = %89
%94 = load i32, i32* %33
%95 = shl i32 %94, 1
store i32 %95, i32* %35
%96 = load i32, i32* %33
%97 = load i32*, i32** %21
%98 = load i32, i32* %97
%99 = sub nsw i32 %96, %98
store i32 %99, i32* %40
%100 = load double*, double** %22
%101 = load i32, i32* %40
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %100, i64 %102
%104 = load double, double* %103
%105 = load double*, double** %23
%106 = load i32, i32* %40
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds double, double* %105, i64 %107
%109 = load double, double* %108
%110 = fsub double %104, %109
store double %110, double* %45
%111 = load double*, double** %22
%112 = load i32, i32* %40
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %111, i64 %113
%115 = load double, double* %114
store double %115, double* %42
%116 = load double*, double** %22
%117 = load i32, i32* %40
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds double, double* %116, i64 %118
%120 = load double, double* %119
%121 = load double*, double** %23
%122 = load i32, i32* %40
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double, double* %121, i64 %123
%125 = load double, double* %124
%126 = fadd double %120, %125
store double %126, double* %51
%127 = load double, double* %51
%128 = load double, double* %42
%129 = fsub double %127, %128
store double %129, double* %52
%130 = load double, double* %45
%131 = call double @llvm.fabs.f64(double %130)
store double %131, double* %31
%132 = load double, double* %51
%133 = call double @llvm.fabs.f64(double %132)
store double %133, double* %32
%134 = load double, double* %31
%135 = load double, double* %32
%136 = call double @llvm.maxnum.f64(double %134, double %135)
store double %136, double* %44
%137 = load double, double* %52
%138 = load double*, double** %20
%139 = load double, double* %138
%140 = load double, double* %44
%141 = fmul double %139, %140
%142 = fcmp olt double %137, %141
br i1 %142, label %143, label %177
143: ; preds = %93
%144 = load i32*, i32** %25
%145 = load i32, i32* %35
%146 = sub nsw i32 %145, 1
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds i32, i32* %144, i64 %147
store i32 -1, i32* %148
%149 = load i32, i32* %33
%150 = load i32, i32* %38
%151 = icmp eq i32 %149, %150
br i1 %151, label %152, label %159
152: ; preds = %143
%153 = load i32, i32* %33
%154 = load i32, i32* %39
%155 = icmp slt i32 %153, %154
br i1 %155, label %156, label %159
156: ; preds = %152
%157 = load i32, i32* %33
%158 = add nsw i32 %157, 1
store i32 %158, i32* %38
br label %159
159: ; preds = %156, %152, %143
%160 = load i32, i32* %48
%161 = load i32, i32* %38
%162 = icmp sge i32 %160, %161
br i1 %162, label %163, label %176
163: ; preds = %159
%164 = load i32, i32* %33
%165 = load i32, i32* %39
%166 = icmp sle i32 %164, %165
br i1 %166, label %167, label %176
167: ; preds = %163
%168 = load i32, i32* %33
%169 = add nsw i32 %168, 1
%170 = load i32*, i32** %25
%171 = load i32, i32* %48
%172 = shl i32 %171, 1
%173 = sub nsw i32 %172, 1
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds i32, i32* %170, i64 %174
store i32 %169, i32* %175
br label %176
176: ; preds = %167, %163, %159
br label %317
177: ; preds = %93
%178 = load i32, i32* %33
store i32 %178, i32* %48
store double 1.000000e+00, double* %41
br label %179
179: ; preds = %229, %177
store i32 0, i32* %43
%180 = load double, double* %45
store double %180, double* %37
%181 = load double*, double** %16
%182 = getelementptr inbounds double, double* %181, i64 1
%183 = load double, double* %182
%184 = load double, double* %37
%185 = fsub double %183, %184
store double %185, double* %53
%186 = load double, double* %53
%187 = fcmp olt double %186, 0.000000e+00
br i1 %187, label %188, label %191
188: ; preds = %179
%189 = load i32, i32* %43
%190 = add nsw i32 %189, 1
store i32 %190, i32* %43
br label %191
191: ; preds = %188, %179
%192 = load i32*, i32** %15
%193 = load i32, i32* %192
store i32 %193, i32* %30
store i32 2, i32* %34
br label %194
194: ; preds = %221, %191
%195 = load i32, i32* %34
%196 = load i32, i32* %30
%197 = icmp sle i32 %195, %196
br i1 %197, label %198, label %224
198: ; preds = %194
%199 = load double*, double** %16
%200 = load i32, i32* %34
%201 = sext i32 %200 to i64
%202 = getelementptr inbounds double, double* %199, i64 %201
%203 = load double, double* %202
%204 = load double, double* %37
%205 = fsub double %203, %204
%206 = load double*, double** %17
%207 = load i32, i32* %34
%208 = sub nsw i32 %207, 1
%209 = sext i32 %208 to i64
%210 = getelementptr inbounds double, double* %206, i64 %209
%211 = load double, double* %210
%212 = load double, double* %53
%213 = fdiv double %211, %212
%214 = fsub double %205, %213
store double %214, double* %53
%215 = load double, double* %53
%216 = fcmp olt double %215, 0.000000e+00
br i1 %216, label %217, label %220
217: ; preds = %198
%218 = load i32, i32* %43
%219 = add nsw i32 %218, 1
store i32 %219, i32* %43
br label %220
220: ; preds = %217, %198
br label %221
221: ; preds = %220
%222 = load i32, i32* %34
%223 = add nsw i32 %222, 1
store i32 %223, i32* %34
br label %194
224: ; preds = %194
%225 = load i32, i32* %43
%226 = load i32, i32* %33
%227 = sub nsw i32 %226, 1
%228 = icmp sgt i32 %225, %227
br i1 %228, label %229, label %241
229: ; preds = %224
%230 = load double*, double** %23
%231 = load i32, i32* %40
%232 = sext i32 %231 to i64
%233 = getelementptr inbounds double, double* %230, i64 %232
%234 = load double, double* %233
%235 = load double, double* %41
%236 = load double, double* %45
%237 = fneg double %234
%238 = call double @llvm.fmuladd.f64(double %237, double %235, double %236)
store double %238, double* %45
%239 = load double, double* %41
%240 = fmul double %239, 2.000000e+00
store double %240, double* %41
br label %179
241: ; preds = %224
store double 1.000000e+00, double* %41
br label %242
242: ; preds = %291, %241
store i32 0, i32* %43
%243 = load double, double* %51
store double %243, double* %37
%244 = load double*, double** %16
%245 = getelementptr inbounds double, double* %244, i64 1
%246 = load double, double* %245
%247 = load double, double* %37
%248 = fsub double %246, %247
store double %248, double* %53
%249 = load double, double* %53
%250 = fcmp olt double %249, 0.000000e+00
br i1 %250, label %251, label %254
251: ; preds = %242
%252 = load i32, i32* %43
%253 = add nsw i32 %252, 1
store i32 %253, i32* %43
br label %254
254: ; preds = %251, %242
%255 = load i32*, i32** %15
%256 = load i32, i32* %255
store i32 %256, i32* %30
store i32 2, i32* %34
br label %257
257: ; preds = %284, %254
%258 = load i32, i32* %34
%259 = load i32, i32* %30
%260 = icmp sle i32 %258, %259
br i1 %260, label %261, label %287
261: ; preds = %257
%262 = load double*, double** %16
%263 = load i32, i32* %34
%264 = sext i32 %263 to i64
%265 = getelementptr inbounds double, double* %262, i64 %264
%266 = load double, double* %265
%267 = load double, double* %37
%268 = fsub double %266, %267
%269 = load double*, double** %17
%270 = load i32, i32* %34
%271 = sub nsw i32 %270, 1
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds double, double* %269, i64 %272
%274 = load double, double* %273
%275 = load double, double* %53
%276 = fdiv double %274, %275
%277 = fsub double %268, %276
store double %277, double* %53
%278 = load double, double* %53
%279 = fcmp olt double %278, 0.000000e+00
br i1 %279, label %280, label %283
280: ; preds = %261
%281 = load i32, i32* %43
%282 = add nsw i32 %281, 1
store i32 %282, i32* %43
br label %283
283: ; preds = %280, %261
br label %284
284: ; preds = %283
%285 = load i32, i32* %34
%286 = add nsw i32 %285, 1
store i32 %286, i32* %34
br label %257
287: ; preds = %257
%288 = load i32, i32* %43
%289 = load i32, i32* %33
%290 = icmp slt i32 %288, %289
br i1 %290, label %291, label %302
291: ; preds = %287
%292 = load double*, double** %23
%293 = load i32, i32* %40
%294 = sext i32 %293 to i64
%295 = getelementptr inbounds double, double* %292, i64 %294
%296 = load double, double* %295
%297 = load double, double* %41
%298 = load double, double* %51
%299 = call double @llvm.fmuladd.f64(double %296, double %297, double %298)
store double %299, double* %51
%300 = load double, double* %41
%301 = fmul double %300, 2.000000e+00
store double %301, double* %41
br label %242
302: ; preds = %287
%303 = load i32, i32* %47
%304 = add nsw i32 %303, 1
store i32 %304, i32* %47
%305 = load i32, i32* %33
%306 = add nsw i32 %305, 1
%307 = load i32*, i32** %25
%308 = load i32, i32* %35
%309 = sub nsw i32 %308, 1
%310 = sext i32 %309 to i64
%311 = getelementptr inbounds i32, i32* %307, i64 %310
store i32 %306, i32* %311
%312 = load i32, i32* %43
%313 = load i32*, i32** %25
%314 = load i32, i32* %35
%315 = sext i32 %314 to i64
%316 = getelementptr inbounds i32, i32* %313, i64 %315
store i32 %312, i32* %316
br label %317
317: ; preds = %302, %176
%318 = load double, double* %45
%319 = load double*, double** %24
%320 = load i32, i32* %35
%321 = sub nsw i32 %320, 1
%322 = sext i32 %321 to i64
%323 = getelementptr inbounds double, double* %319, i64 %322
store double %318, double* %323
%324 = load double, double* %51
%325 = load double*, double** %24
%326 = load i32, i32* %35
%327 = sext i32 %326 to i64
%328 = getelementptr inbounds double, double* %325, i64 %327
store double %324, double* %328
br label %329
329: ; preds = %317
%330 = load i32, i32* %33
%331 = add nsw i32 %330, 1
store i32 %331, i32* %33
br label %89
332: ; preds = %89
%333 = load i32, i32* %38
store i32 %333, i32* %50
store i32 0, i32* %46
br label %334
334: ; preds = %499, %332
%335 = load i32, i32* %38
%336 = sub nsw i32 %335, 1
store i32 %336, i32* %48
%337 = load i32, i32* %38
store i32 %337, i32* %33
%338 = load i32, i32* %47
store i32 %338, i32* %54
%339 = load i32, i32* %54
store i32 %339, i32* %29
store i32 1, i32* %36
br label %340
340: ; preds = %487, %334
%341 = load i32, i32* %36
%342 = load i32, i32* %29
%343 = icmp sle i32 %341, %342
br i1 %343, label %344, label %490
344: ; preds = %340
%345 = load i32, i32* %33
%346 = shl i32 %345, 1
store i32 %346, i32* %35
%347 = load i32, i32* %33
%348 = load i32*, i32** %21
%349 = load i32, i32* %348
%350 = sub nsw i32 %347, %349
store i32 %350, i32* %40
%351 = load i32*, i32** %25
%352 = load i32, i32* %35
%353 = sub nsw i32 %352, 1
%354 = sext i32 %353 to i64
%355 = getelementptr inbounds i32, i32* %351, i64 %354
%356 = load i32, i32* %355
store i32 %356, i32* %49
%357 = load double*, double** %24
%358 = load i32, i32* %35
%359 = sub nsw i32 %358, 1
%360 = sext i32 %359 to i64
%361 = getelementptr inbounds double, double* %357, i64 %360
%362 = load double, double* %361
store double %362, double* %45
%363 = load double*, double** %24
%364 = load i32, i32* %35
%365 = sext i32 %364 to i64
%366 = getelementptr inbounds double, double* %363, i64 %365
%367 = load double, double* %366
store double %367, double* %51
%368 = load double, double* %45
%369 = load double, double* %51
%370 = fadd double %368, %369
%371 = fmul double %370, 5.000000e-01
store double %371, double* %42
%372 = load double, double* %51
%373 = load double, double* %42
%374 = fsub double %372, %373
store double %374, double* %52
%375 = load double, double* %45
%376 = call double @llvm.fabs.f64(double %375)
store double %376, double* %31
%377 = load double, double* %51
%378 = call double @llvm.fabs.f64(double %377)
store double %378, double* %32
%379 = load double, double* %31
%380 = load double, double* %32
%381 = call double @llvm.maxnum.f64(double %379, double %380)
store double %381, double* %44
%382 = load double, double* %52
%383 = load double*, double** %20
%384 = load double, double* %383
%385 = load double, double* %44
%386 = fmul double %384, %385
%387 = fcmp olt double %382, %386
br i1 %387, label %392, label %388
388: ; preds = %344
%389 = load i32, i32* %46
%390 = load i32, i32* %55
%391 = icmp eq i32 %389, %390
br i1 %391, label %392, label %420
392: ; preds = %388, %344
%393 = load i32, i32* %47
%394 = add nsw i32 %393, -1
store i32 %394, i32* %47
%395 = load i32*, i32** %25
%396 = load i32, i32* %35
%397 = sub nsw i32 %396, 1
%398 = sext i32 %397 to i64
%399 = getelementptr inbounds i32, i32* %395, i64 %398
store i32 0, i32* %399
%400 = load i32, i32* %38
%401 = load i32, i32* %33
%402 = icmp eq i32 %400, %401
br i1 %402, label %403, label %405
403: ; preds = %392
%404 = load i32, i32* %49
store i32 %404, i32* %38
br label %418
405: ; preds = %392
%406 = load i32, i32* %48
%407 = load i32, i32* %38
%408 = icmp sge i32 %406, %407
br i1 %408, label %409, label %417
409: ; preds = %405
%410 = load i32, i32* %49
%411 = load i32*, i32** %25
%412 = load i32, i32* %48
%413 = shl i32 %412, 1
%414 = sub nsw i32 %413, 1
%415 = sext i32 %414 to i64
%416 = getelementptr inbounds i32, i32* %411, i64 %415
store i32 %410, i32* %416
br label %417
417: ; preds = %409, %405
br label %418
418: ; preds = %417, %403
%419 = load i32, i32* %49
store i32 %419, i32* %33
br label %486
420: ; preds = %388
%421 = load i32, i32* %33
store i32 %421, i32* %48
store i32 0, i32* %43
%422 = load double, double* %42
store double %422, double* %37
%423 = load double*, double** %16
%424 = getelementptr inbounds double, double* %423, i64 1
%425 = load double, double* %424
%426 = load double, double* %37
%427 = fsub double %425, %426
store double %427, double* %53
%428 = load double, double* %53
%429 = fcmp olt double %428, 0.000000e+00
br i1 %429, label %430, label %433
430: ; preds = %420
%431 = load i32, i32* %43
%432 = add nsw i32 %431, 1
store i32 %432, i32* %43
br label %433
433: ; preds = %430, %420
%434 = load i32*, i32** %15
%435 = load i32, i32* %434
store i32 %435, i32* %30
store i32 2, i32* %34
br label %436
436: ; preds = %463, %433
%437 = load i32, i32* %34
%438 = load i32, i32* %30
%439 = icmp sle i32 %437, %438
br i1 %439, label %440, label %466
440: ; preds = %436
%441 = load double*, double** %16
%442 = load i32, i32* %34
%443 = sext i32 %442 to i64
%444 = getelementptr inbounds double, double* %441, i64 %443
%445 = load double, double* %444
%446 = load double, double* %37
%447 = fsub double %445, %446
%448 = load double*, double** %17
%449 = load i32, i32* %34
%450 = sub nsw i32 %449, 1
%451 = sext i32 %450 to i64
%452 = getelementptr inbounds double, double* %448, i64 %451
%453 = load double, double* %452
%454 = load double, double* %53
%455 = fdiv double %453, %454
%456 = fsub double %447, %455
store double %456, double* %53
%457 = load double, double* %53
%458 = fcmp olt double %457, 0.000000e+00
br i1 %458, label %459, label %462
459: ; preds = %440
%460 = load i32, i32* %43
%461 = add nsw i32 %460, 1
store i32 %461, i32* %43
br label %462
462: ; preds = %459, %440
br label %463
463: ; preds = %462
%464 = load i32, i32* %34
%465 = add nsw i32 %464, 1
store i32 %465, i32* %34
br label %436
466: ; preds = %436
%467 = load i32, i32* %43
%468 = load i32, i32* %33
%469 = sub nsw i32 %468, 1
%470 = icmp sle i32 %467, %469
br i1 %470, label %471, label %478
471: ; preds = %466
%472 = load double, double* %42
%473 = load double*, double** %24
%474 = load i32, i32* %35
%475 = sub nsw i32 %474, 1
%476 = sext i32 %475 to i64
%477 = getelementptr inbounds double, double* %473, i64 %476
store double %472, double* %477
br label %484
478: ; preds = %466
%479 = load double, double* %42
%480 = load double*, double** %24
%481 = load i32, i32* %35
%482 = sext i32 %481 to i64
%483 = getelementptr inbounds double, double* %480, i64 %482
store double %479, double* %483
br label %484
484: ; preds = %478, %471
%485 = load i32, i32* %49
store i32 %485, i32* %33
br label %486
486: ; preds = %484, %418
br label %487
487: ; preds = %486
%488 = load i32, i32* %36
%489 = add nsw i32 %488, 1
store i32 %489, i32* %36
br label %340
490: ; preds = %340
%491 = load i32, i32* %46
%492 = add nsw i32 %491, 1
store i32 %492, i32* %46
%493 = load i32, i32* %47
%494 = icmp sgt i32 %493, 0
br i1 %494, label %495, label %500
495: ; preds = %490
%496 = load i32, i32* %46
%497 = load i32, i32* %55
%498 = icmp sle i32 %496, %497
br i1 %498, label %499, label %500
499: ; preds = %495
br label %334
500: ; preds = %495, %490
%501 = load i32*, i32** %19
%502 = load i32, i32* %501
store i32 %502, i32* %29
%503 = load i32, i32* %50
store i32 %503, i32* %33
br label %504
504: ; preds = %556, %500
%505 = load i32, i32* %33
%506 = load i32, i32* %29
%507 = icmp sle i32 %505, %506
br i1 %507, label %508, label %559
508: ; preds = %504
%509 = load i32, i32* %33
%510 = shl i32 %509, 1
store i32 %510, i32* %35
%511 = load i32, i32* %33
%512 = load i32*, i32** %21
%513 = load i32, i32* %512
%514 = sub nsw i32 %511, %513
store i32 %514, i32* %40
%515 = load i32*, i32** %25
%516 = load i32, i32* %35
%517 = sub nsw i32 %516, 1
%518 = sext i32 %517 to i64
%519 = getelementptr inbounds i32, i32* %515, i64 %518
%520 = load i32, i32* %519
%521 = icmp eq i32 %520, 0
br i1 %521, label %522, label %555
522: ; preds = %508
%523 = load double*, double** %24
%524 = load i32, i32* %35
%525 = sub nsw i32 %524, 1
%526 = sext i32 %525 to i64
%527 = getelementptr inbounds double, double* %523, i64 %526
%528 = load double, double* %527
%529 = load double*, double** %24
%530 = load i32, i32* %35
%531 = sext i32 %530 to i64
%532 = getelementptr inbounds double, double* %529, i64 %531
%533 = load double, double* %532
%534 = fadd double %528, %533
%535 = fmul double %534, 5.000000e-01
%536 = load double*, double** %22
%537 = load i32, i32* %40
%538 = sext i32 %537 to i64
%539 = getelementptr inbounds double, double* %536, i64 %538
store double %535, double* %539
%540 = load double*, double** %24
%541 = load i32, i32* %35
%542 = sext i32 %541 to i64
%543 = getelementptr inbounds double, double* %540, i64 %542
%544 = load double, double* %543
%545 = load double*, double** %22
%546 = load i32, i32* %40
%547 = sext i32 %546 to i64
%548 = getelementptr inbounds double, double* %545, i64 %547
%549 = load double, double* %548
%550 = fsub double %544, %549
%551 = load double*, double** %23
%552 = load i32, i32* %40
%553 = sext i32 %552 to i64
%554 = getelementptr inbounds double, double* %551, i64 %553
store double %550, double* %554
br label %555
555: ; preds = %522, %508
br label %556
556: ; preds = %555
%557 = load i32, i32* %33
%558 = add nsw i32 %557, 1
store i32 %558, i32* %33
br label %504
559: ; preds = %504
ret i32 0
}

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