#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 */
int /*<<< orphan*/ free_matrix(double **, int, int, int, int);
int /*<<< orphan*/ free_vector(double *, int, int);
double **matrix(int, int, int, int);
double *vector(int, int);
void mrqmin(double x[], double y[], double sig[], int ndata, double a[],
            int ia[], int ma, double **covar, double **alpha, double *chisq,
            void (*funcs)(double, double[], double *, double[], int),
            double *alamda) {
  void covsrt(double **covar, int ma, int ia[], int mfit);
  void gaussj(double **a, int n, double **b, int m);
  void mrqcof(double x[], double y[], double sig[], int ndata, double a[],
              int ia[], int ma, double **alpha, double beta[], double *chisq,
              void (*funcs)(double, double[], double *, double[], int));
  int j, k, l;
  static int mfit;
  static double ochisq, *atry, *beta, *da, **oneda;
  if (*alamda < 0.0) {
    atry = vector(1, ma);
    beta = vector(1, ma);
    da = vector(1, ma);
    for (mfit = 0, j = 1; j <= ma; j++)
      if (ia[j])
        mfit++;
    oneda = matrix(1, mfit, 1, 1);
    *alamda = 0.001;
    mrqcof(x, y, sig, ndata, a, ia, ma, alpha, beta, chisq, funcs);
    ochisq = (*chisq);
    for (j = 1; j <= ma; j++)
      atry[j] = a[j];
  }
  for (j = 1; j <= mfit; j++) {
    for (k = 1; k <= mfit; k++)
      covar[j][k] = alpha[j][k];
    covar[j][j] = alpha[j][j] * (1.0 + (*alamda));
    oneda[j][1] = beta[j];
  }
  gaussj(covar, mfit, oneda, 1);
  for (j = 1; j <= mfit; j++)
    da[j] = oneda[j][1];
  if (*alamda == 0.0) {
    covsrt(covar, ma, ia, mfit);
    covsrt(alpha, ma, ia, mfit);
    free_matrix(oneda, 1, mfit, 1, 1);
    free_vector(da, 1, ma);
    free_vector(beta, 1, ma);
    free_vector(atry, 1, ma);
    return;
  }
  for (j = 0, l = 1; l <= ma; l++)
    if (ia[l])
      atry[l] = a[l] + da[++j];
  mrqcof(x, y, sig, ndata, atry, ia, ma, covar, da, chisq, funcs);
  if (*chisq < ochisq) {
    *alamda *= 0.1;
    ochisq = (*chisq);
    for (j = 1; j <= mfit; j++) {
      for (k = 1; k <= mfit; k++)
        alpha[j][k] = covar[j][k];
      beta[j] = da[j];
    }
    for (l = 1; l <= ma; l++)
      a[l] = atry[l];
  } else {
    *alamda *= 10.0;
    *chisq = ochisq;
  }
}


x86 assembly (clang, -O0):
<mrqmin>:
push   %rbp
mov    %rsp,%rbp
push   %r14
push   %rbx
sub    $0x70,%rsp
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %ecx,-0x2c(%rbp)
mov    %r8,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
mov    0x38(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jbe    1310 <mrqmin+0x1a0>
mov    0x10(%rbp),%esi
mov    $0x1,%edi
call   1080 <vector@plt>
mov    %rax,0x2ea2(%rip)
mov    0x10(%rbp),%esi
mov    $0x1,%edi
call   1080 <vector@plt>
mov    %rax,0x2e96(%rip)
mov    0x10(%rbp),%esi
mov    $0x1,%edi
call   1080 <vector@plt>
mov    %rax,0x2e8a(%rip)
movl   $0x0,0x2e60(%rip)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x10(%rbp),%eax
jg     1247 <mrqmin+0xd7>
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1234 <mrqmin+0xc4>
mov    0x2e35(%rip),%eax
add    $0x1,%eax
mov    %eax,0x2e2c(%rip)
jmp    1239 <mrqmin+0xc9>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1207 <mrqmin+0x97>
mov    0x2e13(%rip),%esi
mov    $0x1,%ecx
mov    %ecx,%edi
mov    %ecx,%edx
call   1070 <matrix@plt>
mov    %rax,0x2e26(%rip)
mov    0x38(%rbp),%rax
movsd  0xd92(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
mov    -0x2c(%rbp),%ecx
mov    -0x38(%rbp),%r8
mov    -0x40(%rbp),%r9
mov    0x10(%rbp),%r14d
mov    0x20(%rbp),%rbx
mov    0x2de0(%rip),%r11
mov    0x28(%rbp),%r10
mov    0x30(%rbp),%rax
mov    %r14d,(%rsp)
mov    %rbx,0x8(%rsp)
mov    %r11,0x10(%rsp)
mov    %r10,0x18(%rsp)
mov    %rax,0x20(%rsp)
call   1090 <mrqcof@plt>
mov    0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,0x2d9b(%rip)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x10(%rbp),%eax
jg     130b <mrqmin+0x19b>
mov    -0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2d7c(%rip),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12d4 <mrqmin+0x164>
jmp    1310 <mrqmin+0x1a0>
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x2d40(%rip),%eax
jg     13e4 <mrqmin+0x274>
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    0x2d2a(%rip),%eax
jg     1374 <mrqmin+0x204>
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    132d <mrqmin+0x1bd>
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x38(%rbp),%rax
movsd  0xc83(%rip),%xmm1
addsd  (%rax),%xmm1
mulsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x2cbf(%rip),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2cbf(%rip),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1317 <mrqmin+0x1a7>
mov    0x18(%rbp),%rdi
mov    0x2c72(%rip),%esi
mov    0x2c93(%rip),%rdx
mov    $0x1,%ecx
call   1050 <gaussj@plt>
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x2c51(%rip),%eax
jg     1447 <mrqmin+0x2d7>
mov    0x2c6c(%rip),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  0x8(%rax),%xmm0
mov    0x2c50(%rip),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1406 <mrqmin+0x296>
mov    0x38(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    14ed <mrqmin+0x37d>
jp     14ed <mrqmin+0x37d>
mov    0x18(%rbp),%rdi
mov    0x10(%rbp),%esi
mov    -0x40(%rbp),%rdx
mov    0x2bed(%rip),%ecx
call   1060 <covsrt@plt>
mov    0x20(%rbp),%rdi
mov    0x10(%rbp),%esi
mov    -0x40(%rbp),%rdx
mov    0x2bd7(%rip),%ecx
call   1060 <covsrt@plt>
mov    0x2bf3(%rip),%rdi
mov    0x2bc5(%rip),%edx
mov    $0x1,%r8d
mov    %r8d,%esi
mov    %r8d,%ecx
call   1030 <free_matrix@plt>
mov    0x2bcd(%rip),%rdi
mov    0x10(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
mov    0x2bb1(%rip),%rdi
mov    0x10(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
mov    0x2b95(%rip),%rdi
mov    0x10(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
jmp    16e6 <mrqmin+0x576>
movl   $0x0,-0x44(%rbp)
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    0x10(%rbp),%eax
jg     1561 <mrqmin+0x3f1>
mov    -0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     154e <mrqmin+0x3de>
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2b53(%rip),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x44(%rbp)
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    0x2b2b(%rip),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1553 <mrqmin+0x3e3>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    14fb <mrqmin+0x38b>
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
mov    -0x2c(%rbp),%ecx
mov    0x2af9(%rip),%r8
mov    -0x40(%rbp),%r9
mov    0x10(%rbp),%r14d
mov    0x18(%rbp),%rbx
mov    0x2af6(%rip),%r11
mov    0x28(%rbp),%r10
mov    0x30(%rbp),%rax
mov    %r14d,(%rsp)
mov    %rbx,0x8(%rsp)
mov    %r11,0x10(%rsp)
mov    %r10,0x18(%rsp)
mov    %rax,0x20(%rsp)
call   1090 <mrqcof@plt>
mov    0x28(%rbp),%rax
movsd  (%rax),%xmm1
movsd  0x2aa9(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    16c2 <mrqmin+0x552>
mov    0x38(%rbp),%rax
movsd  0xa3b(%rip),%xmm0
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,0x2a7b(%rip)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x2a63(%rip),%eax
jg     167f <mrqmin+0x50f>
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    0x2a4d(%rip),%eax
jg     1651 <mrqmin+0x4e1>
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    160a <mrqmin+0x49a>
mov    0x2a28(%rip),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2a10(%rip),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    15f4 <mrqmin+0x484>
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    0x10(%rbp),%eax
jg     16bd <mrqmin+0x54d>
mov    0x29d7(%rip),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1686 <mrqmin+0x516>
jmp    16e6 <mrqmin+0x576>
mov    0x38(%rbp),%rax
movsd  0x93a(%rip),%xmm0
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
movsd  0x298a(%rip),%xmm0
mov    0x28(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x70,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

x86 assembly (clang, -O2):
<mrqmin>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r9,0x8(%rsp)
mov    %r8,%r15
mov    %ecx,0x4(%rsp)
mov    %rdx,0x40(%rsp)
mov    %rsi,0x30(%rsp)
mov    %rdi,0x28(%rsp)
mov    0xc8(%rsp),%r13
mov    0xc0(%rsp),%rax
mov    %rax,0x38(%rsp)
mov    0xb8(%rsp),%r12
mov    0xb0(%rsp),%rbx
mov    0xa8(%rsp),%rax
mov    %rax,0x10(%rsp)
mov    0xa0(%rsp),%r14d
xorpd  %xmm0,%xmm0
ucomisd 0x0(%r13),%xmm0
mov    %rbx,0x18(%rsp)
jbe    152b <mrqmin+0x3bb>
mov    $0x1,%edi
mov    %r14d,%esi
call   1080 <vector@plt>
mov    %rax,0x2e74(%rip)
mov    $0x1,%edi
mov    %r14d,%esi
call   1080 <vector@plt>
mov    %rax,0x2e68(%rip)
mov    $0x1,%edi
mov    %r14d,%esi
call   1080 <vector@plt>
mov    %rax,0x2e5c(%rip)
movl   $0x0,0x2e32(%rip)
test   %r14d,%r14d
jle    1344 <mrqmin+0x1d4>
mov    %r14d,%ebp
cmp    $0x1,%r14d
jne    13b9 <mrqmin+0x249>
xor    %esi,%esi
mov    $0x1,%eax
mov    0x8(%rsp),%rdi
test   $0x1,%bpl
je     1265 <mrqmin+0xf5>
cmpl   $0x0,(%rdi,%rax,4)
je     1265 <mrqmin+0xf5>
add    $0x1,%esi
mov    %esi,0x2dfb(%rip)
mov    $0x1,%edi
mov    $0x1,%edx
mov    $0x1,%ecx
call   1070 <matrix@plt>
mov    %rax,0x2e08(%rip)
movabs $0x3f50624dd2f1a9fc,%rax
62 50 3f
mov    %rax,0x0(%r13)
sub    $0x8,%rsp
mov    0x30(%rsp),%rdi
mov    0x38(%rsp),%rsi
mov    0x48(%rsp),%rdx
mov    0xc(%rsp),%ecx
mov    %r15,%r8
mov    0x10(%rsp),%r9
push   0x40(%rsp)
push   %r12
push   0x2dbf(%rip)
push   %rbx
push   %r14
call   1090 <mrqcof@plt>
add    $0x30,%rsp
movsd  (%r12),%xmm0
movsd  %xmm0,0x2d95(%rip)
test   %r14d,%r14d
jle    152b <mrqmin+0x3bb>
add    $0x1,%rbp
mov    0x2d89(%rip),%rax
lea    -0x1(%rbp),%rdx
mov    $0x1,%ecx
cmp    $0x4,%rdx
jb     14bc <mrqmin+0x34c>
lea    0x8(%rax),%rsi
lea    (%r15,%rbp,8),%rdi
cmp    %rdi,%rsi
jae    1318 <mrqmin+0x1a8>
lea    (%rax,%rbp,8),%rsi
lea    0x8(%r15),%rdi
cmp    %rsi,%rdi
jb     14bc <mrqmin+0x34c>
mov    %rdx,%rcx
and    $0xfffffffffffffffc,%rcx
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%ebx
and    $0x3,%ebx
cmp    $0xc,%rdi
jae    1407 <mrqmin+0x297>
xor    %edi,%edi
jmp    1476 <mrqmin+0x306>
mov    $0x1,%edi
xor    %esi,%esi
mov    $0x1,%edx
mov    $0x1,%ecx
call   1070 <matrix@plt>
mov    %rax,0x2d27(%rip)
movabs $0x3f50624dd2f1a9fc,%rax
62 50 3f
mov    %rax,0x0(%r13)
sub    $0x8,%rsp
mov    0x30(%rsp),%rdi
mov    0x38(%rsp),%rsi
mov    0x48(%rsp),%rdx
mov    0xc(%rsp),%ecx
mov    %r15,%r8
mov    0x10(%rsp),%r9
push   0x40(%rsp)
push   %r12
push   0x2cde(%rip)
push   %rbx
push   %r14
call   1090 <mrqcof@plt>
add    $0x30,%rsp
movsd  (%r12),%xmm0
movsd  %xmm0,0x2cb4(%rip)
jmp    152b <mrqmin+0x3bb>
mov    %ebp,%ecx
and    $0xfffffffe,%ecx
neg    %rcx
xor    %esi,%esi
mov    $0x1,%eax
mov    0x8(%rsp),%rdi
jmp    13e6 <mrqmin+0x276>
nop
lea    (%rcx,%rax,1),%rdx
add    $0x2,%rdx
add    $0x2,%rax
cmp    $0x1,%rdx
je     1250 <mrqmin+0xe0>
cmpl   $0x0,(%rdi,%rax,4)
je     13f5 <mrqmin+0x285>
add    $0x1,%esi
mov    %esi,0x2c6b(%rip)
cmpl   $0x0,0x4(%rdi,%rax,4)
je     13d0 <mrqmin+0x260>
add    $0x1,%esi
mov    %esi,0x2c5b(%rip)
jmp    13d0 <mrqmin+0x260>
and    $0xfffffffffffffffc,%rsi
xor    %edi,%edi
nopl   (%rax)
movups 0x8(%r15,%rdi,8),%xmm0
movups 0x18(%r15,%rdi,8),%xmm1
movups %xmm0,0x8(%rax,%rdi,8)
movups %xmm1,0x18(%rax,%rdi,8)
movups 0x28(%r15,%rdi,8),%xmm0
movups 0x38(%r15,%rdi,8),%xmm1
movups %xmm0,0x28(%rax,%rdi,8)
movups %xmm1,0x38(%rax,%rdi,8)
movups 0x48(%r15,%rdi,8),%xmm0
movups 0x58(%r15,%rdi,8),%xmm1
movups %xmm0,0x48(%rax,%rdi,8)
movups %xmm1,0x58(%rax,%rdi,8)
movupd 0x68(%r15,%rdi,8),%xmm0
movupd 0x78(%r15,%rdi,8),%xmm1
movupd %xmm0,0x68(%rax,%rdi,8)
movupd %xmm1,0x78(%rax,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1410 <mrqmin+0x2a0>
test   %rbx,%rbx
je     14ae <mrqmin+0x33e>
nopl   0x0(%rax,%rax,1)
lea    0x0(,%rdi,8),%rsi
or     $0x8,%rsi
movupd (%r15,%rsi,1),%xmm0
movupd 0x10(%r15,%rsi,1),%xmm1
movupd %xmm0,(%rax,%rsi,1)
movupd %xmm1,0x10(%rax,%rsi,1)
add    $0x4,%rdi
add    $0xffffffffffffffff,%rbx
jne    1480 <mrqmin+0x310>
cmp    %rcx,%rdx
mov    0x18(%rsp),%rbx
je     152b <mrqmin+0x3bb>
or     $0x1,%rcx
mov    %ebp,%esi
sub    %ecx,%esi
mov    %rcx,%rdx
not    %rdx
add    %rbp,%rdx
and    $0x3,%rsi
je     14e5 <mrqmin+0x375>
nop
movsd  (%r15,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    14d0 <mrqmin+0x360>
cmp    $0x3,%rdx
jb     152b <mrqmin+0x3bb>
nopl   0x0(%rax,%rax,1)
movsd  (%r15,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x8(%r15,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rax,%rcx,8)
movsd  0x10(%r15,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rax,%rcx,8)
movsd  0x18(%r15,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rax,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rbp
jne    14f0 <mrqmin+0x380>
mov    0x2b2f(%rip),%esi
test   %esi,%esi
jle    17aa <mrqmin+0x63a>
mov    %r15,0x50(%rsp)
mov    %r13,%r15
mov    %r12,0x58(%rsp)
mov    %r14,0x20(%rsp)
mov    0x2b26(%rip),%r14
mov    0x2b2f(%rip),%rdx
lea    0x1(%rsi),%rcx
lea    -0x4(%rsi),%r8
mov    %r8,0x48(%rsp)
shr    $0x2,%r8
add    $0x1,%r8
mov    %esi,%r11d
and    $0xfffffffc,%r11d
lea    0x1(%r11),%rax
mov    %rax,0x60(%rsp)
mov    %r8d,%r10d
and    $0x3,%r10d
and    $0xfffffffffffffffc,%r8
mov    $0x1,%r12d
movsd  0xa69(%rip),%xmm0
jmp    15d3 <mrqmin+0x463>
nopl   0x0(%rax)
movsd  (%r15),%xmm1
addsd  %xmm0,%xmm1
mulsd  0x0(%rbp,%r12,8),%xmm1
movsd  %xmm1,0x0(%r13,%r12,8)
movsd  (%r14,%r12,8),%xmm1
mov    (%rdx,%r12,8),%rax
movsd  %xmm1,0x8(%rax)
add    $0x1,%r12
cmp    %rcx,%r12
je     1796 <mrqmin+0x626>
mov    (%rbx,%r12,8),%rbp
mov    0x10(%rsp),%rax
mov    (%rax,%r12,8),%r13
mov    $0x1,%r9d
cmp    $0x4,%esi
jb     1700 <mrqmin+0x590>
lea    0x8(%r13),%rax
lea    0x0(,%rcx,8),%rdi
add    %rbp,%rdi
cmp    %rdi,%rax
jae    1621 <mrqmin+0x4b1>
lea    0x0(,%rcx,8),%rax
add    %r13,%rax
lea    0x8(%rbp),%rdi
mov    $0x1,%r9d
cmp    %rax,%rdi
jb     1700 <mrqmin+0x590>
cmpq   $0xc,0x48(%rsp)
jae    162d <mrqmin+0x4bd>
xor    %edi,%edi
jmp    16a4 <mrqmin+0x534>
mov    %r8,%rax
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups 0x8(%rbp,%rdi,8),%xmm1
movups 0x18(%rbp,%rdi,8),%xmm2
movups %xmm1,0x8(%r13,%rdi,8)
movups %xmm2,0x18(%r13,%rdi,8)
movups 0x28(%rbp,%rdi,8),%xmm1
movups 0x38(%rbp,%rdi,8),%xmm2
movups %xmm1,0x28(%r13,%rdi,8)
movups %xmm2,0x38(%r13,%rdi,8)
movups 0x48(%rbp,%rdi,8),%xmm1
movups 0x58(%rbp,%rdi,8),%xmm2
movups %xmm1,0x48(%r13,%rdi,8)
movups %xmm2,0x58(%r13,%rdi,8)
movupd 0x68(%rbp,%rdi,8),%xmm1
movups 0x78(%rbp,%rdi,8),%xmm2
movupd %xmm1,0x68(%r13,%rdi,8)
movups %xmm2,0x78(%r13,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rax
jne    1640 <mrqmin+0x4d0>
test   %r10,%r10
je     16de <mrqmin+0x56e>
mov    %r10,%rax
nopl   0x0(%rax)
lea    0x0(,%rdi,8),%rbx
or     $0x8,%rbx
movupd 0x0(%rbp,%rbx,1),%xmm1
movups 0x10(%rbp,%rbx,1),%xmm2
movupd %xmm1,0x0(%r13,%rbx,1)
movups %xmm2,0x10(%r13,%rbx,1)
add    $0x4,%rdi
add    $0xffffffffffffffff,%rax
jne    16b0 <mrqmin+0x540>
mov    0x60(%rsp),%r9
cmp    %r11,%rsi
mov    0x18(%rsp),%rbx
je     15a0 <mrqmin+0x430>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %ecx,%edi
sub    %r9d,%edi
mov    %r9,%rax
not    %rax
add    %rcx,%rax
and    $0x3,%rdi
je     1738 <mrqmin+0x5c8>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  0x0(%rbp,%r9,8),%xmm1
movsd  %xmm1,0x0(%r13,%r9,8)
add    $0x1,%r9
add    $0xffffffffffffffff,%rdi
jne    1720 <mrqmin+0x5b0>
cmp    $0x3,%rax
jb     15a0 <mrqmin+0x430>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%rbp,%r9,8),%xmm1
movsd  %xmm1,0x0(%r13,%r9,8)
movsd  0x8(%rbp,%r9,8),%xmm1
movsd  %xmm1,0x8(%r13,%r9,8)
movsd  0x10(%rbp,%r9,8),%xmm1
movsd  %xmm1,0x10(%r13,%r9,8)
movsd  0x18(%rbp,%r9,8),%xmm1
movsd  %xmm1,0x18(%r13,%r9,8)
add    $0x4,%r9
cmp    %r9,%rcx
jne    1750 <mrqmin+0x5e0>
jmp    15a0 <mrqmin+0x430>
mov    0x20(%rsp),%r14
mov    0x58(%rsp),%r12
mov    %r15,%r13
mov    0x50(%rsp),%r15
jmp    17b1 <mrqmin+0x641>
mov    0x28d7(%rip),%rdx
mov    0x10(%rsp),%rdi
mov    $0x1,%ecx
call   1050 <gaussj@plt>
mov    0x289a(%rip),%ecx
test   %ecx,%ecx
mov    0x8(%rsp),%r9
jle    1887 <mrqmin+0x717>
mov    0x28ae(%rip),%rdx
mov    0x289f(%rip),%rsi
lea    -0x1(%rcx),%rbp
mov    %ecx,%eax
and    $0x3,%eax
mov    $0x1,%edi
cmp    $0x3,%rbp
jb     1852 <mrqmin+0x6e2>
mov    %ecx,%ebp
and    $0xfffffffc,%ebp
xor    %edi,%edi
nopl   0x0(%rax)
mov    0x8(%rdx,%rdi,8),%rbx
movsd  0x8(%rbx),%xmm0
movsd  %xmm0,0x8(%rsi,%rdi,8)
mov    0x10(%rdx,%rdi,8),%rbx
movsd  0x8(%rbx),%xmm0
movsd  %xmm0,0x10(%rsi,%rdi,8)
mov    0x18(%rdx,%rdi,8),%rbx
movsd  0x8(%rbx),%xmm0
movsd  %xmm0,0x18(%rsi,%rdi,8)
mov    0x20(%rdx,%rdi,8),%rbx
movsd  0x8(%rbx),%xmm0
movsd  %xmm0,0x20(%rsi,%rdi,8)
add    $0x4,%rdi
cmp    %rdi,%rbp
jne    1800 <mrqmin+0x690>
add    $0x1,%rdi
mov    0x18(%rsp),%rbx
test   %rax,%rax
je     1887 <mrqmin+0x717>
lea    (%rsi,%rdi,8),%rsi
lea    (%rdx,%rdi,8),%rdx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%rdx,%rdi,8),%rbp
movsd  0x8(%rbp),%xmm0
movsd  %xmm0,(%rsi,%rdi,8)
add    $0x1,%rdi
cmp    %rdi,%rax
jne    1870 <mrqmin+0x700>
movsd  0x0(%r13),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1934 <mrqmin+0x7c4>
jp     1934 <mrqmin+0x7c4>
mov    0x10(%rsp),%rdi
mov    %r14d,%esi
mov    %r9,%rdx
mov    %r9,%rbp
call   1060 <covsrt@plt>
mov    0x27a6(%rip),%ecx
mov    %rbx,%rdi
mov    %r14d,%esi
mov    %rbp,%rdx
call   1060 <covsrt@plt>
mov    0x27b9(%rip),%rdi
mov    0x278b(%rip),%edx
mov    $0x1,%esi
mov    $0x1,%ecx
mov    $0x1,%r8d
call   1030 <free_matrix@plt>
mov    0x278f(%rip),%rdi
mov    $0x1,%esi
mov    %r14d,%edx
call   1040 <free_vector@plt>
mov    0x2773(%rip),%rdi
mov    $0x1,%esi
mov    %r14d,%edx
call   1040 <free_vector@plt>
mov    0x2757(%rip),%rdi
mov    $0x1,%esi
mov    %r14d,%edx
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_vector@plt>
test   %r14d,%r14d
jle    197e <mrqmin+0x80e>
mov    0x2740(%rip),%rax
mov    0x2729(%rip),%r8
mov    %r14d,%ecx
cmp    $0x1,%r14d
jne    1c71 <mrqmin+0xb01>
xor    %esi,%esi
mov    $0x1,%edx
test   $0x1,%cl
je     198c <mrqmin+0x81c>
cmpl   $0x0,(%r9,%rdx,4)
je     198c <mrqmin+0x81c>
movsd  (%r15,%rdx,8),%xmm0
movslq %esi,%rcx
addsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,(%r8,%rdx,8)
jmp    198c <mrqmin+0x81c>
mov    0x26eb(%rip),%r8
mov    0x26f4(%rip),%rax
sub    $0x8,%rsp
mov    0x30(%rsp),%rdi
mov    0x38(%rsp),%rsi
mov    0x48(%rsp),%rdx
mov    0xc(%rsp),%ecx
push   0x40(%rsp)
push   %r12
push   %rax
push   0x30(%rsp)
push   %r14
call   1090 <mrqcof@plt>
add    $0x30,%rsp
movsd  0x26a7(%rip),%xmm1
ucomisd (%r12),%xmm1
movsd  0x0(%r13),%xmm0
jbe    1bec <mrqmin+0xa7c>
mov    %r14,0x20(%rsp)
mulsd  0x630(%rip),%xmm0
movsd  %xmm0,0x0(%r13)
movsd  (%r12),%xmm0
movsd  %xmm0,0x2674(%rip)
mov    0x2666(%rip),%eax
test   %eax,%eax
jle    1c05 <mrqmin+0xa95>
mov    0x2677(%rip),%r12
mov    0x2668(%rip),%r13
lea    0x1(%rax),%r8
lea    -0x4(%rax),%r9
mov    %r9,%r10
shr    $0x2,%r10
add    $0x1,%r10
mov    %eax,%r11d
and    $0xfffffffc,%r11d
lea    0x1(%r11),%rcx
mov    %rcx,0x48(%rsp)
mov    %r10d,%r14d
and    $0x3,%r14d
and    $0xfffffffffffffffc,%r10
mov    $0x1,%ebp
jmp    1a6a <mrqmin+0x8fa>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r12,%rbp,8),%xmm0
movsd  %xmm0,0x0(%r13,%rbp,8)
add    $0x1,%rbp
cmp    %r8,%rbp
je     1c05 <mrqmin+0xa95>
mov    0x10(%rsp),%rcx
mov    (%rcx,%rbp,8),%rbx
mov    0x18(%rsp),%rcx
mov    (%rcx,%rbp,8),%rdi
mov    $0x1,%ecx
cmp    $0x4,%eax
jb     1b70 <mrqmin+0xa00>
lea    0x8(%rdi),%rcx
lea    (%rbx,%r8,8),%rdx
cmp    %rdx,%rcx
jae    1aad <mrqmin+0x93d>
lea    (%rdi,%r8,8),%rdx
lea    0x8(%rbx),%rsi
mov    $0x1,%ecx
cmp    %rdx,%rsi
jb     1b70 <mrqmin+0xa00>
cmp    $0xc,%r9
jae    1ab7 <mrqmin+0x947>
xor    %ecx,%ecx
jmp    1b1e <mrqmin+0x9ae>
mov    %r10,%rdx
xor    %ecx,%ecx
nopl   0x0(%rax)
movups 0x8(%rbx,%rcx,8),%xmm0
movups 0x18(%rbx,%rcx,8),%xmm1
movups %xmm0,0x8(%rdi,%rcx,8)
movups %xmm1,0x18(%rdi,%rcx,8)
movups 0x28(%rbx,%rcx,8),%xmm0
movups 0x38(%rbx,%rcx,8),%xmm1
movups %xmm0,0x28(%rdi,%rcx,8)
movups %xmm1,0x38(%rdi,%rcx,8)
movups 0x48(%rbx,%rcx,8),%xmm0
movups 0x58(%rbx,%rcx,8),%xmm1
movups %xmm0,0x48(%rdi,%rcx,8)
movups %xmm1,0x58(%rdi,%rcx,8)
movupd 0x68(%rbx,%rcx,8),%xmm0
movupd 0x78(%rbx,%rcx,8),%xmm1
movupd %xmm0,0x68(%rdi,%rcx,8)
movupd %xmm1,0x78(%rdi,%rcx,8)
add    $0x10,%rcx
add    $0xfffffffffffffffc,%rdx
jne    1ac0 <mrqmin+0x950>
test   %r14,%r14
je     1b5c <mrqmin+0x9ec>
mov    %r14,%rdx
cs nopw 0x0(%rax,%rax,1)
lea    0x0(,%rcx,8),%rsi
or     $0x8,%rsi
movupd (%rbx,%rsi,1),%xmm0
movupd 0x10(%rbx,%rsi,1),%xmm1
movupd %xmm0,(%rdi,%rsi,1)
movupd %xmm1,0x10(%rdi,%rsi,1)
add    $0x4,%rcx
add    $0xffffffffffffffff,%rdx
jne    1b30 <mrqmin+0x9c0>
mov    0x48(%rsp),%rcx
cmp    %r11,%rax
je     1a50 <mrqmin+0x8e0>
nopw   0x0(%rax,%rax,1)
mov    %r8d,%esi
sub    %ecx,%esi
mov    %rcx,%rdx
not    %rdx
add    %r8,%rdx
and    $0x3,%rsi
je     1ba4 <mrqmin+0xa34>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rbx,%rcx,8),%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1b90 <mrqmin+0xa20>
cmp    $0x3,%rdx
jb     1a50 <mrqmin+0x8e0>
xchg   %ax,%ax
movsd  (%rbx,%rcx,8),%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
movsd  0x8(%rbx,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rdi,%rcx,8)
movsd  0x10(%rbx,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rdi,%rcx,8)
movsd  0x18(%rbx,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rdi,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%r8
jne    1bb0 <mrqmin+0xa40>
jmp    1a50 <mrqmin+0x8e0>
mulsd  0x414(%rip),%xmm0
movsd  %xmm0,0x0(%r13)
movsd  %xmm1,(%r12)
jmp    1e0b <mrqmin+0xc9b>
mov    0x20(%rsp),%rdi
test   %edi,%edi
jle    1e0b <mrqmin+0xc9b>
mov    0x2457(%rip),%rax
mov    %edi,%esi
lea    0x1(%rsi),%rcx
mov    $0x1,%edx
cmp    $0x4,%edi
jb     1d97 <mrqmin+0xc27>
lea    0x8(%r15),%rdi
lea    (%rax,%rcx,8),%rbp
cmp    %rbp,%rdi
jae    1c4b <mrqmin+0xadb>
lea    (%r15,%rcx,8),%rdi
lea    0x8(%rax),%rbp
cmp    %rdi,%rbp
jb     1d97 <mrqmin+0xc27>
mov    %esi,%edx
and    $0xfffffffc,%edx
lea    -0x4(%rdx),%rbp
mov    %rbp,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%ebx
and    $0x3,%ebx
cmp    $0xc,%rbp
jae    1ce9 <mrqmin+0xb79>
xor    %ebp,%ebp
jmp    1d56 <mrqmin+0xbe6>
mov    %ecx,%edi
and    $0xfffffffe,%edi
neg    %rdi
xor    %esi,%esi
mov    $0x1,%edx
jmp    1ca6 <mrqmin+0xb36>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%rdi,%rdx,1),%rbp
add    $0x2,%rbp
add    $0x2,%rdx
cmp    $0x1,%rbp
je     195b <mrqmin+0x7eb>
cmpl   $0x0,(%r9,%rdx,4)
je     1cc5 <mrqmin+0xb55>
movsd  (%r15,%rdx,8),%xmm0
movslq %esi,%rbp
add    $0x1,%esi
addsd  0x8(%rax,%rbp,8),%xmm0
movsd  %xmm0,(%r8,%rdx,8)
cmpl   $0x0,0x4(%r9,%rdx,4)
je     1c90 <mrqmin+0xb20>
movsd  0x8(%r15,%rdx,8),%xmm0
movslq %esi,%rbp
add    $0x1,%esi
addsd  0x8(%rax,%rbp,8),%xmm0
movsd  %xmm0,0x8(%r8,%rdx,8)
jmp    1c90 <mrqmin+0xb20>
and    $0xfffffffffffffffc,%rdi
xor    %ebp,%ebp
nop
movups 0x8(%rax,%rbp,8),%xmm0
movups 0x18(%rax,%rbp,8),%xmm1
movups %xmm0,0x8(%r15,%rbp,8)
movups %xmm1,0x18(%r15,%rbp,8)
movups 0x28(%rax,%rbp,8),%xmm0
movups 0x38(%rax,%rbp,8),%xmm1
movups %xmm0,0x28(%r15,%rbp,8)
movups %xmm1,0x38(%r15,%rbp,8)
movups 0x48(%rax,%rbp,8),%xmm0
movups 0x58(%rax,%rbp,8),%xmm1
movups %xmm0,0x48(%r15,%rbp,8)
movups %xmm1,0x58(%r15,%rbp,8)
movupd 0x68(%rax,%rbp,8),%xmm0
movupd 0x78(%rax,%rbp,8),%xmm1
movupd %xmm0,0x68(%r15,%rbp,8)
movupd %xmm1,0x78(%r15,%rbp,8)
add    $0x10,%rbp
add    $0xfffffffffffffffc,%rdi
jne    1cf0 <mrqmin+0xb80>
test   %rbx,%rbx
je     1d8e <mrqmin+0xc1e>
nopl   0x0(%rax,%rax,1)
lea    0x0(,%rbp,8),%rdi
or     $0x8,%rdi
movupd (%rax,%rdi,1),%xmm0
movupd 0x10(%rax,%rdi,1),%xmm1
movupd %xmm0,(%r15,%rdi,1)
movupd %xmm1,0x10(%r15,%rdi,1)
add    $0x4,%rbp
add    $0xffffffffffffffff,%rbx
jne    1d60 <mrqmin+0xbf0>
cmp    %rdx,%rsi
je     1e0b <mrqmin+0xc9b>
or     $0x1,%rdx
mov    %ecx,%edi
sub    %edx,%edi
mov    %rdx,%rsi
not    %rsi
add    %rcx,%rsi
and    $0x3,%rdi
je     1dc5 <mrqmin+0xc55>
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm0
movsd  %xmm0,(%r15,%rdx,8)
add    $0x1,%rdx
add    $0xffffffffffffffff,%rdi
jne    1db0 <mrqmin+0xc40>
cmp    $0x3,%rsi
jb     1e0b <mrqmin+0xc9b>
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm0
movsd  %xmm0,(%r15,%rdx,8)
movsd  0x8(%rax,%rdx,8),%xmm0
movsd  %xmm0,0x8(%r15,%rdx,8)
movsd  0x10(%rax,%rdx,8),%xmm0
movsd  %xmm0,0x10(%r15,%rdx,8)
movsd  0x18(%rax,%rdx,8),%xmm0
movsd  %xmm0,0x18(%r15,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1dd0 <mrqmin+0xc60>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@mrqmin.mfit = internal unnamed_addr global i32 0
@mrqmin.ochisq = internal unnamed_addr global double 0.000000e+00
@mrqmin.atry = internal unnamed_addr global double* null
@mrqmin.beta = internal unnamed_addr global double* null
@mrqmin.da = internal unnamed_addr global double* null
@mrqmin.oneda = internal unnamed_addr global double** null
define void @mrqmin(double* %0, double* %1, double* %2, i32 %3, double* %4, i32* %5, i32 %6, double** %7, double** %8, double* %9, void (double, double*, double*, double*, i32)* %10, double* nocapture %11) {
%13 = load double, double* %11
%14 = fcmp olt double %13, 0.000000e+00
br i1 %14, label %15, label %204
15: ; preds = %12
%16 = tail call double* @vector(i32 1, i32 %6)
store double* %16, double** @mrqmin.atry
%17 = tail call double* @vector(i32 1, i32 %6)
store double* %17, double** @mrqmin.beta
%18 = tail call double* @vector(i32 1, i32 %6)
store double* %18, double** @mrqmin.da
store i32 0, i32* @mrqmin.mfit
%19 = icmp slt i32 %6, 1
br i1 %19, label %200, label %20
20: ; preds = %15
%21 = add nuw i32 %6, 1
%22 = zext i32 %21 to i64
%23 = add nsw i64 %22, -1
%24 = and i64 %23, 1
%25 = icmp eq i32 %21, 2
br i1 %25, label %50, label %26
26: ; preds = %20
%27 = and i64 %23, -2
br label %28
28: ; preds = %45, %26
%29 = phi i64 [ 1, %26 ], [ %47, %45 ]
%30 = phi i32 [ 0, %26 ], [ %46, %45 ]
%31 = phi i64 [ 0, %26 ], [ %48, %45 ]
%32 = getelementptr inbounds i32, i32* %5, i64 %29
%33 = load i32, i32* %32
%34 = icmp eq i32 %33, 0
br i1 %34, label %37, label %35
35: ; preds = %28
%36 = add nsw i32 %30, 1
store i32 %36, i32* @mrqmin.mfit
br label %37
37: ; preds = %28, %35
%38 = phi i32 [ %30, %28 ], [ %36, %35 ]
%39 = add nuw nsw i64 %29, 1
%40 = getelementptr inbounds i32, i32* %5, i64 %39
%41 = load i32, i32* %40
%42 = icmp eq i32 %41, 0
br i1 %42, label %45, label %43
43: ; preds = %37
%44 = add nsw i32 %38, 1
store i32 %44, i32* @mrqmin.mfit
br label %45
45: ; preds = %43, %37
%46 = phi i32 [ %38, %37 ], [ %44, %43 ]
%47 = add nuw nsw i64 %29, 2
%48 = add i64 %31, 2
%49 = icmp eq i64 %48, %27
br i1 %49, label %50, label %28
50: ; preds = %45, %20
%51 = phi i32 [ undef, %20 ], [ %46, %45 ]
%52 = phi i64 [ 1, %20 ], [ %47, %45 ]
%53 = phi i32 [ 0, %20 ], [ %46, %45 ]
%54 = icmp eq i64 %24, 0
br i1 %54, label %61, label %55
55: ; preds = %50
%56 = getelementptr inbounds i32, i32* %5, i64 %52
%57 = load i32, i32* %56
%58 = icmp eq i32 %57, 0
br i1 %58, label %61, label %59
59: ; preds = %55
%60 = add nsw i32 %53, 1
store i32 %60, i32* @mrqmin.mfit
br label %61
61: ; preds = %55, %59, %50
%62 = phi i32 [ %51, %50 ], [ %53, %55 ], [ %60, %59 ]
%63 = tail call double** @matrix(i32 1, i32 %62, i32 1, i32 1)
store double** %63, double*** @mrqmin.oneda
store double 1.000000e-03, double* %11
%64 = load double*, double** @mrqmin.beta
tail call void @mrqcof(double* %0, double* %1, double* %2, i32 %3, double* %4, i32* %5, i32 %6, double** %8, double* %64, double* %9, void (double, double*, double*, double*, i32)* %10)
%65 = load double, double* %9
store double %65, double* @mrqmin.ochisq
br i1 %19, label %204, label %66
66: ; preds = %61
%67 = load double*, double** @mrqmin.atry
%68 = add nuw i32 %6, 1
%69 = zext i32 %68 to i64
%70 = add nsw i64 %22, -1
%71 = icmp ult i64 %70, 4
br i1 %71, label %162, label %72
72: ; preds = %66
%73 = getelementptr double, double* %67, i64 1
%74 = getelementptr double, double* %67, i64 %22
%75 = getelementptr double, double* %4, i64 1
%76 = getelementptr double, double* %4, i64 %22
%77 = icmp ult double* %73, %76
%78 = icmp ult double* %75, %74
%79 = and i1 %77, %78
br i1 %79, label %162, label %80
80: ; preds = %72
%81 = and i64 %70, -4
%82 = or i64 %81, 1
%83 = add nsw i64 %81, -4
%84 = lshr exact i64 %83, 2
%85 = add nuw nsw i64 %84, 1
%86 = and i64 %85, 3
%87 = icmp ult i64 %83, 12
br i1 %87, label %140, label %88
88: ; preds = %80
%89 = and i64 %85, 9223372036854775804
br label %90
90: ; preds = %90, %88
%91 = phi i64 [ 0, %88 ], [ %137, %90 ]
%92 = phi i64 [ 0, %88 ], [ %138, %90 ]
%93 = or i64 %91, 1
%94 = getelementptr inbounds double, double* %4, i64 %93
%95 = bitcast double* %94 to <2 x double>*
%96 = load <2 x double>, <2 x double>* %95
%97 = getelementptr inbounds double, double* %94, i64 2
%98 = bitcast double* %97 to <2 x double>*
%99 = load <2 x double>, <2 x double>* %98
%100 = getelementptr inbounds double, double* %67, i64 %93
%101 = bitcast double* %100 to <2 x double>*
store <2 x double> %96, <2 x double>* %101
%102 = getelementptr inbounds double, double* %100, i64 2
%103 = bitcast double* %102 to <2 x double>*
store <2 x double> %99, <2 x double>* %103
%104 = or i64 %91, 5
%105 = getelementptr inbounds double, double* %4, i64 %104
%106 = bitcast double* %105 to <2 x double>*
%107 = load <2 x double>, <2 x double>* %106
%108 = getelementptr inbounds double, double* %105, i64 2
%109 = bitcast double* %108 to <2 x double>*
%110 = load <2 x double>, <2 x double>* %109
%111 = getelementptr inbounds double, double* %67, i64 %104
%112 = bitcast double* %111 to <2 x double>*
store <2 x double> %107, <2 x double>* %112
%113 = getelementptr inbounds double, double* %111, i64 2
%114 = bitcast double* %113 to <2 x double>*
store <2 x double> %110, <2 x double>* %114
%115 = or i64 %91, 9
%116 = getelementptr inbounds double, double* %4, i64 %115
%117 = bitcast double* %116 to <2 x double>*
%118 = load <2 x double>, <2 x double>* %117
%119 = getelementptr inbounds double, double* %116, i64 2
%120 = bitcast double* %119 to <2 x double>*
%121 = load <2 x double>, <2 x double>* %120
%122 = getelementptr inbounds double, double* %67, i64 %115
%123 = bitcast double* %122 to <2 x double>*
store <2 x double> %118, <2 x double>* %123
%124 = getelementptr inbounds double, double* %122, i64 2
%125 = bitcast double* %124 to <2 x double>*
store <2 x double> %121, <2 x double>* %125
%126 = or i64 %91, 13
%127 = getelementptr inbounds double, double* %4, i64 %126
%128 = bitcast double* %127 to <2 x double>*
%129 = load <2 x double>, <2 x double>* %128
%130 = getelementptr inbounds double, double* %127, i64 2
%131 = bitcast double* %130 to <2 x double>*
%132 = load <2 x double>, <2 x double>* %131
%133 = getelementptr inbounds double, double* %67, i64 %126
%134 = bitcast double* %133 to <2 x double>*
store <2 x double> %129, <2 x double>* %134
%135 = getelementptr inbounds double, double* %133, i64 2
%136 = bitcast double* %135 to <2 x double>*
store <2 x double> %132, <2 x double>* %136
%137 = add nuw i64 %91, 16
%138 = add nuw i64 %92, 4
%139 = icmp eq i64 %138, %89
br i1 %139, label %140, label %90
140: ; preds = %90, %80
%141 = phi i64 [ 0, %80 ], [ %137, %90 ]
%142 = icmp eq i64 %86, 0
br i1 %142, label %160, label %143
143: ; preds = %140, %143
%144 = phi i64 [ %157, %143 ], [ %141, %140 ]
%145 = phi i64 [ %158, %143 ], [ 0, %140 ]
%146 = or i64 %144, 1
%147 = getelementptr inbounds double, double* %4, i64 %146
%148 = bitcast double* %147 to <2 x double>*
%149 = load <2 x double>, <2 x double>* %148
%150 = getelementptr inbounds double, double* %147, i64 2
%151 = bitcast double* %150 to <2 x double>*
%152 = load <2 x double>, <2 x double>* %151
%153 = getelementptr inbounds double, double* %67, i64 %146
%154 = bitcast double* %153 to <2 x double>*
store <2 x double> %149, <2 x double>* %154
%155 = getelementptr inbounds double, double* %153, i64 2
%156 = bitcast double* %155 to <2 x double>*
store <2 x double> %152, <2 x double>* %156
%157 = add nuw i64 %144, 4
%158 = add i64 %145, 1
%159 = icmp eq i64 %158, %86
br i1 %159, label %160, label %143
160: ; preds = %143, %140
%161 = icmp eq i64 %70, %81
br i1 %161, label %204, label %162
162: ; preds = %72, %66, %160
%163 = phi i64 [ 1, %72 ], [ 1, %66 ], [ %82, %160 ]
%164 = sub nsw i64 %22, %163
%165 = xor i64 %163, -1
%166 = add nsw i64 %165, %22
%167 = and i64 %164, 3
%168 = icmp eq i64 %167, 0
br i1 %168, label %178, label %169
169: ; preds = %162, %169
%170 = phi i64 [ %175, %169 ], [ %163, %162 ]
%171 = phi i64 [ %176, %169 ], [ 0, %162 ]
%172 = getelementptr inbounds double, double* %4, i64 %170
%173 = load double, double* %172
%174 = getelementptr inbounds double, double* %67, i64 %170
store double %173, double* %174
%175 = add nuw nsw i64 %170, 1
%176 = add i64 %171, 1
%177 = icmp eq i64 %176, %167
br i1 %177, label %178, label %169
178: ; preds = %169, %162
%179 = phi i64 [ %163, %162 ], [ %175, %169 ]
%180 = icmp ult i64 %166, 3
br i1 %180, label %204, label %181
181: ; preds = %178, %181
%182 = phi i64 [ %198, %181 ], [ %179, %178 ]
%183 = getelementptr inbounds double, double* %4, i64 %182
%184 = load double, double* %183
%185 = getelementptr inbounds double, double* %67, i64 %182
store double %184, double* %185
%186 = add nuw nsw i64 %182, 1
%187 = getelementptr inbounds double, double* %4, i64 %186
%188 = load double, double* %187
%189 = getelementptr inbounds double, double* %67, i64 %186
store double %188, double* %189
%190 = add nuw nsw i64 %182, 2
%191 = getelementptr inbounds double, double* %4, i64 %190
%192 = load double, double* %191
%193 = getelementptr inbounds double, double* %67, i64 %190
store double %192, double* %193
%194 = add nuw nsw i64 %182, 3
%195 = getelementptr inbounds double, double* %4, i64 %194
%196 = load double, double* %195
%197 = getelementptr inbounds double, double* %67, i64 %194
store double %196, double* %197
%198 = add nuw nsw i64 %182, 4
%199 = icmp eq i64 %198, %69
br i1 %199, label %204, label %181
200: ; preds = %15
%201 = tail call double** @matrix(i32 1, i32 0, i32 1, i32 1)
store double** %201, double*** @mrqmin.oneda
store double 1.000000e-03, double* %11
%202 = load double*, double** @mrqmin.beta
tail call void @mrqcof(double* %0, double* %1, double* %2, i32 %3, double* %4, i32* %5, i32 %6, double** %8, double* %202, double* %9, void (double, double*, double*, double*, i32)* %10)
%203 = load double, double* %9
store double %203, double* @mrqmin.ochisq
br label %204
204: ; preds = %178, %181, %160, %200, %61, %12
%205 = load i32, i32* @mrqmin.mfit
%206 = icmp slt i32 %205, 1
br i1 %206, label %207, label %209
207: ; preds = %204
%208 = load double**, double*** @mrqmin.oneda
br label %363
209: ; preds = %204
%210 = load double*, double** @mrqmin.beta
%211 = load double**, double*** @mrqmin.oneda
%212 = add nuw i32 %205, 1
%213 = zext i32 %212 to i64
%214 = add nsw i64 %213, -1
%215 = add nsw i64 %213, -5
%216 = lshr i64 %215, 2
%217 = add nuw nsw i64 %216, 1
%218 = icmp ult i64 %214, 4
%219 = and i64 %214, -4
%220 = or i64 %219, 1
%221 = and i64 %217, 3
%222 = icmp ult i64 %215, 12
%223 = and i64 %217, 9223372036854775804
%224 = icmp eq i64 %221, 0
%225 = icmp eq i64 %214, %219
br label %226
226: ; preds = %349, %209
%227 = phi i64 [ 1, %209 ], [ %361, %349 ]
%228 = getelementptr inbounds double*, double** %8, i64 %227
%229 = load double*, double** %228
%230 = getelementptr inbounds double*, double** %7, i64 %227
%231 = load double*, double** %230
br i1 %218, label %311, label %232
232: ; preds = %226
%233 = getelementptr double, double* %231, i64 1
%234 = getelementptr double, double* %231, i64 %213
%235 = getelementptr double, double* %229, i64 1
%236 = getelementptr double, double* %229, i64 %213
%237 = icmp ult double* %233, %236
%238 = icmp ult double* %235, %234
%239 = and i1 %237, %238
br i1 %239, label %311, label %240
240: ; preds = %232
br i1 %222, label %291, label %241
241: ; preds = %240, %241
%242 = phi i64 [ %288, %241 ], [ 0, %240 ]
%243 = phi i64 [ %289, %241 ], [ 0, %240 ]
%244 = or i64 %242, 1
%245 = getelementptr inbounds double, double* %229, i64 %244
%246 = bitcast double* %245 to <2 x double>*
%247 = load <2 x double>, <2 x double>* %246
%248 = getelementptr inbounds double, double* %245, i64 2
%249 = bitcast double* %248 to <2 x double>*
%250 = load <2 x double>, <2 x double>* %249
%251 = getelementptr inbounds double, double* %231, i64 %244
%252 = bitcast double* %251 to <2 x double>*
store <2 x double> %247, <2 x double>* %252
%253 = getelementptr inbounds double, double* %251, i64 2
%254 = bitcast double* %253 to <2 x double>*
store <2 x double> %250, <2 x double>* %254
%255 = or i64 %242, 5
%256 = getelementptr inbounds double, double* %229, i64 %255
%257 = bitcast double* %256 to <2 x double>*
%258 = load <2 x double>, <2 x double>* %257
%259 = getelementptr inbounds double, double* %256, i64 2
%260 = bitcast double* %259 to <2 x double>*
%261 = load <2 x double>, <2 x double>* %260
%262 = getelementptr inbounds double, double* %231, i64 %255
%263 = bitcast double* %262 to <2 x double>*
store <2 x double> %258, <2 x double>* %263
%264 = getelementptr inbounds double, double* %262, i64 2
%265 = bitcast double* %264 to <2 x double>*
store <2 x double> %261, <2 x double>* %265
%266 = or i64 %242, 9
%267 = getelementptr inbounds double, double* %229, i64 %266
%268 = bitcast double* %267 to <2 x double>*
%269 = load <2 x double>, <2 x double>* %268
%270 = getelementptr inbounds double, double* %267, i64 2
%271 = bitcast double* %270 to <2 x double>*
%272 = load <2 x double>, <2 x double>* %271
%273 = getelementptr inbounds double, double* %231, i64 %266
%274 = bitcast double* %273 to <2 x double>*
store <2 x double> %269, <2 x double>* %274
%275 = getelementptr inbounds double, double* %273, i64 2
%276 = bitcast double* %275 to <2 x double>*
store <2 x double> %272, <2 x double>* %276
%277 = or i64 %242, 13
%278 = getelementptr inbounds double, double* %229, i64 %277
%279 = bitcast double* %278 to <2 x double>*
%280 = load <2 x double>, <2 x double>* %279
%281 = getelementptr inbounds double, double* %278, i64 2
%282 = bitcast double* %281 to <2 x double>*
%283 = load <2 x double>, <2 x double>* %282
%284 = getelementptr inbounds double, double* %231, i64 %277
%285 = bitcast double* %284 to <2 x double>*
store <2 x double> %280, <2 x double>* %285
%286 = getelementptr inbounds double, double* %284, i64 2
%287 = bitcast double* %286 to <2 x double>*
store <2 x double> %283, <2 x double>* %287
%288 = add nuw i64 %242, 16
%289 = add nuw i64 %243, 4
%290 = icmp eq i64 %289, %223
br i1 %290, label %291, label %241
291: ; preds = %241, %240
%292 = phi i64 [ 0, %240 ], [ %288, %241 ]
br i1 %224, label %310, label %293
293: ; preds = %291, %293
%294 = phi i64 [ %307, %293 ], [ %292, %291 ]
%295 = phi i64 [ %308, %293 ], [ 0, %291 ]
%296 = or i64 %294, 1
%297 = getelementptr inbounds double, double* %229, i64 %296
%298 = bitcast double* %297 to <2 x double>*
%299 = load <2 x double>, <2 x double>* %298
%300 = getelementptr inbounds double, double* %297, i64 2
%301 = bitcast double* %300 to <2 x double>*
%302 = load <2 x double>, <2 x double>* %301
%303 = getelementptr inbounds double, double* %231, i64 %296
%304 = bitcast double* %303 to <2 x double>*
store <2 x double> %299, <2 x double>* %304
%305 = getelementptr inbounds double, double* %303, i64 2
%306 = bitcast double* %305 to <2 x double>*
store <2 x double> %302, <2 x double>* %306
%307 = add nuw i64 %294, 4
%308 = add i64 %295, 1
%309 = icmp eq i64 %308, %221
br i1 %309, label %310, label %293
310: ; preds = %293, %291
br i1 %225, label %349, label %311
311: ; preds = %232, %226, %310
%312 = phi i64 [ 1, %232 ], [ 1, %226 ], [ %220, %310 ]
%313 = sub nsw i64 %213, %312
%314 = xor i64 %312, -1
%315 = add nsw i64 %314, %213
%316 = and i64 %313, 3
%317 = icmp eq i64 %316, 0
br i1 %317, label %327, label %318
318: ; preds = %311, %318
%319 = phi i64 [ %324, %318 ], [ %312, %311 ]
%320 = phi i64 [ %325, %318 ], [ 0, %311 ]
%321 = getelementptr inbounds double, double* %229, i64 %319
%322 = load double, double* %321
%323 = getelementptr inbounds double, double* %231, i64 %319
store double %322, double* %323
%324 = add nuw nsw i64 %319, 1
%325 = add i64 %320, 1
%326 = icmp eq i64 %325, %316
br i1 %326, label %327, label %318
327: ; preds = %318, %311
%328 = phi i64 [ %312, %311 ], [ %324, %318 ]
%329 = icmp ult i64 %315, 3
br i1 %329, label %349, label %330
330: ; preds = %327, %330
%331 = phi i64 [ %347, %330 ], [ %328, %327 ]
%332 = getelementptr inbounds double, double* %229, i64 %331
%333 = load double, double* %332
%334 = getelementptr inbounds double, double* %231, i64 %331
store double %333, double* %334
%335 = add nuw nsw i64 %331, 1
%336 = getelementptr inbounds double, double* %229, i64 %335
%337 = load double, double* %336
%338 = getelementptr inbounds double, double* %231, i64 %335
store double %337, double* %338
%339 = add nuw nsw i64 %331, 2
%340 = getelementptr inbounds double, double* %229, i64 %339
%341 = load double, double* %340
%342 = getelementptr inbounds double, double* %231, i64 %339
store double %341, double* %342
%343 = add nuw nsw i64 %331, 3
%344 = getelementptr inbounds double, double* %229, i64 %343
%345 = load double, double* %344
%346 = getelementptr inbounds double, double* %231, i64 %343
store double %345, double* %346
%347 = add nuw nsw i64 %331, 4
%348 = icmp eq i64 %347, %213
br i1 %348, label %349, label %330
349: ; preds = %327, %330, %310
%350 = getelementptr inbounds double, double* %229, i64 %227
%351 = load double, double* %350
%352 = load double, double* %11
%353 = fadd double %352, 1.000000e+00
%354 = fmul double %351, %353
%355 = getelementptr inbounds double, double* %231, i64 %227
store double %354, double* %355
%356 = getelementptr inbounds double, double* %210, i64 %227
%357 = load double, double* %356
%358 = getelementptr inbounds double*, double** %211, i64 %227
%359 = load double*, double** %358
%360 = getelementptr inbounds double, double* %359, i64 1
store double %357, double* %360
%361 = add nuw nsw i64 %227, 1
%362 = icmp eq i64 %361, %213
br i1 %362, label %363, label %226
363: ; preds = %349, %207
%364 = phi double** [ %208, %207 ], [ %211, %349 ]
tail call void @gaussj(double** %7, i32 %205, double** %364, i32 1)
%365 = load i32, i32* @mrqmin.mfit
%366 = icmp slt i32 %365, 1
br i1 %366, label %421, label %367
367: ; preds = %363
%368 = load double**, double*** @mrqmin.oneda
%369 = load double*, double** @mrqmin.da
%370 = add nuw i32 %365, 1
%371 = zext i32 %370 to i64
%372 = add nsw i64 %371, -1
%373 = add nsw i64 %371, -2
%374 = and i64 %372, 3
%375 = icmp ult i64 %373, 3
br i1 %375, label %407, label %376
376: ; preds = %367
%377 = and i64 %372, -4
br label %378
378: ; preds = %378, %376
%379 = phi i64 [ 1, %376 ], [ %404, %378 ]
%380 = phi i64 [ 0, %376 ], [ %405, %378 ]
%381 = getelementptr inbounds double*, double** %368, i64 %379
%382 = load double*, double** %381
%383 = getelementptr inbounds double, double* %382, i64 1
%384 = load double, double* %383
%385 = getelementptr inbounds double, double* %369, i64 %379
store double %384, double* %385
%386 = add nuw nsw i64 %379, 1
%387 = getelementptr inbounds double*, double** %368, i64 %386
%388 = load double*, double** %387
%389 = getelementptr inbounds double, double* %388, i64 1
%390 = load double, double* %389
%391 = getelementptr inbounds double, double* %369, i64 %386
store double %390, double* %391
%392 = add nuw nsw i64 %379, 2
%393 = getelementptr inbounds double*, double** %368, i64 %392
%394 = load double*, double** %393
%395 = getelementptr inbounds double, double* %394, i64 1
%396 = load double, double* %395
%397 = getelementptr inbounds double, double* %369, i64 %392
store double %396, double* %397
%398 = add nuw nsw i64 %379, 3
%399 = getelementptr inbounds double*, double** %368, i64 %398
%400 = load double*, double** %399
%401 = getelementptr inbounds double, double* %400, i64 1
%402 = load double, double* %401
%403 = getelementptr inbounds double, double* %369, i64 %398
store double %402, double* %403
%404 = add nuw nsw i64 %379, 4
%405 = add i64 %380, 4
%406 = icmp eq i64 %405, %377
br i1 %406, label %407, label %378
407: ; preds = %378, %367
%408 = phi i64 [ 1, %367 ], [ %404, %378 ]
%409 = icmp eq i64 %374, 0
br i1 %409, label %421, label %410
410: ; preds = %407, %410
%411 = phi i64 [ %418, %410 ], [ %408, %407 ]
%412 = phi i64 [ %419, %410 ], [ 0, %407 ]
%413 = getelementptr inbounds double*, double** %368, i64 %411
%414 = load double*, double** %413
%415 = getelementptr inbounds double, double* %414, i64 1
%416 = load double, double* %415
%417 = getelementptr inbounds double, double* %369, i64 %411
store double %416, double* %417
%418 = add nuw nsw i64 %411, 1
%419 = add i64 %412, 1
%420 = icmp eq i64 %419, %374
br i1 %420, label %421, label %410
421: ; preds = %407, %410, %363
%422 = load double, double* %11
%423 = fcmp oeq double %422, 0.000000e+00
br i1 %423, label %437, label %424
424: ; preds = %421
%425 = icmp slt i32 %6, 1
br i1 %425, label %426, label %429
426: ; preds = %424
%427 = load double*, double** @mrqmin.atry
%428 = load double*, double** @mrqmin.da
br label %501
429: ; preds = %424
%430 = load double*, double** @mrqmin.da
%431 = load double*, double** @mrqmin.atry
%432 = zext i32 %6 to i64
%433 = and i64 %432, 1
%434 = icmp eq i32 %6, 1
br i1 %434, label %484, label %435
435: ; preds = %429
%436 = and i64 %432, 4294967294
br label %448
437: ; preds = %421
tail call void @covsrt(double** %7, i32 %6, i32* %5, i32 %365)
%438 = load i32, i32* @mrqmin.mfit
tail call void @covsrt(double** %8, i32 %6, i32* %5, i32 %438)
%439 = load double**, double*** @mrqmin.oneda
%440 = load i32, i32* @mrqmin.mfit
%441 = tail call i32 @free_matrix(double** %439, i32 1, i32 %440, i32 1, i32 1)
%442 = load double*, double** @mrqmin.da
%443 = tail call i32 @free_vector(double* %442, i32 1, i32 %6)
%444 = load double*, double** @mrqmin.beta
%445 = tail call i32 @free_vector(double* %444, i32 1, i32 %6)
%446 = load double*, double** @mrqmin.atry
%447 = tail call i32 @free_vector(double* %446, i32 1, i32 %6)
br label %796
448: ; preds = %479, %435
%449 = phi i64 [ 1, %435 ], [ %481, %479 ]
%450 = phi i32 [ 0, %435 ], [ %480, %479 ]
%451 = phi i64 [ 0, %435 ], [ %482, %479 ]
%452 = getelementptr inbounds i32, i32* %5, i64 %449
%453 = load i32, i32* %452
%454 = icmp eq i32 %453, 0
br i1 %454, label %464, label %455
455: ; preds = %448
%456 = getelementptr inbounds double, double* %4, i64 %449
%457 = load double, double* %456
%458 = add nsw i32 %450, 1
%459 = sext i32 %458 to i64
%460 = getelementptr inbounds double, double* %430, i64 %459
%461 = load double, double* %460
%462 = fadd double %457, %461
%463 = getelementptr inbounds double, double* %431, i64 %449
store double %462, double* %463
br label %464
464: ; preds = %448, %455
%465 = phi i32 [ %458, %455 ], [ %450, %448 ]
%466 = add nuw nsw i64 %449, 1
%467 = getelementptr inbounds i32, i32* %5, i64 %466
%468 = load i32, i32* %467
%469 = icmp eq i32 %468, 0
br i1 %469, label %479, label %470
470: ; preds = %464
%471 = getelementptr inbounds double, double* %4, i64 %466
%472 = load double, double* %471
%473 = add nsw i32 %465, 1
%474 = sext i32 %473 to i64
%475 = getelementptr inbounds double, double* %430, i64 %474
%476 = load double, double* %475
%477 = fadd double %472, %476
%478 = getelementptr inbounds double, double* %431, i64 %466
store double %477, double* %478
br label %479
479: ; preds = %470, %464
%480 = phi i32 [ %473, %470 ], [ %465, %464 ]
%481 = add nuw nsw i64 %449, 2
%482 = add i64 %451, 2
%483 = icmp eq i64 %482, %436
br i1 %483, label %484, label %448
484: ; preds = %479, %429
%485 = phi i64 [ 1, %429 ], [ %481, %479 ]
%486 = phi i32 [ 0, %429 ], [ %480, %479 ]
%487 = icmp eq i64 %433, 0
br i1 %487, label %501, label %488
488: ; preds = %484
%489 = getelementptr inbounds i32, i32* %5, i64 %485
%490 = load i32, i32* %489
%491 = icmp eq i32 %490, 0
br i1 %491, label %501, label %492
492: ; preds = %488
%493 = getelementptr inbounds double, double* %4, i64 %485
%494 = load double, double* %493
%495 = add nsw i32 %486, 1
%496 = sext i32 %495 to i64
%497 = getelementptr inbounds double, double* %430, i64 %496
%498 = load double, double* %497
%499 = fadd double %494, %498
%500 = getelementptr inbounds double, double* %431, i64 %485
store double %499, double* %500
br label %501
501: ; preds = %484, %492, %488, %426
%502 = phi double* [ %428, %426 ], [ %430, %488 ], [ %430, %492 ], [ %430, %484 ]
%503 = phi double* [ %427, %426 ], [ %431, %488 ], [ %431, %492 ], [ %431, %484 ]
tail call void @mrqcof(double* %0, double* %1, double* %2, i32 %3, double* %503, i32* %5, i32 %6, double** %7, double* %502, double* %9, void (double, double*, double*, double*, i32)* %10)
%504 = load double, double* %9
%505 = load double, double* @mrqmin.ochisq
%506 = fcmp olt double %504, %505
%507 = load double, double* %11
br i1 %506, label %508, label %794
508: ; preds = %501
%509 = fmul double %507, 1.000000e-01
store double %509, double* %11
%510 = load double, double* %9
store double %510, double* @mrqmin.ochisq
%511 = load i32, i32* @mrqmin.mfit
%512 = icmp slt i32 %511, 1
br i1 %512, label %634, label %513
513: ; preds = %508
%514 = load double*, double** @mrqmin.da
%515 = load double*, double** @mrqmin.beta
%516 = add nuw i32 %511, 1
%517 = zext i32 %516 to i64
%518 = add nsw i64 %517, -1
%519 = add nsw i64 %517, -5
%520 = lshr i64 %519, 2
%521 = add nuw nsw i64 %520, 1
%522 = icmp ult i64 %518, 4
%523 = and i64 %518, -4
%524 = or i64 %523, 1
%525 = and i64 %521, 3
%526 = icmp ult i64 %519, 12
%527 = and i64 %521, 9223372036854775804
%528 = icmp eq i64 %525, 0
%529 = icmp eq i64 %518, %523
br label %530
530: ; preds = %769, %513
%531 = phi i64 [ 1, %513 ], [ %773, %769 ]
%532 = getelementptr inbounds double*, double** %7, i64 %531
%533 = load double*, double** %532
%534 = getelementptr inbounds double*, double** %8, i64 %531
%535 = load double*, double** %534
br i1 %522, label %615, label %536
536: ; preds = %530
%537 = getelementptr double, double* %535, i64 1
%538 = getelementptr double, double* %535, i64 %517
%539 = getelementptr double, double* %533, i64 1
%540 = getelementptr double, double* %533, i64 %517
%541 = icmp ult double* %537, %540
%542 = icmp ult double* %539, %538
%543 = and i1 %541, %542
br i1 %543, label %615, label %544
544: ; preds = %536
br i1 %526, label %595, label %545
545: ; preds = %544, %545
%546 = phi i64 [ %592, %545 ], [ 0, %544 ]
%547 = phi i64 [ %593, %545 ], [ 0, %544 ]
%548 = or i64 %546, 1
%549 = getelementptr inbounds double, double* %533, i64 %548
%550 = bitcast double* %549 to <2 x double>*
%551 = load <2 x double>, <2 x double>* %550
%552 = getelementptr inbounds double, double* %549, i64 2
%553 = bitcast double* %552 to <2 x double>*
%554 = load <2 x double>, <2 x double>* %553
%555 = getelementptr inbounds double, double* %535, i64 %548
%556 = bitcast double* %555 to <2 x double>*
store <2 x double> %551, <2 x double>* %556
%557 = getelementptr inbounds double, double* %555, i64 2
%558 = bitcast double* %557 to <2 x double>*
store <2 x double> %554, <2 x double>* %558
%559 = or i64 %546, 5
%560 = getelementptr inbounds double, double* %533, i64 %559
%561 = bitcast double* %560 to <2 x double>*
%562 = load <2 x double>, <2 x double>* %561
%563 = getelementptr inbounds double, double* %560, i64 2
%564 = bitcast double* %563 to <2 x double>*
%565 = load <2 x double>, <2 x double>* %564
%566 = getelementptr inbounds double, double* %535, i64 %559
%567 = bitcast double* %566 to <2 x double>*
store <2 x double> %562, <2 x double>* %567
%568 = getelementptr inbounds double, double* %566, i64 2
%569 = bitcast double* %568 to <2 x double>*
store <2 x double> %565, <2 x double>* %569
%570 = or i64 %546, 9
%571 = getelementptr inbounds double, double* %533, i64 %570
%572 = bitcast double* %571 to <2 x double>*
%573 = load <2 x double>, <2 x double>* %572
%574 = getelementptr inbounds double, double* %571, i64 2
%575 = bitcast double* %574 to <2 x double>*
%576 = load <2 x double>, <2 x double>* %575
%577 = getelementptr inbounds double, double* %535, i64 %570
%578 = bitcast double* %577 to <2 x double>*
store <2 x double> %573, <2 x double>* %578
%579 = getelementptr inbounds double, double* %577, i64 2
%580 = bitcast double* %579 to <2 x double>*
store <2 x double> %576, <2 x double>* %580
%581 = or i64 %546, 13
%582 = getelementptr inbounds double, double* %533, i64 %581
%583 = bitcast double* %582 to <2 x double>*
%584 = load <2 x double>, <2 x double>* %583
%585 = getelementptr inbounds double, double* %582, i64 2
%586 = bitcast double* %585 to <2 x double>*
%587 = load <2 x double>, <2 x double>* %586
%588 = getelementptr inbounds double, double* %535, i64 %581
%589 = bitcast double* %588 to <2 x double>*
store <2 x double> %584, <2 x double>* %589
%590 = getelementptr inbounds double, double* %588, i64 2
%591 = bitcast double* %590 to <2 x double>*
store <2 x double> %587, <2 x double>* %591
%592 = add nuw i64 %546, 16
%593 = add nuw i64 %547, 4
%594 = icmp eq i64 %593, %527
br i1 %594, label %595, label %545
595: ; preds = %545, %544
%596 = phi i64 [ 0, %544 ], [ %592, %545 ]
br i1 %528, label %614, label %597
597: ; preds = %595, %597
%598 = phi i64 [ %611, %597 ], [ %596, %595 ]
%599 = phi i64 [ %612, %597 ], [ 0, %595 ]
%600 = or i64 %598, 1
%601 = getelementptr inbounds double, double* %533, i64 %600
%602 = bitcast double* %601 to <2 x double>*
%603 = load <2 x double>, <2 x double>* %602
%604 = getelementptr inbounds double, double* %601, i64 2
%605 = bitcast double* %604 to <2 x double>*
%606 = load <2 x double>, <2 x double>* %605
%607 = getelementptr inbounds double, double* %535, i64 %600
%608 = bitcast double* %607 to <2 x double>*
store <2 x double> %603, <2 x double>* %608
%609 = getelementptr inbounds double, double* %607, i64 2
%610 = bitcast double* %609 to <2 x double>*
store <2 x double> %606, <2 x double>* %610
%611 = add nuw i64 %598, 4
%612 = add i64 %599, 1
%613 = icmp eq i64 %612, %525
br i1 %613, label %614, label %597
614: ; preds = %597, %595
br i1 %529, label %769, label %615
615: ; preds = %536, %530, %614
%616 = phi i64 [ 1, %536 ], [ 1, %530 ], [ %524, %614 ]
%617 = sub nsw i64 %517, %616
%618 = xor i64 %616, -1
%619 = add nsw i64 %618, %517
%620 = and i64 %617, 3
%621 = icmp eq i64 %620, 0
br i1 %621, label %631, label %622
622: ; preds = %615, %622
%623 = phi i64 [ %628, %622 ], [ %616, %615 ]
%624 = phi i64 [ %629, %622 ], [ 0, %615 ]
%625 = getelementptr inbounds double, double* %533, i64 %623
%626 = load double, double* %625
%627 = getelementptr inbounds double, double* %535, i64 %623
store double %626, double* %627
%628 = add nuw nsw i64 %623, 1
%629 = add i64 %624, 1
%630 = icmp eq i64 %629, %620
br i1 %630, label %631, label %622
631: ; preds = %622, %615
%632 = phi i64 [ %616, %615 ], [ %628, %622 ]
%633 = icmp ult i64 %619, 3
br i1 %633, label %769, label %750
634: ; preds = %769, %508
br i1 %425, label %796, label %635
635: ; preds = %634
%636 = load double*, double** @mrqmin.atry
%637 = add nuw i32 %6, 1
%638 = zext i32 %637 to i64
%639 = add nsw i64 %638, -1
%640 = icmp ult i64 %639, 4
br i1 %640, label %731, label %641
641: ; preds = %635
%642 = getelementptr double, double* %4, i64 1
%643 = getelementptr double, double* %4, i64 %638
%644 = getelementptr double, double* %636, i64 1
%645 = getelementptr double, double* %636, i64 %638
%646 = icmp ult double* %642, %645
%647 = icmp ult double* %644, %643
%648 = and i1 %646, %647
br i1 %648, label %731, label %649
649: ; preds = %641
%650 = and i64 %639, -4
%651 = or i64 %650, 1
%652 = add nsw i64 %650, -4
%653 = lshr exact i64 %652, 2
%654 = add nuw nsw i64 %653, 1
%655 = and i64 %654, 3
%656 = icmp ult i64 %652, 12
br i1 %656, label %709, label %657
657: ; preds = %649
%658 = and i64 %654, 9223372036854775804
br label %659
659: ; preds = %659, %657
%660 = phi i64 [ 0, %657 ], [ %706, %659 ]
%661 = phi i64 [ 0, %657 ], [ %707, %659 ]
%662 = or i64 %660, 1
%663 = getelementptr inbounds double, double* %636, i64 %662
%664 = bitcast double* %663 to <2 x double>*
%665 = load <2 x double>, <2 x double>* %664
%666 = getelementptr inbounds double, double* %663, i64 2
%667 = bitcast double* %666 to <2 x double>*
%668 = load <2 x double>, <2 x double>* %667
%669 = getelementptr inbounds double, double* %4, i64 %662
%670 = bitcast double* %669 to <2 x double>*
store <2 x double> %665, <2 x double>* %670
%671 = getelementptr inbounds double, double* %669, i64 2
%672 = bitcast double* %671 to <2 x double>*
store <2 x double> %668, <2 x double>* %672
%673 = or i64 %660, 5
%674 = getelementptr inbounds double, double* %636, i64 %673
%675 = bitcast double* %674 to <2 x double>*
%676 = load <2 x double>, <2 x double>* %675
%677 = getelementptr inbounds double, double* %674, i64 2
%678 = bitcast double* %677 to <2 x double>*
%679 = load <2 x double>, <2 x double>* %678
%680 = getelementptr inbounds double, double* %4, i64 %673
%681 = bitcast double* %680 to <2 x double>*
store <2 x double> %676, <2 x double>* %681
%682 = getelementptr inbounds double, double* %680, i64 2
%683 = bitcast double* %682 to <2 x double>*
store <2 x double> %679, <2 x double>* %683
%684 = or i64 %660, 9
%685 = getelementptr inbounds double, double* %636, i64 %684
%686 = bitcast double* %685 to <2 x double>*
%687 = load <2 x double>, <2 x double>* %686
%688 = getelementptr inbounds double, double* %685, i64 2
%689 = bitcast double* %688 to <2 x double>*
%690 = load <2 x double>, <2 x double>* %689
%691 = getelementptr inbounds double, double* %4, i64 %684
%692 = bitcast double* %691 to <2 x double>*
store <2 x double> %687, <2 x double>* %692
%693 = getelementptr inbounds double, double* %691, i64 2
%694 = bitcast double* %693 to <2 x double>*
store <2 x double> %690, <2 x double>* %694
%695 = or i64 %660, 13
%696 = getelementptr inbounds double, double* %636, i64 %695
%697 = bitcast double* %696 to <2 x double>*
%698 = load <2 x double>, <2 x double>* %697
%699 = getelementptr inbounds double, double* %696, i64 2
%700 = bitcast double* %699 to <2 x double>*
%701 = load <2 x double>, <2 x double>* %700
%702 = getelementptr inbounds double, double* %4, i64 %695
%703 = bitcast double* %702 to <2 x double>*
store <2 x double> %698, <2 x double>* %703
%704 = getelementptr inbounds double, double* %702, i64 2
%705 = bitcast double* %704 to <2 x double>*
store <2 x double> %701, <2 x double>* %705
%706 = add nuw i64 %660, 16
%707 = add nuw i64 %661, 4
%708 = icmp eq i64 %707, %658
br i1 %708, label %709, label %659
709: ; preds = %659, %649
%710 = phi i64 [ 0, %649 ], [ %706, %659 ]
%711 = icmp eq i64 %655, 0
br i1 %711, label %729, label %712
712: ; preds = %709, %712
%713 = phi i64 [ %726, %712 ], [ %710, %709 ]
%714 = phi i64 [ %727, %712 ], [ 0, %709 ]
%715 = or i64 %713, 1
%716 = getelementptr inbounds double, double* %636, i64 %715
%717 = bitcast double* %716 to <2 x double>*
%718 = load <2 x double>, <2 x double>* %717
%719 = getelementptr inbounds double, double* %716, i64 2
%720 = bitcast double* %719 to <2 x double>*
%721 = load <2 x double>, <2 x double>* %720
%722 = getelementptr inbounds double, double* %4, i64 %715
%723 = bitcast double* %722 to <2 x double>*
store <2 x double> %718, <2 x double>* %723
%724 = getelementptr inbounds double, double* %722, i64 2
%725 = bitcast double* %724 to <2 x double>*
store <2 x double> %721, <2 x double>* %725
%726 = add nuw i64 %713, 4
%727 = add i64 %714, 1
%728 = icmp eq i64 %727, %655
br i1 %728, label %729, label %712
729: ; preds = %712, %709
%730 = icmp eq i64 %639, %650
br i1 %730, label %796, label %731
731: ; preds = %641, %635, %729
%732 = phi i64 [ 1, %641 ], [ 1, %635 ], [ %651, %729 ]
%733 = sub nsw i64 %638, %732
%734 = xor i64 %732, -1
%735 = add nsw i64 %734, %638
%736 = and i64 %733, 3
%737 = icmp eq i64 %736, 0
br i1 %737, label %747, label %738
738: ; preds = %731, %738
%739 = phi i64 [ %744, %738 ], [ %732, %731 ]
%740 = phi i64 [ %745, %738 ], [ 0, %731 ]
%741 = getelementptr inbounds double, double* %636, i64 %739
%742 = load double, double* %741
%743 = getelementptr inbounds double, double* %4, i64 %739
store double %742, double* %743
%744 = add nuw nsw i64 %739, 1
%745 = add i64 %740, 1
%746 = icmp eq i64 %745, %736
br i1 %746, label %747, label %738
747: ; preds = %738, %731
%748 = phi i64 [ %732, %731 ], [ %744, %738 ]
%749 = icmp ult i64 %735, 3
br i1 %749, label %796, label %775
750: ; preds = %631, %750
%751 = phi i64 [ %767, %750 ], [ %632, %631 ]
%752 = getelementptr inbounds double, double* %533, i64 %751
%753 = load double, double* %752
%754 = getelementptr inbounds double, double* %535, i64 %751
store double %753, double* %754
%755 = add nuw nsw i64 %751, 1
%756 = getelementptr inbounds double, double* %533, i64 %755
%757 = load double, double* %756
%758 = getelementptr inbounds double, double* %535, i64 %755
store double %757, double* %758
%759 = add nuw nsw i64 %751, 2
%760 = getelementptr inbounds double, double* %533, i64 %759
%761 = load double, double* %760
%762 = getelementptr inbounds double, double* %535, i64 %759
store double %761, double* %762
%763 = add nuw nsw i64 %751, 3
%764 = getelementptr inbounds double, double* %533, i64 %763
%765 = load double, double* %764
%766 = getelementptr inbounds double, double* %535, i64 %763
store double %765, double* %766
%767 = add nuw nsw i64 %751, 4
%768 = icmp eq i64 %767, %517
br i1 %768, label %769, label %750
769: ; preds = %631, %750, %614
%770 = getelementptr inbounds double, double* %514, i64 %531
%771 = load double, double* %770
%772 = getelementptr inbounds double, double* %515, i64 %531
store double %771, double* %772
%773 = add nuw nsw i64 %531, 1
%774 = icmp eq i64 %773, %517
br i1 %774, label %634, label %530
775: ; preds = %747, %775
%776 = phi i64 [ %792, %775 ], [ %748, %747 ]
%777 = getelementptr inbounds double, double* %636, i64 %776
%778 = load double, double* %777
%779 = getelementptr inbounds double, double* %4, i64 %776
store double %778, double* %779
%780 = add nuw nsw i64 %776, 1
%781 = getelementptr inbounds double, double* %636, i64 %780
%782 = load double, double* %781
%783 = getelementptr inbounds double, double* %4, i64 %780
store double %782, double* %783
%784 = add nuw nsw i64 %776, 2
%785 = getelementptr inbounds double, double* %636, i64 %784
%786 = load double, double* %785
%787 = getelementptr inbounds double, double* %4, i64 %784
store double %786, double* %787
%788 = add nuw nsw i64 %776, 3
%789 = getelementptr inbounds double, double* %636, i64 %788
%790 = load double, double* %789
%791 = getelementptr inbounds double, double* %4, i64 %788
store double %790, double* %791
%792 = add nuw nsw i64 %776, 4
%793 = icmp eq i64 %792, %638
br i1 %793, label %796, label %775
794: ; preds = %501
%795 = fmul double %507, 1.000000e+01
store double %795, double* %11
store double %505, double* %9
br label %796
796: ; preds = %747, %775, %729, %634, %794, %437
ret void
}