qrsolv(n, r, ldr, ipvt, diag, qtb, x, sdiag, wa) int n, ldr;
int ipvt[];
double r[], diag[], qtb[], x[], sdiag[], wa[];
{
  int i, ij, ik, kk, j, jp1, k, kp1, l, nsing;
  double cos, cotan, qtbpj, sin, sum, tan, temp;
  static double zero = 0.0;
  static double p25 = 0.25;
  static double p5 = 0.5;
  double fabs(), sqrt();
  kk = 0;
  for (j = 0; j < n; j++) {
    ij = kk;
    ik = kk;
    for (i = j; i < n; i++) {
      r[ij] = r[ik];
      ij += 1;
      ik += ldr;
    }
    x[j] = r[kk];
    wa[j] = qtb[j];
    kk += ldr + 1;
  }
  for (j = 0; j < n; j++) {
    l = ipvt[j];
    if (diag[l] == zero)
      goto L90;
    for (k = j; k < n; k++)
      sdiag[k] = zero;
    sdiag[j] = diag[l];
    qtbpj = zero;
    for (k = j; k < n; k++) {
      if (sdiag[k] == zero)
        continue;
      kk = k + ldr * k;
      if (fabs(r[kk]) < fabs(sdiag[k])) {
        cotan = r[kk] / sdiag[k];
        sin = p5 / sqrt(p25 + p25 * cotan * cotan);
        cos = sin * cotan;
      } else {
        tan = sdiag[k] / r[kk];
        cos = p5 / sqrt(p25 + p25 * tan * tan);
        sin = cos * tan;
      }
      r[kk] = cos * r[kk] + sin * sdiag[k];
      temp = cos * wa[k] + sin * qtbpj;
      qtbpj = -sin * wa[k] + cos * qtbpj;
      wa[k] = temp;
      kp1 = k + 1;
      if (n > kp1) {
        ik = kk + 1;
        for (i = kp1; i < n; i++) {
          temp = cos * r[ik] + sin * sdiag[i];
          sdiag[i] = -sin * r[ik] + cos * sdiag[i];
          r[ik] = temp;
          ik += 1;
        }
      }
    }
  L90:
    kk = j + ldr * j;
    sdiag[j] = r[kk];
    r[kk] = x[j];
  }
  nsing = n;
  for (j = 0; j < n; j++) {
    if ((sdiag[j] == zero) && (nsing == n))
      nsing = j;
    if (nsing < n)
      wa[j] = zero;
  }
  if (nsing < 1)
    goto L150;
  for (k = 0; k < nsing; k++) {
    j = nsing - k - 1;
    sum = zero;
    jp1 = j + 1;
    if (nsing > jp1) {
      ij = jp1 + ldr * j;
      for (i = jp1; i < nsing; i++) {
        sum += r[ij] * wa[i];
        ij += 1;
      }
    }
    wa[j] = (wa[j] - sum) / sdiag[j];
  }
L150:
  for (j = 0; j < n; j++) {
    l = ipvt[j];
    x[l] = wa[j];
  }
}


x86 assembly (clang, -O0):
<qrsolv>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11fd <qrsolv+0xed>
mov    -0x40(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11af <qrsolv+0x9f>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x14(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1169 <qrsolv+0x59>
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    114b <qrsolv+0x3b>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    162c <qrsolv+0x51c>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x54(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd 0x2e0d(%rip),%xmm0
jne    1244 <qrsolv+0x134>
jp     1244 <qrsolv+0x134>
jmp    15db <qrsolv+0x4cb>
mov    -0x44(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1279 <qrsolv+0x169>
movsd  0x2de2(%rip),%xmm0
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    124a <qrsolv+0x13a>
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x2da5(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    15d6 <qrsolv+0x4c6>
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd 0x2d79(%rip),%xmm0
jne    12d8 <qrsolv+0x1c8>
jp     12d8 <qrsolv+0x1c8>
jmp    15c8 <qrsolv+0x4b8>
mov    -0x4c(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xd05(%rip),%xmm0
pand   %xmm0,%xmm1
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xced(%rip),%xmm2
pand   %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    139f <qrsolv+0x28f>
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  0x2ce8(%rip),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0x2cd0(%rip),%xmm1
movsd  0x2cc8(%rip),%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x98(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    1424 <qrsolv+0x314>
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0x2c67(%rip),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  0x2c4f(%rip),%xmm1
movsd  0x2c47(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xa0(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x78(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x78(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
mov    0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0x50(%rbp),%eax
jle    15c3 <qrsolv+0x4b3>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    15be <qrsolv+0x4ae>
movsd  -0x60(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x78(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x60(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1502 <qrsolv+0x3f2>
jmp    15c3 <qrsolv+0x4b3>
jmp    15c8 <qrsolv+0x4b8>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12a6 <qrsolv+0x196>
jmp    15db <qrsolv+0x4cb>
mov    -0x44(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x44(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1204 <qrsolv+0xf4>
mov    -0x8(%rbp),%eax
mov    %eax,-0x58(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    16ac <qrsolv+0x59c>
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd 0x29e6(%rip),%xmm0
jne    1678 <qrsolv+0x568>
jp     1678 <qrsolv+0x568>
mov    -0x58(%rbp),%eax
cmp    -0x8(%rbp),%eax
jne    1678 <qrsolv+0x568>
mov    -0x44(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1699 <qrsolv+0x589>
movsd  0x29b4(%rip),%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    169e <qrsolv+0x58e>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1639 <qrsolv+0x529>
cmpl   $0x1,-0x58(%rbp)
jge    16bb <qrsolv+0x5ab>
jmp    17a4 <qrsolv+0x694>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x58(%rbp),%eax
jge    179f <qrsolv+0x68f>
mov    -0x58(%rbp),%eax
sub    -0x4c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x44(%rbp)
movsd  0x295e(%rip),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x48(%rbp),%eax
jle    1765 <qrsolv+0x655>
mov    -0x48(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x44(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x38(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x58(%rbp),%eax
jge    1760 <qrsolv+0x650>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x80(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1711 <qrsolv+0x601>
jmp    1765 <qrsolv+0x655>
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x80(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    16c2 <qrsolv+0x5b2>
jmp    17a4 <qrsolv+0x694>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    17ed <qrsolv+0x6dd>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x54(%rbp)
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    17ab <qrsolv+0x69b>
mov    -0x4(%rbp),%eax
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qrsolv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xd8,%rsp
mov    %r9,0x40(%rsp)
mov    %r8,0xa8(%rsp)
mov    %rcx,0x70(%rsp)
mov    %rdx,0x90(%rsp)
test   %edi,%edi
jle    1b38 <qrsolv+0xa18>
mov    %rsi,%r13
mov    0x90(%rsp),%rax
lea    0x1(%rax),%ecx
mov    %ecx,0x64(%rsp)
movslq %ecx,%rcx
movslq %eax,%rbx
mov    %edi,%r11d
mov    %rcx,0x38(%rsp)
lea    0x0(,%rcx,8),%rax
mov    %rax,0x10(%rsp)
lea    (%rbx,%rbx,2),%rax
mov    %rax,0x30(%rsp)
lea    0x0(,%rbx,4),%rsi
lea    (%rbx,%rbx,1),%rax
mov    %rax,0x28(%rsp)
xor    %r15d,%r15d
mov    %edi,%r14d
mov    %r13,%rdx
xor    %eax,%eax
mov    %rax,0x50(%rsp)
xor    %r12d,%r12d
xor    %r8d,%r8d
mov    %r11,0x20(%rsp)
mov    %edi,0xc(%rsp)
jmp    1223 <qrsolv+0x103>
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%r8,8),%xmm0
mov    0x110(%rsp),%rax
movsd  %xmm0,(%rax,%r12,8)
mov    0x40(%rsp),%rax
movsd  (%rax,%r12,8),%xmm0
mov    0x120(%rsp),%r9
movsd  %xmm0,(%r9,%r12,8)
mov    0x38(%rsp),%rax
add    %rax,%r8
add    $0x1,%r12
sub    %rax,0x50(%rsp)
add    0x10(%rsp),%rdx
mov    0x18(%rsp),%r15
add    $0xffffffffffffffff,%r15
add    $0x3,%r14b
mov    0x20(%rsp),%r11
cmp    %r11,%r12
mov    0xc(%rsp),%edi
je     130c <qrsolv+0x1ec>
movzbl %r14b,%r14d
and    $0x3,%r14d
mov    %edi,%r10d
sub    %r12d,%r10d
mov    %r12d,%r9d
not    %r9d
add    %edi,%r9d
mov    %r8,%rax
mov    %r8,%rbp
mov    %r12d,%ecx
test   $0x3,%r10b
je     1277 <qrsolv+0x157>
xor    %ecx,%ecx
xor    %eax,%eax
nopl   (%rax)
movsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
add    %rbx,%rax
add    $0x1,%rcx
cmp    %ecx,%r14d
jne    1250 <qrsolv+0x130>
mov    0x50(%rsp),%r10
sub    %r10,%rax
mov    %rcx,%rbp
sub    %r10,%rbp
sub    %r15d,%ecx
mov    %r15,0x18(%rsp)
cmp    $0x3,%r9d
jb     11c0 <qrsolv+0xa0>
mov    0x30(%rsp),%rdi
lea    (%rdi,%rax,1),%r10
mov    0x28(%rsp),%rdi
lea    (%rdi,%rax,1),%r11
lea    (%rbx,%rax,1),%r9
add    $0x3,%rbp
mov    0xc(%rsp),%r15d
sub    %ecx,%r15d
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rax,8),%xmm0
movsd  %xmm0,-0x18(%r13,%rbp,8)
add    %rbx,%rax
movsd  0x0(%r13,%r9,8),%xmm0
movsd  %xmm0,-0x10(%r13,%rbp,8)
add    %rbx,%rax
movsd  0x0(%r13,%r11,8),%xmm0
movsd  %xmm0,-0x8(%r13,%rbp,8)
add    %rbx,%rax
movsd  0x0(%r13,%r10,8),%xmm0
movsd  %xmm0,0x0(%r13,%rbp,8)
add    %rbx,%rax
add    %rsi,%r10
add    %rsi,%r11
add    %rsi,%r9
add    $0x4,%rbp
add    $0xfffffffc,%r15d
jne    12b0 <qrsolv+0x190>
jmp    11c0 <qrsolv+0xa0>
test   %edi,%edi
jle    1b38 <qrsolv+0xa18>
lea    0x8(%r13),%rax
mov    %rax,0xb8(%rsp)
lea    0x10(%r13),%rax
mov    %rax,0xd0(%rsp)
mov    0x118(%rsp),%rsi
lea    0x8(%rsi),%rax
mov    %rax,0xa0(%rsp)
xor    %eax,%eax
mov    %rax,0x78(%rsp)
mov    $0x1,%eax
mov    %rax,0x98(%rsp)
movl   $0x2,0x60(%rsp)
movl   $0xfffffffe,0x5c(%rsp)
ff
xorpd  %xmm8,%xmm8
xor    %edx,%edx
mov    %rbx,0x48(%rsp)
jmp    13f0 <qrsolv+0x2d0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x88(%rsp),%rdx
mov    %rdx,%rax
imul   %rbx,%rax
add    %rdx,%rax
movsd  0x0(%r13,%rax,8),%xmm0
mov    0xb0(%rsp),%rcx
movsd  %xmm0,(%rcx)
mov    0x110(%rsp),%rcx
movsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,0x0(%r13,%rax,8)
add    $0x1,%rdx
mov    0x78(%rsp),%rax
add    0x64(%rsp),%eax
mov    %rax,0x78(%rsp)
addq   $0x1,0x98(%rsp)
addl   $0x1,0x60(%rsp)
addl   $0xffffffff,0x5c(%rsp)
addq   $0x8,0xa0(%rsp)
cmp    %r11,%rdx
je     18f7 <qrsolv+0x7d7>
lea    (%rsi,%rdx,8),%rax
mov    %rax,0xb0(%rsp)
mov    0x70(%rsp),%rax
mov    %rdx,0x88(%rsp)
movslq (%rax,%rdx,4),%rbp
mov    0xa8(%rsp),%rax
movsd  (%rax,%rbp,8),%xmm0
ucomisd %xmm8,%xmm0
jne    1427 <qrsolv+0x307>
jnp    1380 <qrsolv+0x260>
mov    0x88(%rsp),%r15
lea    0x2(%r15),%r12
mov    %r15d,%eax
not    %eax
add    %edi,%eax
lea    0x8(,%rax,8),%rdx
mov    0xb0(%rsp),%rbx
mov    %rbx,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    %r12,%r10
xorpd  %xmm8,%xmm8
mov    0x118(%rsp),%rsi
mov    0xc(%rsp),%edi
mov    0x20(%rsp),%r11
mov    0x120(%rsp),%r9
mov    0xa8(%rsp),%rax
movsd  (%rax,%rbp,8),%xmm10
movsd  %xmm10,(%rbx)
xorpd  %xmm9,%xmm9
mov    0xa0(%rsp),%rbp
mov    0x5c(%rsp),%r8d
mov    0x60(%rsp),%ecx
xor    %edx,%edx
mov    0x98(%rsp),%rax
mov    %rax,0x18(%rsp)
mov    0x78(%rsp),%rax
mov    %eax,0x50(%rsp)
mov    %r15,%rax
mov    %r12,0x68(%rsp)
xchg   %ax,%ax
cmp    %ecx,%edi
mov    %ecx,%r15d
cmovg  %edi,%r15d
lea    (%r10,%rdx,1),%r14d
cmp    %r14d,%edi
cmovg  %edi,%r14d
lea    (%r10,%rdx,1),%rbx
cmp    %ebx,%edi
mov    %ebx,%r12d
cmovg  %edi,%r12d
ucomisd %xmm8,%xmm10
jne    1510 <qrsolv+0x3f0>
jp     1510 <qrsolv+0x3f0>
add    $0x1,%rax
mov    0x48(%rsp),%rbx
cmp    %r11,%rax
jne    1820 <qrsolv+0x700>
jmp    1380 <qrsolv+0x260>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbx,0x10(%rsp)
mov    %r8d,0x28(%rsp)
mov    %rbp,0x30(%rsp)
mov    %rax,%rbp
imul   0x90(%rsp),%eax
movslq %eax,%rbx
mov    %rbp,%rax
add    %rbp,%rbx
movsd  0x0(%r13,%rbx,8),%xmm11
movapd %xmm11,%xmm0
movapd 0xab9(%rip),%xmm2
andpd  %xmm2,%xmm0
movapd %xmm10,%xmm1
andpd  %xmm2,%xmm1
ucomisd %xmm0,%xmm1
mov    %ecx,0x40(%rsp)
mov    %rdx,0x38(%rsp)
jbe    1592 <qrsolv+0x472>
divsd  %xmm10,%xmm11
movapd %xmm11,%xmm0
movsd  0xaab(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm11,%xmm0
addsd  %xmm1,%xmm0
ucomisd %xmm8,%xmm0
jb     15c1 <qrsolv+0x4a1>
sqrtsd %xmm0,%xmm0
jmp    161d <qrsolv+0x4fd>
divsd  %xmm11,%xmm10
movapd %xmm10,%xmm0
movsd  0xa7c(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm10,%xmm0
addsd  %xmm1,%xmm0
ucomisd %xmm8,%xmm0
jb     1637 <qrsolv+0x517>
sqrtsd %xmm0,%xmm0
jmp    1693 <qrsolv+0x573>
movsd  %xmm9,0x80(%rsp)
movapd %xmm11,0xc0(%rsp)
call   1040 <sqrt@plt>
movapd 0xc0(%rsp),%xmm11
mov    %rbp,%rax
mov    0x38(%rsp),%rdx
mov    0x40(%rsp),%ecx
movsd  0x80(%rsp),%xmm9
mov    0x68(%rsp),%r10
xorpd  %xmm8,%xmm8
mov    0x118(%rsp),%rsi
mov    0xc(%rsp),%edi
mov    0x20(%rsp),%r11
mov    0x120(%rsp),%r9
mov    0x10(%rsp),%rbp
movsd  0x9fd(%rip),%xmm10
divsd  %xmm0,%xmm10
mulsd  %xmm10,%xmm11
jmp    16ab <qrsolv+0x58b>
movsd  %xmm9,0x80(%rsp)
movapd %xmm10,0xc0(%rsp)
call   1040 <sqrt@plt>
movapd 0xc0(%rsp),%xmm10
mov    %rbp,%rax
mov    0x38(%rsp),%rdx
mov    0x40(%rsp),%ecx
movsd  0x80(%rsp),%xmm9
mov    0x68(%rsp),%r10
xorpd  %xmm8,%xmm8
mov    0x118(%rsp),%rsi
mov    0xc(%rsp),%edi
mov    0x20(%rsp),%r11
mov    0x120(%rsp),%r9
mov    0x10(%rsp),%rbp
movsd  0x987(%rip),%xmm11
divsd  %xmm0,%xmm11
mulsd  %xmm11,%xmm10
movsd  (%rsi,%rax,8),%xmm0
mulsd  %xmm10,%xmm0
movsd  0x0(%r13,%rbx,8),%xmm1
mulsd  %xmm11,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rbx,8)
movsd  (%r9,%rax,8),%xmm0
movapd %xmm9,%xmm1
mulsd  %xmm10,%xmm1
movapd %xmm11,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm11,%xmm9
mulsd  %xmm10,%xmm0
subsd  %xmm0,%xmm9
movsd  %xmm2,(%r9,%rax,8)
add    $0x1,%rax
cmp    %r11,%rax
jge    1782 <qrsolv+0x662>
mov    %rax,0x10(%rsp)
mov    0x88(%rsp),%rax
add    %edx,%eax
mov    $0xfffffffe,%ecx
sub    %eax,%ecx
movslq 0x50(%rsp),%r11
movapd %xmm10,%xmm0
xorpd  0x8e1(%rip),%xmm0
add    %ecx,%r14d
je     179f <qrsolv+0x67f>
add    %ecx,%r12d
lea    (%r11,%r12,1),%rbx
mov    0xb8(%rsp),%rax
lea    (%rax,%r11,8),%rdx
mov    0x18(%rsp),%rax
lea    (%rsi,%rax,8),%rax
mov    0xd0(%rsp),%rcx
lea    (%rcx,%rbx,8),%rcx
cmp    %rcx,%rax
jae    1848 <qrsolv+0x728>
lea    (%rsi,%rbp,8),%rcx
lea    (%rcx,%r12,8),%rcx
cmp    %rcx,%rdx
jae    1848 <qrsolv+0x728>
mov    0x18(%rsp),%rcx
mov    0x68(%rsp),%r10
jmp    17a4 <qrsolv+0x684>
mov    0x30(%rsp),%rbp
mov    0x48(%rsp),%rbx
mov    0x28(%rsp),%r8d
cmp    %r11,%rax
jne    1820 <qrsolv+0x700>
jmp    1380 <qrsolv+0x260>
mov    0x18(%rsp),%rcx
mov    0x30(%rsp),%rbp
mov    0x48(%rsp),%rbx
mov    0x28(%rsp),%r8d
mov    0xb8(%rsp),%rax
lea    (%rax,%r11,8),%rax
nop
movsd  (%rax),%xmm1
movsd  (%rsi,%rcx,8),%xmm2
movapd %xmm10,%xmm3
mulsd  %xmm2,%xmm3
movapd %xmm11,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm11,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
movsd  %xmm4,(%rax)
add    $0x1,%rcx
add    $0x8,%rax
cmp    %edi,%ecx
jl     17c0 <qrsolv+0x6a0>
mov    0x40(%rsp),%ecx
mov    0x38(%rsp),%rdx
mov    0x20(%rsp),%r11
mov    0x10(%rsp),%rax
cmp    %r11,%rax
je     1380 <qrsolv+0x260>
nopl   (%rax)
addq   $0x1,0x18(%rsp)
mov    0x64(%rsp),%ebx
add    %ebx,0x50(%rsp)
movsd  (%rsi,%rax,8),%xmm10
add    $0x1,%rdx
add    $0x1,%ecx
add    $0xffffffff,%r8d
add    $0x8,%rbp
jmp    14c0 <qrsolv+0x3a0>
mov    0x28(%rsp),%r8d
add    %r8d,%r15d
add    $0x1,%r15
and    $0xfffffffffffffffe,%r15
add    $0x1,%r14
mov    %r14,%r10
and    $0xfffffffffffffffe,%r10
mov    0x18(%rsp),%rax
lea    (%rax,%r10,1),%rcx
add    %r10,%r11
movapd %xmm10,%xmm1
unpcklpd %xmm10,%xmm1
movapd %xmm11,%xmm2
unpcklpd %xmm11,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
xor    %ebx,%ebx
mov    0x30(%rsp),%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rdx,%rbx,8),%xmm4
movupd 0x0(%rbp,%rbx,8),%xmm5
movapd %xmm1,%xmm6
mulpd  %xmm5,%xmm6
movapd %xmm2,%xmm7
mulpd  %xmm4,%xmm7
addpd  %xmm6,%xmm7
mulpd  %xmm2,%xmm5
mulpd  %xmm3,%xmm4
addpd  %xmm5,%xmm4
movupd %xmm4,0x0(%rbp,%rbx,8)
movupd %xmm7,(%rdx,%rbx,8)
add    $0x2,%rbx
cmp    %rbx,%r15
jne    18a0 <qrsolv+0x780>
cmp    %r10,%r14
mov    0x68(%rsp),%r10
mov    0x48(%rsp),%rbx
jne    17b3 <qrsolv+0x693>
jmp    1801 <qrsolv+0x6e1>
test   %edi,%edi
jle    1b38 <qrsolv+0xa18>
mov    %r11,%r8
add    $0xffffffffffffffff,%r8
je     197e <qrsolv+0x85e>
mov    %r11d,%ebp
and    $0xfffffffe,%ebp
xor    %edx,%edx
xorpd  %xmm0,%xmm0
mov    %edi,%esi
jmp    1929 <qrsolv+0x809>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdx
cmp    %rbp,%rdx
je     1978 <qrsolv+0x858>
mov    %rdx,%rcx
cmp    %edi,%esi
mov    %ecx,%eax
cmovne %esi,%eax
mov    0x118(%rsp),%rbx
ucomisd (%rbx,%rdx,8),%xmm0
cmovne %esi,%eax
cmovp  %esi,%eax
cmp    %edi,%eax
jge    1952 <qrsolv+0x832>
movq   $0x0,(%r9,%rcx,8)
lea    0x1(%rcx),%rdx
cmp    %edi,%eax
mov    %edx,%esi
cmovne %eax,%esi
ucomisd 0x8(%rbx,%rcx,8),%xmm0
cmovne %eax,%esi
cmovp  %eax,%esi
cmp    %edi,%esi
jge    1920 <qrsolv+0x800>
movq   $0x0,0x8(%r9,%rcx,8)
jmp    1920 <qrsolv+0x800>
add    $0x2,%rcx
jmp    1982 <qrsolv+0x862>
xor    %ecx,%ecx
mov    %edi,%esi
mov    %esi,%r14d
test   $0x1,%r11b
je     19ba <qrsolv+0x89a>
cmp    %edi,%esi
mov    %ecx,%r14d
cmovne %esi,%r14d
xorpd  %xmm0,%xmm0
mov    0x118(%rsp),%rax
ucomisd (%rax,%rcx,8),%xmm0
cmovne %esi,%r14d
cmovp  %esi,%r14d
cmp    %edi,%r14d
jge    19ba <qrsolv+0x89a>
movq   $0x0,(%r9,%rcx,8)
mov    %r11,%rcx
test   %r14d,%r14d
jle    1a6d <qrsolv+0x94d>
mov    %r14d,%edx
lea    -0x1(%r14),%esi
mov    0x90(%rsp),%r11
imul   %r11d,%esi
add    %r14d,%esi
not    %r11d
xor    %ebp,%ebp
jmp    1a21 <qrsolv+0x901>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorpd  %xmm0,%xmm0
movsd  (%r9,%rdi,8),%xmm1
subsd  %xmm0,%xmm1
mov    0x118(%rsp),%rax
divsd  (%rax,%rdi,8),%xmm1
movsd  %xmm1,(%r9,%rdi,8)
add    $0x1,%rbp
add    $0xffffffff,%r14d
add    %r11d,%esi
cmp    %rdx,%rbp
je     1a6d <qrsolv+0x94d>
mov    %rbp,%rdi
not    %rdi
add    %rdx,%rdi
test   %rbp,%rbp
je     19f0 <qrsolv+0x8d0>
movslq %esi,%rax
lea    0x0(,%rax,8),%rbx
add    %r13,%rbx
movslq %r14d,%rax
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rbx),%xmm1
mulsd  (%r9,%rax,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rax
add    $0x8,%rbx
cmp    %rdx,%rax
jl     1a50 <qrsolv+0x930>
jmp    19f4 <qrsolv+0x8d4>
cmpl   $0x0,0xc(%rsp)
jle    1b38 <qrsolv+0xa18>
mov    %ecx,%ebx
and    $0x3,%ebx
cmp    $0x3,%r8
jae    1a8f <qrsolv+0x96f>
xor    %eax,%eax
mov    0x110(%rsp),%rsi
jmp    1afb <qrsolv+0x9db>
and    $0xfffffffc,%ecx
xor    %eax,%eax
mov    0x110(%rsp),%rsi
mov    0x70(%rsp),%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movslq (%rdi,%rax,4),%rdx
movsd  (%r9,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
movslq 0x4(%rdi,%rax,4),%rdx
movsd  0x8(%r9,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
movslq 0x8(%rdi,%rax,4),%rdx
movsd  0x10(%r9,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
movslq 0xc(%rdi,%rax,4),%rdx
movsd  0x18(%r9,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
add    $0x4,%rax
cmp    %rax,%rcx
jne    1ab0 <qrsolv+0x990>
test   %rbx,%rbx
je     1b38 <qrsolv+0xa18>
mov    %rsi,%rbp
lea    (%r9,%rax,8),%rdx
mov    0x70(%rsp),%rsi
lea    (%rsi,%rax,4),%rax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq (%rax,%rsi,4),%rdi
movsd  (%rdx,%rsi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdi,8)
add    $0x1,%rsi
cmp    %rsi,%rbx
jne    1b20 <qrsolv+0xa00>
add    $0xd8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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