double hinit(unsigned n, double x, double *y, double *pars, double posneg,
             double *f0, double *f1, double *yy1, int iord, double hmax,
             double *atoler, double *rtoler, int itoler) {
  double dnf, dny, atoli, rtoli, sk, h, h1, der2, der12, sqr;
  unsigned i;
  dnf = 0.0;
  dny = 0.0;
  atoli = atoler[0];
  rtoli = rtoler[0];
  if (!itoler)
    for (i = 0; i < n; i++) {
      sk = atoli + rtoli * fabs(y[i]);
      sqr = f0[i] / sk;
      dnf += sqr * sqr;
      sqr = y[i] / sk;
      dny += sqr * sqr;
    }
  else
    for (i = 0; i < n; i++) {
      sk = atoler[i] + rtoler[i] * fabs(y[i]);
      sqr = f0[i] / sk;
      dnf += sqr * sqr;
      sqr = y[i] / sk;
      dny += sqr * sqr;
    }
  if ((dnf <= 1.0E-10) || (dny <= 1.0E-10))
    h = 1.0E-6;
  else
    h = sqrt(dny / dnf) * 0.01;
  h = min_d(h, hmax);
  h = sign(h, posneg);
  for (i = 0; i < n; i++)
    yy1[i] = y[i] + h * f0[i];
  safeOdeFunc(n, x + h, yy1, pars, f1);
  der2 = 0.0;
  if (!itoler)
    for (i = 0; i < n; i++) {
      sk = atoli + rtoli * fabs(y[i]);
      sqr = (f1[i] - f0[i]) / sk;
      der2 += sqr * sqr;
    }
  else
    for (i = 0; i < n; i++) {
      sk = atoler[i] + rtoler[i] * fabs(y[i]);
      sqr = (f1[i] - f0[i]) / sk;
      der2 += sqr * sqr;
    }
  der2 = sqrt(der2) / h;
  der12 = max_d(fabs(der2), sqrt(dnf));
  if (der12 <= 1.0E-15)
    h1 = max_d(1.0E-6, fabs(h) * 1.0E-3);
  else
    h1 = pow(0.01 / der12, 1.0 / (double)iord);
  h = min_d(100.0 * h, min_d(h1, hmax));
  if (isnan(h)) {
    __asm__("int $3");
  };
  return sign(h, posneg);
}


x86 assembly (clang, -O0):
<hinit>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
movsd  %xmm2,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
cmpl   $0x0,0x28(%rbp)
jne    12c2 <hinit+0x152>
movl   $0x0,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    12bd <hinit+0x14d>
movsd  -0x60(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
movaps 0xde4(%rip),%xmm3
pand   %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
jmp    11ed <hinit+0x7d>
jmp    13b5 <hinit+0x245>
movl   $0x0,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    13b0 <hinit+0x240>
mov    0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm1
mov    0x20(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
movaps 0xcf1(%rip),%xmm3
pand   %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
jmp    12cc <hinit+0x15c>
jmp    13b5 <hinit+0x245>
movsd  0xc53(%rip),%xmm0
ucomisd -0x50(%rbp),%xmm0
jae    13db <hinit+0x26b>
movsd  0xc40(%rip),%xmm0
ucomisd -0x58(%rbp),%xmm0
jb     13ed <hinit+0x27d>
movsd  0xc3d(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    140d <hinit+0x29d>
movsd  -0x58(%rbp),%xmm0
divsd  -0x50(%rbp),%xmm0
call   1090 <sqrt@plt>
movsd  0xc14(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mov    $0x2,%al
call   1030 <min_d@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
mov    $0x2,%al
call   1050 <sign@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    14a8 <hinit+0x338>
mov    -0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x78(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
jmp    144b <hinit+0x2db>
mov    -0x4(%rbp),%edi
movsd  -0x10(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
mov    -0x40(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x38(%rbp),%rcx
mov    $0x1,%al
call   1040 <safeOdeFunc@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
cmpl   $0x0,0x28(%rbp)
jne    1593 <hinit+0x423>
movl   $0x0,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    158e <hinit+0x41e>
movsd  -0x60(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
movaps 0xaea(%rip),%xmm3
pand   %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x9c(%rbp),%ecx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
jmp    14e7 <hinit+0x377>
jmp    165d <hinit+0x4ed>
movl   $0x0,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    1658 <hinit+0x4e8>
mov    0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm1
mov    0x20(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
movaps 0xa20(%rip),%xmm3
pand   %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x9c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x9c(%rbp),%ecx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
jmp    159d <hinit+0x42d>
jmp    165d <hinit+0x4ed>
movsd  -0x88(%rbp),%xmm0
ff
call   1090 <sqrt@plt>
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movaps 0x97a(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
call   1090 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xa8(%rbp),%xmm0
ff
mov    $0x2,%al
call   1070 <max_d@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0x966(%rip),%xmm0
ucomisd -0x90(%rbp),%xmm0
ff
jb     1709 <hinit+0x599>
movsd  -0x78(%rbp),%xmm1
movaps 0x924(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0x950(%rip),%xmm0
mulsd  %xmm0,%xmm1
movsd  0x92c(%rip),%xmm0
mov    $0x2,%al
call   1070 <max_d@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1734 <hinit+0x5c4>
movsd  0x907(%rip),%xmm0
divsd  -0x90(%rbp),%xmm0
ff
cvtsi2sdl 0x10(%rbp),%xmm2
movsd  0x90a(%rip),%xmm1
divsd  %xmm2,%xmm1
call   1080 <pow@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  0x904(%rip),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mov    $0x2,%al
call   1030 <min_d@plt>
movsd  -0xb0(%rbp),%xmm0
ff
mov    %eax,%edi
mov    $0x1,%al
call   1030 <min_d@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
mov    $0x1,%al
call   1060 <isnan@plt>
cmp    $0x0,%eax
je     178a <hinit+0x61a>
int3
movsd  -0x78(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
mov    $0x2,%al
call   1050 <sign@plt>
cvtsi2sd %eax,%xmm0
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<hinit>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
movapd %xmm2,%xmm8
mov    %r9,%rbp
mov    %r8,%r14
mov    %rcx,%r15
movsd  %xmm1,0x40(%rsp)
mov    %rsi,%rbx
movapd %xmm0,%xmm9
mov    0xa0(%rsp),%r8
mov    0x98(%rsp),%rsi
movsd  (%rsi),%xmm7
movsd  (%r8),%xmm6
mov    %edi,%r13d
cmpl   $0x0,0xa8(%rsp)
movsd  %xmm2,0x38(%rsp)
movsd  %xmm7,0x50(%rsp)
movsd  %xmm6,0x48(%rsp)
mov    %rdx,0x18(%rsp)
movsd  %xmm0,0x30(%rsp)
je     1249 <hinit+0xd9>
test   %edi,%edi
je     12e4 <hinit+0x174>
lea    -0x1(%r13),%rax
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
movapd 0xe39(%rip),%xmm2
movapd %xmm6,%xmm4
movapd %xmm7,%xmm3
nop
movsd  (%rbx,%rcx,8),%xmm0
movapd %xmm0,%xmm5
andpd  %xmm2,%xmm5
mulsd  %xmm4,%xmm5
movhpd (%r15,%rcx,8),%xmm0
addsd  %xmm3,%xmm5
unpcklpd %xmm5,%xmm5
divpd  %xmm5,%xmm0
mulpd  %xmm0,%xmm0
cmp    %rcx,%rax
je     1243 <hinit+0xd3>
addpd  %xmm0,%xmm1
movsd  0x8(%rsi,%rcx,8),%xmm3
movsd  0x8(%r8,%rcx,8),%xmm4
add    $0x1,%rcx
jmp    1200 <hinit+0x90>
addpd  %xmm1,%xmm0
jmp    12a0 <hinit+0x130>
test   %edi,%edi
je     12e4 <hinit+0x174>
xorpd  %xmm3,%xmm3
xor    %eax,%eax
movapd 0xdd1(%rip),%xmm1
nop
movapd %xmm3,%xmm2
movsd  (%rbx,%rax,8),%xmm0
movapd %xmm0,%xmm3
andpd  %xmm1,%xmm3
mulsd  %xmm6,%xmm3
addsd  %xmm7,%xmm3
movhpd (%r15,%rax,8),%xmm0
unpcklpd %xmm3,%xmm3
divpd  %xmm3,%xmm0
mulpd  %xmm0,%xmm0
movapd %xmm0,%xmm3
addpd  %xmm2,%xmm3
add    $0x1,%rax
cmp    %rax,%r13
jne    1260 <hinit+0xf0>
addpd  %xmm2,%xmm0
movapd %xmm0,%xmm3
unpckhpd %xmm0,%xmm3
movsd  0xd58(%rip),%xmm2
ucomisd %xmm3,%xmm2
movsd  0xd44(%rip),%xmm1
movapd %xmm3,0x20(%rsp)
jae    12fb <hinit+0x18b>
ucomisd %xmm0,%xmm2
jae    12fb <hinit+0x18b>
divsd  %xmm3,%xmm0
xorpd  %xmm1,%xmm1
mov    %edi,%r12d
ucomisd %xmm1,%xmm0
jb     1300 <hinit+0x190>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
jmp    1310 <hinit+0x1a0>
mov    %edi,%r12d
movsd  0xd11(%rip),%xmm1
xorpd  %xmm0,%xmm0
movapd %xmm0,0x20(%rsp)
jmp    1318 <hinit+0x1a8>
mov    %edi,%r12d
jmp    1318 <hinit+0x1a8>
call   1090 <sqrt@plt>
movsd  0x38(%rsp),%xmm8
movapd %xmm0,%xmm1
mulsd  0xcf8(%rip),%xmm1
movapd %xmm1,%xmm0
movapd %xmm8,%xmm1
mov    $0x2,%al
call   1030 <min_d@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x40(%rsp),%xmm1
mov    $0x2,%al
call   1050 <sign@plt>
xorps  %xmm5,%xmm5
cvtsi2sd %eax,%xmm5
mov    %r12d,%edi
test   %r12d,%r12d
movapd %xmm5,(%rsp)
je     13da <hinit+0x26a>
cmp    $0x4,%edi
movsd  0x30(%rsp),%xmm0
jae    13ff <hinit+0x28f>
xor    %eax,%eax
mov    0x18(%rsp),%rdx
mov    0xa8(%rsp),%r12d
mov    %rax,%rcx
not    %rcx
test   $0x1,%r13b
je     1397 <hinit+0x227>
movsd  (%r15,%rax,8),%xmm1
mulsd  %xmm5,%xmm1
addsd  (%rbx,%rax,8),%xmm1
movsd  %xmm1,0x0(%rbp,%rax,8)
or     $0x1,%rax
add    %r13,%rcx
je     1537 <hinit+0x3c7>
movsd  (%r15,%rax,8),%xmm1
mulsd  %xmm5,%xmm1
addsd  (%rbx,%rax,8),%xmm1
movsd  %xmm1,0x0(%rbp,%rax,8)
movsd  0x8(%r15,%rax,8),%xmm1
mulsd  %xmm5,%xmm1
addsd  0x8(%rbx,%rax,8),%xmm1
movsd  %xmm1,0x8(%rbp,%rax,8)
add    $0x2,%rax
cmp    %rax,%r13
jne    13a0 <hinit+0x230>
jmp    1537 <hinit+0x3c7>
movsd  0x30(%rsp),%xmm0
addsd  %xmm5,%xmm0
mov    %rbp,%rsi
mov    0x18(%rsp),%rdx
mov    %r14,%rcx
mov    $0x1,%al
call   1040 <safeOdeFunc@plt>
xorpd  %xmm0,%xmm0
jmp    1602 <hinit+0x492>
lea    0x0(,%r13,8),%rax
add    %rbp,%rax
lea    (%rbx,%r13,8),%rcx
cmp    %rbp,%rcx
seta   %r8b
lea    (%r15,%r13,8),%rcx
cmp    %rbx,%rax
seta   %dl
cmp    %rbp,%rcx
seta   %cl
cmp    %r15,%rax
seta   %sil
xor    %eax,%eax
test   %dl,%r8b
jne    1365 <hinit+0x1f5>
and    %sil,%cl
mov    0x18(%rsp),%rdx
mov    0xa8(%rsp),%r12d
jne    1372 <hinit+0x202>
mov    %r13d,%eax
and    $0xfffffffc,%eax
movapd %xmm5,%xmm6
unpcklpd %xmm5,%xmm6
lea    -0x4(%rax),%rcx
mov    %rcx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rcx,%rcx
je     173a <hinit+0x5ca>
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ecx,%ecx
nopl   0x0(%rax)
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movupd (%r15,%rcx,8),%xmm3
movupd 0x10(%r15,%rcx,8),%xmm4
mulpd  %xmm6,%xmm3
addpd  %xmm1,%xmm3
mulpd  %xmm6,%xmm4
addpd  %xmm2,%xmm4
movupd %xmm3,0x0(%rbp,%rcx,8)
movupd %xmm4,0x10(%rbp,%rcx,8)
movupd 0x20(%rbx,%rcx,8),%xmm1
movupd 0x30(%rbx,%rcx,8),%xmm2
movupd 0x20(%r15,%rcx,8),%xmm3
movupd 0x30(%r15,%rcx,8),%xmm4
mulpd  %xmm6,%xmm3
addpd  %xmm1,%xmm3
mulpd  %xmm6,%xmm4
addpd  %xmm2,%xmm4
movupd %xmm3,0x20(%rbp,%rcx,8)
movupd %xmm4,0x30(%rbp,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1480 <hinit+0x310>
test   $0x1,%r8b
je     152e <hinit+0x3be>
movupd (%rbx,%rcx,8),%xmm1
movupd 0x10(%rbx,%rcx,8),%xmm2
movupd (%r15,%rcx,8),%xmm3
movupd 0x10(%r15,%rcx,8),%xmm4
mulpd  %xmm6,%xmm3
addpd  %xmm1,%xmm3
mulpd  %xmm6,%xmm4
addpd  %xmm2,%xmm4
movupd %xmm3,0x0(%rbp,%rcx,8)
movupd %xmm4,0x10(%rbp,%rcx,8)
cmp    %r13,%rax
jne    1372 <hinit+0x202>
addsd  %xmm5,%xmm0
mov    %rbp,%rsi
mov    %r14,%rcx
mov    $0x1,%al
mov    %edi,%ebp
call   1040 <safeOdeFunc@plt>
xorpd  %xmm0,%xmm0
test   %r12d,%r12d
je     15b6 <hinit+0x446>
test   %ebp,%ebp
je     1602 <hinit+0x492>
xor    %eax,%eax
movapd 0xacb(%rip),%xmm1
mov    0xa0(%rsp),%rcx
mov    0x98(%rsp),%rdx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rbx,%rax,8),%xmm2
andpd  %xmm1,%xmm2
mulsd  (%rcx,%rax,8),%xmm2
addsd  (%rdx,%rax,8),%xmm2
movsd  (%r14,%rax,8),%xmm3
subsd  (%r15,%rax,8),%xmm3
divsd  %xmm2,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm3,%xmm0
add    $0x1,%rax
cmp    %rax,%r13
jne    1580 <hinit+0x410>
jmp    1602 <hinit+0x492>
test   %ebp,%ebp
je     1602 <hinit+0x492>
xor    %eax,%eax
movapd 0xa6c(%rip),%xmm1
movsd  0x50(%rsp),%xmm4
movsd  0x48(%rsp),%xmm5
movsd  (%rbx,%rax,8),%xmm2
andpd  %xmm1,%xmm2
mulsd  %xmm5,%xmm2
addsd  %xmm4,%xmm2
movsd  (%r14,%rax,8),%xmm3
subsd  (%r15,%rax,8),%xmm3
divsd  %xmm2,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm3,%xmm0
add    $0x1,%rax
cmp    %rax,%r13
jne    15d0 <hinit+0x460>
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
jb     1615 <hinit+0x4a5>
xorps  %xmm2,%xmm2
sqrtsd %xmm0,%xmm2
jmp    1622 <hinit+0x4b2>
call   1090 <sqrt@plt>
xorpd  %xmm3,%xmm3
movapd %xmm0,%xmm2
movapd (%rsp),%xmm1
movapd 0x20(%rsp),%xmm0
divsd  %xmm1,%xmm2
andpd  0x9f7(%rip),%xmm2
ucomisd %xmm3,%xmm0
jb     1648 <hinit+0x4d8>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
jmp    165d <hinit+0x4ed>
movapd %xmm2,0x20(%rsp)
call   1090 <sqrt@plt>
movapd 0x20(%rsp),%xmm2
movapd %xmm0,%xmm1
movapd %xmm2,%xmm0
mov    $0x2,%al
call   1070 <max_d@plt>
test   %eax,%eax
jle    16a3 <hinit+0x533>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x995(%rip),%xmm0
divsd  %xmm1,%xmm0
xorps  %xmm2,%xmm2
cvtsi2sdl 0x90(%rsp),%xmm2
movsd  0x985(%rip),%xmm1
divsd  %xmm2,%xmm1
call   1080 <pow@plt>
movapd (%rsp),%xmm1
jmp    16d3 <hinit+0x563>
movapd 0x985(%rip),%xmm1
andpd  (%rsp),%xmm1
mulsd  0x968(%rip),%xmm1
movsd  0x940(%rip),%xmm0
mov    $0x2,%al
call   1070 <max_d@plt>
movapd (%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0x94d(%rip),%xmm1
movapd %xmm1,(%rsp)
movsd  0x38(%rsp),%xmm1
mov    $0x2,%al
call   1030 <min_d@plt>
movaps (%rsp),%xmm0
mov    %eax,%edi
mov    $0x1,%al
call   1030 <min_d@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rsp)
mov    $0x1,%al
call   1060 <isnan@plt>
test   %eax,%eax
je     1712 <hinit+0x5a2>
int3
movsd  (%rsp),%xmm0
movsd  0x40(%rsp),%xmm1
mov    $0x2,%al
call   1050 <sign@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ecx,%ecx
test   $0x1,%r8b
jne    14fa <hinit+0x38a>
jmp    152e <hinit+0x3be>

LLVM IR (-O0):
define double @hinit(i32 %0, double %1, double* %2, double* %3, double %4, double* %5, double* %6, double* %7, i32 %8, double %9, double* %10, double* %11, i32 %12) {
%14 = alloca i32
%15 = alloca double
%16 = alloca double*
%17 = alloca double*
%18 = alloca double
%19 = alloca double*
%20 = alloca double*
%21 = alloca double*
%22 = alloca i32
%23 = alloca double
%24 = alloca double*
%25 = alloca double*
%26 = alloca i32
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca i32
store i32 %0, i32* %14
store double %1, double* %15
store double* %2, double** %16
store double* %3, double** %17
store double %4, double* %18
store double* %5, double** %19
store double* %6, double** %20
store double* %7, double** %21
store i32 %8, i32* %22
store double %9, double* %23
store double* %10, double** %24
store double* %11, double** %25
store i32 %12, i32* %26
store double 0.000000e+00, double* %27
store double 0.000000e+00, double* %28
%38 = load double*, double** %24
%39 = getelementptr inbounds double, double* %38, i64 0
%40 = load double, double* %39
store double %40, double* %29
%41 = load double*, double** %25
%42 = getelementptr inbounds double, double* %41, i64 0
%43 = load double, double* %42
store double %43, double* %30
%44 = load i32, i32* %26
%45 = icmp ne i32 %44, 0
br i1 %45, label %87, label %46
46: ; preds = %13
store i32 0, i32* %37
br label %47
47: ; preds = %83, %46
%48 = load i32, i32* %37
%49 = load i32, i32* %14
%50 = icmp ult i32 %48, %49
br i1 %50, label %51, label %86
51: ; preds = %47
%52 = load double, double* %29
%53 = load double, double* %30
%54 = load double*, double** %16
%55 = load i32, i32* %37
%56 = zext i32 %55 to i64
%57 = getelementptr inbounds double, double* %54, i64 %56
%58 = load double, double* %57
%59 = call double @llvm.fabs.f64(double %58)
%60 = call double @llvm.fmuladd.f64(double %53, double %59, double %52)
store double %60, double* %31
%61 = load double*, double** %19
%62 = load i32, i32* %37
%63 = zext i32 %62 to i64
%64 = getelementptr inbounds double, double* %61, i64 %63
%65 = load double, double* %64
%66 = load double, double* %31
%67 = fdiv double %65, %66
store double %67, double* %36
%68 = load double, double* %36
%69 = load double, double* %36
%70 = load double, double* %27
%71 = call double @llvm.fmuladd.f64(double %68, double %69, double %70)
store double %71, double* %27
%72 = load double*, double** %16
%73 = load i32, i32* %37
%74 = zext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
%76 = load double, double* %75
%77 = load double, double* %31
%78 = fdiv double %76, %77
store double %78, double* %36
%79 = load double, double* %36
%80 = load double, double* %36
%81 = load double, double* %28
%82 = call double @llvm.fmuladd.f64(double %79, double %80, double %81)
store double %82, double* %28
br label %83
83: ; preds = %51
%84 = load i32, i32* %37
%85 = add i32 %84, 1
store i32 %85, i32* %37
br label %47
86: ; preds = %47
br label %136
87: ; preds = %13
store i32 0, i32* %37
br label %88
88: ; preds = %132, %87
%89 = load i32, i32* %37
%90 = load i32, i32* %14
%91 = icmp ult i32 %89, %90
br i1 %91, label %92, label %135
92: ; preds = %88
%93 = load double*, double** %24
%94 = load i32, i32* %37
%95 = zext i32 %94 to i64
%96 = getelementptr inbounds double, double* %93, i64 %95
%97 = load double, double* %96
%98 = load double*, double** %25
%99 = load i32, i32* %37
%100 = zext i32 %99 to i64
%101 = getelementptr inbounds double, double* %98, i64 %100
%102 = load double, double* %101
%103 = load double*, double** %16
%104 = load i32, i32* %37
%105 = zext i32 %104 to i64
%106 = getelementptr inbounds double, double* %103, i64 %105
%107 = load double, double* %106
%108 = call double @llvm.fabs.f64(double %107)
%109 = call double @llvm.fmuladd.f64(double %102, double %108, double %97)
store double %109, double* %31
%110 = load double*, double** %19
%111 = load i32, i32* %37
%112 = zext i32 %111 to i64
%113 = getelementptr inbounds double, double* %110, i64 %112
%114 = load double, double* %113
%115 = load double, double* %31
%116 = fdiv double %114, %115
store double %116, double* %36
%117 = load double, double* %36
%118 = load double, double* %36
%119 = load double, double* %27
%120 = call double @llvm.fmuladd.f64(double %117, double %118, double %119)
store double %120, double* %27
%121 = load double*, double** %16
%122 = load i32, i32* %37
%123 = zext i32 %122 to i64
%124 = getelementptr inbounds double, double* %121, i64 %123
%125 = load double, double* %124
%126 = load double, double* %31
%127 = fdiv double %125, %126
store double %127, double* %36
%128 = load double, double* %36
%129 = load double, double* %36
%130 = load double, double* %28
%131 = call double @llvm.fmuladd.f64(double %128, double %129, double %130)
store double %131, double* %28
br label %132
132: ; preds = %92
%133 = load i32, i32* %37
%134 = add i32 %133, 1
store i32 %134, i32* %37
br label %88
135: ; preds = %88
br label %136
136: ; preds = %135, %86
%137 = load double, double* %27
%138 = fcmp ole double %137, 1.000000e-10
br i1 %138, label %142, label %139
139: ; preds = %136
%140 = load double, double* %28
%141 = fcmp ole double %140, 1.000000e-10
br i1 %141, label %142, label %143
142: ; preds = %139, %136
store double 0x3EB0C6F7A0B5ED8D, double* %32
br label %149
143: ; preds = %139
%144 = load double, double* %28
%145 = load double, double* %27
%146 = fdiv double %144, %145
%147 = call double @sqrt(double %146)
%148 = fmul double %147, 1.000000e-02
store double %148, double* %32
br label %149
149: ; preds = %143, %142
%150 = load double, double* %32
%151 = load double, double* %23
%152 = call i32 (double, double, ...) bitcast (i32 (...)* @min_d to i32 (double, double, ...)*)(double %150, double %151)
%153 = sitofp i32 %152 to double
store double %153, double* %32
%154 = load double, double* %32
%155 = load double, double* %18
%156 = call i32 (double, double, ...) bitcast (i32 (...)* @sign to i32 (double, double, ...)*)(double %154, double %155)
%157 = sitofp i32 %156 to double
store double %157, double* %32
store i32 0, i32* %37
br label %158
158: ; preds = %179, %149
%159 = load i32, i32* %37
%160 = load i32, i32* %14
%161 = icmp ult i32 %159, %160
br i1 %161, label %162, label %182
162: ; preds = %158
%163 = load double*, double** %16
%164 = load i32, i32* %37
%165 = zext i32 %164 to i64
%166 = getelementptr inbounds double, double* %163, i64 %165
%167 = load double, double* %166
%168 = load double, double* %32
%169 = load double*, double** %19
%170 = load i32, i32* %37
%171 = zext i32 %170 to i64
%172 = getelementptr inbounds double, double* %169, i64 %171
%173 = load double, double* %172
%174 = call double @llvm.fmuladd.f64(double %168, double %173, double %167)
%175 = load double*, double** %21
%176 = load i32, i32* %37
%177 = zext i32 %176 to i64
%178 = getelementptr inbounds double, double* %175, i64 %177
store double %174, double* %178
br label %179
179: ; preds = %162
%180 = load i32, i32* %37
%181 = add i32 %180, 1
store i32 %181, i32* %37
br label %158
182: ; preds = %158
%183 = load i32, i32* %14
%184 = load double, double* %15
%185 = load double, double* %32
%186 = fadd double %184, %185
%187 = load double*, double** %21
%188 = load double*, double** %17
%189 = load double*, double** %20
%190 = call i32 (i32, double, double*, double*, double*, ...) bitcast (i32 (...)* @safeOdeFunc to i32 (i32, double, double*, double*, double*, ...)*)(i32 %183, double %186, double* %187, double* %188, double* %189)
store double 0.000000e+00, double* %34
%191 = load i32, i32* %26
%192 = icmp ne i32 %191, 0
br i1 %192, label %229, label %193
193: ; preds = %182
store i32 0, i32* %37
br label %194
194: ; preds = %225, %193
%195 = load i32, i32* %37
%196 = load i32, i32* %14
%197 = icmp ult i32 %195, %196
br i1 %197, label %198, label %228
198: ; preds = %194
%199 = load double, double* %29
%200 = load double, double* %30
%201 = load double*, double** %16
%202 = load i32, i32* %37
%203 = zext i32 %202 to i64
%204 = getelementptr inbounds double, double* %201, i64 %203
%205 = load double, double* %204
%206 = call double @llvm.fabs.f64(double %205)
%207 = call double @llvm.fmuladd.f64(double %200, double %206, double %199)
store double %207, double* %31
%208 = load double*, double** %20
%209 = load i32, i32* %37
%210 = zext i32 %209 to i64
%211 = getelementptr inbounds double, double* %208, i64 %210
%212 = load double, double* %211
%213 = load double*, double** %19
%214 = load i32, i32* %37
%215 = zext i32 %214 to i64
%216 = getelementptr inbounds double, double* %213, i64 %215
%217 = load double, double* %216
%218 = fsub double %212, %217
%219 = load double, double* %31
%220 = fdiv double %218, %219
store double %220, double* %36
%221 = load double, double* %36
%222 = load double, double* %36
%223 = load double, double* %34
%224 = call double @llvm.fmuladd.f64(double %221, double %222, double %223)
store double %224, double* %34
br label %225
225: ; preds = %198
%226 = load i32, i32* %37
%227 = add i32 %226, 1
store i32 %227, i32* %37
br label %194
228: ; preds = %194
br label %273
229: ; preds = %182
store i32 0, i32* %37
br label %230
230: ; preds = %269, %229
%231 = load i32, i32* %37
%232 = load i32, i32* %14
%233 = icmp ult i32 %231, %232
br i1 %233, label %234, label %272
234: ; preds = %230
%235 = load double*, double** %24
%236 = load i32, i32* %37
%237 = zext i32 %236 to i64
%238 = getelementptr inbounds double, double* %235, i64 %237
%239 = load double, double* %238
%240 = load double*, double** %25
%241 = load i32, i32* %37
%242 = zext i32 %241 to i64
%243 = getelementptr inbounds double, double* %240, i64 %242
%244 = load double, double* %243
%245 = load double*, double** %16
%246 = load i32, i32* %37
%247 = zext i32 %246 to i64
%248 = getelementptr inbounds double, double* %245, i64 %247
%249 = load double, double* %248
%250 = call double @llvm.fabs.f64(double %249)
%251 = call double @llvm.fmuladd.f64(double %244, double %250, double %239)
store double %251, double* %31
%252 = load double*, double** %20
%253 = load i32, i32* %37
%254 = zext i32 %253 to i64
%255 = getelementptr inbounds double, double* %252, i64 %254
%256 = load double, double* %255
%257 = load double*, double** %19
%258 = load i32, i32* %37
%259 = zext i32 %258 to i64
%260 = getelementptr inbounds double, double* %257, i64 %259
%261 = load double, double* %260
%262 = fsub double %256, %261
%263 = load double, double* %31
%264 = fdiv double %262, %263
store double %264, double* %36
%265 = load double, double* %36
%266 = load double, double* %36
%267 = load double, double* %34
%268 = call double @llvm.fmuladd.f64(double %265, double %266, double %267)
store double %268, double* %34
br label %269
269: ; preds = %234
%270 = load i32, i32* %37
%271 = add i32 %270, 1
store i32 %271, i32* %37
br label %230
272: ; preds = %230
br label %273
273: ; preds = %272, %228
%274 = load double, double* %34
%275 = call double @sqrt(double %274)
%276 = load double, double* %32
%277 = fdiv double %275, %276
store double %277, double* %34
%278 = load double, double* %34
%279 = call double @llvm.fabs.f64(double %278)
%280 = load double, double* %27
%281 = call double @sqrt(double %280)
%282 = call i32 (double, double, ...) bitcast (i32 (...)* @max_d to i32 (double, double, ...)*)(double %279, double %281)
%283 = sitofp i32 %282 to double
store double %283, double* %35
%284 = load double, double* %35
%285 = fcmp ole double %284, 1.000000e-15
br i1 %285, label %286, label %292
286: ; preds = %273
%287 = load double, double* %32
%288 = call double @llvm.fabs.f64(double %287)
%289 = fmul double %288, 1.000000e-03
%290 = call i32 (double, double, ...) bitcast (i32 (...)* @max_d to i32 (double, double, ...)*)(double 0x3EB0C6F7A0B5ED8D, double %289)
%291 = sitofp i32 %290 to double
store double %291, double* %33
br label %299
292: ; preds = %273
%293 = load double, double* %35
%294 = fdiv double 1.000000e-02, %293
%295 = load i32, i32* %22
%296 = sitofp i32 %295 to double
%297 = fdiv double 1.000000e+00, %296
%298 = call double @pow(double %294, double %297)
store double %298, double* %33
br label %299
299: ; preds = %292, %286
%300 = load double, double* %32
%301 = fmul double 1.000000e+02, %300
%302 = load double, double* %33
%303 = load double, double* %23
%304 = call i32 (double, double, ...) bitcast (i32 (...)* @min_d to i32 (double, double, ...)*)(double %302, double %303)
%305 = call i32 (double, i32, ...) bitcast (i32 (...)* @min_d to i32 (double, i32, ...)*)(double %301, i32 %304)
%306 = sitofp i32 %305 to double
store double %306, double* %32
%307 = load double, double* %32
%308 = call i32 (double, ...) bitcast (i32 (...)* @isnan to i32 (double, ...)*)(double %307)
%309 = icmp ne i32 %308, 0
br i1 %309, label %310, label %311
310: ; preds = %299
br label %311
311: ; preds = %310, %299
%312 = load double, double* %32
%313 = load double, double* %18
%314 = call i32 (double, double, ...) bitcast (i32 (...)* @sign to i32 (double, double, ...)*)(double %312, double %313)
%315 = sitofp i32 %314 to double
ret double %315
}

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