#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 */
scalar_t__ fabs(double);
double sqrt(double);
void qrsolv(int n, double *r, int ldr, const int *ipvt, const double *diag,
            const double *qtb, double *x, double *sdiag, double *wa) {
  int i, j, k, l;
  double cos, sin, sum, temp;
  int nsing;
  double qtbpj;
  for (j = 0; j < n; ++j) {
    for (i = j; i < n; ++i) {
      r[i + j * ldr] = r[j + i * ldr];
    }
    x[j] = r[j + j * ldr];
    wa[j] = qtb[j];
  }
  for (j = 0; j < n; ++j) {
    l = ipvt[j] - 1;
    if (diag[l] != 0.) {
      for (k = j; k < n; ++k) {
        sdiag[k] = 0.;
      }
      sdiag[j] = diag[l];
      qtbpj = 0.;
      for (k = j; k < n; ++k) {
        if (sdiag[k] != 0.) {
          if (fabs(r[k + k * ldr]) < fabs(sdiag[k])) {
            double cotan;
            cotan = r[k + k * ldr] / sdiag[k];
            sin = .5 / sqrt(.25 + .25 * (cotan * cotan));
            cos = sin * cotan;
          } else {
            double tan;
            tan = sdiag[k] / r[k + k * ldr];
            cos = .5 / sqrt(.25 + .25 * (tan * tan));
            sin = cos * tan;
          }
          temp = cos * wa[k] + sin * qtbpj;
          qtbpj = -sin * wa[k] + cos * qtbpj;
          wa[k] = temp;
          r[k + k * ldr] = cos * r[k + k * ldr] + sin * sdiag[k];
          if (n > k + 1) {
            for (i = k + 1; i < n; ++i) {
              temp = cos * r[i + k * ldr] + sin * sdiag[i];
              sdiag[i] = -sin * r[i + k * ldr] + cos * sdiag[i];
              r[i + k * ldr] = temp;
            }
          }
        }
      }
    }
    sdiag[j] = r[j + j * ldr];
    r[j + j * ldr] = x[j];
  }
  nsing = n;
  for (j = 0; j < n; ++j) {
    if (sdiag[j] == 0. && nsing == n) {
      nsing = j;
    }
    if (nsing < n) {
      wa[j] = 0.;
    }
  }
  if (nsing >= 1) {
    for (k = 1; k <= nsing; ++k) {
      j = nsing - k;
      sum = 0.;
      if (nsing > j + 1) {
        for (i = j + 1; i < nsing; ++i) {
          sum += r[i + j * ldr] * wa[i];
        }
      }
      wa[j] = (wa[j] - sum) / sdiag[j];
    }
  }
  for (j = 0; j < n; ++j) {
    l = ipvt[j] - 1;
    x[l] = wa[j];
  }
  return;
}


x86 assembly (clang, -O0):
<qrsolv>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11fd <qrsolv+0xdd>
mov    -0x38(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11b0 <qrsolv+0x90>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x34(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1166 <qrsolv+0x46>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1154 <qrsolv+0x34>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1621 <qrsolv+0x501>
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
sub    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1246 <qrsolv+0x126>
jp     1246 <qrsolv+0x126>
jmp    15c9 <qrsolv+0x4a9>
mov    -0x38(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1276 <qrsolv+0x156>
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    124c <qrsolv+0x12c>
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15c4 <qrsolv+0x4a4>
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12cf <qrsolv+0x1af>
jp     12cf <qrsolv+0x1af>
jmp    15b1 <qrsolv+0x491>
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
mov    %rax,-0x88(%rbp)
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
mov    %rax,%rcx
mov    -0x88(%rbp),%rax
cmp    %rcx,%rax
jge    1389 <qrsolv+0x269>
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  0xcb4(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xc94(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    13f5 <qrsolv+0x2d5>
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
movsd  0xc43(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xc23(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%eax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jle    15ac <qrsolv+0x48c>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15a7 <qrsolv+0x487>
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%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  -0x60(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    14d9 <qrsolv+0x3b9>
jmp    15ac <qrsolv+0x48c>
jmp    15b1 <qrsolv+0x491>
jmp    15b6 <qrsolv+0x496>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    129e <qrsolv+0x17e>
jmp    15c9 <qrsolv+0x4a9>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1204 <qrsolv+0xe4>
mov    -0x4(%rbp),%eax
mov    %eax,-0x64(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    169b <qrsolv+0x57b>
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    166c <qrsolv+0x54c>
jp     166c <qrsolv+0x54c>
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jne    166c <qrsolv+0x54c>
mov    -0x38(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1688 <qrsolv+0x568>
mov    0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    168d <qrsolv+0x56d>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    162e <qrsolv+0x50e>
cmpl   $0x1,-0x64(%rbp)
jl     1778 <qrsolv+0x658>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x64(%rbp),%eax
jg     1773 <qrsolv+0x653>
mov    -0x64(%rbp),%eax
sub    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x64(%rbp),%eax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jle    1739 <qrsolv+0x619>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x64(%rbp),%eax
jge    1734 <qrsolv+0x614>
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    16e3 <qrsolv+0x5c3>
jmp    1739 <qrsolv+0x619>
mov    0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x58(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    16ac <qrsolv+0x58c>
jmp    1778 <qrsolv+0x658>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    17c4 <qrsolv+0x6a4>
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
sub    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    177f <qrsolv+0x65f>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qrsolv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xe8,%rsp
mov    %r8,0x80(%rsp)
mov    %edx,0x14(%rsp)
mov    %rsi,0x8(%rsp)
test   %edi,%edi
jle    1a6f <qrsolv+0x93f>
mov    0x128(%rsp),%rbx
movslq 0x14(%rsp),%r15
mov    %edi,0x10(%rsp)
mov    %edi,%edx
lea    0x8(,%r15,8),%r8
xor    %eax,%eax
mov    %rdx,%r10
mov    0x8(%rsp),%rsi
nop
xor    %edi,%edi
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rsi,%rbp,8)
add    $0x1,%rbp
add    %r15,%rdi
cmp    %rbp,%rdx
jne    1190 <qrsolv+0x60>
mov    %eax,%edi
imul   0x14(%rsp),%edi
movslq %edi,%rdi
add    %rax,%rdi
mov    0x8(%rsp),%rbp
movsd  0x0(%rbp,%rdi,8),%xmm0
mov    0x120(%rsp),%rdi
movsd  %xmm0,(%rdi,%rax,8)
movsd  (%r9,%rax,8),%xmm0
mov    0x130(%rsp),%rdi
movsd  %xmm0,(%rdi,%rax,8)
add    $0x1,%rax
add    %r8,%rsi
add    $0xffffffffffffffff,%rdx
cmp    %r10,%rax
jne    1180 <qrsolv+0x50>
cmpl   $0x0,0x10(%rsp)
jle    1a6f <qrsolv+0x93f>
mov    0x8(%rsp),%rax
add    $0x8,%rax
lea    0x8(%rbx),%rdx
mov    %rdx,0xc0(%rsp)
mov    $0x1,%edx
mov    %rdx,0x78(%rsp)
mov    $0x2,%edx
mov    %rdx,0x70(%rsp)
movq   $0xffffffffffffffff,0x68(%rsp)
ff ff
xor    %edx,%edx
mov    %rax,0xb8(%rsp)
mov    %rax,0x60(%rsp)
mov    %r10,%rbp
mov    %rcx,0x58(%rsp)
mov    %r15,0x90(%rsp)
mov    %r10,0x30(%rsp)
jmp    12c1 <qrsolv+0x191>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x50(%rsp),%rdx
mov    %rdx,%rax
imul   %r15,%rax
add    %rdx,%rax
mov    0x8(%rsp),%rsi
movsd  (%rsi,%rax,8),%xmm0
mov    0x88(%rsp),%rcx
movsd  %xmm0,(%rcx)
mov    0x120(%rsp),%rcx
movsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
add    $0x1,%rdx
addq   $0x1,0x78(%rsp)
addq   $0x1,0x70(%rsp)
addq   $0xffffffffffffffff,0x68(%rsp)
addq   $0x8,0x60(%rsp)
cmp    %rbp,%rdx
mov    0x58(%rsp),%rcx
je     17f5 <qrsolv+0x6c5>
lea    (%rbx,%rdx,8),%rax
mov    %rax,0x88(%rsp)
mov    %rdx,0x50(%rsp)
movslq (%rcx,%rdx,4),%r13
mov    0x80(%rsp),%rax
movsd  -0x8(%rax,%r13,8),%xmm0
ucomisd 0xd23(%rip),%xmm0
jne    12f5 <qrsolv+0x1c5>
jnp    1260 <qrsolv+0x130>
mov    0x50(%rsp),%r12
lea    0x2(%r12),%rax
mov    %rax,0xd8(%rsp)
mov    %r12,%rax
not    %rax
mov    %rax,0xc8(%rsp)
mov    %r12d,%eax
not    %eax
add    0x10(%rsp),%eax
lea    0x8(,%rax,8),%rdx
mov    %rbp,%r14
mov    0x88(%rsp),%rbp
mov    %rbp,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x80(%rsp),%rax
movsd  -0x8(%rax,%r13,8),%xmm0
movsd  %xmm0,0x0(%rbp)
mov    %r14,%rbp
xorpd  %xmm8,%xmm8
mov    0x60(%rsp),%r9
mov    0x68(%rsp),%r11
mov    0x70(%rsp),%rdi
xor    %r8d,%r8d
mov    0x78(%rsp),%rax
mov    %rax,0x28(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmp    %rbp,%rdi
mov    %rbp,%rax
cmova  %rdi,%rax
mov    %rax,0xe0(%rsp)
mov    0xd8(%rsp),%rax
lea    (%rax,%r8,1),%r13
cmp    %rbp,%r13
mov    %rbp,%rax
cmova  %r13,%rax
mov    %rbp,%rcx
cmova  %r13,%rcx
cmovbe %rbp,%r13
ucomisd 0xc55(%rip),%xmm0
jne    13c3 <qrsolv+0x293>
jnp    1460 <qrsolv+0x330>
mov    %rcx,0x98(%rsp)
mov    %rax,0xd0(%rsp)
mov    %r8,0x18(%rsp)
mov    %rdi,0x20(%rsp)
mov    %r11,0x38(%rsp)
mov    %r9,0x40(%rsp)
movsd  %xmm8,0x48(%rsp)
mov    %r12d,%eax
imul   0x14(%rsp),%eax
movslq %eax,%r15
add    %r12,%r15
mov    0x8(%rsp),%rbp
movsd  0x0(%rbp,%r15,8),%xmm0
call   1040 <fabs@plt>
mov    %rax,%r14
mov    0x128(%rsp),%rbx
movsd  (%rbx,%r12,8),%xmm0
call   1040 <fabs@plt>
cmp    %rax,%r14
jge    1480 <qrsolv+0x350>
movsd  0x0(%rbp,%r15,8),%xmm9
divsd  (%rbx,%r12,8),%xmm9
movapd %xmm9,%xmm0
mulsd  %xmm9,%xmm0
movsd  0xbb9(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xbb9(%rip),%xmm0
jb     14b7 <qrsolv+0x387>
sqrtsd %xmm0,%xmm0
jmp    14d8 <qrsolv+0x3a8>
nop
add    $0x1,%r12
mov    0x128(%rsp),%rbx
jmp    167b <qrsolv+0x54b>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%r12,8),%xmm10
divsd  0x0(%rbp,%r15,8),%xmm10
movapd %xmm10,%xmm0
mulsd  %xmm10,%xmm0
movsd  0xb61(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xb61(%rip),%xmm0
jb     150d <qrsolv+0x3dd>
sqrtsd %xmm0,%xmm0
jmp    152e <qrsolv+0x3fe>
movapd %xmm9,0xa0(%rsp)
call   1050 <sqrt@plt>
movapd 0xa0(%rsp),%xmm9
mov    0x128(%rsp),%rbx
mov    0x30(%rsp),%rbp
movsd  0x48(%rsp),%xmm8
mov    0x40(%rsp),%r9
mov    0x38(%rsp),%r11
mov    0x20(%rsp),%rdi
mov    0x18(%rsp),%r8
movsd  0xb07(%rip),%xmm10
divsd  %xmm0,%xmm10
mulsd  %xmm10,%xmm9
jmp    1561 <qrsolv+0x431>
movapd %xmm10,0xa0(%rsp)
call   1050 <sqrt@plt>
movapd 0xa0(%rsp),%xmm10
mov    0x128(%rsp),%rbx
mov    0x30(%rsp),%rbp
movsd  0x48(%rsp),%xmm8
mov    0x40(%rsp),%r9
mov    0x38(%rsp),%r11
mov    0x20(%rsp),%rdi
mov    0x18(%rsp),%r8
movsd  0xab1(%rip),%xmm9
divsd  %xmm0,%xmm9
mulsd  %xmm9,%xmm10
mov    0x130(%rsp),%rax
movsd  (%rax,%r12,8),%xmm0
movapd %xmm8,%xmm1
mulsd  %xmm10,%xmm1
movapd %xmm9,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm9,%xmm8
mulsd  %xmm10,%xmm0
subsd  %xmm0,%xmm8
movsd  %xmm2,(%rax,%r12,8)
movsd  (%rbx,%r12,8),%xmm0
mulsd  %xmm10,%xmm0
mov    0x8(%rsp),%rax
movsd  (%rax,%r15,8),%xmm1
mulsd  %xmm9,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%r15,8)
lea    0x1(%r12),%r10
cmp    %rbp,%r10
jge    1670 <qrsolv+0x540>
mov    0xc8(%rsp),%rcx
sub    %r8,%rcx
mov    0xd0(%rsp),%rax
add    %rcx,%rax
mov    0x90(%rsp),%r15
imul   %r15,%r12
movapd %xmm10,%xmm0
xorpd  0xa23(%rip),%xmm0
mov    0x28(%rsp),%r14
mov    %r14,%rcx
cmp    $0x2,%rax
jae    16ad <qrsolv+0x57d>
mov    0x8(%rsp),%rax
lea    (%rax,%r12,8),%rdx
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%rcx,8),%xmm1
mov    %ecx,%esi
movsd  (%rbx,%rsi,8),%xmm2
movapd %xmm10,%xmm3
mulsd  %xmm2,%xmm3
movapd %xmm9,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm9,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
movsd  %xmm4,(%rdx,%rcx,8)
add    $0x1,%rcx
cmp    %rbp,%rcx
jb     1620 <qrsolv+0x4f0>
mov    %r10,%r12
jmp    167b <qrsolv+0x54b>
nopw   0x0(%rax,%rax,1)
mov    %r10,%r12
mov    0x90(%rsp),%r15
mov    0x28(%rsp),%r14
cmp    %rbp,%r12
je     1260 <qrsolv+0x130>
add    $0x1,%r14
mov    %r14,0x28(%rsp)
movsd  (%rbx,%r12,8),%xmm0
add    $0x1,%r8
add    $0x1,%rdi
add    $0xffffffffffffffff,%r11
add    $0x8,%r9
jmp    1380 <qrsolv+0x250>
mov    0x50(%rsp),%rcx
lea    (%rcx,%r8,1),%rdx
mov    $0xfffffffffffffffe,%rcx
sub    %rdx,%rcx
add    %rcx,%r13
mov    %r13d,%esi
mov    %r14,%rcx
add    %r14d,%esi
jb     160f <qrsolv+0x4df>
shr    $0x20,%r13
mov    %r14,%rcx
jne    160f <qrsolv+0x4df>
neg    %rdx
mov    %r14d,%ecx
mov    0x98(%rsp),%rsi
add    %rsi,%rdx
add    $0xfffffffffffffffe,%rdx
add    %rcx,%rdx
mov    0xc0(%rsp),%rsi
lea    (%rsi,%rdx,8),%rdx
lea    (%r14,%r12,1),%rsi
mov    %rax,%r13
mov    0x8(%rsp),%rax
lea    (%rax,%rsi,8),%rsi
mov    %r13,%rax
cmp    %rdx,%rsi
jae    1746 <qrsolv+0x616>
lea    (%rbx,%rcx,8),%rdx
mov    0x98(%rsp),%rsi
add    %r12,%rsi
mov    0xb8(%rsp),%rcx
lea    (%rcx,%rsi,8),%rsi
add    $0xfffffffffffffff8,%rsi
mov    %r14,%rcx
cmp    %rsi,%rdx
jb     160f <qrsolv+0x4df>
mov    %rax,%r8
mov    0xe0(%rsp),%rax
add    %r11,%rax
and    $0xfffffffffffffffe,%rax
and    $0xfffffffffffffffe,%r8
lea    (%r14,%r8,1),%rcx
movapd %xmm10,%xmm1
unpcklpd %xmm10,%xmm1
movapd %xmm9,%xmm2
unpcklpd %xmm9,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
lea    (%r9,%r12,8),%rsi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%r14,%rdi,1),%edx
movupd (%rsi,%rdi,8),%xmm4
movupd (%rbx,%rdx,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,(%rbx,%rdx,8)
movupd %xmm7,(%rsi,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rax
jne    1790 <qrsolv+0x660>
cmp    %r8,%r13
mov    0x20(%rsp),%rdi
mov    0x18(%rsp),%r8
jne    160f <qrsolv+0x4df>
mov    %r10,%r12
cmp    %rbp,%r12
jne    1689 <qrsolv+0x559>
jmp    1260 <qrsolv+0x130>
mov    %rbp,%r14
mov    0x10(%rsp),%ebp
test   %ebp,%ebp
jle    1a6f <qrsolv+0x93f>
mov    %r14,%r8
add    $0xffffffffffffffff,%r8
je     1886 <qrsolv+0x756>
mov    %r14,%rcx
and    $0xfffffffe,%ecx
xor    %esi,%esi
xorpd  %xmm0,%xmm0
mov    %ebp,%edx
jmp    1829 <qrsolv+0x6f9>
nopl   (%rax)
add    $0x1,%rsi
cmp    %rcx,%rsi
je     1880 <qrsolv+0x750>
mov    %rsi,%rax
cmp    %ebp,%edx
mov    %eax,%edi
cmovne %edx,%edi
ucomisd (%rbx,%rsi,8),%xmm0
cmovne %edx,%edi
cmovp  %edx,%edi
cmp    %ebp,%edi
jge    1852 <qrsolv+0x722>
mov    0x130(%rsp),%rdx
movq   $0x0,(%rdx,%rax,8)
lea    0x1(%rax),%rsi
cmp    %ebp,%edi
mov    %esi,%edx
cmovne %edi,%edx
ucomisd 0x8(%rbx,%rax,8),%xmm0
cmovne %edi,%edx
cmovp  %edi,%edx
cmp    %ebp,%edx
jge    1820 <qrsolv+0x6f0>
mov    0x130(%rsp),%rdi
movq   $0x0,0x8(%rdi,%rax,8)
jmp    1820 <qrsolv+0x6f0>
add    $0x2,%rax
jmp    188a <qrsolv+0x75a>
xor    %eax,%eax
mov    %ebp,%edx
mov    %edx,%r11d
test   $0x1,%r14b
je     18c2 <qrsolv+0x792>
cmp    %ebp,%edx
mov    %eax,%r11d
cmovne %edx,%r11d
xorpd  %xmm0,%xmm0
ucomisd (%rbx,%rax,8),%xmm0
cmovne %edx,%r11d
cmovp  %edx,%r11d
cmp    %ebp,%r11d
jge    18c2 <qrsolv+0x792>
mov    0x130(%rsp),%rcx
movq   $0x0,(%rcx,%rax,8)
test   %r11d,%r11d
jle    19f3 <qrsolv+0x8c3>
mov    %r11d,%edx
lea    0x1(%r11),%r9d
lea    -0x1(%rdx),%rax
imul   %r15,%rax
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
shl    $0x3,%r15
neg    %r15
mov    %r15,%r13
mov    $0x1,%ebp
xor    %r15d,%r15d
jmp    1933 <qrsolv+0x803>
nopw   0x0(%rax,%rax,1)
mov    0x130(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
subsd  %xmm0,%xmm1
divsd  (%rbx,%rax,8),%xmm1
movsd  %xmm1,(%rcx,%rax,8)
add    $0x1,%rbp
add    $0xffffffff,%r11d
add    $0x1,%r15d
add    %r13,%rdi
cmp    %r9,%rbp
je     19f3 <qrsolv+0x8c3>
mov    %rdx,%rax
sub    %rbp,%rax
lea    0x1(%rax),%rcx
xorpd  %xmm0,%xmm0
cmp    %rdx,%rcx
jge    1900 <qrsolv+0x7d0>
lea    -0x1(%r15),%r10d
movslq %r11d,%rsi
test   $0x3,%r15b
je     1981 <qrsolv+0x851>
mov    %r15d,%ecx
and    $0x3,%ecx
xorpd  %xmm0,%xmm0
nopl   (%rax)
movsd  (%rdi,%rsi,8),%xmm1
mov    0x130(%rsp),%rbx
mulsd  (%rbx,%rsi,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rsi
add    $0xffffffff,%ecx
jne    1960 <qrsolv+0x830>
jmp    1985 <qrsolv+0x855>
xorpd  %xmm0,%xmm0
cmp    $0x3,%r10d
mov    0x128(%rsp),%rbx
jb     1900 <qrsolv+0x7d0>
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rsi,8),%xmm1
movsd  0x8(%rdi,%rsi,8),%xmm2
mov    0x130(%rsp),%rcx
mulsd  (%rcx,%rsi,8),%xmm1
mulsd  0x8(%rcx,%rsi,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x10(%rdi,%rsi,8),%xmm3
mulsd  0x10(%rcx,%rsi,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x18(%rdi,%rsi,8),%xmm0
mulsd  0x18(%rcx,%rsi,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rsi
cmp    %esi,%edx
jne    19a0 <qrsolv+0x870>
jmp    1900 <qrsolv+0x7d0>
cmpl   $0x0,0x10(%rsp)
mov    0x58(%rsp),%rsi
mov    0x120(%rsp),%rdi
jle    1a6f <qrsolv+0x93f>
test   %r8,%r8
je     1a81 <qrsolv+0x951>
mov    %r14d,%ecx
and    $0xfffffffe,%ecx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movslq (%rsi,%rax,4),%rdx
mov    0x130(%rsp),%rbp
movsd  0x0(%rbp,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rdi,%rdx,8)
movslq 0x4(%rsi,%rax,4),%rdx
movsd  0x8(%rbp,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rdi,%rdx,8)
add    $0x2,%rax
cmp    %rax,%rcx
jne    1a20 <qrsolv+0x8f0>
test   $0x1,%r14b
je     1a6f <qrsolv+0x93f>
movslq (%rsi,%rax,4),%rcx
mov    0x130(%rsp),%rdx
movsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rdi,%rcx,8)
add    $0xe8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %eax,%eax
test   $0x1,%r14b
jne    1a58 <qrsolv+0x928>
jmp    1a6f <qrsolv+0x93f>

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

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