#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 const sqrt(double const);
int qrsolv_(int n, double *r, int ldr, int *ipvt, double *diag, double *qtb,
            double *x, double *sdiag, double *wa) {
  const double p5 = .5;
  const double p25 = .25;
  double tanx, cosx, sinx, cotan;
  double sum, temp, qtbpj;
  int i, j, k, l;
  int nsing;
  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];
    if (diag[l] == 0.0) {
      goto store;
    }
    for (k = j; k < n; ++k) {
      sdiag[k] = 0.0;
    }
    sdiag[j] = diag[l];
    qtbpj = 0.0;
    for (k = j; k < n; ++k) {
      if (sdiag[k] == 0.0) {
        continue;
      }
      if (fabs(r[k + k * ldr]) < fabs(sdiag[k])) {
        cotan = r[k + k * ldr] / sdiag[k];
        sinx = p5 / sqrt(p25 + p25 * (cotan * cotan));
        cosx = sinx * cotan;
      } else {
        tanx = sdiag[k] / r[k + k * ldr];
        cosx = p5 / sqrt(p25 + p25 * (tanx * tanx));
        sinx = cosx * tanx;
      }
      r[k + k * ldr] = cosx * r[k + k * ldr] + sinx * sdiag[k];
      temp = cosx * wa[k] + sinx * qtbpj;
      qtbpj = -sinx * wa[k] + cosx * qtbpj;
      wa[k] = temp;
      for (i = k + 1; i < n; ++i) {
        temp = cosx * r[i + k * ldr] + sinx * sdiag[i];
        sdiag[i] = -sinx * r[i + k * ldr] + cosx * sdiag[i];
        r[i + k * ldr] = temp;
      }
    }
  store:
    sdiag[j] = r[j + j * ldr];
    r[j + j * ldr] = x[j];
  }
  nsing = n;
  for (j = 0; j < n; ++j) {
    if (sdiag[j] == 0.0 && nsing == n) {
      nsing = j;
    }
    if (nsing < n) {
      wa[j] = 0.0;
    }
  }
  if (nsing > 0) {
    for (k = 0; k < nsing; ++k) {
      j = nsing - k - 1;
      sum = 0.0;
      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];
    x[l] = wa[j];
  }
  return 0;
}


x86 assembly (clang, -O0):
<qrsolv_>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  0xeb3(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe9e(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x80(%rbp)
mov    -0x80(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1217 <qrsolv_+0xf7>
mov    -0x80(%rbp),%eax
mov    %eax,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11ca <qrsolv_+0xaa>
mov    -0x10(%rbp),%rax
mov    -0x80(%rbp),%ecx
mov    -0x7c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    -0x80(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    1180 <qrsolv_+0x60>
mov    -0x10(%rbp),%rax
mov    -0x80(%rbp),%ecx
mov    -0x80(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    116e <qrsolv_+0x4e>
movl   $0x0,-0x80(%rbp)
mov    -0x80(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1683 <qrsolv_+0x563>
mov    -0x20(%rbp),%rax
movslq -0x80(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x88(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x88(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1263 <qrsolv_+0x143>
jp     1263 <qrsolv_+0x143>
jmp    162b <qrsolv_+0x50b>
mov    -0x80(%rbp),%eax
mov    %eax,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12a2 <qrsolv_+0x182>
mov    0x18(%rbp),%rax
movslq -0x84(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    126c <qrsolv_+0x14c>
mov    -0x28(%rbp),%rax
movslq -0x88(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x80(%rbp),%eax
mov    %eax,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1626 <qrsolv_+0x506>
mov    0x18(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1307 <qrsolv_+0x1e7>
jp     1307 <qrsolv_+0x1e7>
jmp    1612 <qrsolv_+0x4f2>
mov    -0x10(%rbp),%rax
mov    -0x84(%rbp),%ecx
mov    -0x84(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
mov    %rax,-0x98(%rbp)
mov    0x18(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
mov    %rax,%rcx
mov    -0x98(%rbp),%rax
cmp    %rcx,%rax
jge    13d3 <qrsolv_+0x2b3>
mov    -0x10(%rbp),%rax
mov    -0x84(%rbp),%ecx
mov    -0x84(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x84(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  0xc62(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xc52(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1448 <qrsolv_+0x328>
mov    0x18(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x84(%rbp),%ecx
mov    -0x84(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  0xbe8(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xbd8(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x84(%rbp),%ecx
mov    -0x84(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x84(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x84(%rbp),%ecx
mov    -0x84(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x58(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x20(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x70(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    160d <qrsolv_+0x4ed>
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    -0x84(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mov    0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x58(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    -0x84(%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 -0x7c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    -0x84(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    1536 <qrsolv_+0x416>
jmp    1612 <qrsolv_+0x4f2>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    12d0 <qrsolv_+0x1b0>
jmp    162b <qrsolv_+0x50b>
mov    -0x10(%rbp),%rax
mov    -0x80(%rbp),%ecx
mov    -0x80(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x80(%rbp),%ecx
mov    -0x80(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    121e <qrsolv_+0xfe>
mov    -0x4(%rbp),%eax
mov    %eax,-0x8c(%rbp)
movl   $0x0,-0x80(%rbp)
mov    -0x80(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1709 <qrsolv_+0x5e9>
mov    0x18(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    16d7 <qrsolv_+0x5b7>
jp     16d7 <qrsolv_+0x5b7>
mov    -0x8c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jne    16d7 <qrsolv_+0x5b7>
mov    -0x80(%rbp),%eax
mov    %eax,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    16f6 <qrsolv_+0x5d6>
mov    0x20(%rbp),%rax
movslq -0x80(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    16fb <qrsolv_+0x5db>
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    1693 <qrsolv_+0x573>
cmpl   $0x0,-0x8c(%rbp)
jle    17ee <qrsolv_+0x6ce>
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x8c(%rbp),%eax
jge    17e9 <qrsolv_+0x6c9>
mov    -0x8c(%rbp),%eax
sub    -0x84(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x80(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    -0x8c(%rbp),%eax
jge    17a9 <qrsolv_+0x689>
mov    -0x10(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    -0x80(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    1755 <qrsolv_+0x635>
mov    0x20(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x68(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x80(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    1720 <qrsolv_+0x600>
jmp    17ee <qrsolv_+0x6ce>
movl   $0x0,-0x80(%rbp)
mov    -0x80(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    183d <qrsolv_+0x71d>
mov    -0x20(%rbp),%rax
movslq -0x80(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x88(%rbp)
mov    0x20(%rbp),%rax
movslq -0x80(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x88(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    17f5 <qrsolv_+0x6d5>
xor    %eax,%eax
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qrsolv_>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xe8,%rsp
mov    %r8,0x98(%rsp)
mov    %rcx,0x50(%rsp)
mov    %edx,0x14(%rsp)
mov    %rsi,0x8(%rsp)
test   %edi,%edi
jle    1a07 <qrsolv_+0x8d7>
mov    0x130(%rsp),%r8
mov    0x128(%rsp),%rax
mov    0x120(%rsp),%rbx
movslq 0x14(%rsp),%r15
mov    %edi,0x4(%rsp)
mov    %edi,%r12d
lea    0x8(,%r15,8),%r10
xor    %ecx,%ecx
mov    %r12,%rdx
mov    0x8(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
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    11b0 <qrsolv_+0x80>
mov    %ecx,%edi
imul   0x14(%rsp),%edi
movslq %edi,%rdi
add    %rcx,%rdi
mov    0x8(%rsp),%rbp
movsd  0x0(%rbp,%rdi,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
movsd  (%r9,%rcx,8),%xmm0
movsd  %xmm0,(%r8,%rcx,8)
add    $0x1,%rcx
add    %r10,%rsi
add    $0xffffffffffffffff,%rdx
cmp    %r12,%rcx
jne    11a0 <qrsolv_+0x70>
cmpl   $0x0,0x4(%rsp)
jle    1a07 <qrsolv_+0x8d7>
mov    0x8(%rsp),%rcx
add    $0x8,%rcx
lea    0x8(%rax),%rdx
mov    %rdx,0xd0(%rsp)
mov    $0x1,%edx
mov    %rdx,0x90(%rsp)
mov    $0x2,%edx
mov    %rdx,0x88(%rsp)
movq   $0xffffffffffffffff,0x80(%rsp)
xor    %r9d,%r9d
xorpd  %xmm8,%xmm8
mov    %rcx,0xc8(%rsp)
mov    %rcx,0x78(%rsp)
mov    %r15,0x30(%rsp)
mov    %r12,0x28(%rsp)
jmp    12d1 <qrsolv_+0x1a1>
cs nopw 0x0(%rax,%rax,1)
mov    %r9,%rcx
imul   %r15,%rcx
add    %r9,%rcx
mov    0x8(%rsp),%rsi
movsd  (%rsi,%rcx,8),%xmm0
mov    0xa0(%rsp),%rdx
movsd  %xmm0,(%rdx)
mov    0x120(%rsp),%rbx
movsd  (%rbx,%r9,8),%xmm0
movsd  %xmm0,(%rsi,%rcx,8)
add    $0x1,%r9
addq   $0x1,0x90(%rsp)
addq   $0x1,0x88(%rsp)
addq   $0xffffffffffffffff,0x80(%rsp)
addq   $0x8,0x78(%rsp)
cmp    %r12,%r9
je     181d <qrsolv_+0x6ed>
lea    (%rax,%r9,8),%rcx
mov    %rcx,0xa0(%rsp)
mov    0x50(%rsp),%rcx
movslq (%rcx,%r9,4),%rbx
mov    0x98(%rsp),%rcx
movsd  (%rcx,%rbx,8),%xmm0
ucomisd %xmm8,%xmm0
jne    1300 <qrsolv_+0x1d0>
jnp    1270 <qrsolv_+0x140>
lea    0x2(%r9),%r14
mov    %r9,%rax
not    %rax
mov    %rax,0xd8(%rsp)
mov    %r9d,%eax
not    %eax
add    0x4(%rsp),%eax
lea    0x8(,%rax,8),%rdx
mov    0xa0(%rsp),%r13
mov    %r13,%rdi
xor    %esi,%esi
mov    %r9,%rbp
call   1030 <memset@plt>
mov    %rbp,%r9
mov    %r14,%rcx
xorpd  %xmm8,%xmm8
mov    0x130(%rsp),%r8
mov    0x98(%rsp),%rax
movsd  (%rax,%rbx,8),%xmm0
movsd  %xmm0,0x0(%r13)
xorpd  %xmm9,%xmm9
mov    0x78(%rsp),%r10
mov    0x80(%rsp),%rdx
mov    0x88(%rsp),%rsi
xor    %r11d,%r11d
mov    %rbp,%rbx
mov    0x90(%rsp),%rax
mov    %rax,0x20(%rsp)
mov    %rbp,0x60(%rsp)
mov    %r14,0x58(%rsp)
jmp    13d8 <qrsolv_+0x2a8>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rbx
mov    0x128(%rsp),%rax
mov    0x20(%rsp),%r14
cmp    %r12,%rbx
je     1270 <qrsolv_+0x140>
add    $0x1,%r14
mov    %r14,0x20(%rsp)
movsd  (%rax,%rbx,8),%xmm0
add    $0x1,%r11
add    $0x1,%rsi
add    $0xffffffffffffffff,%rdx
add    $0x8,%r10
cmp    %r12,%rsi
mov    %r12,%rax
cmova  %rsi,%rax
lea    (%rcx,%r11,1),%r13
cmp    %r12,%r13
mov    %r12,%rbp
cmova  %r13,%rbp
mov    %r12,%rdi
cmova  %r13,%rdi
cmovbe %r12,%r13
ucomisd %xmm8,%xmm0
jne    1404 <qrsolv_+0x2d4>
jnp    13a0 <qrsolv_+0x270>
mov    %rdi,0xa8(%rsp)
mov    %rbp,0x68(%rsp)
mov    %rax,0xe0(%rsp)
mov    %r11,0x18(%rsp)
mov    %rsi,0x70(%rsp)
mov    %rdx,0x48(%rsp)
mov    %r10,0x38(%rsp)
movsd  %xmm9,0x40(%rsp)
mov    %ebx,%eax
imul   0x14(%rsp),%eax
movslq %eax,%rbp
add    %rbx,%rbp
mov    0x8(%rsp),%r12
movsd  (%r12,%rbp,8),%xmm0
call   1040 <fabs@plt>
mov    %rax,%r14
mov    0x128(%rsp),%r15
movsd  (%r15,%rbx,8),%xmm0
call   1040 <fabs@plt>
cmp    %rax,%r14
jge    14a7 <qrsolv_+0x377>
movsd  (%r12,%rbp,8),%xmm11
divsd  (%r15,%rbx,8),%xmm11
movapd %xmm11,%xmm0
mulsd  %xmm11,%xmm0
movsd  0xb76(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm8,%xmm8
ucomisd %xmm8,%xmm0
jb     14e5 <qrsolv_+0x3b5>
mov    %r15,%rax
sqrtsd %xmm0,%xmm0
jmp    150b <qrsolv_+0x3db>
movsd  (%r15,%rbx,8),%xmm10
divsd  (%r12,%rbp,8),%xmm10
movapd %xmm10,%xmm0
mulsd  %xmm10,%xmm0
movsd  0xb3b(%rip),%xmm1
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm8,%xmm8
ucomisd %xmm8,%xmm0
jb     1543 <qrsolv_+0x413>
mov    %r15,%rax
sqrtsd %xmm0,%xmm0
jmp    1569 <qrsolv_+0x439>
movapd %xmm11,0xb0(%rsp)
call   1050 <sqrt@plt>
movapd 0xb0(%rsp),%xmm11
xorpd  %xmm8,%xmm8
mov    0x128(%rsp),%rax
mov    0x130(%rsp),%r8
mov    0x30(%rsp),%r15
mov    0x28(%rsp),%r12
movsd  0x40(%rsp),%xmm9
mov    0x38(%rsp),%r10
mov    0x18(%rsp),%r11
movsd  0xad1(%rip),%xmm10
divsd  %xmm0,%xmm10
mulsd  %xmm10,%xmm11
jmp    159f <qrsolv_+0x46f>
movapd %xmm10,0xb0(%rsp)
call   1050 <sqrt@plt>
movapd 0xb0(%rsp),%xmm10
xorpd  %xmm8,%xmm8
mov    0x128(%rsp),%rax
mov    0x130(%rsp),%r8
mov    0x30(%rsp),%r15
mov    0x28(%rsp),%r12
movsd  0x40(%rsp),%xmm9
mov    0x38(%rsp),%r10
mov    0x18(%rsp),%r11
movsd  0xa73(%rip),%xmm11
divsd  %xmm0,%xmm11
mulsd  %xmm11,%xmm10
movsd  (%rax,%rbx,8),%xmm0
mulsd  %xmm10,%xmm0
mov    0x8(%rsp),%rcx
movsd  (%rcx,%rbp,8),%xmm1
mulsd  %xmm11,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx,%rbp,8)
movsd  (%r8,%rbx,8),%xmm0
movapd %xmm9,%xmm1
mulsd  %xmm10,%xmm1
movapd %xmm11,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm11,%xmm9
mulsd  %xmm10,%xmm0
subsd  %xmm0,%xmm9
movsd  %xmm2,(%r8,%rbx,8)
lea    0x1(%rbx),%r9
cmp    %r12,%r9
jge    16c2 <qrsolv_+0x592>
mov    0xd8(%rsp),%rcx
sub    %r11,%rcx
mov    0x68(%rsp),%rbp
add    %rcx,%rbp
imul   %r15,%rbx
movapd %xmm10,%xmm0
xorpd  0x9ec(%rip),%xmm0
mov    0x20(%rsp),%r14
mov    %r14,%rcx
cmp    $0x2,%rbp
mov    0x60(%rsp),%r8
jae    16de <qrsolv_+0x5ae>
mov    0x8(%rsp),%rdx
lea    (%rdx,%rbx,8),%rdx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rdx,%rcx,8),%xmm1
mov    %ecx,%esi
movsd  (%rax,%rsi,8),%xmm2
movapd %xmm10,%xmm3
mulsd  %xmm2,%xmm3
movapd %xmm11,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm11,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rax,%rsi,8)
movsd  %xmm4,(%rdx,%rcx,8)
add    $0x1,%rcx
cmp    %r12,%rcx
jb     1650 <qrsolv_+0x520>
mov    %r9,%rbx
mov    0x58(%rsp),%rcx
mov    0x48(%rsp),%rdx
mov    0x70(%rsp),%rsi
mov    0x20(%rsp),%r14
mov    %r8,%r9
mov    0x130(%rsp),%r8
cmp    %r12,%rbx
jne    13ba <qrsolv_+0x28a>
jmp    1270 <qrsolv_+0x140>
mov    %r9,%rbx
mov    0x60(%rsp),%r9
mov    0x58(%rsp),%rcx
mov    0x48(%rsp),%rdx
mov    0x70(%rsp),%rsi
jmp    13ac <qrsolv_+0x27c>
lea    (%r8,%r11,1),%rdx
mov    $0xfffffffffffffffe,%rcx
sub    %rdx,%rcx
add    %rcx,%r13
mov    %r13d,%esi
mov    %r14,%rcx
add    %r14d,%esi
jb     163b <qrsolv_+0x50b>
shr    $0x20,%r13
mov    %r14,%rcx
jne    163b <qrsolv_+0x50b>
neg    %rdx
mov    %r14d,%ecx
mov    0xa8(%rsp),%rsi
add    %rsi,%rdx
add    $0xfffffffffffffffe,%rdx
add    %rcx,%rdx
mov    0xd0(%rsp),%rsi
lea    (%rsi,%rdx,8),%rdx
lea    (%r14,%rbx,1),%rsi
mov    0x8(%rsp),%rdi
lea    (%rdi,%rsi,8),%rsi
cmp    %rdx,%rsi
jae    176c <qrsolv_+0x63c>
lea    (%rax,%rcx,8),%rdx
mov    0xa8(%rsp),%rsi
add    %rbx,%rsi
mov    0xc8(%rsp),%rcx
lea    (%rcx,%rsi,8),%rsi
add    $0xfffffffffffffff8,%rsi
mov    %r14,%rcx
cmp    %rsi,%rdx
jb     163b <qrsolv_+0x50b>
mov    0xe0(%rsp),%rdx
add    0x48(%rsp),%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rbp,0x68(%rsp)
mov    %rbp,%r13
and    $0xfffffffffffffffe,%r13
lea    (%r14,%r13,1),%rcx
movapd %xmm10,%xmm1
unpcklpd %xmm10,%xmm1
movapd %xmm11,%xmm2
unpcklpd %xmm11,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
lea    (%r10,%rbx,8),%rsi
xor    %edi,%edi
nop
lea    (%r14,%rdi,1),%ebp
movupd (%rsi,%rdi,8),%xmm4
movupd (%rax,%rbp,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,(%rax,%rbp,8)
movupd %xmm7,(%rsi,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rdx
jne    17b0 <qrsolv_+0x680>
cmp    %r13,0x68(%rsp)
mov    0x60(%rsp),%r8
mov    0x18(%rsp),%r11
jne    163b <qrsolv_+0x50b>
mov    %r9,%rbx
mov    0x58(%rsp),%rcx
mov    0x48(%rsp),%rdx
mov    0x70(%rsp),%rsi
jmp    16a9 <qrsolv_+0x579>
cmpl   $0x0,0x4(%rsp)
jle    1a07 <qrsolv_+0x8d7>
mov    %r12,%r9
add    $0xffffffffffffffff,%r9
je     18a6 <qrsolv_+0x776>
mov    %r12d,%r10d
and    $0xfffffffe,%r10d
xor    %edi,%edi
xorpd  %xmm0,%xmm0
mov    0x4(%rsp),%esi
mov    0x4(%rsp),%ecx
jmp    1859 <qrsolv_+0x729>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %r10,%rdi
je     18a0 <qrsolv_+0x770>
mov    %rdi,%rdx
cmp    %ecx,%esi
mov    %edx,%ebp
cmovne %esi,%ebp
ucomisd (%rax,%rdi,8),%xmm0
cmovne %esi,%ebp
cmovp  %esi,%ebp
cmp    %ecx,%ebp
jge    187a <qrsolv_+0x74a>
movq   $0x0,(%r8,%rdx,8)
lea    0x1(%rdx),%rdi
cmp    %ecx,%ebp
mov    %edi,%esi
cmovne %ebp,%esi
ucomisd 0x8(%rax,%rdx,8),%xmm0
cmovne %ebp,%esi
cmovp  %ebp,%esi
cmp    %ecx,%esi
jge    1850 <qrsolv_+0x720>
movq   $0x0,0x8(%r8,%rdx,8)
jmp    1850 <qrsolv_+0x720>
add    $0x2,%rdx
jmp    18ac <qrsolv_+0x77c>
xor    %edx,%edx
mov    0x4(%rsp),%esi
mov    %esi,%ecx
test   $0x1,%r12b
je     18da <qrsolv_+0x7aa>
mov    0x4(%rsp),%edi
cmp    %edi,%esi
mov    %edx,%ecx
cmovne %esi,%ecx
xorpd  %xmm0,%xmm0
ucomisd (%rax,%rdx,8),%xmm0
cmovne %esi,%ecx
cmovp  %esi,%ecx
cmp    %edi,%ecx
jge    18da <qrsolv_+0x7aa>
movq   $0x0,(%r8,%rdx,8)
test   %ecx,%ecx
jle    196a <qrsolv_+0x83a>
mov    %ecx,%edx
lea    -0x1(%rdx),%rsi
imul   %r15,%rsi
mov    0x8(%rsp),%rdi
lea    (%rdi,%rsi,8),%rsi
shl    $0x3,%r15
neg    %r15
xor    %ebp,%ebp
jmp    1930 <qrsolv_+0x800>
xorpd  %xmm0,%xmm0
movsd  (%r8,%rdi,8),%xmm1
subsd  %xmm0,%xmm1
divsd  (%rax,%rdi,8),%xmm1
movsd  %xmm1,(%r8,%rdi,8)
add    $0x1,%rbp
add    $0xffffffff,%ecx
add    %r15,%rsi
cmp    %rdx,%rbp
mov    0x120(%rsp),%rbx
je     196a <qrsolv_+0x83a>
mov    %rbp,%rdi
not    %rdi
add    %rdx,%rdi
test   %rbp,%rbp
je     1900 <qrsolv_+0x7d0>
movslq %ecx,%rbx
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rsi,%rbx,8),%xmm1
mulsd  (%r8,%rbx,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rbx
cmp    %rdx,%rbx
jl     1950 <qrsolv_+0x820>
jmp    1904 <qrsolv_+0x7d4>
cmpl   $0x0,0x4(%rsp)
jle    1a07 <qrsolv_+0x8d7>
mov    %r12d,%ecx
and    $0x3,%ecx
cmp    $0x3,%r9
jae    1985 <qrsolv_+0x855>
xor    %eax,%eax
jmp    19db <qrsolv_+0x8ab>
and    $0xfffffffc,%r12d
xor    %eax,%eax
mov    0x50(%rsp),%rsi
movslq (%rsi,%rax,4),%rdx
movsd  (%r8,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
movslq 0x4(%rsi,%rax,4),%rdx
movsd  0x8(%r8,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
movslq 0x8(%rsi,%rax,4),%rdx
movsd  0x10(%r8,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
movslq 0xc(%rsi,%rax,4),%rdx
movsd  0x18(%r8,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
add    $0x4,%rax
cmp    %rax,%r12
jne    1990 <qrsolv_+0x860>
test   %rcx,%rcx
je     1a07 <qrsolv_+0x8d7>
lea    (%r8,%rax,8),%rdx
mov    0x50(%rsp),%rsi
lea    (%rsi,%rax,4),%rax
xor    %esi,%esi
nop
movslq (%rax,%rsi,4),%rdi
movsd  (%rdx,%rsi,8),%xmm0
movsd  %xmm0,(%rbx,%rdi,8)
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    19f0 <qrsolv_+0x8c0>
xor    %eax,%eax
add    $0xe8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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