#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ erhand(char *);
double fabs(double);
double sqrt(double);
void tqli(double *d, double *e, int n, double **z) {
  int m, l, iter, i, k;
  double s, r, p, g, f, dd, c, b;
  for (i = 1; i < n; i++)
    e[i - 1] = e[i];
  e[n - 1] = 0.0;
  for (l = 0; l < n; l++) {
    iter = 0;
    do {
      for (m = l; m < n - 1; m++) {
        dd = fabs(d[m]) + fabs(d[m + 1]);
        if (fabs(e[m]) + dd == dd)
          break;
      }
      if (m != l) {
        if (iter++ == 30)
          erhand("No convergence in TLQI.");
        g = (d[l + 1] - d[l]) / (2.0 * e[l]);
        r = sqrt((g * g) + 1.0);
        g = d[m] - d[l] + e[l] / (g + ((g) < 0 ? -fabs(r) : fabs(r)));
        s = c = 1.0;
        p = 0.0;
        for (i = m - 1; i >= l; i--) {
          f = s * e[i];
          b = c * e[i];
          if (fabs(f) >= fabs(g)) {
            c = g / f;
            r = sqrt((c * c) + 1.0);
            e[i + 1] = f * r;
            c *= (s = 1.0 / r);
          } else {
            s = f / g;
            r = sqrt((s * s) + 1.0);
            e[i + 1] = g * r;
            s *= (c = 1.0 / r);
          }
          g = d[i + 1] - p;
          r = (d[i] - g) * s + 2.0 * c * b;
          p = s * r;
          d[i + 1] = g + p;
          g = c * r - b;
          for (k = 0; k < n; k++) {
            f = z[k][i + 1];
            z[k][i + 1] = s * z[k][i] + c * f;
            z[k][i] = c * z[k][i] - s * f;
          }
        }
        d[l] = d[l] - p;
        e[l] = g;
        e[m] = 0.0;
      }
    } while (m != l);
  }
}


x86 assembly (clang, -O0):
<tqli>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    117a <tqli+0x5a>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1141 <tqli+0x21>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    170a <tqli+0x5ea>
movl   $0x0,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1244 <tqli+0x124>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe2c(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xe0f(%rip),%xmm2
pand   %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xdee(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0x68(%rbp),%xmm0
ucomisd -0x68(%rbp),%xmm0
jne    1231 <tqli+0x111>
jp     1231 <tqli+0x111>
jmp    1244 <tqli+0x124>
jmp    1236 <tqli+0x116>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11af <tqli+0x8f>
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
je     16e6 <tqli+0x5c6>
mov    -0x2c(%rbp),%eax
mov    %eax,%ecx
add    $0x1,%ecx
mov    %ecx,-0x2c(%rbp)
cmp    $0x1e,%eax
jne    1270 <tqli+0x150>
lea    0xdb5(%rip),%rdi
call   1040 <erhand@plt>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  0xd79(%rip),%xmm1
mulsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd4d(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x58(%rbp),%xmm0
jbe    1354 <tqli+0x234>
movsd  -0x48(%rbp),%xmm0
movaps 0xcd4(%rip),%xmm1
pand   %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
jmp    136c <tqli+0x24c>
movsd  -0x48(%rbp),%xmm0
movaps 0xca0(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x80(%rbp),%xmm2
movsd  -0x98(%rbp),%xmm3
ff
addsd  %xmm3,%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xc6e(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0xc61(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x28(%rbp),%eax
jl     16a5 <tqli+0x585>
movsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x70(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x60(%rbp),%xmm0
movaps 0xbf5(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  -0x58(%rbp),%xmm1
movaps 0xbe5(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jb     1499 <tqli+0x379>
movsd  -0x58(%rbp),%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xbc2(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  -0x60(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0xb90(%rip),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mulsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    1504 <tqli+0x3e4>
movsd  -0x60(%rbp),%xmm0
divsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xb52(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  -0x58(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0xb20(%rip),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  0xad9(%rip),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x58(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x78(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1692 <tqli+0x572>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x70(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    15bb <tqli+0x49b>
jmp    1697 <tqli+0x577>
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    13c5 <tqli+0x2a5>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    16eb <tqli+0x5cb>
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jne    11a9 <tqli+0x89>
jmp    16fc <tqli+0x5dc>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1196 <tqli+0x76>
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tqli>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %rcx,%rbx
mov    %rsi,%r15
mov    %rdi,%r12
cmp    $0x2,%edx
mov    %rdx,(%rsp)
jl     117d <tqli+0x4d>
lea    0x8(%r15),%rsi
lea    -0x1(%rdx),%ebp
lea    0x0(,%rbp,8),%rdx
mov    %r15,%rdi
call   1030 <memmove@plt>
mov    (%rsp),%rdx
mov    %rbp,0x40(%rsp)
movq   $0x0,(%r15,%rbp,8)
jmp    1199 <tqli+0x69>
movslq %edx,%rax
movq   $0x0,-0x8(%r15,%rax,8)
xor    %eax,%eax
mov    %rax,0x40(%rsp)
cmp    $0x1,%edx
jne    167c <tqli+0x54c>
mov    %edx,%ebp
mov    %rbp,0x48(%rsp)
and    $0xfffffffe,%ebp
mov    $0xffffffffffffffff,%r14
xor    %esi,%esi
movsd  0xe6c(%rip),%xmm5
jmp    11d2 <tqli+0xa2>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r14
mov    %r13,%rsi
cmp    0x48(%rsp),%r13
je     167c <tqli+0x54c>
mov    0x40(%rsp),%rax
cmp    %eax,%esi
cmovg  %esi,%eax
mov    %rax,0x88(%rsp)
neg    %rax
mov    %rax,0x90(%rsp)
lea    0x1(%rsi),%r13
movl   $0x0,0xc(%rsp)
mov    %rsi,0x10(%rsp)
mov    %r14,0x80(%rsp)
mov    %r13,0x78(%rsp)
jmp    1241 <tqli+0x111>
nop
movsd  (%r12,%rsi,8),%xmm1
mov    0x80(%rsp),%r14
mov    0x78(%rsp),%r13
subsd  %xmm8,%xmm1
addl   $0x1,0xc(%rsp)
movsd  %xmm1,(%r12,%rsi,8)
movsd  %xmm6,(%r15,%rsi,8)
movq   $0x0,(%r15,%rdi,8)
mov    %r14,%rax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x90(%rsp),%rcx
add    %rax,%rcx
cmp    $0xffffffffffffffff,%rcx
je     12b0 <tqli+0x180>
movdqu 0x8(%r12,%rax,8),%xmm0
movdqa 0xd90(%rip),%xmm2
pand   %xmm2,%xmm0
pshufd $0xee,%xmm0,%xmm1
addsd  %xmm0,%xmm1
movq   0x8(%r15,%rax,8),%xmm0
pand   %xmm2,%xmm0
addsd  %xmm1,%xmm0
add    $0x1,%rax
ucomisd %xmm1,%xmm0
jne    1250 <tqli+0x120>
jp     1250 <tqli+0x120>
mov    %eax,%edi
cmp    %rsi,%rdi
jne    12c1 <tqli+0x191>
jmp    11c0 <tqli+0x90>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x88(%rsp),%rdi
cmp    %rsi,%rdi
je     11c0 <tqli+0x90>
cmpl   $0x1e,0xc(%rsp)
mov    %rdi,0x18(%rsp)
jne    12ef <tqli+0x1bf>
lea    0xd5c(%rip),%rdi
call   1050 <erhand@plt>
mov    0x18(%rsp),%rdi
mov    0x10(%rsp),%rsi
movsd  0xd35(%rip),%xmm5
mov    (%rsp),%rdx
movsd  (%r12,%r13,8),%xmm7
subsd  (%r12,%rsi,8),%xmm7
movsd  (%r15,%rsi,8),%xmm0
addsd  %xmm0,%xmm0
divsd  %xmm0,%xmm7
movapd %xmm7,%xmm0
mulsd  %xmm7,%xmm0
addsd  %xmm5,%xmm0
ucomisd 0xd0b(%rip),%xmm0
jb     1330 <tqli+0x200>
sqrtsd %xmm0,%xmm0
jmp    1357 <tqli+0x227>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm7,0x30(%rsp)
call   1040 <sqrt@plt>
movsd  0x30(%rsp),%xmm7
mov    0x18(%rsp),%rdi
mov    0x10(%rsp),%rsi
movsd  0xccd(%rip),%xmm5
mov    (%rsp),%rdx
movsd  (%r12,%rdi,8),%xmm2
movsd  (%r12,%rsi,8),%xmm1
subsd  %xmm1,%xmm2
movsd  (%r15,%rsi,8),%xmm6
movapd %xmm0,%xmm3
andpd  0xc87(%rip),%xmm3
orpd   0xc8f(%rip),%xmm0
xorpd  %xmm8,%xmm8
movapd %xmm7,%xmm4
cmpltsd %xmm8,%xmm4
andpd  %xmm4,%xmm0
andnpd %xmm3,%xmm4
orpd   %xmm0,%xmm4
addsd  %xmm7,%xmm4
divsd  %xmm4,%xmm6
addsd  %xmm2,%xmm6
cmp    %rsi,%rdi
jbe    1223 <tqli+0xf3>
xorpd  %xmm8,%xmm8
mov    %rdi,%r13
movapd %xmm5,%xmm9
movapd %xmm5,%xmm7
jmp    13e7 <tqli+0x2b7>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mulsd  %xmm1,%xmm6
subsd  %xmm9,%xmm6
movapd %xmm1,%xmm9
cmp    %rsi,%r13
jle    1210 <tqli+0xe0>
mov    %r13,%r14
add    $0xffffffffffffffff,%r13
movsd  -0x8(%r15,%r14,8),%xmm3
mulsd  %xmm3,%xmm7
movapd %xmm7,%xmm0
movapd 0xbfb(%rip),%xmm2
andpd  %xmm2,%xmm0
movapd %xmm6,%xmm1
andpd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jae    1440 <tqli+0x310>
divsd  %xmm6,%xmm7
movapd %xmm7,%xmm0
mulsd  %xmm7,%xmm0
addsd  %xmm5,%xmm0
ucomisd 0xbf9(%rip),%xmm0
jb     1470 <tqli+0x340>
sqrtsd %xmm0,%xmm0
jmp    14cb <tqli+0x39b>
nopw   0x0(%rax,%rax,1)
divsd  %xmm7,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm6,%xmm0
addsd  %xmm5,%xmm0
ucomisd 0xbd0(%rip),%xmm0
jb     14f0 <tqli+0x3c0>
sqrtsd %xmm0,%xmm0
jmp    154b <tqli+0x41b>
nopw   0x0(%rax,%rax,1)
movapd %xmm6,0x30(%rsp)
movapd %xmm8,0x60(%rsp)
movsd  %xmm9,0x28(%rsp)
movapd %xmm7,0x50(%rsp)
movsd  %xmm3,0x20(%rsp)
call   1040 <sqrt@plt>
movsd  0x20(%rsp),%xmm3
movapd 0x50(%rsp),%xmm7
movsd  0x28(%rsp),%xmm9
movapd 0x60(%rsp),%xmm8
movapd 0x30(%rsp),%xmm6
mov    0x18(%rsp),%rdi
mov    0x10(%rsp),%rsi
movsd  0xb59(%rip),%xmm5
mov    (%rsp),%rdx
mulsd  %xmm0,%xmm6
movapd %xmm5,%xmm1
divsd  %xmm0,%xmm1
mulsd  %xmm1,%xmm7
unpcklpd %xmm7,%xmm1
movapd %xmm6,%xmm7
jmp    1563 <tqli+0x433>
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm6,0x30(%rsp)
movapd %xmm8,0x60(%rsp)
movsd  %xmm9,0x28(%rsp)
movapd %xmm7,0x50(%rsp)
movsd  %xmm3,0x20(%rsp)
call   1040 <sqrt@plt>
movsd  0x20(%rsp),%xmm3
movapd 0x50(%rsp),%xmm7
movsd  0x28(%rsp),%xmm9
movapd 0x60(%rsp),%xmm8
movapd 0x30(%rsp),%xmm6
mov    0x18(%rsp),%rdi
mov    0x10(%rsp),%rsi
movsd  0xad9(%rip),%xmm5
mov    (%rsp),%rdx
mulsd  %xmm0,%xmm7
movapd %xmm5,%xmm1
divsd  %xmm0,%xmm1
mulsd  %xmm1,%xmm6
unpcklpd %xmm1,%xmm6
movapd %xmm6,%xmm1
mulsd  %xmm3,%xmm9
movsd  %xmm7,(%r15,%r14,8)
movsd  (%r12,%r14,8),%xmm0
subsd  %xmm8,%xmm0
movsd  (%r12,%r13,8),%xmm6
subsd  %xmm0,%xmm6
movapd %xmm1,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm9,%xmm2
movapd %xmm1,%xmm7
unpckhpd %xmm1,%xmm7
mulsd  %xmm7,%xmm6
addsd  %xmm2,%xmm6
movapd %xmm7,%xmm8
mulsd  %xmm6,%xmm8
addsd  %xmm8,%xmm0
movsd  %xmm0,(%r12,%r14,8)
movapd %xmm7,%xmm0
xorpd  0xa4f(%rip),%xmm0
cmp    $0x1,%edx
jne    15d0 <tqli+0x4a0>
xor    %eax,%eax
jmp    163e <tqli+0x50e>
nopw   0x0(%rax,%rax,1)
movapd %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
mov    (%rbx,%rax,8),%rcx
movsd  (%rcx,%r14,8),%xmm3
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
movsd  (%rcx,%r13,8),%xmm4
unpcklpd %xmm4,%xmm4
mulpd  %xmm1,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%rcx,%r13,8)
mov    0x8(%rbx,%rax,8),%rcx
movsd  (%rcx,%r14,8),%xmm3
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
movsd  (%rcx,%r13,8),%xmm4
unpcklpd %xmm4,%xmm4
mulpd  %xmm1,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,(%rcx,%r13,8)
add    $0x2,%rax
cmp    %rax,%rbp
jne    15e0 <tqli+0x4b0>
testb  $0x1,0x48(%rsp)
je     13d0 <tqli+0x2a0>
mov    (%rbx,%rax,8),%rax
movsd  (%rax,%r14,8),%xmm2
unpcklpd %xmm2,%xmm2
unpcklpd %xmm1,%xmm0
mulpd  %xmm2,%xmm0
movsd  (%rax,%r13,8),%xmm2
unpcklpd %xmm2,%xmm2
mulpd  %xmm1,%xmm2
addpd  %xmm0,%xmm2
movupd %xmm2,(%rax,%r13,8)
jmp    13d0 <tqli+0x2a0>
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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