#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);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
double sqrt(double);
void qrsolv(double *r, int ldr, int N, int *ipvt, double *diag, double *qtb,
            double *x, double *sdiag) {
  int i, j, jp1, k, kp1, l, nsing, t;
  double cos, cotan, p5, p25, qtbpj, sin, sum, tan, temp, zero;
  double *wa;
  wa = (double *)malloc(sizeof(double) * N);
  p5 = 5.0e-1;
  p25 = 2.5e-1;
  zero = 0.0;
  for (j = 0; j < N; ++j) {
    for (i = j; i < N; ++i) {
      r[i * N + j] = r[j * N + i];
    }
    x[j] = r[j * N + j];
    wa[j] = qtb[j];
  }
  for (j = 0; j < N; ++j) {
    l = ipvt[j];
    if (diag[l] != zero) {
      for (k = j; k < N; ++k) {
        sdiag[k] = zero;
      }
      sdiag[j] = diag[l];
      qtbpj = zero;
      for (k = j; k < N; ++k) {
        t = k * N;
        if (sdiag[k] != zero) {
          if (fabs(r[t + k]) < fabs(sdiag[k])) {
            cotan = r[t + k] / sdiag[k];
            sin = p5 / sqrt(p25 + p25 * cotan * cotan);
            cos = sin * cotan;
          } else {
            tan = sdiag[k] / r[t + k];
            cos = p5 / sqrt(p25 + p25 * tan * tan);
            sin = cos * tan;
          }
          r[t + k] = cos * r[t + k] + sin * sdiag[k];
          temp = cos * wa[k] + sin * qtbpj;
          qtbpj = -sin * wa[k] + cos * qtbpj;
          wa[k] = temp;
          kp1 = k + 1;
          if (N >= kp1 + 1) {
            for (i = kp1; i < N; ++i) {
              temp = cos * r[i * N + k] + sin * sdiag[i];
              sdiag[i] = -sin * r[i * N + k] + cos * sdiag[i];
              r[i * N + k] = temp;
            }
          }
        }
      }
    }
    sdiag[j] = r[j * N + j];
    r[j * N + j] = x[j];
  }
  nsing = N;
  for (j = 1; j <= N; ++j) {
    if (sdiag[j - 1] == zero && nsing == N) {
      nsing = j - 1;
    }
    if (nsing < N) {
      wa[j - 1] = zero;
    }
  }
  if (nsing >= 1) {
    for (k = 1; k <= nsing; ++k) {
      j = nsing - k + 1;
      sum = zero;
      jp1 = j + 1;
      if (nsing >= jp1) {
        for (i = jp1; i <= nsing; ++i) {
          sum = sum + r[(i - 1) * N + j - 1] * wa[i - 1];
        }
      }
      wa[j - 1] = (wa[j - 1] - sum) / sdiag[j - 1];
    }
  }
  for (j = 0; j < N; ++j) {
    l = ipvt[j];
    x[l] = wa[j];
  }
  free(wa);
}


x86 assembly (clang, -O0):
<qrsolv>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movslq -0x10(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0xa0(%rbp)
movsd  0xe81(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xe6c(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1251 <qrsolv+0x111>
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1203 <qrsolv+0xc3>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11bd <qrsolv+0x7d>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xa0(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11ab <qrsolv+0x6b>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    16a5 <qrsolv+0x565>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x98(%rbp),%xmm0
ff
jne    1298 <qrsolv+0x158>
jp     1298 <qrsolv+0x158>
jmp    1651 <qrsolv+0x511>
mov    -0x30(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12cd <qrsolv+0x18d>
movsd  -0x98(%rbp),%xmm0
ff
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    129e <qrsolv+0x15e>
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    164c <qrsolv+0x50c>
mov    -0x38(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x98(%rbp),%xmm0
ff
jne    1336 <qrsolv+0x1f6>
jp     1336 <qrsolv+0x1f6>
jmp    1639 <qrsolv+0x4f9>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <fabs@plt>
mov    %rax,-0xa8(%rbp)
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <fabs@plt>
mov    %rax,%rcx
mov    -0xa8(%rbp),%rax
cmp    %rcx,%rax
jge    13f3 <qrsolv+0x2b3>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x68(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1060 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xb0(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1474 <qrsolv+0x334>
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x68(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1060 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x78(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x78(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0xa0(%rbp),%rax
movslq -0x38(%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    -0xa0(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
mov    -0xa0(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jl     1634 <qrsolv+0x4f4>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    162f <qrsolv+0x4ef>
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x78(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x2c(%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    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1561 <qrsolv+0x421>
jmp    1634 <qrsolv+0x4f4>
jmp    1639 <qrsolv+0x4f9>
jmp    163e <qrsolv+0x4fe>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12fa <qrsolv+0x1ba>
jmp    1651 <qrsolv+0x511>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1258 <qrsolv+0x118>
mov    -0x10(%rbp),%eax
mov    %eax,-0x44(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1735 <qrsolv+0x5f5>
mov    0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x98(%rbp),%xmm0
ff
jne    16f9 <qrsolv+0x5b9>
jp     16f9 <qrsolv+0x5b9>
mov    -0x44(%rbp),%eax
cmp    -0x10(%rbp),%eax
jne    16f9 <qrsolv+0x5b9>
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1722 <qrsolv+0x5e2>
movsd  -0x98(%rbp),%xmm0
ff
mov    -0xa0(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1727 <qrsolv+0x5e7>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    16b2 <qrsolv+0x572>
cmpl   $0x1,-0x44(%rbp)
jl     183c <qrsolv+0x6fc>
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     1837 <qrsolv+0x6f7>
mov    -0x44(%rbp),%eax
sub    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
movsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x80(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x34(%rbp),%eax
jl     17e8 <qrsolv+0x6a8>
mov    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     17e3 <qrsolv+0x6a3>
movsd  -0x80(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xa0(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1786 <qrsolv+0x646>
jmp    17e8 <qrsolv+0x6a8>
mov    -0xa0(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x80(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0xa0(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1746 <qrsolv+0x606>
jmp    183c <qrsolv+0x6fc>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1888 <qrsolv+0x748>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0xa0(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1843 <qrsolv+0x703>
mov    -0xa0(%rbp),%rdi
call   1030 <free@plt>
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qrsolv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9,%r14
mov    %r8,0x78(%rsp)
mov    %rcx,%r12
mov    %edx,%ebx
mov    %rdi,%r13
mov    0xc8(%rsp),%rbp
mov    0xc0(%rsp),%r15
lea    0x0(,%rbx,8),%edi
call   1060 <malloc@plt>
test   %ebx,%ebx
jle    1848 <qrsolv+0x6f8>
mov    %r13,%r10
mov    %rbx,(%rsp)
mov    %ebx,%r13d
lea    0x0(,%r13,8),%rdx
lea    0x8(,%r13,8),%r8
xor    %r9d,%r9d
mov    %r13,%rsi
mov    %r10,%rcx
nopl   0x0(%rax,%rax,1)
xor    %ebx,%ebx
mov    %rcx,%rdi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx,%rbx,8),%xmm0
movsd  %xmm0,(%rdi)
add    %rdx,%rdi
add    $0x1,%rbx
cmp    %rbx,%rsi
jne    11d0 <qrsolv+0x80>
mov    %r9d,%edi
imul   (%rsp),%edi
movslq %edi,%rdi
add    %r9,%rdi
movsd  (%r10,%rdi,8),%xmm0
movsd  %xmm0,(%r15,%r9,8)
movsd  (%r14,%r9,8),%xmm0
movsd  %xmm0,(%rax,%r9,8)
add    $0x1,%r9
add    %r8,%rcx
add    $0xffffffffffffffff,%rsi
cmp    %r13,%r9
jne    11c0 <qrsolv+0x70>
mov    (%rsp),%r14
test   %r14d,%r14d
jle    1848 <qrsolv+0x6f8>
mov    %r12,0x58(%rsp)
mov    %rax,0x8(%rsp)
mov    %r15,0x60(%rsp)
lea    0x0(,%r13,8),%r15
mov    $0x1,%eax
mov    %rax,0x70(%rsp)
xor    %edx,%edx
mov    %r13,0x68(%rsp)
mov    %r10,0x28(%rsp)
mov    %r13,0x30(%rsp)
jmp    12a9 <qrsolv+0x159>
nopl   0x0(%rax,%rax,1)
mov    0x38(%rsp),%rdx
mov    %rdx,%rax
imul   %r13,%rax
add    %rdx,%rax
movsd  (%r10,%rax,8),%xmm0
mov    0x80(%rsp),%rcx
movsd  %xmm0,(%rcx)
mov    0x60(%rsp),%rcx
movsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%r10,%rax,8)
add    $0x1,%rdx
addq   $0x1,0x70(%rsp)
add    %r13,0x68(%rsp)
cmp    %r13,%rdx
je     1581 <qrsolv+0x431>
lea    0x0(,%rdx,8),%rax
add    %rbp,%rax
mov    %rax,0x80(%rsp)
mov    0x58(%rsp),%rax
mov    %rdx,0x38(%rsp)
movslq (%rax,%rdx,4),%r12
mov    0x78(%rsp),%rax
movsd  (%rax,%r12,8),%xmm0
ucomisd 0xd33(%rip),%xmm0
jne    12e5 <qrsolv+0x195>
jnp    1260 <qrsolv+0x110>
mov    0x38(%rsp),%rax
not    %eax
add    %r14d,%eax
lea    0x8(,%rax,8),%rdx
mov    0x80(%rsp),%rbx
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
mov    0x28(%rsp),%r10
mov    0x78(%rsp),%rax
movsd  (%rax,%r12,8),%xmm0
movsd  %xmm0,(%rbx)
xorpd  %xmm5,%xmm5
mov    0x68(%rsp),%r8
mov    0x38(%rsp),%r12
mov    0x70(%rsp),%rdi
ucomisd 0xcd8(%rip),%xmm0
jne    1340 <qrsolv+0x1f0>
jnp    13c0 <qrsolv+0x270>
mov    %rdi,0x10(%rsp)
mov    %r8,0x18(%rsp)
movsd  %xmm5,0x20(%rsp)
mov    %r12d,%eax
imul   %r14d,%eax
movslq %eax,%rbx
add    %r12,%rbx
movsd  (%r10,%rbx,8),%xmm0
mov    %r10,%r13
call   1050 <fabs@plt>
mov    %rax,%r14
movsd  0x0(%rbp,%r12,8),%xmm0
call   1050 <fabs@plt>
cmp    %rax,%r14
jge    13e0 <qrsolv+0x290>
movsd  0x0(%r13,%rbx,8),%xmm6
divsd  0x0(%rbp,%r12,8),%xmm6
movapd %xmm6,%xmm0
movsd  0xc67(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm6,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xc63(%rip),%xmm0
jb     1419 <qrsolv+0x2c9>
mov    %r13,%r10
sqrtsd %xmm0,%xmm0
jmp    142f <qrsolv+0x2df>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r12
cmp    %r13,%r12
jne    156e <qrsolv+0x41e>
jmp    1260 <qrsolv+0x110>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%rbp,%r12,8),%xmm7
divsd  0x0(%r13,%rbx,8),%xmm7
movapd %xmm7,%xmm0
movsd  0xc06(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm7,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xc02(%rip),%xmm0
jb     145a <qrsolv+0x30a>
mov    %r13,%r10
sqrtsd %xmm0,%xmm0
jmp    1470 <qrsolv+0x320>
movsd  %xmm6,0x40(%rsp)
call   1070 <sqrt@plt>
movsd  0x40(%rsp),%xmm6
mov    0x28(%rsp),%r10
mov    (%rsp),%r14
mov    0x30(%rsp),%r13
movsd  0x20(%rsp),%xmm5
mov    0x18(%rsp),%r8
mov    0x10(%rsp),%rdi
movsd  0xbb8(%rip),%xmm7
divsd  %xmm0,%xmm7
mulsd  %xmm7,%xmm6
jmp    1499 <qrsolv+0x349>
movapd %xmm7,0x40(%rsp)
call   1070 <sqrt@plt>
movapd 0x40(%rsp),%xmm7
mov    0x28(%rsp),%r10
mov    (%rsp),%r14
mov    0x30(%rsp),%r13
movsd  0x20(%rsp),%xmm5
mov    0x18(%rsp),%r8
mov    0x10(%rsp),%rdi
movsd  0xb77(%rip),%xmm6
divsd  %xmm0,%xmm6
mulsd  %xmm6,%xmm7
movsd  0x0(%rbp,%r12,8),%xmm0
mulsd  %xmm7,%xmm0
movsd  (%r10,%rbx,8),%xmm1
mulsd  %xmm6,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r10,%rbx,8)
mov    0x8(%rsp),%rax
movsd  (%rax,%r12,8),%xmm0
movapd %xmm5,%xmm1
mulsd  %xmm7,%xmm1
movapd %xmm6,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm6,%xmm5
mulsd  %xmm7,%xmm0
subsd  %xmm0,%xmm5
movsd  %xmm2,(%rax,%r12,8)
lea    0x2(%r12),%rcx
lea    0x1(%r12),%rax
cmp    %r13,%rcx
jg     1562 <qrsolv+0x412>
cmp    %r13,%rax
jge    1562 <qrsolv+0x412>
movapd %xmm7,%xmm0
xorpd  0xb17(%rip),%xmm0
add    %r8,%r12
lea    (%r10,%r12,8),%rcx
mov    %rdi,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rcx),%xmm1
mov    %edx,%esi
movsd  0x0(%rbp,%rsi,8),%xmm2
movapd %xmm7,%xmm3
mulsd  %xmm2,%xmm3
movapd %xmm6,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm6,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x0(%rbp,%rsi,8)
movsd  %xmm4,(%rcx)
add    $0x1,%rdx
add    %r15,%rcx
cmp    %r13,%rdx
jb     1520 <qrsolv+0x3d0>
mov    %rax,%r12
cmp    %r13,%r12
je     1260 <qrsolv+0x110>
add    $0x1,%rdi
movsd  0x0(%rbp,%r12,8),%xmm0
add    %r13,%r8
jmp    1330 <qrsolv+0x1e0>
test   %r14d,%r14d
jle    1727 <qrsolv+0x5d7>
mov    (%rsp),%rbx
cmp    $0x1,%ebx
jne    1757 <qrsolv+0x607>
xor    %eax,%eax
mov    %ebx,%ecx
mov    0x8(%rsp),%rdi
mov    %ecx,%edx
test   $0x1,%r13b
je     15cb <qrsolv+0x47b>
cmp    %ebx,%ecx
mov    %eax,%edx
cmovne %ecx,%edx
xorpd  %xmm0,%xmm0
ucomisd 0x0(%rbp,%rax,8),%xmm0
cmovne %ecx,%edx
cmovp  %ecx,%edx
cmp    %ebx,%edx
jge    15cb <qrsolv+0x47b>
movq   $0x0,(%rdi,%rax,8)
test   %edx,%edx
jle    1727 <qrsolv+0x5d7>
lea    0x1(%rdx),%ecx
mov    %edx,%eax
mov    %rax,0x40(%rsp)
lea    -0x1(%rdx),%eax
mov    %eax,0x10(%rsp)
mov    (%rsp),%rax
lea    (%rax,%rax,1),%edi
mov    $0x1,%r12d
xor    %r15d,%r15d
xorpd  %xmm0,%xmm0
mov    %rcx,0x18(%rsp)
mov    %ecx,%r14d
jmp    1650 <qrsolv+0x500>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x8(%rsp),%rax
movsd  (%rax,%r11,8),%xmm2
subsd  %xmm1,%xmm2
divsd  0x0(%rbp,%r11,8),%xmm2
movsd  %xmm2,(%rax,%r11,8)
add    $0x1,%r12
add    $0xffffffff,%r14d
add    $0x1,%r15d
addl   $0xffffffff,0x10(%rsp)
mov    0x20(%rsp),%rdx
add    $0xffffffff,%edx
cmp    0x18(%rsp),%r12
je     1727 <qrsolv+0x5d7>
mov    %rdx,0x20(%rsp)
mov    0x40(%rsp),%rax
mov    %rax,%r11
sub    %r12,%r11
lea    0x2(%r11),%rcx
xorpd  %xmm1,%xmm1
cmp    %rax,%rcx
jg     1610 <qrsolv+0x4c0>
movslq %r14d,%rbx
xorpd  %xmm1,%xmm1
test   $0x1,%r15b
je     16a1 <qrsolv+0x551>
lea    -0x1(%rbx),%rcx
imul   (%rsp),%ecx
add    %r11d,%ecx
movslq %ecx,%rcx
movsd  (%r10,%rcx,8),%xmm1
mov    0x8(%rsp),%rax
mulsd  -0x8(%rax,%rbx,8),%xmm1
addsd  %xmm0,%xmm1
add    $0x1,%rbx
cmp    $0x1,%r15d
je     1610 <qrsolv+0x4c0>
mov    0x18(%rsp),%rsi
sub    %ebx,%esi
mov    (%rsp),%r8
mov    %r8d,%edx
imul   %ebx,%edx
add    0x10(%rsp),%edx
mov    0x8(%rsp),%rax
lea    (%rax,%rbx,8),%rcx
add    $0xffffffff,%ebx
imul   %r8d,%ebx
add    0x20(%rsp),%ebx
mov    $0xfffffffe,%r8d
xor    %r9d,%r9d
nopl   (%rax)
lea    (%r8,%rbx,1),%r10d
add    $0x1,%r10d
movslq %r10d,%rax
mov    0x28(%rsp),%r10
movsd  (%r10,%rax,8),%xmm2
mulsd  -0x8(%rcx,%r9,8),%xmm2
addsd  %xmm1,%xmm2
movslq %edx,%rdx
movsd  (%r10,%rdx,8),%xmm1
mulsd  (%rcx,%r9,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x2,%r9
add    %edi,%edx
add    %edi,%r8d
cmp    %r9d,%esi
jne    16e0 <qrsolv+0x590>
jmp    1610 <qrsolv+0x4c0>
cmpl   $0x0,(%rsp)
mov    0x60(%rsp),%rbp
mov    0x8(%rsp),%rax
mov    0x58(%rsp),%rsi
jle    1848 <qrsolv+0x6f8>
lea    -0x1(%r13),%rcx
mov    %r13d,%ebx
and    $0x3,%ebx
cmp    $0x3,%rcx
jae    17bf <qrsolv+0x66f>
xor    %ecx,%ecx
jmp    181b <qrsolv+0x6cb>
mov    %r13d,%edx
and    $0xfffffffe,%edx
xor    %eax,%eax
xorpd  %xmm0,%xmm0
mov    %ebx,%ecx
mov    0x8(%rsp),%rdi
jmp    177d <qrsolv+0x62d>
nopl   0x0(%rax)
add    $0x2,%rax
cmp    %rax,%rdx
je     15a0 <qrsolv+0x450>
cmp    %ebx,%ecx
mov    %eax,%esi
cmovne %ecx,%esi
ucomisd 0x0(%rbp,%rax,8),%xmm0
cmovne %ecx,%esi
cmovp  %ecx,%esi
cmp    %ebx,%esi
jge    179c <qrsolv+0x64c>
movq   $0x0,(%rdi,%rax,8)
lea    0x1(%rax),%ecx
cmp    %ebx,%esi
cmovne %esi,%ecx
ucomisd 0x8(%rbp,%rax,8),%xmm0
cmovne %esi,%ecx
cmovp  %esi,%ecx
cmp    %ebx,%ecx
jge    1770 <qrsolv+0x620>
movq   $0x0,0x8(%rdi,%rax,8)
jmp    1770 <qrsolv+0x620>
and    $0xfffffffc,%r13d
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movslq (%rsi,%rcx,4),%rdx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
movslq 0x4(%rsi,%rcx,4),%rdx
movsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
movslq 0x8(%rsi,%rcx,4),%rdx
movsd  0x10(%rax,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
movslq 0xc(%rsi,%rcx,4),%rdx
movsd  0x18(%rax,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
add    $0x4,%rcx
cmp    %rcx,%r13
jne    17d0 <qrsolv+0x680>
test   %rbx,%rbx
je     1848 <qrsolv+0x6f8>
lea    (%rax,%rcx,8),%rdx
lea    (%rsi,%rcx,4),%rcx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movslq (%rcx,%rsi,4),%rdi
movsd  (%rdx,%rsi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdi,8)
add    $0x1,%rsi
cmp    %rsi,%rbx
jne    1830 <qrsolv+0x6e0>
mov    %rax,%rdi
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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