void qreig(n, d, e2) int n;
double *d;
double *e2;
{
  int i, k, m;
  double b, b2, f, g, h, p2, r2, s2;
  f = 0.0;
  b2 = 0.0;
  b = 0.0;
  e2[n] = 0.0;
  for (k = 1; k <= n; k++) {
    h = 1.0e-6 * 1.0e-6 * (d[k] * d[k] + e2[k]);
    if (b2 < h) {
      b = sqrt(h);
      b2 = h;
    };
    for (m = k; m <= n; m++) {
      if (e2[m] <= b2)
        goto cont1;
    };
  cont1:
    if (m == k)
      goto root;
  nextit:
    g = d[k];
    p2 = sqrt(e2[k]);
    h = (d[k + 1] - g) / (2.0 * p2);
    r2 = sqrt(h * h + 1.0);
    if (h < 0.0) {
      h = p2 / (h - r2);
    } else {
      h = p2 / (h + r2);
    };
    d[k] = h;
    h = g - h;
    f = f + h;
    for (i = k + 1; i <= n; i++) {
      d[i] = d[i] - h;
    };
    g = d[m];
    if (g == 0.0)
      g = b;
    h = g;
    s2 = 0.0;
    for (i = m - 1; i >= k; i = i - 1) {
      p2 = g * h;
      r2 = p2 + e2[i];
      e2[i + 1] = s2 * r2;
      s2 = e2[i] / r2;
      d[i + 1] = h + s2 * (h + d[i]);
      g = d[i] - e2[i] / g;
      if (g == 0.0)
        g = b;
      h = g * p2 / r2;
    };
    e2[k] = s2 * g * h;
    d[k] = h;
    if (e2[k] > b2)
      goto nextit;
  root:
    h = d[k] + f;
    for (i = k; i >= 2; i = i - 1) {
      if (h < d[i - 1]) {
        d[i] = d[i - 1];
      } else {
        goto cont2;
      };
    };
    i = 1;
  cont2:
    d[i] = h;
  };
}


x86 assembly (clang, -O0):
<qreig>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x4(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     15ac <qreig+0x49c>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe6b(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11cb <qreig+0xbb>
movsd  -0x50(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x30(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1211 <qreig+0x101>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     11fe <qreig+0xee>
jmp    1216 <qreig+0x106>
jmp    1203 <qreig+0xf3>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11d1 <qreig+0xc1>
jmp    1216 <qreig+0x106>
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jne    1227 <qreig+0x117>
jmp    1501 <qreig+0x3f1>
jmp    122c <qreig+0x11c>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  0xd9c(%rip),%xmm1
mulsd  -0x58(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd70(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x50(%rbp),%xmm0
jbe    12d1 <qreig+0x1c1>
movsd  -0x58(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
subsd  -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    12e9 <qreig+0x1d9>
movsd  -0x58(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
addsd  -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     135b <qreig+0x24b>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1322 <qreig+0x212>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    138f <qreig+0x27f>
jp     138f <qreig+0x27f>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jl     14b1 <qreig+0x3a1>
movsd  -0x48(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x68(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
divsd  -0x48(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    148f <qreig+0x37f>
jp     148f <qreig+0x37f>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    13aa <qreig+0x29a>
movsd  -0x68(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    14fc <qreig+0x3ec>
jmp    122c <qreig+0x11c>
jmp    1501 <qreig+0x3f1>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x2,-0x1c(%rbp)
jl     1585 <qreig+0x475>
movsd  -0x50(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    156d <qreig+0x45d>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1572 <qreig+0x462>
jmp    158c <qreig+0x47c>
jmp    1577 <qreig+0x467>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    151e <qreig+0x40e>
movl   $0x1,-0x1c(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1152 <qreig+0x42>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qreig>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xc8,%rsp
movslq %edi,%rcx
movq   $0x0,(%rdx,%rcx,8)
mov    %rcx,0xa0(%rsp)
test   %ecx,%ecx
jle    1688 <qreig+0x578>
mov    %rdx,%r14
mov    %rsi,%rbx
mov    %rdi,%rax
lea    0x1(%rdi),%ebp
lea    0x2(%rdi),%ecx
mov    %rcx,0x60(%rsp)
add    $0xfffffffe,%eax
lea    0x40(%rsi),%r13
xorpd  %xmm9,%xmm9
mov    $0x1,%r15d
mov    $0x2,%ecx
mov    %rcx,0x8(%rsp)
xor    %r12d,%r12d
mov    %rax,0x68(%rsp)
xorpd  %xmm8,%xmm8
xorpd  %xmm10,%xmm10
xorpd  %xmm11,%xmm11
jmp    11ca <qreig+0xba>
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x30(%rsp),%r12
mov    %eax,%eax
movsd  %xmm1,(%rbx,%rax,8)
add    $0x1,%r15
mov    0x8(%rsp),%rax
add    $0x1,%rax
add    $0x1,%r12d
add    $0x8,%r13
mov    0x2c(%rsp),%ecx
add    $0xffffffff,%ecx
mov    %rax,0x8(%rsp)
cmp    0x60(%rsp),%rax
mov    %ecx,%eax
je     1688 <qreig+0x578>
mov    %eax,0x2c(%rsp)
mov    %eax,%ecx
add    $0x1,%rcx
and    $0xfffffffffffffffc,%rcx
add    $0xfffffffffffffffc,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    0x68(%rsp),%rax
sub    %r12d,%eax
mov    %rax,0x88(%rsp)
add    $0x1,%rax
mov    %rax,0x80(%rsp)
mov    %rax,%rdx
movabs $0x1fffffffc,%rax
and    %rax,%rdx
mov    %rdx,%rax
mov    %rdx,0x38(%rsp)
lea    -0x4(%rdx),%rax
mov    %rax,0x78(%rsp)
movsd  (%rbx,%r15,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  (%r14,%r15,8),%xmm0
shr    $0x2,%rax
mulsd  0xdc3(%rip),%xmm0
ucomisd %xmm10,%xmm0
jbe    129e <qreig+0x18e>
ucomisd %xmm9,%xmm0
jb     1260 <qreig+0x150>
xorps  %xmm11,%xmm11
sqrtsd %xmm0,%xmm11
movapd %xmm0,%xmm10
jmp    129e <qreig+0x18e>
nopl   0x0(%rax,%rax,1)
movsd  %xmm0,0x50(%rsp)
movsd  %xmm8,0x18(%rsp)
mov    %rcx,0x10(%rsp)
mov    %rax,0x20(%rsp)
call   1030 <sqrt@plt>
mov    0x20(%rsp),%rax
mov    0x10(%rsp),%rcx
movsd  0x18(%rsp),%xmm8
xorpd  %xmm9,%xmm9
movapd %xmm0,%xmm11
movsd  0x50(%rsp),%xmm10
and    $0xfffffffffffffffe,%rcx
add    $0x1,%rax
mov    %rax,0x20(%rsp)
mov    %r15,%rax
xchg   %ax,%ax
ucomisd (%r14,%rax,8),%xmm10
jae    12c3 <qreig+0x1b3>
add    $0x1,%rax
cmp    %rax,%rbp
jne    12b0 <qreig+0x1a0>
mov    %ebp,%eax
mov    %eax,%esi
cmp    %rsi,%r15
mov    %r12,0x30(%rsp)
jne    1330 <qreig+0x220>
movsd  (%rbx,%r15,8),%xmm1
addsd  %xmm8,%xmm1
mov    $0x1,%eax
cmp    $0x2,%r15
jb     1190 <qreig+0x80>
mov    %r15,%rcx
mov    0x30(%rsp),%r12
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    -0x1(%rcx),%rdx
mov    %edx,%esi
movsd  (%rbx,%rsi,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1680 <qreig+0x570>
movsd  %xmm0,(%rbx,%rcx,8)
lea    0x1(%rdx),%rsi
mov    %rdx,%rcx
cmp    $0x2,%rsi
jg     1300 <qreig+0x1f0>
jmp    1195 <qreig+0x85>
nopl   0x0(%rax)
lea    0x1(%r15),%eax
mov    %rax,0xa8(%rsp)
movslq %esi,%r12
movsd  (%rbx,%r15,8),%xmm7
movsd  (%r14,%r15,8),%xmm0
mov    0x8(%rsp),%rax
mov    0x38(%rsp),%rdx
add    %rdx,%rax
mov    %rax,0x70(%rsp)
movapd %xmm11,0xb0(%rsp)
movsd  %xmm10,0x48(%rsp)
mov    %rcx,0x10(%rsp)
mov    %rsi,0x40(%rsp)
jmp    13b2 <qreig+0x2a2>
nopw   0x0(%rax,%rax,1)
movapd %xmm1,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm7,%xmm8
mulsd  %xmm1,%xmm2
movsd  %xmm2,(%r14,%r15,8)
movsd  %xmm1,(%rbx,%r15,8)
movsd  (%r14,%r15,8),%xmm0
ucomisd %xmm10,%xmm0
movapd %xmm1,%xmm7
jbe    12d5 <qreig+0x1c5>
ucomisd %xmm9,%xmm0
movsd  %xmm8,0x18(%rsp)
movapd %xmm7,0x50(%rsp)
jb     13d0 <qreig+0x2c0>
xorps  %xmm3,%xmm3
sqrtsd %xmm0,%xmm3
jmp    13e4 <qreig+0x2d4>
nop
call   1030 <sqrt@plt>
movapd 0x50(%rsp),%xmm7
xorpd  %xmm9,%xmm9
movapd %xmm0,%xmm3
mov    0xa8(%rsp),%rax
movsd  (%rbx,%rax,8),%xmm4
subsd  %xmm7,%xmm4
movapd %xmm3,%xmm0
addsd  %xmm3,%xmm0
divsd  %xmm0,%xmm4
movapd %xmm4,%xmm0
mulsd  %xmm4,%xmm0
addsd  0xbf7(%rip),%xmm0
ucomisd %xmm9,%xmm0
jb     1420 <qreig+0x310>
sqrtsd %xmm0,%xmm0
jmp    1454 <qreig+0x344>
xchg   %ax,%ax
movsd  %xmm3,0x98(%rsp)
movsd  %xmm4,0x90(%rsp)
call   1030 <sqrt@plt>
movsd  0x90(%rsp),%xmm4
movsd  0x98(%rsp),%xmm3
movapd 0x50(%rsp),%xmm7
xorpd  %xmm9,%xmm9
movsd  0x18(%rsp),%xmm8
movapd 0xb0(%rsp),%xmm11
movsd  0x48(%rsp),%xmm10
mov    0x10(%rsp),%rcx
mov    0x40(%rsp),%rdx
movapd %xmm0,%xmm1
xorpd  0xb8e(%rip),%xmm1
movapd %xmm4,%xmm2
cmpltsd %xmm9,%xmm2
andpd  %xmm2,%xmm1
andnpd %xmm0,%xmm2
orpd   %xmm1,%xmm2
addsd  %xmm4,%xmm2
divsd  %xmm2,%xmm3
movsd  %xmm3,(%rbx,%r15,8)
subsd  %xmm3,%xmm7
cmp    0xa0(%rsp),%r15
jge    1596 <qreig+0x486>
mov    0x8(%rsp),%rax
cmpl   $0x3,0x88(%rsp)
03
jb     1580 <qreig+0x470>
movapd %xmm7,%xmm0
unpcklpd %xmm7,%xmm0
cmpq   $0x0,0x78(%rsp)
je     165f <qreig+0x54f>
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movupd -0x30(%r13,%rax,8),%xmm1
movupd -0x20(%r13,%rax,8),%xmm2
movupd -0x10(%r13,%rax,8),%xmm3
movupd 0x0(%r13,%rax,8),%xmm4
subpd  %xmm0,%xmm1
subpd  %xmm0,%xmm2
movupd %xmm1,-0x30(%r13,%rax,8)
movupd %xmm2,-0x20(%r13,%rax,8)
subpd  %xmm0,%xmm3
subpd  %xmm0,%xmm4
movupd %xmm3,-0x10(%r13,%rax,8)
movupd %xmm4,0x0(%r13,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rcx
jne    14f0 <qreig+0x3e0>
testb  $0x1,0x20(%rsp)
je     156c <qreig+0x45c>
add    0x8(%rsp),%rax
movupd (%rbx,%rax,8),%xmm1
movupd 0x10(%rbx,%rax,8),%xmm2
subpd  %xmm0,%xmm1
subpd  %xmm0,%xmm2
movupd %xmm1,(%rbx,%rax,8)
movupd %xmm2,0x10(%rbx,%rax,8)
mov    0x70(%rsp),%rax
mov    0x38(%rsp),%rcx
cmp    %rcx,0x80(%rsp)
je     1596 <qreig+0x486>
movsd  (%rbx,%rax,8),%xmm0
subsd  %xmm7,%xmm0
movsd  %xmm0,(%rbx,%rax,8)
add    $0x1,%rax
cmp    %eax,%ebp
jne    1580 <qreig+0x470>
movsd  (%rbx,%rdx,8),%xmm2
xorpd  %xmm0,%xmm0
movapd %xmm2,%xmm1
cmpeqsd %xmm0,%xmm1
movapd %xmm1,%xmm3
andnpd %xmm2,%xmm3
andpd  %xmm11,%xmm1
orpd   %xmm3,%xmm1
cmp    %r12,%r15
jge    1380 <qreig+0x270>
xorpd  %xmm5,%xmm5
mov    %r12,%rax
movapd %xmm1,%xmm2
nopl   (%rax)
movapd %xmm2,%xmm3
mulsd  %xmm1,%xmm3
movsd  -0x8(%r14,%rax,8),%xmm0
movapd %xmm3,%xmm4
addsd  %xmm0,%xmm4
mulsd  %xmm4,%xmm5
movsd  %xmm5,(%r14,%rax,8)
divsd  %xmm4,%xmm0
movsd  -0x8(%rbx,%rax,8),%xmm5
movapd %xmm1,%xmm6
addsd  %xmm5,%xmm6
mulsd  %xmm0,%xmm6
addsd  %xmm1,%xmm6
movsd  %xmm6,(%rbx,%rax,8)
movsd  -0x8(%r14,%rax,8),%xmm1
lea    -0x1(%rax),%rcx
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm5
movapd %xmm5,%xmm2
cmpeqsd %xmm9,%xmm2
movapd %xmm2,%xmm1
andnpd %xmm5,%xmm1
andpd  %xmm11,%xmm2
orpd   %xmm1,%xmm2
movapd %xmm2,%xmm1
mulsd  %xmm3,%xmm1
divsd  %xmm4,%xmm1
mov    %rcx,%rax
movapd %xmm0,%xmm5
cmp    %r15,%rcx
jg     15d0 <qreig+0x4c0>
jmp    1384 <qreig+0x274>
xor    %eax,%eax
testb  $0x1,0x20(%rsp)
jne    1549 <qreig+0x439>
jmp    156c <qreig+0x45c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rcx,%rax
jmp    1195 <qreig+0x85>
add    $0xc8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @qreig(i32 %0, double* nocapture %1, double* nocapture %2) {
%4 = sext i32 %0 to i64
%5 = getelementptr inbounds double, double* %2, i64 %4
store double 0.000000e+00, double* %5
%6 = icmp slt i32 %0, 1
br i1 %6, label %215, label %7
7: ; preds = %3
%8 = add nuw i32 %0, 1
%9 = add nuw i32 %0, 2
%10 = zext i32 %9 to i64
%11 = zext i32 %8 to i64
%12 = add i32 %0, -2
br label %13
13: ; preds = %7, %207
%14 = phi i32 [ 0, %7 ], [ %214, %207 ]
%15 = phi i64 [ 2, %7 ], [ %212, %207 ]
%16 = phi i64 [ 1, %7 ], [ %211, %207 ]
%17 = phi double [ 0.000000e+00, %7 ], [ %194, %207 ]
%18 = phi double [ 0.000000e+00, %7 ], [ %41, %207 ]
%19 = phi double [ 0.000000e+00, %7 ], [ %40, %207 ]
%20 = sub i32 %12, %14
%21 = zext i32 %20 to i64
%22 = add nuw nsw i64 %21, 1
%23 = and i64 %22, 8589934588
%24 = add nsw i64 %23, -4
%25 = lshr exact i64 %24, 2
%26 = add nuw nsw i64 %25, 1
%27 = sub i32 %12, %14
%28 = zext i32 %27 to i64
%29 = add nuw nsw i64 %28, 1
%30 = getelementptr inbounds double, double* %1, i64 %16
%31 = load double, double* %30
%32 = getelementptr inbounds double, double* %2, i64 %16
%33 = load double, double* %32
%34 = tail call double @llvm.fmuladd.f64(double %31, double %31, double %33)
%35 = fmul double %34, 0x3D719799812DEA11
%36 = fcmp olt double %18, %35
br i1 %36, label %37, label %39
37: ; preds = %13
%38 = tail call double @sqrt(double %35)
br label %39
39: ; preds = %13, %37
%40 = phi double [ %38, %37 ], [ %19, %13 ]
%41 = phi double [ %35, %37 ], [ %18, %13 ]
br label %42
42: ; preds = %39, %47
%43 = phi i64 [ %16, %39 ], [ %48, %47 ]
%44 = getelementptr inbounds double, double* %2, i64 %43
%45 = load double, double* %44
%46 = fcmp ugt double %45, %41
br i1 %46, label %47, label %50
47: ; preds = %42
%48 = add nuw nsw i64 %43, 1
%49 = icmp eq i64 %48, %11
br i1 %49, label %52, label %42
50: ; preds = %42
%51 = trunc i64 %43 to i32
br label %52
52: ; preds = %47, %50
%53 = phi i32 [ %51, %50 ], [ %8, %47 ]
%54 = zext i32 %53 to i64
%55 = icmp eq i64 %16, %54
br i1 %55, label %56, label %58
56: ; preds = %52
%57 = load double, double* %30
br label %192
58: ; preds = %52
%59 = add nuw i64 %16, 1
%60 = and i64 %59, 4294967295
%61 = getelementptr inbounds double, double* %1, i64 %60
%62 = icmp slt i64 %16, %4
%63 = getelementptr inbounds double, double* %1, i64 %54
%64 = sext i32 %53 to i64
%65 = icmp slt i64 %16, %64
%66 = load double, double* %30
%67 = load double, double* %32
%68 = icmp ult i32 %27, 3
%69 = and i64 %29, 8589934588
%70 = add nuw i64 %15, %69
%71 = and i64 %26, 1
%72 = icmp eq i64 %24, 0
%73 = and i64 %26, 9223372036854775806
%74 = icmp eq i64 %71, 0
%75 = icmp eq i64 %29, %69
br label %76
76: ; preds = %58, %184
%77 = phi double [ %190, %184 ], [ %67, %58 ]
%78 = phi double [ %186, %184 ], [ %66, %58 ]
%79 = phi double [ %93, %184 ], [ %17, %58 ]
%80 = tail call double @sqrt(double %77)
%81 = load double, double* %61
%82 = fsub double %81, %78
%83 = fmul double %80, 2.000000e+00
%84 = fdiv double %82, %83
%85 = tail call double @llvm.fmuladd.f64(double %84, double %84, double 1.000000e+00)
%86 = tail call double @sqrt(double %85)
%87 = fcmp olt double %84, 0.000000e+00
%88 = fneg double %86
%89 = select i1 %87, double %88, double %86
%90 = fadd double %84, %89
%91 = fdiv double %80, %90
store double %91, double* %30
%92 = fsub double %78, %91
%93 = fadd double %79, %92
br i1 %62, label %94, label %154
94: ; preds = %76
br i1 %68, label %144, label %95
95: ; preds = %94
%96 = insertelement <2 x double> poison, double %92, i64 0
%97 = shufflevector <2 x double> %96, <2 x double> poison, <2 x i32> zeroinitializer
%98 = insertelement <2 x double> poison, double %92, i64 0
%99 = shufflevector <2 x double> %98, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %72, label %129, label %100
100: ; preds = %95, %100
%101 = phi i64 [ %126, %100 ], [ 0, %95 ]
%102 = phi i64 [ %127, %100 ], [ 0, %95 ]
%103 = add i64 %15, %101
%104 = getelementptr inbounds double, double* %1, i64 %103
%105 = bitcast double* %104 to <2 x double>*
%106 = load <2 x double>, <2 x double>* %105
%107 = getelementptr inbounds double, double* %104, i64 2
%108 = bitcast double* %107 to <2 x double>*
%109 = load <2 x double>, <2 x double>* %108
%110 = fsub <2 x double> %106, %97
%111 = fsub <2 x double> %109, %99
%112 = bitcast double* %104 to <2 x double>*
store <2 x double> %110, <2 x double>* %112
%113 = bitcast double* %107 to <2 x double>*
store <2 x double> %111, <2 x double>* %113
%114 = or i64 %101, 4
%115 = add i64 %15, %114
%116 = getelementptr inbounds double, double* %1, i64 %115
%117 = bitcast double* %116 to <2 x double>*
%118 = load <2 x double>, <2 x double>* %117
%119 = getelementptr inbounds double, double* %116, i64 2
%120 = bitcast double* %119 to <2 x double>*
%121 = load <2 x double>, <2 x double>* %120
%122 = fsub <2 x double> %118, %97
%123 = fsub <2 x double> %121, %99
%124 = bitcast double* %116 to <2 x double>*
store <2 x double> %122, <2 x double>* %124
%125 = bitcast double* %119 to <2 x double>*
store <2 x double> %123, <2 x double>* %125
%126 = add nuw i64 %101, 8
%127 = add i64 %102, 2
%128 = icmp eq i64 %127, %73
br i1 %128, label %129, label %100
129: ; preds = %100, %95
%130 = phi i64 [ 0, %95 ], [ %126, %100 ]
br i1 %74, label %143, label %131
131: ; preds = %129
%132 = add i64 %15, %130
%133 = getelementptr inbounds double, double* %1, i64 %132
%134 = bitcast double* %133 to <2 x double>*
%135 = load <2 x double>, <2 x double>* %134
%136 = getelementptr inbounds double, double* %133, i64 2
%137 = bitcast double* %136 to <2 x double>*
%138 = load <2 x double>, <2 x double>* %137
%139 = fsub <2 x double> %135, %97
%140 = fsub <2 x double> %138, %99
%141 = bitcast double* %133 to <2 x double>*
store <2 x double> %139, <2 x double>* %141
%142 = bitcast double* %136 to <2 x double>*
store <2 x double> %140, <2 x double>* %142
br label %143
143: ; preds = %129, %131
br i1 %75, label %154, label %144
144: ; preds = %94, %143
%145 = phi i64 [ %15, %94 ], [ %70, %143 ]
br label %146
146: ; preds = %144, %146
%147 = phi i64 [ %151, %146 ], [ %145, %144 ]
%148 = getelementptr inbounds double, double* %1, i64 %147
%149 = load double, double* %148
%150 = fsub double %149, %92
store double %150, double* %148
%151 = add nuw nsw i64 %147, 1
%152 = trunc i64 %151 to i32
%153 = icmp eq i32 %8, %152
br i1 %153, label %154, label %146
154: ; preds = %146, %143, %76
%155 = load double, double* %63
%156 = fcmp oeq double %155, 0.000000e+00
%157 = select i1 %156, double %40, double %155
br i1 %65, label %158, label %184
158: ; preds = %154, %158
%159 = phi i64 [ %163, %158 ], [ %64, %154 ]
%160 = phi double [ %170, %158 ], [ 0.000000e+00, %154 ]
%161 = phi double [ %182, %158 ], [ %157, %154 ]
%162 = phi double [ %180, %158 ], [ %157, %154 ]
%163 = add nsw i64 %159, -1
%164 = fmul double %162, %161
%165 = getelementptr inbounds double, double* %2, i64 %163
%166 = load double, double* %165
%167 = fadd double %164, %166
%168 = fmul double %160, %167
%169 = getelementptr inbounds double, double* %2, i64 %159
store double %168, double* %169
%170 = fdiv double %166, %167
%171 = getelementptr inbounds double, double* %1, i64 %163
%172 = load double, double* %171
%173 = fadd double %161, %172
%174 = tail call double @llvm.fmuladd.f64(double %170, double %173, double %161)
%175 = getelementptr inbounds double, double* %1, i64 %159
store double %174, double* %175
%176 = load double, double* %165
%177 = fdiv double %176, %162
%178 = fsub double %172, %177
%179 = fcmp oeq double %178, 0.000000e+00
%180 = select i1 %179, double %40, double %178
%181 = fmul double %164, %180
%182 = fdiv double %181, %167
%183 = icmp sgt i64 %163, %16
br i1 %183, label %158, label %184
184: ; preds = %158, %154
%185 = phi double [ %157, %154 ], [ %180, %158 ]
%186 = phi double [ %157, %154 ], [ %182, %158 ]
%187 = phi double [ 0.000000e+00, %154 ], [ %170, %158 ]
%188 = fmul double %185, %187
%189 = fmul double %186, %188
store double %189, double* %32
store double %186, double* %30
%190 = load double, double* %32
%191 = fcmp ogt double %190, %41
br i1 %191, label %76, label %192
192: ; preds = %184, %56
%193 = phi double [ %57, %56 ], [ %186, %184 ]
%194 = phi double [ %17, %56 ], [ %93, %184 ]
%195 = fadd double %194, %193
%196 = icmp ugt i64 %16, 1
br i1 %196, label %197, label %207
197: ; preds = %192, %204
%198 = phi i64 [ %199, %204 ], [ %16, %192 ]
%199 = add nsw i64 %198, -1
%200 = and i64 %199, 4294967295
%201 = getelementptr inbounds double, double* %1, i64 %200
%202 = load double, double* %201
%203 = fcmp olt double %195, %202
br i1 %203, label %204, label %207
204: ; preds = %197
%205 = getelementptr inbounds double, double* %1, i64 %198
store double %202, double* %205
%206 = icmp sgt i64 %198, 2
br i1 %206, label %197, label %207
207: ; preds = %197, %204, %192
%208 = phi i64 [ 1, %192 ], [ 1, %204 ], [ %198, %197 ]
%209 = and i64 %208, 4294967295
%210 = getelementptr inbounds double, double* %1, i64 %209
store double %195, double* %210
%211 = add nuw nsw i64 %16, 1
%212 = add nuw nsw i64 %15, 1
%213 = icmp eq i64 %212, %10
%214 = add i32 %14, 1
br i1 %213, label %215, label %13
215: ; preds = %207, %3
ret void
}