#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 */
double DBL_EPSILON;
double enorm_(int, double *);
double sqrt(double);
int qrfac_(int m, int n, double *a, int lda, int *ipvt, double *rdiag,
           double *acnorm, double *wa) {
  const double p05 = .05;
  const double epsmch = DBL_EPSILON;
  double ajnorm, d, sum, temp;
  int kmax, minmn, jp1;
  int i, j, k;
  for (j = 0; j < n; ++j) {
    acnorm[j] = enorm_(m, a + j * lda);
    rdiag[j] = acnorm[j];
    wa[j] = rdiag[j];
    ipvt[j] = j;
  }
  minmn = ((m) <= (n) ? (m) : (n));
  for (j = 0; j < minmn; ++j) {
    kmax = j;
    for (k = j; k < n; ++k) {
      if (rdiag[k] > rdiag[kmax]) {
        kmax = k;
      }
    }
    if (kmax != j) {
      for (i = 0; i < m; ++i) {
        temp = a[i + j * lda];
        a[i + j * lda] = a[i + kmax * lda];
        a[i + kmax * lda] = temp;
      }
      rdiag[kmax] = rdiag[j];
      wa[kmax] = wa[j];
      k = ipvt[j];
      ipvt[j] = ipvt[kmax];
      ipvt[kmax] = k;
    }
    ajnorm = enorm_(m - j, &a[j + j * lda]);
    if (ajnorm == 0.0) {
      rdiag[j] = -ajnorm;
      continue;
    }
    if (a[j + j * lda] < 0.0) {
      ajnorm = -ajnorm;
    }
    for (i = j; i < m; ++i) {
      a[i + j * lda] /= ajnorm;
    }
    a[j + j * lda] += 1.0;
    jp1 = j + 1;
    for (k = jp1; k < n; ++k) {
      sum = 0.0;
      for (i = j; i < m; ++i) {
        sum += a[i + j * lda] * a[i + k * lda];
      }
      temp = sum / a[j + j * lda];
      for (i = j; i < m; ++i) {
        a[i + k * lda] -= temp * a[i + j * lda];
      }
      if (rdiag[k] != 0.0) {
        temp = a[j + k * lda] / rdiag[k];
        d = 1.0 - temp * temp;
        if (d > 0) {
          rdiag[k] *= sqrt(d);
        } else {
          rdiag[k] = 0.0;
        }
        d = rdiag[k] / wa[k];
        if (p05 * (d * d) <= epsmch) {
          rdiag[k] = enorm_(m - jp1, &a[jp1 + k * lda]);
          wa[k] = rdiag[k];
        }
      }
    }
    rdiag[j] = -ajnorm;
  }
  return 0;
}


x86 assembly (clang, -O0):
<qrfac_>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movsd  0xeb0(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    0x2e8c(%rip),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11f1 <qrfac_+0xd1>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x6c(%rbp),%eax
imul   -0x14(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
call   1030 <enorm_@plt>
mov    0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%edx
mov    -0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    116c <qrfac_+0x4c>
mov    -0x4(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     1208 <qrfac_+0xe8>
mov    -0x4(%rbp),%eax
mov    %eax,-0x74(%rbp)
jmp    120e <qrfac_+0xee>
mov    -0x8(%rbp),%eax
mov    %eax,-0x74(%rbp)
mov    -0x74(%rbp),%eax
mov    %eax,-0x60(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x60(%rbp),%eax
jge    1757 <qrfac_+0x637>
mov    -0x6c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1278 <qrfac_+0x158>
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1265 <qrfac_+0x145>
mov    -0x70(%rbp),%eax
mov    %eax,-0x5c(%rbp)
jmp    126a <qrfac_+0x14a>
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    1233 <qrfac_+0x113>
mov    -0x5c(%rbp),%eax
cmp    -0x6c(%rbp),%eax
je     1375 <qrfac_+0x255>
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    130f <qrfac_+0x1ef>
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x5c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x5c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    128b <qrfac_+0x16b>
mov    -0x28(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x70(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x70(%rbp),%edx
mov    -0x20(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x4(%rbp),%edi
sub    -0x6c(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x6c(%rbp),%eax
mov    -0x6c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
call   1030 <enorm_@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    13e4 <qrfac_+0x2c4>
jp     13e4 <qrfac_+0x2c4>
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1749 <qrfac_+0x629>
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    1426 <qrfac_+0x306>
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    146c <qrfac_+0x34c>
movsd  -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    142c <qrfac_+0x30c>
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xb81(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1720 <qrfac_+0x600>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1516 <qrfac_+0x3f6>
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x70(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    14ba <qrfac_+0x39a>
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15b1 <qrfac_+0x491>
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x6c(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x70(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    153e <qrfac_+0x41e>
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    15d6 <qrfac_+0x4b6>
jp     15d6 <qrfac_+0x4b6>
jmp    170d <qrfac_+0x5ed>
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    -0x70(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0x9db(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1669 <qrfac_+0x549>
movsd  -0x48(%rbp),%xmm0
call   1040 <sqrt@plt>
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1679 <qrfac_+0x559>
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x70(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  0x956(%rip),%xmm1
mulsd  %xmm0,%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     1708 <qrfac_+0x5e8>
mov    -0x4(%rbp),%edi
sub    -0x64(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x64(%rbp),%eax
mov    -0x70(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
call   1030 <enorm_@plt>
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x70(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x70(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    170d <qrfac_+0x5ed>
jmp    1712 <qrfac_+0x5f2>
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    14a0 <qrfac_+0x380>
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    121b <qrfac_+0xfb>
xor    %eax,%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qrfac_>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x168,%rsp
mov    %r8,%r12
mov    %rdx,0x38(%rsp)
mov    %edi,%r13d
mov    0x1a8(%rsp),%r11
mov    0x2e9d(%rip),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,0x138(%rsp)
mov    %ecx,0x60(%rsp)
movslq %ecx,%rax
mov    %rax,0x30(%rsp)
mov    %esi,%eax
mov    %esi,0x48(%rsp)
test   %esi,%esi
mov    %rax,%rsi
mov    %r9,0x20(%rsp)
mov    %rax,0x28(%rsp)
jle    11df <qrfac_+0xbf>
mov    0x1a0(%rsp),%r14
mov    0x30(%rsp),%rax
lea    0x0(,%rax,8),%r15
xor    %ebx,%ebx
mov    0x38(%rsp),%rbp
nopw   0x0(%rax,%rax,1)
mov    %r13d,%edi
mov    %rbp,%rsi
call   1030 <enorm_@plt>
mov    0x28(%rsp),%rsi
mov    0x1a8(%rsp),%r11
mov    0x20(%rsp),%r9
movsd  %xmm0,(%r14,%rbx,8)
movsd  %xmm0,(%r9,%rbx,8)
movsd  %xmm0,(%r11,%rbx,8)
mov    %ebx,(%r12,%rbx,4)
add    $0x1,%rbx
add    %r15,%rbp
cmp    %rbx,%rsi
jne    11a0 <qrfac_+0x80>
mov    0x48(%rsp),%edx
cmp    %edx,%r13d
mov    %r13d,%eax
cmovg  %edx,%eax
test   %eax,%eax
mov    0x38(%rsp),%r10
mov    %r13d,0xc(%rsp)
mov    0x30(%rsp),%rcx
jle    1c57 <qrfac_+0xb37>
movslq %edx,%rdx
mov    %rdx,0xf0(%rsp)
mov    %eax,%r15d
mov    0xc(%rsp),%ebp
lea    (%rcx,%rbp,1),%rax
mov    %rax,0xe0(%rsp)
lea    0x1(%rcx),%rax
mov    %rax,0x70(%rsp)
lea    (%r10,%rbp,8),%r13
mov    %ebp,%eax
and    $0xfffffffc,%eax
mov    %rax,0xd8(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,0xd0(%rsp)
mov    %rax,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %rdx,%rax
mov    %rdx,0xc8(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0xc0(%rsp)
mov    %rbp,%rdi
neg    %rdi
lea    0x30(%r10),%rax
mov    %rax,0xb8(%rsp)
lea    0x8(%r10),%rax
mov    %rax,0x78(%rsp)
lea    0x10(%r10),%rax
mov    %rax,0x68(%rsp)
lea    0x8(%r10,%rcx,8),%rax
mov    %rax,0x98(%rsp)
lea    0x0(,%rcx,8),%rax
mov    %rax,0xa8(%rsp)
lea    0x10(%r10,%rcx,8),%rax
mov    %rax,0x90(%rsp)
mov    $0x1,%eax
mov    %rax,0x88(%rsp)
xor    %eax,%eax
mov    %rax,0x80(%rsp)
xor    %edx,%edx
mov    %r15,0xe8(%rsp)
mov    %r13,0x58(%rsp)
mov    %r12,0xf8(%rsp)
mov    %rdi,0x18(%rsp)
jmp    1357 <qrfac_+0x237>
nopl   0x0(%rax)
movapd 0x150(%rsp),%xmm0
xorpd  0xcff(%rip),%xmm0
movlpd %xmm0,(%r9,%rdx,8)
mov    0x50(%rsp),%rcx
mov    %rcx,%rdx
mov    0xf8(%rsp),%r12
mov    0xe8(%rsp),%r15
mov    0x58(%rsp),%r13
addq   $0x1,0x88(%rsp)
addq   $0x1,0x80(%rsp)
mov    0xa8(%rsp),%rax
add    %rax,0x98(%rsp)
add    %rax,0x90(%rsp)
cmp    %r15,%rcx
je     1c57 <qrfac_+0xb37>
mov    %esi,%eax
sub    %edx,%eax
mov    %rdx,%rcx
not    %rcx
mov    %rcx,0x148(%rsp)
add    %rcx,%rsi
mov    %rdx,%rcx
mov    %rdx,0x10(%rsp)
and    $0x3,%rax
je     13c0 <qrfac_+0x2a0>
mov    0x10(%rsp),%rdx
mov    %rdx,%rcx
mov    %edx,%edi
jmp    139c <qrfac_+0x27c>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rcx
mov    %edx,%edi
add    $0xffffffffffffffff,%rax
je     13c0 <qrfac_+0x2a0>
movsd  (%r9,%rcx,8),%xmm0
movslq %edi,%rdx
ucomisd (%r9,%rdx,8),%xmm0
mov    %ecx,%edx
ja     1390 <qrfac_+0x270>
mov    %edi,%edx
jmp    1390 <qrfac_+0x270>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %edx,%eax
cmp    $0x3,%rsi
mov    0x28(%rsp),%rsi
jae    14df <qrfac_+0x3bf>
mov    %eax,%ecx
mov    0x10(%rsp),%rbx
cmp    %rcx,%rbx
je     165f <qrfac_+0x53f>
cmpl   $0x0,0xc(%rsp)
mov    0x18(%rsp),%rdi
jle    1630 <qrfac_+0x510>
mov    0x10(%rsp),%rcx
mov    0x30(%rsp),%rsi
imul   %rsi,%rcx
mov    %eax,%edx
imul   %esi,%edx
movslq %edx,%r8
cmpl   $0x4,0xc(%rsp)
jb     143e <qrfac_+0x31e>
lea    (%r10,%rcx,8),%rdx
lea    0x0(,%r8,8),%rsi
add    %r13,%rsi
cmp    %rsi,%rdx
jae    1535 <qrfac_+0x415>
lea    0x0(,%rcx,8),%rdx
add    %r13,%rdx
lea    (%r10,%r8,8),%rsi
cmp    %rdx,%rsi
jae    1535 <qrfac_+0x415>
xor    %edx,%edx
mov    %rdx,%rsi
test   $0x1,%bpl
je     1475 <qrfac_+0x355>
lea    (%rdx,%rcx,1),%rsi
movsd  (%r10,%rsi,8),%xmm0
lea    (%rdx,%r8,1),%rdi
movsd  (%r10,%rdi,8),%xmm1
movsd  %xmm1,(%r10,%rsi,8)
movsd  %xmm0,(%r10,%rdi,8)
mov    0x18(%rsp),%rdi
mov    %rdx,%rsi
or     $0x1,%rsi
not    %rdx
cmp    %rdi,%rdx
je     1630 <qrfac_+0x510>
mov    0x78(%rsp),%rdi
lea    (%rdi,%r8,8),%rdx
lea    (%rdi,%rcx,8),%rcx
xchg   %ax,%ax
movsd  -0x8(%rcx,%rsi,8),%xmm0
movsd  -0x8(%rdx,%rsi,8),%xmm1
movsd  %xmm1,-0x8(%rcx,%rsi,8)
movsd  %xmm0,-0x8(%rdx,%rsi,8)
movsd  (%rcx,%rsi,8),%xmm0
movsd  (%rdx,%rsi,8),%xmm1
movsd  %xmm1,(%rcx,%rsi,8)
movsd  %xmm0,(%rdx,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%rbp
jne    1490 <qrfac_+0x370>
jmp    1630 <qrfac_+0x510>
nopw   0x0(%rax,%rax,1)
add    $0x4,%rcx
mov    %eax,%edx
cmp    %rsi,%rcx
je     13d1 <qrfac_+0x2b1>
movslq %edx,%rax
movsd  (%r9,%rcx,8),%xmm0
ucomisd (%r9,%rax,8),%xmm0
mov    %ecx,%eax
ja     14f4 <qrfac_+0x3d4>
mov    %edx,%eax
movsd  0x8(%r9,%rcx,8),%xmm0
movslq %eax,%rdx
ucomisd (%r9,%rdx,8),%xmm0
jbe    1509 <qrfac_+0x3e9>
lea    0x1(%rcx),%eax
movsd  0x10(%r9,%rcx,8),%xmm0
movslq %eax,%rdx
ucomisd (%r9,%rdx,8),%xmm0
jbe    151e <qrfac_+0x3fe>
lea    0x2(%rcx),%eax
movsd  0x18(%r9,%rcx,8),%xmm0
movslq %eax,%rdx
ucomisd (%r9,%rdx,8),%xmm0
jbe    14d0 <qrfac_+0x3b0>
lea    0x3(%rcx),%eax
jmp    14d0 <qrfac_+0x3b0>
cmpq   $0x0,0xd0(%rsp)
je     1c42 <qrfac_+0xb22>
mov    0xb8(%rsp),%rdx
lea    (%rdx,%rcx,8),%rdi
lea    (%rdx,%r8,8),%rbx
mov    0xc0(%rsp),%rdx
xor    %esi,%esi
xchg   %ax,%ax
movups -0x30(%rdi,%rsi,8),%xmm0
movups -0x20(%rdi,%rsi,8),%xmm1
movups -0x30(%rbx,%rsi,8),%xmm2
movups -0x20(%rbx,%rsi,8),%xmm3
movups %xmm2,-0x30(%rdi,%rsi,8)
movups %xmm3,-0x20(%rdi,%rsi,8)
movups %xmm0,-0x30(%rbx,%rsi,8)
movups %xmm1,-0x20(%rbx,%rsi,8)
movupd -0x10(%rdi,%rsi,8),%xmm0
movupd (%rdi,%rsi,8),%xmm1
movupd -0x10(%rbx,%rsi,8),%xmm2
movupd (%rbx,%rsi,8),%xmm3
movupd %xmm2,-0x10(%rdi,%rsi,8)
movupd %xmm3,(%rdi,%rsi,8)
movupd %xmm0,-0x10(%rbx,%rsi,8)
movupd %xmm1,(%rbx,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdx
jne    1560 <qrfac_+0x440>
testb  $0x1,0xc8(%rsp)
01
je     1603 <qrfac_+0x4e3>
lea    (%rsi,%rcx,1),%rdx
movupd (%r10,%rdx,8),%xmm0
movupd 0x10(%r10,%rdx,8),%xmm1
add    %r8,%rsi
movupd (%r10,%rsi,8),%xmm2
movupd 0x10(%r10,%rsi,8),%xmm3
movupd %xmm2,(%r10,%rdx,8)
movupd %xmm3,0x10(%r10,%rdx,8)
movupd %xmm0,(%r10,%rsi,8)
movupd %xmm1,0x10(%r10,%rsi,8)
mov    0xd8(%rsp),%rsi
mov    %rsi,%rdx
cmp    %rbp,%rsi
mov    0x58(%rsp),%r13
mov    0x18(%rsp),%rdi
jne    1440 <qrfac_+0x320>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x10(%rsp),%rbx
movsd  (%r9,%rbx,8),%xmm0
cltq
movsd  %xmm0,(%r9,%rax,8)
movsd  (%r11,%rbx,8),%xmm0
movsd  %xmm0,(%r11,%rax,8)
mov    (%r12,%rbx,4),%ecx
mov    (%r12,%rax,4),%edx
mov    %edx,(%r12,%rbx,4)
mov    %ecx,(%r12,%rax,4)
mov    %ebx,%eax
imul   0x60(%rsp),%eax
cltq
add    %rbx,%rax
lea    (%r10,%rax,8),%r14
mov    0xc(%rsp),%edi
sub    %ebx,%edi
mov    %r14,%rsi
call   1030 <enorm_@plt>
xorpd  %xmm6,%xmm6
ucomisd %xmm6,%xmm0
jne    16c0 <qrfac_+0x5a0>
jp     16c0 <qrfac_+0x5a0>
xorpd  0x96f(%rip),%xmm0
mov    0x20(%rsp),%r9
movlpd %xmm0,(%r9,%rbx,8)
add    $0x1,%rbx
mov    %rbx,%rcx
mov    0x38(%rsp),%r10
mov    0x1a8(%rsp),%r11
mov    0x28(%rsp),%rsi
mov    %rbx,%rdx
jmp    1324 <qrfac_+0x204>
nopl   (%rax)
mov    %rbx,%r15
mov    0x30(%rsp),%r8
imul   %r8,%r15
movsd  (%r14),%xmm2
movapd %xmm0,%xmm1
xorpd  0x923(%rip),%xmm1
cmpltsd %xmm6,%xmm2
andpd  %xmm2,%xmm1
andnpd %xmm0,%xmm2
orpd   %xmm1,%xmm2
movapd %xmm2,%xmm1
mov    %rbp,%rcx
sub    %rbx,%rcx
mov    %rbx,%rax
mov    %rcx,%rdx
mov    %rcx,0x40(%rsp)
cmp    $0x2,%rcx
mov    0x38(%rsp),%r10
mov    0x1a8(%rsp),%r11
mov    0x28(%rsp),%rsi
mov    0x18(%rsp),%rdi
mov    %rbx,%rdx
movapd %xmm2,0x150(%rsp)
jb     17e7 <qrfac_+0x6c7>
mov    0x40(%rsp),%rax
and    $0xfffffffffffffffe,%rax
movapd %xmm1,%xmm0
unpcklpd %xmm1,%xmm0
lea    -0x2(%rax),%rdx
mov    %rdx,%rcx
shr    %rcx
add    $0x1,%rcx
test   %rdx,%rdx
je     1c2d <qrfac_+0xb0d>
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
mov    0x70(%rsp),%rdx
imul   0x10(%rsp),%rdx
mov    0x68(%rsp),%rdi
lea    (%rdi,%rdx,8),%rdi
xor    %edx,%edx
mov    0x20(%rsp),%r9
nopw   0x0(%rax,%rax,1)
movupd -0x10(%rdi,%rdx,8),%xmm1
movupd (%rdi,%rdx,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,-0x10(%rdi,%rdx,8)
divpd  %xmm0,%xmm2
movupd %xmm2,(%rdi,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1780 <qrfac_+0x660>
test   $0x1,%cl
je     17c5 <qrfac_+0x6a5>
add    0x10(%rsp),%rdx
add    %r15,%rdx
movupd (%r10,%rdx,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%r10,%rdx,8)
cmp    %rax,0x40(%rsp)
mov    0x28(%rsp),%rsi
mov    0x18(%rsp),%rdi
mov    0x10(%rsp),%rdx
movapd 0x150(%rsp),%xmm1
je     1807 <qrfac_+0x6e7>
add    %rdx,%rax
lea    (%r10,%r15,8),%rcx
mov    0x20(%rsp),%r9
movsd  (%rcx,%rax,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rax
cmp    %rax,%rbp
jne    17f0 <qrfac_+0x6d0>
mov    %r15,0xa0(%rsp)
movsd  (%r14),%xmm0
addsd  0x7f4(%rip),%xmm0
mov    %r14,0x140(%rsp)
movsd  %xmm0,(%r14)
lea    0x1(%rdx),%rcx
mov    %rcx,%rax
mov    %rcx,0x50(%rsp)
cmp    0xf0(%rsp),%rcx
jge    12f0 <qrfac_+0x1d0>
mov    0x80(%rsp),%rcx
imul   %r8,%rcx
mov    0xe0(%rsp),%rax
add    %rcx,%rax
mov    %rax,0x118(%rsp)
add    %r8,%rcx
mov    0xc(%rsp),%eax
sub    0x50(%rsp),%eax
mov    %eax,0x64(%rsp)
add    %rdx,%rcx
mov    %rcx,0x130(%rsp)
mov    0x70(%rsp),%rax
imul   %rdx,%rax
lea    (%r10,%rax,8),%rcx
mov    %rcx,0x110(%rsp)
mov    0x58(%rsp),%rcx
mov    0xa0(%rsp),%rbx
lea    (%rcx,%rbx,8),%rcx
mov    %rcx,0x128(%rsp)
lea    (%rdx,%rbx,1),%rcx
mov    %rcx,0x120(%rsp)
mov    0x40(%rsp),%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x108(%rsp)
add    %rdx,%rcx
mov    %rcx,0x100(%rsp)
mov    0x78(%rsp),%rcx
lea    (%rcx,%rbx,8),%r14
mov    0x68(%rsp),%rcx
lea    (%rcx,%rax,8),%r12
mov    0x90(%rsp),%rax
lea    (%rax,%rdx,8),%r15
mov    0x98(%rsp),%r13
xor    %eax,%eax
mov    %rax,0xb0(%rsp)
mov    0x88(%rsp),%rbx
jmp    1934 <qrfac_+0x814>
nopl   0x0(%rax)
add    $0x1,%rbx
addq   $0x1,0xb0(%rsp)
mov    0xa8(%rsp),%rax
add    %rax,%r13
add    %rax,%r15
cmp    %rsi,%rbx
je     12f0 <qrfac_+0x1d0>
mov    %rbx,%rax
imul   %r8,%rax
mov    %rax,0x48(%rsp)
xorpd  %xmm1,%xmm1
mov    %rdx,%rax
testb  $0x1,0x40(%rsp)
je     1977 <qrfac_+0x857>
mov    0x120(%rsp),%rax
movsd  (%r10,%rax,8),%xmm0
mov    0x48(%rsp),%rax
add    %rdx,%rax
mulsd  (%r10,%rax,8),%xmm0
addsd  %xmm6,%xmm0
mov    0x50(%rsp),%rax
movapd %xmm0,%xmm1
cmp    %rdi,0x148(%rsp)
je     19c0 <qrfac_+0x8a0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%r14,%rax,8),%xmm2
movsd  (%r14,%rax,8),%xmm0
mulsd  -0x8(%r13,%rax,8),%xmm2
mulsd  0x0(%r13,%rax,8),%xmm0
addsd  %xmm1,%xmm2
addsd  %xmm2,%xmm0
add    $0x2,%rax
movapd %xmm0,%xmm1
cmp    %rax,%rbp
jne    1990 <qrfac_+0x870>
xorpd  0x638(%rip),%xmm0
mov    0x140(%rsp),%rax
divsd  (%rax),%xmm0
mov    %rdx,%rax
cmpq   $0x4,0x40(%rsp)
jb     1aa0 <qrfac_+0x980>
mov    0xb0(%rsp),%rax
imul   %r8,%rax
mov    0x130(%rsp),%rcx
add    %rax,%rcx
lea    (%r10,%rcx,8),%rcx
cmp    0x128(%rsp),%rcx
jae    1a21 <qrfac_+0x901>
add    0x118(%rsp),%rax
lea    (%r10,%rax,8),%rcx
mov    %rdx,%rax
cmp    %rcx,0x110(%rsp)
jb     1aa0 <qrfac_+0x980>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %eax,%eax
mov    0x108(%rsp),%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x10(%r12,%rax,8),%xmm2
movupd (%r12,%rax,8),%xmm3
movupd -0x10(%r15,%rax,8),%xmm4
movupd (%r15,%rax,8),%xmm5
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,-0x10(%r15,%rax,8)
movupd %xmm3,(%r15,%rax,8)
add    $0x4,%rax
cmp    %rax,%rcx
jne    1a40 <qrfac_+0x920>
mov    0x100(%rsp),%rax
cmp    %rcx,0x40(%rsp)
je     1b2a <qrfac_+0xa0a>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %ebp,%edx
sub    %eax,%edx
mov    %rax,%rcx
test   $0x1,%dl
je     1ad9 <qrfac_+0x9b9>
mov    0xa0(%rsp),%rcx
add    %rax,%rcx
mov    0x48(%rsp),%rdx
add    %rax,%rdx
movsd  (%r10,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%r10,%rdx,8),%xmm1
movsd  %xmm1,(%r10,%rdx,8)
lea    0x1(%rax),%rcx
not    %rax
cmp    %rdi,%rax
mov    0x10(%rsp),%rdx
je     1b2a <qrfac_+0xa0a>
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%r14,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  -0x8(%r13,%rcx,8),%xmm1
movsd  %xmm1,-0x8(%r13,%rcx,8)
movsd  (%r14,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x0(%r13,%rcx,8),%xmm1
movsd  %xmm1,0x0(%r13,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rbp
jne    1af0 <qrfac_+0x9d0>
movsd  (%r9,%rbx,8),%xmm0
ucomisd %xmm6,%xmm0
jne    1b3c <qrfac_+0xa1c>
jnp    1910 <qrfac_+0x7f0>
mov    0x48(%rsp),%rax
add    %rdx,%rax
movsd  (%r10,%rax,8),%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm2,%xmm2
movsd  0x4b6(%rip),%xmm1
subsd  %xmm2,%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jbe    1bb0 <qrfac_+0xa90>
ucomisd %xmm6,%xmm1
jb     1b77 <qrfac_+0xa57>
xorps  %xmm0,%xmm0
sqrtsd %xmm1,%xmm0
jmp    1baa <qrfac_+0xa8a>
movapd %xmm1,%xmm0
call   1040 <sqrt@plt>
mov    0x10(%rsp),%rdx
xorpd  %xmm6,%xmm6
mov    0x18(%rsp),%rdi
mov    0x28(%rsp),%rsi
mov    0x30(%rsp),%r8
mov    0x1a8(%rsp),%r11
mov    0x38(%rsp),%r10
mov    0x20(%rsp),%r9
mulsd  (%r9,%rbx,8),%xmm0
movsd  %xmm0,(%r9,%rbx,8)
divsd  (%r11,%rbx,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  0x450(%rip),%xmm0
movsd  0x138(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jb     1910 <qrfac_+0x7f0>
mov    0x48(%rsp),%rax
add    0x50(%rsp),%rax
lea    (%r10,%rax,8),%rsi
mov    0x64(%rsp),%edi
call   1030 <enorm_@plt>
mov    0x10(%rsp),%rdx
xorpd  %xmm6,%xmm6
mov    0x18(%rsp),%rdi
mov    0x28(%rsp),%rsi
mov    0x30(%rsp),%r8
mov    0x1a8(%rsp),%r11
mov    0x38(%rsp),%r10
mov    0x20(%rsp),%r9
movsd  %xmm0,(%r9,%rbx,8)
movsd  %xmm0,(%r11,%rbx,8)
jmp    1910 <qrfac_+0x7f0>
xor    %edx,%edx
mov    0x20(%rsp),%r9
test   $0x1,%cl
jne    17ad <qrfac_+0x68d>
jmp    17c5 <qrfac_+0x6a5>
xor    %esi,%esi
testb  $0x1,0xc8(%rsp)
01
jne    15c8 <qrfac_+0x4a8>
jmp    1603 <qrfac_+0x4e3>
xor    %eax,%eax
add    $0x168,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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