int simq(A, B, X, n, flag, IPS)
double A[], B[], X[];
int n, flag;
int IPS[];
{
  int i, j, ij, ip, ipj, ipk, ipn;
  int idxpiv, iback;
  int k, kp, kp1, kpk, kpn;
  int nip, nkp, nm1;
  double em, q, rownrm, big, size, pivot, sum;
  nm1 = n - 1;
  if (flag < 0)
    goto solve;
  ij = 0;
  for (i = 0; i < n; i++) {
    IPS[i] = i;
    rownrm = 0.0;
    for (j = 0; j < n; j++) {
      q = ((A[ij]) < 0 ? -(A[ij]) : (A[ij]));
      if (rownrm < q)
        rownrm = q;
      ++ij;
    }
    if (rownrm == 0.0) {
      puts("SIMQ ROWNRM=0");
      return (1);
    }
    X[i] = 1.0 / rownrm;
  }
  for (k = 0; k < nm1; k++) {
    big = 0.0;
    idxpiv = 0;
    for (i = k; i < n; i++) {
      ip = IPS[i];
      ipk = n * ip + k;
      size = ((A[ipk]) < 0 ? -(A[ipk]) : (A[ipk])) * X[ip];
      if (size > big) {
        big = size;
        idxpiv = i;
      }
    }
    if (big == 0.0) {
      puts("SIMQ BIG=0");
      return (2);
    }
    if (idxpiv != k) {
      j = IPS[k];
      IPS[k] = IPS[idxpiv];
      IPS[idxpiv] = j;
    }
    kp = IPS[k];
    kpk = n * kp + k;
    pivot = A[kpk];
    kp1 = k + 1;
    for (i = kp1; i < n; i++) {
      ip = IPS[i];
      ipk = n * ip + k;
      em = -A[ipk] / pivot;
      A[ipk] = -em;
      nip = n * ip;
      nkp = n * kp;
      for (j = kp1; j < n; j++) {
        ipj = nip + j;
        A[ipj] = A[ipj] + em * A[nkp + j];
      }
    }
  }
  kpn = n * IPS[n - 1] + n - 1;
  if (A[kpn] == 0.0) {
    puts("SIMQ A[kpn]=0");
    return (3);
  }
solve:
  ip = IPS[0];
  X[0] = B[ip];
  for (i = 1; i < n; i++) {
    ip = IPS[i];
    ipj = n * ip;
    sum = 0.0;
    for (j = 0; j < i; j++) {
      sum += A[ipj] * X[j];
      ++ipj;
    }
    X[i] = B[ip] - sum;
  }
  ipn = n * IPS[n - 1] + n - 1;
  X[n - 1] = X[n - 1] / A[ipn];
  for (iback = 1; iback < n; iback++) {
    i = nm1 - iback;
    ip = IPS[i];
    nip = n * ip;
    sum = 0.0;
    for (j = i + 1; j < n; j++)
      sum += A[nip + j] * X[j];
    X[i] = (X[i] - sum) / A[nip + i];
  }
  return (0);
}


x86 assembly (clang, -O0):
<simq>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8d,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x74(%rbp)
cmpl   $0x0,-0x28(%rbp)
jge    114a <simq+0x3a>
jmp    15d2 <simq+0x4c2>
movl   $0x0,-0x3c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    129d <simq+0x18d>
mov    -0x34(%rbp),%edx
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    123d <simq+0x12d>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    11d7 <simq+0xc7>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
jmp    11ec <simq+0xdc>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x88(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
jbe    1226 <simq+0x116>
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1184 <simq+0x74>
movsd  -0x90(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1272 <simq+0x162>
jp     1272 <simq+0x162>
lea    0xda9(%rip),%rdi
mov    $0x0,%al
call   1030 <puts@plt>
movl   $0x1,-0x4(%rbp)
jmp    17e2 <simq+0x6d2>
movsd  0xd86(%rip),%xmm0
divsd  -0x90(%rbp),%xmm0
ff
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1158 <simq+0x48>
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x74(%rbp),%eax
jge    1576 <simq+0x466>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x0,-0x50(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    13a7 <simq+0x297>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x40(%rbp),%eax
add    -0x58(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    1336 <simq+0x226>
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
jmp    134b <simq+0x23b>
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
ucomisd -0x98(%rbp),%xmm0
ff
jbe    1394 <simq+0x284>
movsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x34(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    1399 <simq+0x289>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12c8 <simq+0x1b8>
movsd  -0x98(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    13dc <simq+0x2cc>
jp     13dc <simq+0x2cc>
lea    0xc4d(%rip),%rdi
mov    $0x0,%al
call   1030 <puts@plt>
movl   $0x2,-0x4(%rbp)
jmp    17e2 <simq+0x6d2>
mov    -0x50(%rbp),%eax
cmp    -0x58(%rbp),%eax
je     141a <simq+0x30a>
mov    -0x30(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x30(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x38(%rbp),%edx
mov    -0x30(%rbp),%rax
movslq -0x50(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x30(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x5c(%rbp),%eax
add    -0x58(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1563 <simq+0x453>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x40(%rbp),%eax
add    -0x58(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
imul   -0x40(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x5c(%rbp),%eax
mov    %eax,-0x70(%rbp)
mov    -0x60(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1550 <simq+0x440>
mov    -0x6c(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x80(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x70(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    14f4 <simq+0x3e4>
jmp    1555 <simq+0x445>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1459 <simq+0x349>
jmp    1568 <simq+0x458>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    12a4 <simq+0x194>
mov    -0x24(%rbp),%eax
mov    -0x30(%rbp),%rcx
mov    -0x24(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
imul   (%rcx,%rdx,4),%eax
add    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    15cd <simq+0x4bd>
jp     15cd <simq+0x4bd>
lea    0xa67(%rip),%rdi
mov    $0x0,%al
call   1030 <puts@plt>
movl   $0x3,-0x4(%rbp)
jmp    17e2 <simq+0x6d2>
jmp    15d2 <simq+0x4c2>
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    16b2 <simq+0x5a2>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x40(%rbp),%eax
mov    %eax,-0x44(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1682 <simq+0x572>
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xb0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    162d <simq+0x51d>
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0xb0(%rbp),%xmm0
ff
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    15f7 <simq+0x4e7>
mov    -0x24(%rbp),%eax
mov    -0x30(%rbp),%rcx
mov    -0x24(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
imul   (%rcx,%rdx,4),%eax
add    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    17db <simq+0x6cb>
mov    -0x74(%rbp),%eax
sub    -0x54(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x40(%rbp),%eax
mov    %eax,-0x6c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1799 <simq+0x689>
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xb0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1748 <simq+0x638>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0xb0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1707 <simq+0x5f7>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<simq>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9,%r12
mov    %rsi,0x30(%rsp)
lea    -0x1(%rcx),%r13d
movslq %ecx,%rax
mov    %rax,0x10(%rsp)
test   %r8d,%r8d
js     15ab <simq+0x49b>
test   %ecx,%ecx
jle    157a <simq+0x46a>
mov    %r13d,%eax
mov    %rax,0x40(%rsp)
lea    0x1(%rax),%rbx
mov    %ecx,%r11d
mov    %ecx,%r8d
and    $0xfffffffe,%r8d
xor    %r14d,%r14d
lea    0xeb2(%rip),%r9
xorpd  %xmm2,%xmm2
movsd  0xe9d(%rip),%xmm8
movapd 0xe84(%rip),%xmm9
movapd 0xe7c(%rip),%xmm1
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
mov    %r14d,(%r12,%r14,4)
movslq %esi,%rsi
mov    %rsi,%rax
xorpd  %xmm4,%xmm4
test   %r13d,%r13d
je     1206 <simq+0xf6>
xorpd  %xmm4,%xmm4
mov    %r8d,%ebp
mov    %rsi,%rax
nopl   (%rax)
movsd  (%rdi,%rax,8),%xmm3
movsd  0x8(%rdi,%rax,8),%xmm6
movapd %xmm3,%xmm7
xorpd  %xmm1,%xmm7
movapd %xmm3,%xmm5
cmpltsd %xmm2,%xmm5
andpd  %xmm5,%xmm7
andnpd %xmm3,%xmm5
orpd   %xmm7,%xmm5
maxsd  %xmm4,%xmm5
movapd %xmm6,%xmm3
xorpd  %xmm1,%xmm3
movapd %xmm6,%xmm4
cmpltsd %xmm2,%xmm4
andpd  %xmm4,%xmm3
andnpd %xmm6,%xmm4
orpd   %xmm3,%xmm4
maxsd  %xmm5,%xmm4
add    $0x2,%rax
add    $0xfffffffe,%ebp
jne    11b0 <simq+0xa0>
test   $0x1,%cl
je     1236 <simq+0x126>
movsd  (%rdi,%rax,8),%xmm3
movapd %xmm3,%xmm5
xorpd  %xmm9,%xmm5
movapd %xmm3,%xmm6
cmpltsd %xmm2,%xmm6
andpd  %xmm6,%xmm5
andnpd %xmm3,%xmm6
orpd   %xmm5,%xmm6
maxsd  %xmm4,%xmm6
movapd %xmm6,%xmm4
ucomisd %xmm2,%xmm4
jne    1242 <simq+0x132>
jnp    181f <simq+0x70f>
add    %rbx,%rsi
movapd %xmm8,%xmm3
divsd  %xmm4,%xmm3
movsd  %xmm3,(%rdx,%r14,8)
add    $0x1,%r14
cmp    %r11,%r14
jne    1190 <simq+0x80>
cmp    $0x2,%ecx
jl     157a <simq+0x46a>
mov    %r11,%rax
neg    %rax
mov    %rax,0x80(%rsp)
lea    -0x1(%r11),%r8
lea    0x18(%rdi),%rax
mov    %rax,(%rsp)
lea    0x8(%rdi),%rax
mov    %rax,0x28(%rsp)
mov    $0x1,%r9d
xor    %r14d,%r14d
xorpd  %xmm8,%xmm8
mov    %r13d,0x1c(%rsp)
mov    %rcx,0x38(%rsp)
jmp    12df <simq+0x1cf>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
mov    0x50(%rsp),%r8
add    $0xffffffffffffffff,%r8
addq   $0x8,(%rsp)
mov    0x48(%rsp),%r14
cmp    0x40(%rsp),%r14
mov    %rbx,%r12
mov    0x1c(%rsp),%r13d
mov    0x38(%rsp),%rcx
je     157a <simq+0x46a>
cmp    %r11,%r14
je     1826 <simq+0x716>
mov    %r8,%r15
and    $0xfffffffffffffffc,%r15
mov    %r14,%rax
not    %rax
add    %r11,%rax
mov    %rax,0x8(%rsp)
lea    (%rdi,%rax,8),%rax
mov    %rax,0x20(%rsp)
xor    %eax,%eax
xorpd  %xmm3,%xmm3
mov    %r14,%rsi
nop
movapd %xmm3,%xmm4
movslq (%r12,%rsi,4),%rbp
mov    %ebp,%ebx
imul   %ecx,%ebx
add    %r14d,%ebx
movslq %ebx,%rbx
movsd  (%rdi,%rbx,8),%xmm5
movapd %xmm5,%xmm6
xorpd  %xmm9,%xmm6
movapd %xmm5,%xmm3
cmpltsd %xmm8,%xmm3
andpd  %xmm3,%xmm6
andnpd %xmm5,%xmm3
orpd   %xmm6,%xmm3
mulsd  (%rdx,%rbp,8),%xmm3
ucomisd %xmm4,%xmm3
cmova  %esi,%eax
maxsd  %xmm4,%xmm3
add    $0x1,%rsi
cmp    %rsi,%r11
jne    1310 <simq+0x200>
ucomisd %xmm8,%xmm3
jne    136d <simq+0x25d>
jnp    1826 <simq+0x716>
mov    %r8,0x50(%rsp)
mov    %eax,%ecx
cmp    %rcx,%r14
je     138b <simq+0x27b>
mov    (%r12,%r14,4),%ecx
cltq
mov    (%r12,%rax,4),%esi
mov    %esi,(%r12,%r14,4)
mov    %ecx,(%r12,%rax,4)
mov    %r12,%rbx
lea    0x1(%r14),%rax
mov    %rax,0x48(%rsp)
cmp    %r11,%rax
jae    12b0 <simq+0x1a0>
movslq (%rbx,%r14,4),%rsi
imul   0x10(%rsp),%rsi
lea    (%rsi,%r14,1),%eax
cltq
movsd  (%rdi,%rax,8),%xmm3
lea    (%r9,%rsi,1),%rax
lea    (%rdi,%rax,8),%rcx
mov    %rcx,0x68(%rsp)
mov    0x20(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
mov    %rax,0x78(%rsp)
mov    0x8(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x60(%rsp)
add    %r9,%rax
mov    %rax,0x58(%rsp)
mov    (%rsp),%rax
lea    (%rax,%rsi,8),%r12
mov    0x28(%rsp),%rax
mov    %rsi,0x70(%rsp)
lea    (%rax,%rsi,8),%r13
mov    %r9,%r8
jmp    141d <simq+0x30d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
cmp    %r11,%r8
je     12b0 <simq+0x1a0>
movslq (%rbx,%r8,4),%rax
imul   0x10(%rsp),%rax
lea    (%rax,%r14,1),%ecx
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm4
xorpd  %xmm1,%xmm4
divsd  %xmm3,%xmm4
movapd %xmm4,%xmm5
xorpd  %xmm1,%xmm5
movlpd %xmm5,(%rdi,%rcx,8)
mov    %r9,%r10
cmpq   $0x4,0x8(%rsp)
jb     14f0 <simq+0x3e0>
lea    (%r9,%rax,1),%rsi
lea    (%rdi,%rsi,8),%rcx
cmp    0x78(%rsp),%rcx
jae    1479 <simq+0x369>
mov    0x20(%rsp),%rcx
lea    (%rcx,%rsi,8),%rcx
mov    %r9,%r10
cmp    %rcx,0x68(%rsp)
jb     14f0 <simq+0x3e0>
movapd %xmm4,%xmm5
unpcklpd %xmm4,%xmm5
mov    (%rsp),%rcx
lea    (%rcx,%rax,8),%rbp
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movupd -0x10(%rbp,%rsi,8),%xmm6
movupd 0x0(%rbp,%rsi,8),%xmm7
movupd -0x10(%r12,%rsi,8),%xmm0
movupd (%r12,%rsi,8),%xmm2
mulpd  %xmm5,%xmm0
addpd  %xmm6,%xmm0
mulpd  %xmm5,%xmm2
addpd  %xmm7,%xmm2
movupd %xmm0,-0x10(%rbp,%rsi,8)
movupd %xmm2,0x0(%rbp,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%r15
jne    1490 <simq+0x380>
mov    0x58(%rsp),%r10
mov    0x60(%rsp),%rcx
cmp    %rcx,0x8(%rsp)
je     1410 <simq+0x300>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r11d,%ecx
sub    %r10d,%ecx
mov    %r10,%rbp
test   $0x1,%cl
je     1521 <simq+0x411>
lea    (%r10,%rax,1),%rcx
mov    0x70(%rsp),%rsi
add    %r10,%rsi
movsd  (%rdi,%rsi,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
lea    0x1(%r10),%rbp
not    %r10
cmp    0x80(%rsp),%r10
je     1410 <simq+0x300>
mov    0x28(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%r13,%rbp,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  -0x8(%rax,%rbp,8),%xmm0
movsd  %xmm0,-0x8(%rax,%rbp,8)
movsd  0x0(%r13,%rbp,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  (%rax,%rbp,8),%xmm0
movsd  %xmm0,(%rax,%rbp,8)
add    $0x2,%rbp
cmp    %rbp,%r11
jne    1540 <simq+0x430>
jmp    1410 <simq+0x300>
movslq %r13d,%rax
mov    (%r12,%rax,4),%eax
imul   %ecx,%eax
add    %r13d,%eax
cltq
movsd  (%rdi,%rax,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    15ab <simq+0x49b>
jp     15ab <simq+0x49b>
mov    $0x3,%ebx
lea    0xa8b(%rip),%r9
jmp    1832 <simq+0x722>
movslq (%r12),%rax
mov    0x30(%rsp),%rsi
movsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx)
mov    %rcx,%r9
cmp    $0x2,%ecx
jl     16dc <simq+0x5cc>
mov    %r9d,%r11d
lea    0x18(%rdi),%r8
mov    $0x1,%esi
xor    %r10d,%r10d
jmp    1605 <simq+0x4f5>
nopw   0x0(%rax,%rax,1)
mov    0x30(%rsp),%rax
movsd  (%rax,%r14,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
add    $0x1,%rsi
add    $0x1,%r10
cmp    %r11,%rsi
je     16dc <simq+0x5cc>
movslq (%r12,%rsi,4),%r14
mov    %r14d,%eax
imul   %r9d,%eax
movslq %eax,%r15
cmp    $0x3,%r10
jae    1630 <simq+0x520>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
test   $0x3,%sil
jne    16a5 <simq+0x595>
jmp    15e0 <simq+0x4d0>
nopl   0x0(%rax,%rax,1)
mov    %rsi,%rbx
and    $0xfffffffffffffffc,%rbx
lea    (%r8,%r15,8),%rbp
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x18(%rbp,%rax,8),%xmm1
movsd  -0x10(%rbp,%rax,8),%xmm2
mulsd  (%rdx,%rax,8),%xmm1
mulsd  0x8(%rdx,%rax,8),%xmm2
addsd  %xmm0,%xmm1
movsd  -0x8(%rbp,%rax,8),%xmm3
mulsd  0x10(%rdx,%rax,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x0(%rbp,%rax,8),%xmm0
mulsd  0x18(%rdx,%rax,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rax
cmp    %rax,%rbx
jne    1650 <simq+0x540>
add    %rax,%r15
test   $0x3,%sil
je     15e0 <simq+0x4d0>
mov    %esi,%ecx
and    $0x3,%ecx
lea    (%rdi,%r15,8),%rbx
lea    (%rdx,%rax,8),%rax
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rbx,%rbp,8),%xmm1
mulsd  (%rax,%rbp,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rbp
cmp    %rbp,%rcx
jne    16c0 <simq+0x5b0>
jmp    15e0 <simq+0x4d0>
movslq %r13d,%r11
mov    (%r12,%r11,4),%esi
mov    %r9,%rcx
imul   %ecx,%esi
add    %r13d,%esi
xor    %eax,%eax
movsd  (%rdx,%r11,8),%xmm0
movslq %esi,%rsi
divsd  (%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%r11,8)
cmp    $0x2,%ecx
jl     183e <simq+0x72e>
mov    %ecx,%ecx
lea    0x18(%rdi),%r8
mov    $0x1,%r10d
xor    %r14d,%r14d
jmp    1750 <simq+0x640>
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
movsd  (%rdx,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
add    %rsi,%r15
divsd  (%rdi,%r15,8),%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
add    $0x1,%r10
add    $0xffffffff,%r13d
add    $0x1,%r14d
cmp    %rcx,%r10
je     181b <simq+0x70b>
mov    %r11,%rsi
sub    %r10,%rsi
movslq (%r12,%rsi,4),%r15
imul   0x10(%rsp),%r15
lea    0x1(%rsi),%rax
cmp    %rcx,%rax
jge    1720 <simq+0x610>
lea    0x1(%r14),%eax
movslq %r13d,%rbp
movslq %r15d,%r15
test   $0x3,%al
je     17b3 <simq+0x6a3>
mov    %r10d,%ebx
and    $0x3,%ebx
lea    (%rdi,%r15,8),%rax
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rax,%rbp,8),%xmm1
mulsd  (%rdx,%rbp,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rbp
add    $0xffffffff,%ebx
jne    1790 <simq+0x680>
cmp    $0x3,%r14d
jb     1724 <simq+0x614>
jmp    17c1 <simq+0x6b1>
xorpd  %xmm0,%xmm0
cmp    $0x3,%r14d
jb     1724 <simq+0x614>
lea    (%r8,%r15,8),%rbx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x18(%rbx,%rbp,8),%xmm1
movsd  -0x10(%rbx,%rbp,8),%xmm2
mulsd  (%rdx,%rbp,8),%xmm1
mulsd  0x8(%rdx,%rbp,8),%xmm2
addsd  %xmm0,%xmm1
movsd  -0x8(%rbx,%rbp,8),%xmm3
mulsd  0x10(%rdx,%rbp,8),%xmm3
addsd  %xmm1,%xmm2
movsd  (%rbx,%rbp,8),%xmm0
mulsd  0x18(%rdx,%rbp,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rbp
cmp    %ebp,%ecx
jne    17d0 <simq+0x6c0>
jmp    1724 <simq+0x614>
xor    %eax,%eax
jmp    183e <simq+0x72e>
mov    $0x1,%ebx
jmp    1832 <simq+0x722>
mov    $0x2,%ebx
lea    0x7f4(%rip),%r9
mov    %r9,%rdi
xor    %eax,%eax
call   1030 <puts@plt>
mov    %ebx,%eax
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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