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;
  double fabs();
  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 = fabs(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 = fabs(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    $0xb0,%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    153a <simq+0x42a>
movl   $0x0,-0x3c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1251 <simq+0x141>
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    11f1 <simq+0xe1>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe5c(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x88(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
jbe    11da <simq+0xca>
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    1226 <simq+0x116>
jp     1226 <simq+0x116>
lea    0xe05(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x1,-0x4(%rbp)
jmp    174a <simq+0x63a>
movsd  0xde2(%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    14de <simq+0x3ce>
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    130f <simq+0x1ff>
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
movaps 0xd49(%rip),%xmm1
pand   %xmm1,%xmm0
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    12fc <simq+0x1ec>
movsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x34(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    1301 <simq+0x1f1>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    127c <simq+0x16c>
movsd  -0x98(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1344 <simq+0x234>
jp     1344 <simq+0x234>
lea    0xcf5(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x2,-0x4(%rbp)
jmp    174a <simq+0x63a>
mov    -0x50(%rbp),%eax
cmp    -0x58(%rbp),%eax
je     1382 <simq+0x272>
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    14cb <simq+0x3bb>
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    14b8 <simq+0x3a8>
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    145c <simq+0x34c>
jmp    14bd <simq+0x3ad>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    13c1 <simq+0x2b1>
jmp    14d0 <simq+0x3c0>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1258 <simq+0x148>
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    1535 <simq+0x425>
jp     1535 <simq+0x425>
lea    0xb0f(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x3,-0x4(%rbp)
jmp    174a <simq+0x63a>
jmp    153a <simq+0x42a>
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    161a <simq+0x50a>
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    15ea <simq+0x4da>
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    1595 <simq+0x485>
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    155f <simq+0x44f>
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    1743 <simq+0x633>
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    1701 <simq+0x5f1>
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    16b0 <simq+0x5a0>
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    166f <simq+0x55f>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<simq>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
lea    -0x1(%rcx),%r13d
movslq %ecx,%rax
mov    %rax,0x20(%rsp)
test   %r8d,%r8d
js     1581 <simq+0x471>
mov    %rsi,0x40(%rsp)
test   %ecx,%ecx
jle    154a <simq+0x43a>
mov    %r13d,%eax
mov    %rax,0x58(%rsp)
lea    0x1(%rax),%rbx
mov    %ecx,%r11d
mov    %ecx,%r8d
and    $0xfffffffe,%r8d
xor    %r10d,%r10d
lea    0xec5(%rip),%r15
xorpd  %xmm1,%xmm1
movsd  0xeb1(%rip),%xmm2
movapd 0xe88(%rip),%xmm9
movapd 0xe80(%rip),%xmm3
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r10d,(%r9,%r10,4)
cltq
mov    %rax,%rsi
xorpd  %xmm4,%xmm4
test   %r13d,%r13d
je     11d8 <simq+0xc8>
xorpd  %xmm5,%xmm5
mov    %r8d,%ebp
mov    %rax,%rsi
nopl   0x0(%rax)
movsd  (%rdi,%rsi,8),%xmm6
movsd  0x8(%rdi,%rsi,8),%xmm4
andpd  %xmm3,%xmm6
maxsd  %xmm5,%xmm6
andpd  %xmm3,%xmm4
maxsd  %xmm6,%xmm4
add    $0x2,%rsi
movapd %xmm4,%xmm5
add    $0xfffffffe,%ebp
jne    11b0 <simq+0xa0>
test   $0x1,%cl
je     11ef <simq+0xdf>
movsd  (%rdi,%rsi,8),%xmm5
andpd  %xmm9,%xmm5
maxsd  %xmm4,%xmm5
movapd %xmm5,%xmm4
ucomisd %xmm1,%xmm4
jne    11fb <simq+0xeb>
jnp    1800 <simq+0x6f0>
add    %rbx,%rax
movapd %xmm2,%xmm5
divsd  %xmm4,%xmm5
movsd  %xmm5,(%rdx,%r10,8)
add    $0x1,%r10
cmp    %r11,%r10
jne    1190 <simq+0x80>
cmp    $0x2,%ecx
jl     154a <simq+0x43a>
mov    %r11,%rax
neg    %rax
mov    %rax,0x48(%rsp)
lea    -0x1(%r11),%r8
lea    0x18(%rdi),%rax
mov    %rax,0x10(%rsp)
lea    0x8(%rdi),%rax
mov    %rax,0x38(%rsp)
mov    $0x1,%eax
mov    %rax,0x8(%rsp)
xor    %r10d,%r10d
xorpd  %xmm8,%xmm8
movapd 0xdb3(%rip),%xmm2
mov    %r13d,0x2c(%rsp)
mov    %rcx,0x50(%rsp)
jmp    12a2 <simq+0x192>
nopl   0x0(%rax)
addq   $0x1,0x8(%rsp)
mov    0x68(%rsp),%r8
add    $0xffffffffffffffff,%r8
addq   $0x8,0x10(%rsp)
mov    0x60(%rsp),%r10
cmp    0x58(%rsp),%r10
mov    %rbx,%r9
mov    0x2c(%rsp),%r13d
mov    0x50(%rsp),%rcx
je     154a <simq+0x43a>
mov    $0x2,%r14d
cmp    %r11,%r10
jae    180b <simq+0x6fb>
mov    %r8,%r15
and    $0xfffffffffffffffc,%r15
mov    %r10,%rax
not    %rax
add    %r11,%rax
mov    %rax,0x18(%rsp)
lea    (%rdi,%rax,8),%rax
mov    %rax,0x30(%rsp)
xor    %eax,%eax
xorpd  %xmm3,%xmm3
mov    %r10,%rbp
nopl   0x0(%rax,%rax,1)
movapd %xmm3,%xmm4
movslq (%r9,%rbp,4),%rsi
mov    %esi,%ebx
imul   %ecx,%ebx
add    %r10d,%ebx
movslq %ebx,%rbx
movsd  (%rdi,%rbx,8),%xmm3
andpd  %xmm9,%xmm3
mulsd  (%rdx,%rsi,8),%xmm3
ucomisd %xmm4,%xmm3
cmova  %ebp,%eax
maxsd  %xmm4,%xmm3
add    $0x1,%rbp
cmp    %rbp,%r11
jne    12e0 <simq+0x1d0>
ucomisd %xmm8,%xmm3
jne    1323 <simq+0x213>
jnp    180b <simq+0x6fb>
mov    %r8,0x68(%rsp)
mov    %eax,%ecx
cmp    %rcx,%r10
je     1341 <simq+0x231>
mov    (%r9,%r10,4),%ecx
cltq
mov    (%r9,%rax,4),%esi
mov    %esi,(%r9,%r10,4)
mov    %ecx,(%r9,%rax,4)
mov    %r9,%rbx
lea    0x1(%r10),%rax
mov    %rax,0x60(%rsp)
cmp    %r11,%rax
mov    0x48(%rsp),%r9
jae    1270 <simq+0x160>
movslq (%rbx,%r10,4),%rsi
imul   0x20(%rsp),%rsi
lea    (%rsi,%r10,1),%eax
cltq
movsd  (%rdi,%rax,8),%xmm3
mov    0x8(%rsp),%r8
lea    (%r8,%rsi,1),%rax
lea    (%rdi,%rax,8),%rcx
mov    %rcx,0x80(%rsp)
mov    0x30(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
mov    %rax,0x90(%rsp)
mov    0x18(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x78(%rsp)
add    %r8,%rax
mov    %rax,0x70(%rsp)
mov    0x10(%rsp),%rax
lea    (%rax,%rsi,8),%r12
mov    0x38(%rsp),%rax
mov    %rsi,0x88(%rsp)
lea    (%rax,%rsi,8),%r13
jmp    13dd <simq+0x2cd>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
cmp    %r11,%r8
je     1270 <simq+0x160>
movslq (%rbx,%r8,4),%rax
imul   0x20(%rsp),%rax
lea    (%rax,%r10,1),%ecx
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm4
xorpd  %xmm2,%xmm4
divsd  %xmm3,%xmm4
movapd %xmm4,%xmm5
xorpd  %xmm2,%xmm5
movlpd %xmm5,(%rdi,%rcx,8)
mov    0x8(%rsp),%r14
cmpq   $0x4,0x18(%rsp)
jb     14c0 <simq+0x3b0>
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,1),%rsi
lea    (%rdi,%rsi,8),%rcx
cmp    0x90(%rsp),%rcx
jae    1448 <simq+0x338>
mov    0x30(%rsp),%rcx
lea    (%rcx,%rsi,8),%rcx
mov    0x8(%rsp),%r14
cmp    %rcx,0x80(%rsp)
jb     14c0 <simq+0x3b0>
movapd %xmm4,%xmm5
unpcklpd %xmm4,%xmm5
mov    0x10(%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),%xmm1
movupd (%r12,%rsi,8),%xmm0
mulpd  %xmm5,%xmm1
addpd  %xmm6,%xmm1
mulpd  %xmm5,%xmm0
addpd  %xmm7,%xmm0
movupd %xmm1,-0x10(%rbp,%rsi,8)
movupd %xmm0,0x0(%rbp,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%r15
jne    1460 <simq+0x350>
mov    0x70(%rsp),%r14
mov    0x78(%rsp),%rcx
cmp    %rcx,0x18(%rsp)
je     13d0 <simq+0x2c0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r11d,%ecx
sub    %r14d,%ecx
mov    %r14,%rbp
test   $0x1,%cl
je     14f4 <simq+0x3e4>
lea    (%r14,%rax,1),%rcx
mov    0x88(%rsp),%rsi
add    %r14,%rsi
movsd  (%rdi,%rsi,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
lea    0x1(%r14),%rbp
not    %r14
cmp    %r9,%r14
je     13d0 <simq+0x2c0>
mov    0x38(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
nopl   0x0(%rax)
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    1510 <simq+0x400>
jmp    13d0 <simq+0x2c0>
movslq %r13d,%rax
mov    (%r9,%rax,4),%eax
imul   %ecx,%eax
add    %r13d,%eax
cltq
movsd  (%rdi,%rax,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
mov    0x40(%rsp),%rsi
jne    1581 <simq+0x471>
jp     1581 <simq+0x471>
mov    $0x3,%r14d
lea    0xac5(%rip),%rdi
jmp    1812 <simq+0x702>
mov    %r9,%r11
movslq (%r9),%rax
movsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx)
mov    %rcx,%r8
cmp    $0x2,%ecx
jl     16bc <simq+0x5ac>
mov    %r8d,%r9d
lea    0x18(%rdi),%rax
mov    %rax,0x8(%rsp)
mov    $0x1,%r10d
xor    %r12d,%r12d
jmp    15e1 <simq+0x4d1>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rsi,%r14,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%r10,8)
add    $0x1,%r10
add    $0x1,%r12
cmp    %r9,%r10
je     16bc <simq+0x5ac>
movslq (%r11,%r10,4),%r14
mov    %r14d,%eax
imul   %r8d,%eax
movslq %eax,%r15
cmp    $0x3,%r12
jae    1610 <simq+0x500>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
test   $0x3,%r10b
jne    1685 <simq+0x575>
jmp    15c0 <simq+0x4b0>
nopw   0x0(%rax,%rax,1)
mov    %r10,%rbx
and    $0xfffffffffffffffc,%rbx
mov    0x8(%rsp),%rax
lea    (%rax,%r15,8),%rbp
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 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    1630 <simq+0x520>
add    %rax,%r15
test   $0x3,%r10b
je     15c0 <simq+0x4b0>
mov    %r10d,%ebp
and    $0x3,%ebp
lea    (%rdi,%r15,8),%rcx
lea    (%rdx,%rax,8),%rax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx,%rbx,8),%xmm1
mulsd  (%rax,%rbx,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rbx
cmp    %rbx,%rbp
jne    16a0 <simq+0x590>
jmp    15c0 <simq+0x4b0>
movslq %r13d,%r12
mov    (%r11,%r12,4),%eax
imul   %r8d,%eax
add    %r13d,%eax
xor    %r14d,%r14d
movsd  (%rdx,%r12,8),%xmm0
cltq
divsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%r12,8)
cmp    $0x2,%r8d
jl     1819 <simq+0x709>
mov    %r8d,%ecx
lea    0x18(%rdi),%r8
mov    $0x1,%r10d
xor    %r14d,%r14d
jmp    1730 <simq+0x620>
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     17fb <simq+0x6eb>
mov    %r12,%rsi
sub    %r10,%rsi
movslq (%r11,%rsi,4),%r15
imul   0x20(%rsp),%r15
lea    0x1(%rsi),%rax
cmp    %rcx,%rax
jge    1700 <simq+0x5f0>
lea    0x1(%r14),%eax
movslq %r13d,%rbp
movslq %r15d,%r15
test   $0x3,%al
je     1793 <simq+0x683>
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    1770 <simq+0x660>
cmp    $0x3,%r14d
jb     1704 <simq+0x5f4>
jmp    17a1 <simq+0x691>
xorpd  %xmm0,%xmm0
cmp    $0x3,%r14d
jb     1704 <simq+0x5f4>
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    17b0 <simq+0x6a0>
jmp    1704 <simq+0x5f4>
xor    %r14d,%r14d
jmp    1819 <simq+0x709>
mov    $0x1,%r14d
mov    %r15,%rdi
jmp    1812 <simq+0x702>
lea    0x824(%rip),%rdi
xor    %eax,%eax
call   1030 <printf@plt>
mov    %r14d,%eax
add    $0x98,%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 %277
43: ; preds = %6
store i32 0, i32* %16
store i32 0, i32* %14
br label %44
44: ; preds = %88, %43
%45 = load i32, i32* %14
%46 = load i32, i32* %11
%47 = icmp slt i32 %45, %46
br i1 %47, label %48, label %91
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 = %73, %48
%55 = load i32, i32* %15
%56 = load i32, i32* %11
%57 = icmp slt i32 %55, %56
br i1 %57, label %58, label %76
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 = call double @llvm.fabs.f64(double %63)
store double %64, double* %32
%65 = load double, double* %33
%66 = load double, double* %32
%67 = fcmp olt double %65, %66
br i1 %67, label %68, label %70
68: ; preds = %58
%69 = load double, double* %32
store double %69, double* %33
br label %70
70: ; preds = %68, %58
%71 = load i32, i32* %16
%72 = add nsw i32 %71, 1
store i32 %72, i32* %16
br label %73
73: ; preds = %70
%74 = load i32, i32* %15
%75 = add nsw i32 %74, 1
store i32 %75, i32* %15
br label %54
76: ; preds = %54
%77 = load double, double* %33
%78 = fcmp oeq double %77, 0.000000e+00
br i1 %78, label %79, label %81
79: ; preds = %76
%80 = 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 %431
81: ; preds = %76
%82 = load double, double* %33
%83 = fdiv double 1.000000e+00, %82
%84 = load double*, double** %10
%85 = load i32, i32* %14
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
store double %83, double* %87
br label %88
88: ; preds = %81
%89 = load i32, i32* %14
%90 = add nsw i32 %89, 1
store i32 %90, i32* %14
br label %44
91: ; preds = %44
store i32 0, i32* %23
br label %92
92: ; preds = %253, %91
%93 = load i32, i32* %23
%94 = load i32, i32* %30
%95 = icmp slt i32 %93, %94
br i1 %95, label %96, label %256
96: ; preds = %92
store double 0.000000e+00, double* %34
store i32 0, i32* %21
%97 = load i32, i32* %23
store i32 %97, i32* %14
br label %98
98: ; preds = %132, %96
%99 = load i32, i32* %14
%100 = load i32, i32* %11
%101 = icmp slt i32 %99, %100
br i1 %101, label %102, label %135
102: ; preds = %98
%103 = load i32*, i32** %13
%104 = load i32, i32* %14
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds i32, i32* %103, i64 %105
%107 = load i32, i32* %106
store i32 %107, i32* %17
%108 = load i32, i32* %11
%109 = load i32, i32* %17
%110 = mul nsw i32 %108, %109
%111 = load i32, i32* %23
%112 = add nsw i32 %110, %111
store i32 %112, i32* %19
%113 = load double*, double** %8
%114 = load i32, i32* %19
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds double, double* %113, i64 %115
%117 = load double, double* %116
%118 = call double @llvm.fabs.f64(double %117)
%119 = load double*, double** %10
%120 = load i32, i32* %17
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %119, i64 %121
%123 = load double, double* %122
%124 = fmul double %118, %123
store double %124, double* %35
%125 = load double, double* %35
%126 = load double, double* %34
%127 = fcmp ogt double %125, %126
br i1 %127, label %128, label %131
128: ; preds = %102
%129 = load double, double* %35
store double %129, double* %34
%130 = load i32, i32* %14
store i32 %130, i32* %21
br label %131
131: ; preds = %128, %102
br label %132
132: ; preds = %131
%133 = load i32, i32* %14
%134 = add nsw i32 %133, 1
store i32 %134, i32* %14
br label %98
135: ; preds = %98
%136 = load double, double* %34
%137 = fcmp oeq double %136, 0.000000e+00
br i1 %137, label %138, label %140
138: ; preds = %135
%139 = 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 %431
140: ; preds = %135
%141 = load i32, i32* %21
%142 = load i32, i32* %23
%143 = icmp ne i32 %141, %142
br i1 %143, label %144, label %164
144: ; preds = %140
%145 = load i32*, i32** %13
%146 = load i32, i32* %23
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds i32, i32* %145, i64 %147
%149 = load i32, i32* %148
store i32 %149, i32* %15
%150 = load i32*, i32** %13
%151 = load i32, i32* %21
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds i32, i32* %150, i64 %152
%154 = load i32, i32* %153
%155 = load i32*, i32** %13
%156 = load i32, i32* %23
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds i32, i32* %155, i64 %157
store i32 %154, i32* %158
%159 = load i32, i32* %15
%160 = load i32*, i32** %13
%161 = load i32, i32* %21
%162 = sext i32 %161 to i64
%163 = getelementptr inbounds i32, i32* %160, i64 %162
store i32 %159, i32* %163
br label %164
164: ; preds = %144, %140
%165 = load i32*, i32** %13
%166 = load i32, i32* %23
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds i32, i32* %165, i64 %167
%169 = load i32, i32* %168
store i32 %169, i32* %24
%170 = load i32, i32* %11
%171 = load i32, i32* %24
%172 = mul nsw i32 %170, %171
%173 = load i32, i32* %23
%174 = add nsw i32 %172, %173
store i32 %174, i32* %26
%175 = load double*, double** %8
%176 = load i32, i32* %26
%177 = sext i32 %176 to i64
%178 = getelementptr inbounds double, double* %175, i64 %177
%179 = load double, double* %178
store double %179, double* %36
%180 = load i32, i32* %23
%181 = add nsw i32 %180, 1
store i32 %181, i32* %25
%182 = load i32, i32* %25
store i32 %182, i32* %14
br label %183
183: ; preds = %249, %164
%184 = load i32, i32* %14
%185 = load i32, i32* %11
%186 = icmp slt i32 %184, %185
br i1 %186, label %187, label %252
187: ; preds = %183
%188 = load i32*, i32** %13
%189 = load i32, i32* %14
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds i32, i32* %188, i64 %190
%192 = load i32, i32* %191
store i32 %192, i32* %17
%193 = load i32, i32* %11
%194 = load i32, i32* %17
%195 = mul nsw i32 %193, %194
%196 = load i32, i32* %23
%197 = add nsw i32 %195, %196
store i32 %197, i32* %19
%198 = load double*, double** %8
%199 = load i32, i32* %19
%200 = sext i32 %199 to i64
%201 = getelementptr inbounds double, double* %198, i64 %200
%202 = load double, double* %201
%203 = fneg double %202
%204 = load double, double* %36
%205 = fdiv double %203, %204
store double %205, double* %31
%206 = load double, double* %31
%207 = fneg double %206
%208 = load double*, double** %8
%209 = load i32, i32* %19
%210 = sext i32 %209 to i64
%211 = getelementptr inbounds double, double* %208, i64 %210
store double %207, double* %211
%212 = load i32, i32* %11
%213 = load i32, i32* %17
%214 = mul nsw i32 %212, %213
store i32 %214, i32* %28
%215 = load i32, i32* %11
%216 = load i32, i32* %24
%217 = mul nsw i32 %215, %216
store i32 %217, i32* %29
%218 = load i32, i32* %25
store i32 %218, i32* %15
br label %219
219: ; preds = %245, %187
%220 = load i32, i32* %15
%221 = load i32, i32* %11
%222 = icmp slt i32 %220, %221
br i1 %222, label %223, label %248
223: ; preds = %219
%224 = load i32, i32* %28
%225 = load i32, i32* %15
%226 = add nsw i32 %224, %225
store i32 %226, i32* %18
%227 = load double*, double** %8
%228 = load i32, i32* %18
%229 = sext i32 %228 to i64
%230 = getelementptr inbounds double, double* %227, i64 %229
%231 = load double, double* %230
%232 = load double, double* %31
%233 = load double*, double** %8
%234 = load i32, i32* %29
%235 = load i32, i32* %15
%236 = add nsw i32 %234, %235
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds double, double* %233, i64 %237
%239 = load double, double* %238
%240 = call double @llvm.fmuladd.f64(double %232, double %239, double %231)
%241 = load double*, double** %8
%242 = load i32, i32* %18
%243 = sext i32 %242 to i64
%244 = getelementptr inbounds double, double* %241, i64 %243
store double %240, double* %244
br label %245
245: ; preds = %223
%246 = load i32, i32* %15
%247 = add nsw i32 %246, 1
store i32 %247, i32* %15
br label %219
248: ; preds = %219
br label %249
249: ; preds = %248
%250 = load i32, i32* %14
%251 = add nsw i32 %250, 1
store i32 %251, i32* %14
br label %183
252: ; preds = %183
br label %253
253: ; preds = %252
%254 = load i32, i32* %23
%255 = add nsw i32 %254, 1
store i32 %255, i32* %23
br label %92
256: ; preds = %92
%257 = load i32, i32* %11
%258 = load i32*, i32** %13
%259 = load i32, i32* %11
%260 = sub nsw i32 %259, 1
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds i32, i32* %258, i64 %261
%263 = load i32, i32* %262
%264 = mul nsw i32 %257, %263
%265 = load i32, i32* %11
%266 = add nsw i32 %264, %265
%267 = sub nsw i32 %266, 1
store i32 %267, i32* %27
%268 = load double*, double** %8
%269 = load i32, i32* %27
%270 = sext i32 %269 to i64
%271 = getelementptr inbounds double, double* %268, i64 %270
%272 = load double, double* %271
%273 = fcmp oeq double %272, 0.000000e+00
br i1 %273, label %274, label %276
274: ; preds = %256
%275 = 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 %431
276: ; preds = %256
br label %277
277: ; preds = %276, %42
%278 = load i32*, i32** %13
%279 = getelementptr inbounds i32, i32* %278, i64 0
%280 = load i32, i32* %279
store i32 %280, i32* %17
%281 = load double*, double** %9
%282 = load i32, i32* %17
%283 = sext i32 %282 to i64
%284 = getelementptr inbounds double, double* %281, i64 %283
%285 = load double, double* %284
%286 = load double*, double** %10
%287 = getelementptr inbounds double, double* %286, i64 0
store double %285, double* %287
store i32 1, i32* %14
br label %288
288: ; preds = %335, %277
%289 = load i32, i32* %14
%290 = load i32, i32* %11
%291 = icmp slt i32 %289, %290
br i1 %291, label %292, label %338
292: ; preds = %288
%293 = load i32*, i32** %13
%294 = load i32, i32* %14
%295 = sext i32 %294 to i64
%296 = getelementptr inbounds i32, i32* %293, i64 %295
%297 = load i32, i32* %296
store i32 %297, i32* %17
%298 = load i32, i32* %11
%299 = load i32, i32* %17
%300 = mul nsw i32 %298, %299
store i32 %300, i32* %18
store double 0.000000e+00, double* %37
store i32 0, i32* %15
br label %301
301: ; preds = %320, %292
%302 = load i32, i32* %15
%303 = load i32, i32* %14
%304 = icmp slt i32 %302, %303
br i1 %304, label %305, label %323
305: ; preds = %301
%306 = load double*, double** %8
%307 = load i32, i32* %18
%308 = sext i32 %307 to i64
%309 = getelementptr inbounds double, double* %306, i64 %308
%310 = load double, double* %309
%311 = load double*, double** %10
%312 = load i32, i32* %15
%313 = sext i32 %312 to i64
%314 = getelementptr inbounds double, double* %311, i64 %313
%315 = load double, double* %314
%316 = load double, double* %37
%317 = call double @llvm.fmuladd.f64(double %310, double %315, double %316)
store double %317, double* %37
%318 = load i32, i32* %18
%319 = add nsw i32 %318, 1
store i32 %319, i32* %18
br label %320
320: ; preds = %305
%321 = load i32, i32* %15
%322 = add nsw i32 %321, 1
store i32 %322, i32* %15
br label %301
323: ; preds = %301
%324 = load double*, double** %9
%325 = load i32, i32* %17
%326 = sext i32 %325 to i64
%327 = getelementptr inbounds double, double* %324, i64 %326
%328 = load double, double* %327
%329 = load double, double* %37
%330 = fsub double %328, %329
%331 = load double*, double** %10
%332 = load i32, i32* %14
%333 = sext i32 %332 to i64
%334 = getelementptr inbounds double, double* %331, i64 %333
store double %330, double* %334
br label %335
335: ; preds = %323
%336 = load i32, i32* %14
%337 = add nsw i32 %336, 1
store i32 %337, i32* %14
br label %288
338: ; preds = %288
%339 = load i32, i32* %11
%340 = load i32*, i32** %13
%341 = load i32, i32* %11
%342 = sub nsw i32 %341, 1
%343 = sext i32 %342 to i64
%344 = getelementptr inbounds i32, i32* %340, i64 %343
%345 = load i32, i32* %344
%346 = mul nsw i32 %339, %345
%347 = load i32, i32* %11
%348 = add nsw i32 %346, %347
%349 = sub nsw i32 %348, 1
store i32 %349, i32* %20
%350 = load double*, double** %10
%351 = load i32, i32* %11
%352 = sub nsw i32 %351, 1
%353 = sext i32 %352 to i64
%354 = getelementptr inbounds double, double* %350, i64 %353
%355 = load double, double* %354
%356 = load double*, double** %8
%357 = load i32, i32* %20
%358 = sext i32 %357 to i64
%359 = getelementptr inbounds double, double* %356, i64 %358
%360 = load double, double* %359
%361 = fdiv double %355, %360
%362 = load double*, double** %10
%363 = load i32, i32* %11
%364 = sub nsw i32 %363, 1
%365 = sext i32 %364 to i64
%366 = getelementptr inbounds double, double* %362, i64 %365
store double %361, double* %366
store i32 1, i32* %22
br label %367
367: ; preds = %427, %338
%368 = load i32, i32* %22
%369 = load i32, i32* %11
%370 = icmp slt i32 %368, %369
br i1 %370, label %371, label %430
371: ; preds = %367
%372 = load i32, i32* %30
%373 = load i32, i32* %22
%374 = sub nsw i32 %372, %373
store i32 %374, i32* %14
%375 = load i32*, i32** %13
%376 = load i32, i32* %14
%377 = sext i32 %376 to i64
%378 = getelementptr inbounds i32, i32* %375, i64 %377
%379 = load i32, i32* %378
store i32 %379, i32* %17
%380 = load i32, i32* %11
%381 = load i32, i32* %17
%382 = mul nsw i32 %380, %381
store i32 %382, i32* %28
store double 0.000000e+00, double* %37
%383 = load i32, i32* %14
%384 = add nsw i32 %383, 1
store i32 %384, i32* %15
br label %385
385: ; preds = %404, %371
%386 = load i32, i32* %15
%387 = load i32, i32* %11
%388 = icmp slt i32 %386, %387
br i1 %388, label %389, label %407
389: ; preds = %385
%390 = load double*, double** %8
%391 = load i32, i32* %28
%392 = load i32, i32* %15
%393 = add nsw i32 %391, %392
%394 = sext i32 %393 to i64
%395 = getelementptr inbounds double, double* %390, i64 %394
%396 = load double, double* %395
%397 = load double*, double** %10
%398 = load i32, i32* %15
%399 = sext i32 %398 to i64
%400 = getelementptr inbounds double, double* %397, i64 %399
%401 = load double, double* %400
%402 = load double, double* %37
%403 = call double @llvm.fmuladd.f64(double %396, double %401, double %402)
store double %403, double* %37
br label %404
404: ; preds = %389
%405 = load i32, i32* %15
%406 = add nsw i32 %405, 1
store i32 %406, i32* %15
br label %385
407: ; preds = %385
%408 = load double*, double** %10
%409 = load i32, i32* %14
%410 = sext i32 %409 to i64
%411 = getelementptr inbounds double, double* %408, i64 %410
%412 = load double, double* %411
%413 = load double, double* %37
%414 = fsub double %412, %413
%415 = load double*, double** %8
%416 = load i32, i32* %28
%417 = load i32, i32* %14
%418 = add nsw i32 %416, %417
%419 = sext i32 %418 to i64
%420 = getelementptr inbounds double, double* %415, i64 %419
%421 = load double, double* %420
%422 = fdiv double %414, %421
%423 = load double*, double** %10
%424 = load i32, i32* %14
%425 = sext i32 %424 to i64
%426 = getelementptr inbounds double, double* %423, i64 %425
store double %422, double* %426
br label %427
427: ; preds = %407
%428 = load i32, i32* %22
%429 = add nsw i32 %428, 1
store i32 %429, i32* %22
br label %367
430: ; preds = %367
store i32 0, i32* %7
br label %431
431: ; preds = %430, %274, %138, %79
%432 = load i32, i32* %7
ret i32 %432
}

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