int emrqmin(x, y, sig, ndata, a, ma, lista, mfit, covar, alpha, chisq, funcs,
            alamda)
float x[], y[], sig[], a[], **covar, **alpha, *chisq, *alamda;
int ndata, ma, lista[], mfit;
void (*funcs)();
{
  int k, kk, j, ihit;
  static float *da, *atry, **oneda, *beta, ochisq;
  float *vector(), **matrix();
  float newa;
  void mrqcof(), covsrt(), nrerror(), free_matrix(), free_vector();
  int gaussj();
  if (*alamda < 0.0) {
    oneda = matrix(1, mfit, 1, 1);
    atry = vector(1, ma);
    da = vector(1, ma);
    beta = vector(1, ma);
    kk = mfit + 1;
    for (j = 1; j <= ma; j++) {
      ihit = 0;
      for (k = 1; k <= mfit; k++)
        if (lista[k] == j)
          ihit++;
      if (ihit == 0) {
        lista[kk++] = j;
      } else if (ihit > 1) {
        nrerror("Bad LISTA permutation in MRQMIN-1");
        return (0);
      }
    }
    if (kk != ma + 1) {
      nrerror("Bad LISTA permutation in MRQMIN-2");
      return (0);
    }
    *alamda = 0.001;
    mrqcof(x, y, sig, ndata, a, ma, lista, mfit, alpha, beta, chisq, funcs);
    ochisq = (*chisq);
  }
  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];
  }
  if (gaussj(covar, mfit, oneda, 1) == 0)
    return (0);
  for (j = 1; j <= mfit; j++)
    da[j] = oneda[j][1];
  if (*alamda == 0.0) {
    covsrt(covar, ma, lista, mfit);
    free_vector(beta, 1, ma);
    free_vector(da, 1, ma);
    free_vector(atry, 1, ma);
    free_matrix(oneda, 1, mfit, 1, 1);
    return (1);
  }
  for (j = 1; j <= mfit; j++) {
    newa = a[lista[j]] + da[j];
    if ((j > 1) & (j % 2) & (newa <= 0.0))
      atry[lista[j]] = atry[lista[j]] / 10.0;
    else
      atry[lista[j]] = newa;
  }
  mrqcof(x, y, sig, ndata, atry, ma, lista, mfit, 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];
      a[lista[j]] = atry[lista[j]];
    }
  } else {
    *alamda *= 10.0;
    *chisq = ochisq;
  }
  return (1);
}


x86 assembly (clang, -O0):
<emrqmin>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %rbx
sub    $0x88,%rsp
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %rdi,-0x28(%rbp)
mov    %rsi,-0x30(%rbp)
mov    %rdx,-0x38(%rbp)
mov    %ecx,-0x3c(%rbp)
mov    %r8,-0x48(%rbp)
mov    %r9d,-0x4c(%rbp)
mov    0x40(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    13a1 <emrqmin+0x221>
mov    0x18(%rbp),%esi
mov    $0x1,%ecx
mov    %ecx,%edi
mov    %ecx,%edx
mov    $0x0,%al
call   1080 <matrix@plt>
mov    %rax,0x2e83(%rip)
mov    -0x4c(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1090 <vector@plt>
mov    %rax,0x2e65(%rip)
mov    -0x4c(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1090 <vector@plt>
mov    %rax,0x2e47(%rip)
mov    -0x4c(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1090 <vector@plt>
mov    %rax,0x2e49(%rip)
mov    0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jg     1300 <emrqmin+0x180>
movl   $0x0,-0x5c(%rbp)
movl   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     129d <emrqmin+0x11d>
mov    0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x58(%rbp),%eax
jne    128a <emrqmin+0x10a>
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    128f <emrqmin+0x10f>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1261 <emrqmin+0xe1>
cmpl   $0x0,-0x5c(%rbp)
jne    12c4 <emrqmin+0x144>
mov    -0x58(%rbp),%edx
mov    0x10(%rbp),%rax
mov    -0x54(%rbp),%ecx
mov    %ecx,%esi
add    $0x1,%esi
mov    %esi,-0x54(%rbp)
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
jmp    12ed <emrqmin+0x16d>
cmpl   $0x1,-0x5c(%rbp)
jle    12e8 <emrqmin+0x168>
lea    0xd4b(%rip),%rdi
mov    $0x0,%al
call   1040 <nrerror@plt>
movl   $0x0,-0x1c(%rbp)
jmp    1822 <emrqmin+0x6a2>
jmp    12ed <emrqmin+0x16d>
jmp    12f2 <emrqmin+0x172>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1247 <emrqmin+0xc7>
mov    -0x54(%rbp),%eax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
je     132b <emrqmin+0x1ab>
lea    0xd2a(%rip),%rdi
mov    $0x0,%al
call   1040 <nrerror@plt>
movl   $0x0,-0x1c(%rbp)
jmp    1822 <emrqmin+0x6a2>
mov    0x40(%rbp),%rax
movss  0xcc9(%rip),%xmm0
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rdi
mov    -0x30(%rbp),%rsi
mov    -0x38(%rbp),%rdx
mov    -0x3c(%rbp),%ecx
mov    -0x48(%rbp),%r8
mov    -0x4c(%rbp),%r9d
mov    0x10(%rbp),%r15
mov    0x18(%rbp),%r14d
mov    0x28(%rbp),%rbx
mov    0x2d1b(%rip),%r11
mov    0x30(%rbp),%r10
mov    0x38(%rbp),%rax
mov    %r15,(%rsp)
mov    %r14d,0x8(%rsp)
mov    %rbx,0x10(%rsp)
mov    %r11,0x18(%rsp)
mov    %r10,0x20(%rsp)
mov    %rax,0x28(%rsp)
mov    $0x0,%al
call   10a0 <mrqcof@plt>
mov    0x30(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,0x2ce7(%rip)
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     147f <emrqmin+0x2ff>
movl   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     13ff <emrqmin+0x27f>
mov    0x28(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x20(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    13bb <emrqmin+0x23b>
mov    0x28(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x58(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    0x40(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xbec(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    0x20(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x58(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    0x2c2c(%rip),%rax
movslq -0x58(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x2c14(%rip),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movss  %xmm0,0x4(%rax)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    13a8 <emrqmin+0x228>
mov    0x20(%rbp),%rdi
mov    0x18(%rbp),%esi
mov    0x2beb(%rip),%rdx
mov    $0x1,%ecx
mov    $0x0,%al
call   1060 <gaussj@plt>
cmp    $0x0,%eax
jne    14ae <emrqmin+0x32e>
movl   $0x0,-0x1c(%rbp)
jmp    1822 <emrqmin+0x6a2>
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     14f3 <emrqmin+0x373>
mov    0x2bb0(%rip),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movss  0x4(%rax),%xmm0
mov    0x2b8c(%rip),%rax
movslq -0x58(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    14b5 <emrqmin+0x335>
mov    0x40(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1592 <emrqmin+0x412>
jp     1592 <emrqmin+0x412>
mov    0x20(%rbp),%rdi
mov    -0x4c(%rbp),%esi
mov    0x10(%rbp),%rdx
mov    0x18(%rbp),%ecx
mov    $0x0,%al
call   1070 <covsrt@plt>
mov    0x2b52(%rip),%rdi
mov    -0x4c(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1050 <free_vector@plt>
mov    0x2b24(%rip),%rdi
mov    -0x4c(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1050 <free_vector@plt>
mov    0x2b16(%rip),%rdi
mov    -0x4c(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1050 <free_vector@plt>
mov    0x2b08(%rip),%rdi
mov    0x18(%rbp),%edx
mov    $0x1,%r8d
mov    %r8d,%esi
mov    %r8d,%ecx
mov    $0x0,%al
call   1030 <free_matrix@plt>
movl   $0x1,-0x1c(%rbp)
jmp    1822 <emrqmin+0x6a2>
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     168b <emrqmin+0x50b>
mov    -0x48(%rbp),%rax
mov    0x10(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x2aa7(%rip),%rax
movslq -0x58(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x60(%rbp)
cmpl   $0x1,-0x58(%rbp)
setg   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x64(%rbp)
mov    -0x58(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    -0x64(%rbp),%eax
and    %edx,%eax
movss  -0x60(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
setae  %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     165b <emrqmin+0x4db>
mov    0x2a57(%rip),%rax
mov    0x10(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0x9d2(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    0x2a2b(%rip),%rax
mov    0x10(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1678 <emrqmin+0x4f8>
movss  -0x60(%rbp),%xmm0
mov    0x2a09(%rip),%rax
mov    0x10(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    167d <emrqmin+0x4fd>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1599 <emrqmin+0x419>
mov    -0x28(%rbp),%rdi
mov    -0x30(%rbp),%rsi
mov    -0x38(%rbp),%rdx
mov    -0x3c(%rbp),%ecx
mov    0x29cf(%rip),%r8
mov    -0x4c(%rbp),%r9d
mov    0x10(%rbp),%r15
mov    0x18(%rbp),%r14d
mov    0x20(%rbp),%rbx
mov    0x29b0(%rip),%r11
mov    0x30(%rbp),%r10
mov    0x38(%rbp),%rax
mov    %r15,(%rsp)
mov    %r14d,0x8(%rsp)
mov    %rbx,0x10(%rsp)
mov    %r11,0x18(%rsp)
mov    %r10,0x20(%rsp)
mov    %rax,0x28(%rsp)
mov    $0x0,%al
call   10a0 <mrqcof@plt>
mov    0x30(%rbp),%rax
movss  (%rax),%xmm1
movss  0x2994(%rip),%xmm0
ucomiss %xmm1,%xmm0
jbe    17eb <emrqmin+0x66b>
mov    0x40(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0x8ff(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    0x30(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,0x295b(%rip)
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     17e6 <emrqmin+0x666>
movl   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     178b <emrqmin+0x60b>
mov    0x20(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x28(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1747 <emrqmin+0x5c7>
mov    0x28d6(%rip),%rax
movslq -0x58(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x28de(%rip),%rax
movslq -0x58(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    0x28be(%rip),%rax
mov    0x10(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
mov    0x10(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1734 <emrqmin+0x5b4>
jmp    181b <emrqmin+0x69b>
mov    0x40(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0x809(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax)
movss  0x2875(%rip),%xmm0
mov    0x30(%rbp),%rax
movss  %xmm0,(%rax)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x88,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

x86 assembly (clang, -O2):
<emrqmin>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9d,%r14d
mov    %r8,0x18(%rsp)
mov    %ecx,0x10(%rsp)
mov    %rdx,0x40(%rsp)
mov    %rsi,0x38(%rsp)
mov    %rdi,0x30(%rsp)
mov    0xb0(%rsp),%rax
mov    0x98(%rsp),%rbx
mov    0x88(%rsp),%r13d
mov    0x80(%rsp),%r12
xorps  %xmm0,%xmm0
ucomiss (%rax),%xmm0
mov    %r14,0x8(%rsp)
jbe    1432 <emrqmin+0x2b2>
mov    $0x1,%ebp
mov    $0x1,%edi
mov    %r13d,%esi
mov    $0x1,%edx
mov    $0x1,%ecx
xor    %eax,%eax
call   1080 <matrix@plt>
mov    %rax,0x2e79(%rip)
mov    $0x1,%edi
mov    %r14d,%esi
xor    %eax,%eax
call   1090 <vector@plt>
mov    %rax,0x2e5b(%rip)
mov    $0x1,%edi
mov    %r14d,%esi
xor    %eax,%eax
call   1090 <vector@plt>
mov    %rax,0x2e3d(%rip)
mov    $0x1,%edi
mov    %r14d,%esi
xor    %eax,%eax
call   1090 <vector@plt>
mov    %rax,%r10
mov    %rax,0x2e3c(%rip)
lea    0x1(%r13),%eax
test   %r14d,%r14d
mov    %r14,%rsi
jle    13c7 <emrqmin+0x247>
mov    %eax,%ecx
lea    -0x1(%rcx),%rbx
lea    -0x9(%rcx),%r8
mov    %r8,%r15
shr    $0x3,%r15
add    $0x1,%r15
mov    %rbx,%r14
and    $0xfffffffffffffff8,%r14
lea    0x1(%r14),%r9
movabs $0x3ffffffffffffffc,%r11
ff ff 3f
or     $0x2,%r11
and    %r15,%r11
jmp    12b4 <emrqmin+0x134>
nopw   0x0(%rax,%rax,1)
test   %esi,%esi
mov    0x8(%rsp),%rsi
jne    18b3 <emrqmin+0x733>
movslq %eax,%rdx
add    $0x1,%eax
mov    %ebp,(%r12,%rdx,4)
lea    0x1(%rbp),%edx
cmp    %esi,%ebp
mov    %edx,%ebp
je     13c7 <emrqmin+0x247>
test   %r13d,%r13d
jle    129d <emrqmin+0x11d>
cmp    $0x8,%rbx
jae    12d0 <emrqmin+0x150>
xor    %esi,%esi
mov    $0x1,%edi
jmp    13a0 <emrqmin+0x220>
nopl   0x0(%rax,%rax,1)
movd   %ebp,%xmm0
pshufd $0x0,%xmm0,%xmm1
pxor   %xmm0,%xmm0
cmp    $0x8,%r8
jae    12eb <emrqmin+0x16b>
xor    %esi,%esi
pxor   %xmm2,%xmm2
jmp    1346 <emrqmin+0x1c6>
mov    %r11,%rdi
xor    %esi,%esi
pxor   %xmm2,%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movdqu 0x4(%r12,%rsi,4),%xmm3
movdqu 0x14(%r12,%rsi,4),%xmm4
movdqu 0x24(%r12,%rsi,4),%xmm5
movdqu 0x34(%r12,%rsi,4),%xmm6
pcmpeqd %xmm1,%xmm3
psubd  %xmm3,%xmm0
pcmpeqd %xmm1,%xmm4
psubd  %xmm4,%xmm2
pcmpeqd %xmm1,%xmm5
psubd  %xmm5,%xmm0
pcmpeqd %xmm1,%xmm6
psubd  %xmm6,%xmm2
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1300 <emrqmin+0x180>
test   $0x1,%r15b
je     1371 <emrqmin+0x1f1>
shl    $0x2,%rsi
or     $0x4,%rsi
movdqu (%r12,%rsi,1),%xmm3
movdqu 0x10(%r12,%rsi,1),%xmm4
pcmpeqd %xmm1,%xmm3
psubd  %xmm3,%xmm0
pcmpeqd %xmm1,%xmm4
psubd  %xmm4,%xmm2
paddd  %xmm2,%xmm0
pshufd $0xee,%xmm0,%xmm1
paddd  %xmm0,%xmm1
pshufd $0x55,%xmm1,%xmm0
paddd  %xmm1,%xmm0
movd   %xmm0,%esi
mov    %r9,%rdi
cmp    %r14,%rbx
je     13b4 <emrqmin+0x234>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %edx,%edx
cmp    %ebp,(%r12,%rdi,4)
sete   %dl
add    %edx,%esi
add    $0x1,%rdi
cmp    %rdi,%rcx
jne    13a0 <emrqmin+0x220>
cmp    $0x1,%esi
jne    1290 <emrqmin+0x110>
mov    0x8(%rsp),%rsi
jmp    12a7 <emrqmin+0x127>
lea    0x1(%rsi),%ecx
cmp    %ecx,%eax
jne    1821 <emrqmin+0x6a1>
mov    0xb0(%rsp),%rax
movl   $0x3a83126f,(%rax)
mov    %rsi,%r9
mov    0x30(%rsp),%rdi
mov    0x38(%rsp),%rsi
mov    0x40(%rsp),%rdx
mov    0x10(%rsp),%ecx
mov    0x18(%rsp),%r8
xor    %eax,%eax
push   0xa8(%rsp)
mov    0xa8(%rsp),%rbp
push   %rbp
push   %r10
mov    0xb0(%rsp),%rbx
push   %rbx
push   %r13
push   %r12
call   10a0 <mrqcof@plt>
add    $0x30,%rsp
movss  0x0(%rbp),%xmm0
movss  %xmm0,0x2c56(%rip)
mov    0x90(%rsp),%rdi
test   %r13d,%r13d
jle    1759 <emrqmin+0x5d9>
mov    %r13,%rax
mov    0x2c33(%rip),%r13
mov    0x2c24(%rip),%rdx
mov    %eax,%r15d
lea    0x1(%r15),%rbp
lea    -0x8(%r15),%r9
mov    %r9,0x28(%rsp)
shr    $0x3,%r9
add    $0x1,%r9
mov    %r15d,%r11d
and    $0xfffffff8,%r11d
lea    0x1(%r11),%rax
mov    %rax,0x20(%rsp)
mov    %r9d,%r10d
and    $0x3,%r10d
and    $0xfffffffffffffffc,%r9
mov    $0x1,%ecx
movsd  0xb6c(%rip),%xmm0
jmp    14eb <emrqmin+0x36b>
cs nopw 0x0(%rax,%rax,1)
movss  (%r14,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mov    0xb0(%rsp),%rax
movss  (%rax),%xmm2
cvtss2sd %xmm2,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
movss  %xmm1,(%rdi,%rcx,4)
movd   0x0(%r13,%rcx,4),%xmm1
mov    (%rdx,%rcx,8),%rax
movd   %xmm1,0x4(%rax)
add    $0x1,%rcx
cmp    %rbp,%rcx
je     1684 <emrqmin+0x504>
mov    (%rbx,%rcx,8),%r14
mov    0x90(%rsp),%rax
mov    (%rax,%rcx,8),%rdi
mov    $0x1,%r8d
cmp    $0x8,%r15d
jb     1600 <emrqmin+0x480>
lea    0x4(%rdi),%rax
lea    (%r14,%rbp,4),%rsi
cmp    %rsi,%rax
jae    152f <emrqmin+0x3af>
lea    (%rdi,%rbp,4),%rax
lea    0x4(%r14),%rsi
mov    $0x1,%r8d
cmp    %rax,%rsi
jb     1600 <emrqmin+0x480>
cmpq   $0x18,0x28(%rsp)
jae    153b <emrqmin+0x3bb>
xor    %esi,%esi
jmp    15a6 <emrqmin+0x426>
mov    %r9,%rax
xor    %esi,%esi
movups 0x4(%r14,%rsi,4),%xmm1
movups 0x14(%r14,%rsi,4),%xmm2
movups %xmm1,0x4(%rdi,%rsi,4)
movups %xmm2,0x14(%rdi,%rsi,4)
movups 0x24(%r14,%rsi,4),%xmm1
movups 0x34(%r14,%rsi,4),%xmm2
movups %xmm1,0x24(%rdi,%rsi,4)
movups %xmm2,0x34(%rdi,%rsi,4)
movups 0x44(%r14,%rsi,4),%xmm1
movups 0x54(%r14,%rsi,4),%xmm2
movups %xmm1,0x44(%rdi,%rsi,4)
movups %xmm2,0x54(%rdi,%rsi,4)
movdqu 0x64(%r14,%rsi,4),%xmm1
movdqu 0x74(%r14,%rsi,4),%xmm2
movdqu %xmm1,0x64(%rdi,%rsi,4)
movdqu %xmm2,0x74(%rdi,%rsi,4)
add    $0x20,%rsi
add    $0xfffffffffffffffc,%rax
jne    1540 <emrqmin+0x3c0>
test   %r10,%r10
je     15de <emrqmin+0x45e>
mov    %r10,%rax
xchg   %ax,%ax
lea    0x0(,%rsi,4),%rbx
or     $0x4,%rbx
movdqu (%r14,%rbx,1),%xmm1
movdqu 0x10(%r14,%rbx,1),%xmm2
movdqu %xmm1,(%rdi,%rbx,1)
movdqu %xmm2,0x10(%rdi,%rbx,1)
add    $0x8,%rsi
add    $0xffffffffffffffff,%rax
jne    15b0 <emrqmin+0x430>
mov    0x20(%rsp),%r8
cmp    %r11,%r15
mov    0x98(%rsp),%rbx
je     14a0 <emrqmin+0x320>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %ebp,%esi
sub    %r8d,%esi
mov    %r8,%rax
not    %rax
add    %rbp,%rax
and    $0x3,%rsi
je     1636 <emrqmin+0x4b6>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movd   (%r14,%r8,4),%xmm1
movd   %xmm1,(%rdi,%r8,4)
add    $0x1,%r8
add    $0xffffffffffffffff,%rsi
jne    1620 <emrqmin+0x4a0>
cmp    $0x3,%rax
jb     14a0 <emrqmin+0x320>
movss  (%r14,%r8,4),%xmm1
movss  %xmm1,(%rdi,%r8,4)
movss  0x4(%r14,%r8,4),%xmm1
movss  %xmm1,0x4(%rdi,%r8,4)
movss  0x8(%r14,%r8,4),%xmm1
movss  %xmm1,0x8(%rdi,%r8,4)
movd   0xc(%r14,%r8,4),%xmm1
movd   %xmm1,0xc(%rdi,%r8,4)
add    $0x4,%r8
cmp    %r8,%rbp
jne    1640 <emrqmin+0x4c0>
jmp    14a0 <emrqmin+0x320>
xor    %r14d,%r14d
mov    0x90(%rsp),%rdi
mov    0x88(%rsp),%r13d
mov    %r13d,%esi
mov    $0x1,%ecx
xor    %eax,%eax
call   1060 <gaussj@plt>
test   %eax,%eax
je     1b94 <emrqmin+0xa14>
test   %r13d,%r13d
jle    177a <emrqmin+0x5fa>
mov    0x29ba(%rip),%rcx
mov    0x29a3(%rip),%rdx
add    $0xfffffffffffffffe,%rbp
mov    %r15d,%eax
and    $0x3,%eax
mov    $0x1,%esi
cmp    $0x3,%rbp
jb     172d <emrqmin+0x5ad>
and    $0xfffffffc,%r15d
xor    %esi,%esi
mov    0x8(%rcx,%rsi,8),%rdi
movss  0x4(%rdi),%xmm0
movss  %xmm0,0x4(%rdx,%rsi,4)
mov    0x10(%rcx,%rsi,8),%rdi
movss  0x4(%rdi),%xmm0
movss  %xmm0,0x8(%rdx,%rsi,4)
mov    0x18(%rcx,%rsi,8),%rdi
movss  0x4(%rdi),%xmm0
movss  %xmm0,0xc(%rdx,%rsi,4)
mov    0x20(%rcx,%rsi,8),%rdi
movss  0x4(%rdi),%xmm0
movss  %xmm0,0x10(%rdx,%rsi,4)
add    $0x4,%rsi
cmp    %rsi,%r15
jne    16e0 <emrqmin+0x560>
add    $0x1,%rsi
test   %rax,%rax
je     177a <emrqmin+0x5fa>
lea    (%rdx,%rsi,4),%rdx
lea    (%rcx,%rsi,8),%rcx
xor    %esi,%esi
nopl   0x0(%rax)
mov    (%rcx,%rsi,8),%rdi
movss  0x4(%rdi),%xmm0
movss  %xmm0,(%rdx,%rsi,4)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1740 <emrqmin+0x5c0>
jmp    177a <emrqmin+0x5fa>
mov    0x2918(%rip),%rdx
xor    %r14d,%r14d
mov    %r13d,%esi
mov    $0x1,%ecx
xor    %eax,%eax
call   1060 <gaussj@plt>
test   %eax,%eax
je     1b94 <emrqmin+0xa14>
mov    0xb0(%rsp),%r15
movss  (%r15),%xmm0
pxor   %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    182d <emrqmin+0x6ad>
jp     182d <emrqmin+0x6ad>
mov    0x90(%rsp),%rdi
mov    0x8(%rsp),%rbx
mov    %ebx,%esi
mov    %r12,%rdx
mov    %r13d,%ecx
xor    %eax,%eax
call   1070 <covsrt@plt>
mov    0x28c3(%rip),%rdi
mov    $0x1,%r14d
mov    $0x1,%esi
mov    %ebx,%edx
xor    %eax,%eax
call   1050 <free_vector@plt>
mov    0x2890(%rip),%rdi
mov    $0x1,%esi
mov    %ebx,%edx
xor    %eax,%eax
call   1050 <free_vector@plt>
mov    0x2883(%rip),%rdi
mov    $0x1,%esi
mov    %ebx,%edx
xor    %eax,%eax
call   1050 <free_vector@plt>
mov    0x2876(%rip),%rdi
mov    $0x1,%esi
mov    %r13d,%edx
mov    $0x1,%ecx
mov    $0x1,%r8d
xor    %eax,%eax
call   1030 <free_matrix@plt>
jmp    1b94 <emrqmin+0xa14>
lea    0x80e(%rip),%rdi
jmp    18ba <emrqmin+0x73a>
test   %r13d,%r13d
jle    18c9 <emrqmin+0x749>
mov    0x282b(%rip),%rbp
mov    0x282c(%rip),%r8
lea    0x1(%r13),%eax
mov    $0x1,%ecx
xorps  %xmm0,%xmm0
movss  0x7b8(%rip),%xmm1
mov    0x18(%rsp),%rbx
mov    0xa0(%rsp),%r14
jmp    187f <emrqmin+0x6ff>
nopw   0x0(%rax,%rax,1)
movss  %xmm2,(%r8,%rdx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
je     18df <emrqmin+0x75f>
movslq (%r12,%rcx,4),%rdx
movss  (%rbx,%rdx,4),%xmm2
addss  0x0(%rbp,%rcx,4),%xmm2
xor    %esi,%esi
cmp    $0x2,%rcx
setae  %sil
and    %ecx,%esi
xor    %edi,%edi
ucomiss %xmm2,%xmm0
setae  %dil
test   %edi,%esi
je     1870 <emrqmin+0x6f0>
movss  (%r8,%rdx,4),%xmm2
divss  %xmm1,%xmm2
jmp    1870 <emrqmin+0x6f0>
lea    0x75a(%rip),%rdi
xor    %r14d,%r14d
xor    %eax,%eax
call   1040 <nrerror@plt>
jmp    1b94 <emrqmin+0xa14>
mov    0x27a0(%rip),%r8
mov    0x2791(%rip),%rbp
mov    0xa0(%rsp),%r14
mov    0x30(%rsp),%rdi
mov    0x38(%rsp),%rsi
mov    0x40(%rsp),%rdx
mov    0x10(%rsp),%ecx
mov    0x8(%rsp),%r9
xor    %eax,%eax
push   0xa8(%rsp)
push   %r14
push   %rbp
push   0xa8(%rsp)
push   %r13
push   %r12
call   10a0 <mrqcof@plt>
add    $0x30,%rsp
movss  0x2769(%rip),%xmm0
ucomiss (%r14),%xmm0
movss  (%r15),%xmm1
jbe    1b7c <emrqmin+0x9fc>
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
mulsd  0x6cb(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r15)
movss  (%r14),%xmm0
movss  %xmm0,0x2735(%rip)
mov    $0x1,%r14d
test   %r13d,%r13d
mov    0x18(%rsp),%r9
jle    1b94 <emrqmin+0xa14>
mov    0x26fa(%rip),%r15
mov    0x270b(%rip),%r10
mov    %r13,%rax
mov    0x26f1(%rip),%r13
mov    %eax,%r11d
lea    0x1(%r11),%r8
lea    -0x8(%r11),%rax
mov    %rax,0x8(%rsp)
shr    $0x3,%rax
add    $0x1,%rax
mov    %r11d,%ecx
and    $0xfffffff8,%ecx
mov    %rcx,0x28(%rsp)
add    $0x1,%rcx
mov    %rcx,0x20(%rsp)
mov    %eax,%r14d
and    $0x3,%r14d
and    $0xfffffffffffffffc,%rax
mov    %rax,0x10(%rsp)
mov    $0x1,%ebp
mov    0x90(%rsp),%rbx
jmp    19fa <emrqmin+0x87a>
nopw   0x0(%rax,%rax,1)
movss  (%r15,%rbp,4),%xmm0
movss  %xmm0,(%r10,%rbp,4)
movslq (%r12,%rbp,4),%rax
movss  0x0(%r13,%rax,4),%xmm0
movss  %xmm0,(%r9,%rax,4)
add    $0x1,%rbp
cmp    %r8,%rbp
je     1b8e <emrqmin+0xa0e>
mov    (%rbx,%rbp,8),%rax
mov    0x98(%rsp),%rcx
mov    (%rcx,%rbp,8),%rcx
mov    $0x1,%edx
cmp    $0x8,%r11d
jb     1b00 <emrqmin+0x980>
lea    0x4(%rcx),%rdx
lea    (%rax,%r8,4),%rsi
cmp    %rsi,%rdx
jae    1a3c <emrqmin+0x8bc>
lea    (%rcx,%r8,4),%rsi
lea    0x4(%rax),%rdi
mov    $0x1,%edx
cmp    %rsi,%rdi
jb     1b00 <emrqmin+0x980>
cmpq   $0x18,0x8(%rsp)
jae    1a48 <emrqmin+0x8c8>
xor    %edx,%edx
jmp    1aaa <emrqmin+0x92a>
mov    0x10(%rsp),%r9
xor    %edx,%edx
nop
movups 0x4(%rax,%rdx,4),%xmm0
movups 0x14(%rax,%rdx,4),%xmm1
movups %xmm0,0x4(%rcx,%rdx,4)
movups %xmm1,0x14(%rcx,%rdx,4)
movups 0x24(%rax,%rdx,4),%xmm0
movups 0x34(%rax,%rdx,4),%xmm1
movups %xmm0,0x24(%rcx,%rdx,4)
movups %xmm1,0x34(%rcx,%rdx,4)
movups 0x44(%rax,%rdx,4),%xmm0
movups 0x54(%rax,%rdx,4),%xmm1
movups %xmm0,0x44(%rcx,%rdx,4)
movups %xmm1,0x54(%rcx,%rdx,4)
movups 0x64(%rax,%rdx,4),%xmm0
movups 0x74(%rax,%rdx,4),%xmm1
movups %xmm0,0x64(%rcx,%rdx,4)
movups %xmm1,0x74(%rcx,%rdx,4)
add    $0x20,%rdx
add    $0xfffffffffffffffc,%r9
jne    1a50 <emrqmin+0x8d0>
test   %r14,%r14
mov    0x18(%rsp),%r9
mov    0x90(%rsp),%rbx
je     1ae8 <emrqmin+0x968>
mov    %r14,%rsi
nop
lea    0x0(,%rdx,4),%rdi
or     $0x4,%rdi
movups (%rax,%rdi,1),%xmm0
movups 0x10(%rax,%rdi,1),%xmm1
movups %xmm0,(%rcx,%rdi,1)
movups %xmm1,0x10(%rcx,%rdi,1)
add    $0x8,%rdx
add    $0xffffffffffffffff,%rsi
jne    1ac0 <emrqmin+0x940>
mov    0x20(%rsp),%rdx
cmp    0x28(%rsp),%r11
je     19d0 <emrqmin+0x850>
nopl   0x0(%rax,%rax,1)
mov    %r8d,%edi
sub    %edx,%edi
mov    %rdx,%rsi
not    %rsi
add    %r8,%rsi
and    $0x3,%rdi
je     1b34 <emrqmin+0x9b4>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rax,%rdx,4),%xmm0
movss  %xmm0,(%rcx,%rdx,4)
add    $0x1,%rdx
add    $0xffffffffffffffff,%rdi
jne    1b20 <emrqmin+0x9a0>
cmp    $0x3,%rsi
jb     19d0 <emrqmin+0x850>
xchg   %ax,%ax
movss  (%rax,%rdx,4),%xmm0
movss  %xmm0,(%rcx,%rdx,4)
movss  0x4(%rax,%rdx,4),%xmm0
movss  %xmm0,0x4(%rcx,%rdx,4)
movss  0x8(%rax,%rdx,4),%xmm0
movss  %xmm0,0x8(%rcx,%rdx,4)
movss  0xc(%rax,%rdx,4),%xmm0
movss  %xmm0,0xc(%rcx,%rdx,4)
add    $0x4,%rdx
cmp    %rdx,%r8
jne    1b40 <emrqmin+0x9c0>
jmp    19d0 <emrqmin+0x850>
mulss  0x48c(%rip),%xmm1
movss  %xmm1,(%r15)
movss  %xmm0,(%r14)
mov    $0x1,%r14d
mov    %r14d,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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