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