void sor_1_iter(float *grad, float *err, float *dflow, float *flow,
                float *prev_flow, int nx, int ny, float om, float *scale,
                float l1, float l2, float l3, float *s1, float *s2, float *s3) {
  int index, i, j;
  float new_flow, ts2;
  float psi2(), psi();
  for (i = 1; i < ny - 1; i++) {
    for (j = 1 + (i % 2); j < nx - 1; j = j + 2) {
      index = (i * nx) + j;
      new_flow = flow[index] + dflow[index];
      ts2 = s2[index];
      ts2 = 2.0 * ts2 * ts2;
      dflow[index] =
          dflow[index] -
          ((om * scale[index]) *
           ((l1 * grad[index] * psi(err[index], s1[index])) +
            (l2 *
             (psi2(new_flow - (flow[index - nx] + dflow[index - nx]), ts2) +
              psi2(new_flow - (flow[index - 1] + dflow[index - 1]), ts2) +
              psi2(new_flow - (flow[index + nx] + dflow[index + nx]), ts2) +
              psi2(new_flow - (flow[index + 1] + dflow[index + 1]), ts2))) +
            (l3 * psi(new_flow - prev_flow[index], s3[index]))));
    }
  }
  for (i = 1; i < ny - 1; i++) {
    for (j = 2 + (i % 2); j < nx - 1; j = j + 2) {
      index = (i * nx) + j;
      new_flow = flow[index] + dflow[index];
      ts2 = s2[index];
      ts2 = 2.0 * ts2 * ts2;
      dflow[index] =
          dflow[index] -
          ((om * scale[index]) *
           ((l1 * grad[index] * psi(err[index], s1[index])) +
            (l2 *
             (psi2(new_flow - (flow[index - nx] + dflow[index - nx]), ts2) +
              psi2(new_flow - (flow[index - 1] + dflow[index - 1]), ts2) +
              psi2(new_flow - (flow[index + nx] + dflow[index + nx]), ts2) +
              psi2(new_flow - (flow[index + 1] + dflow[index + 1]), ts2))) +
            (l3 * psi(new_flow - prev_flow[index], s3[index]))));
    }
  }
  fix_flow(dflow, nx, ny);
}


x86 assembly (clang, -O0):
<sor_1_iter>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movss  %xmm0,-0x30(%rbp)
movss  %xmm1,-0x34(%rbp)
movss  %xmm2,-0x38(%rbp)
movss  %xmm3,-0x3c(%rbp)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1486 <sor_1_iter+0x356>
mov    -0x44(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    $0x1,%edx
mov    %edx,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1473 <sor_1_iter+0x343>
mov    -0x44(%rbp),%eax
imul   -0x2c(%rbp),%eax
add    -0x48(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdfd(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x30(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x34(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x64(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1030 <psi@plt>
movss  %xmm0,-0x68(%rbp)
movss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x6c(%rbp)
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movss  %xmm0,-0x78(%rbp)
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movaps %xmm0,%xmm1
movss  -0x78(%rbp),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x74(%rbp)
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movaps %xmm0,%xmm1
movss  -0x74(%rbp),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x70(%rbp)
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movss  -0x70(%rbp),%xmm3
movss  -0x6c(%rbp),%xmm1
movss  -0x68(%rbp),%xmm2
movaps %xmm0,%xmm4
movss  -0x64(%rbp),%xmm0
addss  %xmm4,%xmm3
mulss  %xmm3,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x5c(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x4c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1030 <psi@plt>
movss  -0x60(%rbp),%xmm2
movss  -0x5c(%rbp),%xmm3
movss  -0x58(%rbp),%xmm1
movaps %xmm0,%xmm4
movss  -0x54(%rbp),%xmm0
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x48(%rbp)
jmp    11a3 <sor_1_iter+0x73>
jmp    1478 <sor_1_iter+0x348>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1181 <sor_1_iter+0x51>
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    17c2 <sor_1_iter+0x692>
mov    -0x44(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    $0x2,%edx
mov    %edx,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    17af <sor_1_iter+0x67f>
mov    -0x44(%rbp),%eax
imul   -0x2c(%rbp),%eax
add    -0x48(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xaf1(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x80(%rbp)
movss  -0x30(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x7c(%rbp)
movss  -0x34(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1030 <psi@plt>
movss  %xmm0,-0x90(%rbp)
ff
movss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x94(%rbp)
ff
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movss  %xmm0,-0xa0(%rbp)
ff
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movaps %xmm0,%xmm1
movss  -0xa0(%rbp),%xmm0
ff
addss  %xmm1,%xmm0
movss  %xmm0,-0x9c(%rbp)
ff
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movaps %xmm0,%xmm1
movss  -0x9c(%rbp),%xmm0
ff
addss  %xmm1,%xmm0
movss  %xmm0,-0x98(%rbp)
ff
movss  -0x4c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x50(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
movss  -0x98(%rbp),%xmm3
ff
movss  -0x94(%rbp),%xmm1
ff
movss  -0x90(%rbp),%xmm2
ff
movaps %xmm0,%xmm4
movss  -0x8c(%rbp),%xmm0
ff
addss  %xmm4,%xmm3
mulss  %xmm3,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x84(%rbp)
ff
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x88(%rbp)
ff
movss  -0x4c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1030 <psi@plt>
movss  -0x88(%rbp),%xmm2
ff
movss  -0x84(%rbp),%xmm3
ff
movss  -0x80(%rbp),%xmm1
movaps %xmm0,%xmm4
movss  -0x7c(%rbp),%xmm0
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x48(%rbp)
jmp    14af <sor_1_iter+0x37f>
jmp    17b4 <sor_1_iter+0x684>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    148d <sor_1_iter+0x35d>
mov    -0x18(%rbp),%rdi
mov    -0x2c(%rbp),%esi
mov    0x10(%rbp),%edx
mov    $0x0,%al
call   1040 <fix_flow@plt>
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sor_1_iter>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
movss  %xmm2,0x74(%rsp)
mov    %r9,0xb0(%rsp)
mov    %rdx,0x80(%rsp)
mov    0x140(%rsp),%eax
mov    %rax,0x88(%rsp)
cmp    $0x3,%eax
jl     18cc <sor_1_iter+0x79c>
movaps %xmm0,%xmm4
mov    %rcx,%r10
mov    0x160(%rsp),%rbp
mov    0x158(%rsp),%r15
mov    0x150(%rsp),%r14
mov    0x148(%rsp),%rax
mov    %rax,0xa8(%rsp)
mov    0x88(%rsp),%rax
add    $0xffffffff,%eax
mov    0xb0(%rsp),%rcx
lea    -0x1(%rcx),%edx
movslq %ecx,%rbx
mov    %edx,0x64(%rsp)
movslq %edx,%rcx
mov    %rcx,0xc8(%rsp)
mov    %eax,%eax
mov    %rax,0xc0(%rsp)
mov    %rbp,0xd0(%rsp)
lea    0x0(%rbp,%rbx,4),%r9
lea    0x0(,%rbx,4),%rax
mov    %rax,0xb8(%rsp)
mov    %r8,0xf0(%rsp)
lea    (%r8,%rbx,4),%r8
mov    %r14,0xd8(%rsp)
lea    (%r14,%rbx,4),%r11
lea    (%r10,%rbx,8),%r14
mov    0x80(%rsp),%r13
lea    0x0(%r13,%rbx,8),%r12
lea    0x4(%r10,%rbx,4),%rcx
lea    0x0(%r13,%rbx,4),%rbp
mov    %r15,0xe8(%rsp)
lea    (%r15,%rbx,4),%rax
mov    0xa8(%rsp),%rdx
lea    (%rdx,%rbx,4),%rdx
mov    %rdi,0xf8(%rsp)
lea    (%rdi,%rbx,4),%rdi
mov    %rsi,0x100(%rsp)
mov    %rbx,0xe0(%rsp)
lea    (%rsi,%rbx,4),%rsi
mov    %r9,%rbx
mov    $0x1,%r15d
mov    %r10,0x78(%rsp)
mov    %r10,%r9
movss  %xmm3,0x70(%rsp)
movss  %xmm1,0x6c(%rsp)
movss  %xmm0,0x68(%rsp)
jmp    12f3 <sor_1_iter+0x1c3>
nopl   0x0(%rax)
mov    0x90(%rsp),%r15
add    $0x1,%r15
mov    0xb8(%rsp),%rbx
mov    0x48(%rsp),%r8
add    %rbx,%r8
mov    0x38(%rsp),%r12
add    %rbx,%r12
add    %rbx,%r11
mov    0x40(%rsp),%r14
add    %rbx,%r14
mov    0x30(%rsp),%r13
add    %rbx,%r13
mov    0x50(%rsp),%r9
add    %rbx,%r9
mov    0x58(%rsp),%r10
add    %rbx,%r10
add    %rbx,%rcx
add    %rbx,%rbp
add    %rbx,%rax
add    %rbx,%rdx
add    %rbx,%rdi
add    %rbx,%rsi
cmp    0xc0(%rsp),%r15
mov    %r8,%rbx
mov    %r12,%r8
mov    %r13,%r12
mov    %r10,%r13
je     1554 <sor_1_iter+0x424>
mov    %r12,0x30(%rsp)
mov    %r8,0x38(%rsp)
mov    %r14,0x40(%rsp)
mov    %rbx,0x48(%rsp)
mov    %r9,0x50(%rsp)
mov    %r13,0x58(%rsp)
mov    %r15,0x90(%rsp)
mov    %r15d,%ebx
and    $0x1,%ebx
add    $0x1,%ebx
cmp    0x64(%rsp),%ebx
jge    1280 <sor_1_iter+0x150>
mov    %r11,0xa0(%rsp)
mov    %rax,0x98(%rsp)
nopl   0x0(%rax)
movss  0x0(%rbp,%rbx,4),%xmm0
movss  %xmm0,0x2c(%rsp)
movss  -0x4(%rcx,%rbx,4),%xmm2
addss  %xmm0,%xmm2
movss  %xmm2,0x10(%rsp)
movss  (%rax,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movaps %xmm0,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsd2ss %xmm2,%xmm0
movss  %xmm0,0x18(%rsp)
movss  (%rdi,%rbx,4),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,0x24(%rsp)
movss  (%rsi,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movss  (%r11,%rbx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movss  (%rdx,%rbx,4),%xmm2
mulss  %xmm4,%xmm2
movss  %xmm2,0x28(%rsp)
mov    $0x2,%al
mov    %rdx,%r15
mov    %rdi,%r13
mov    %rsi,%r12
mov    %rcx,%r14
call   1030 <psi@plt>
mov    0x50(%rsp),%rax
movss  (%rax,%rbx,4),%xmm1
mov    0x58(%rsp),%rax
addss  (%rax,%rbx,4),%xmm1
movss  %xmm0,0x14(%rsp)
movss  0x10(%rsp),%xmm0
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  0x18(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  %xmm1,0x18(%rsp)
mov    $0x2,%al
call   1050 <psi2@plt>
movss  %xmm0,0xc(%rsp)
movss  -0x8(%r14,%rbx,4),%xmm0
addss  -0x4(%rbp,%rbx,4),%xmm0
movss  0x10(%rsp),%xmm1
subss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movsd  0x18(%rsp),%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
addss  0xc(%rsp),%xmm0
movss  %xmm0,0xc(%rsp)
mov    0x40(%rsp),%rax
movss  (%rax,%rbx,4),%xmm0
mov    0x30(%rsp),%rax
addss  (%rax,%rbx,4),%xmm0
movss  0x10(%rsp),%xmm1
subss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movsd  0x18(%rsp),%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
addss  0xc(%rsp),%xmm0
movss  %xmm0,0xc(%rsp)
movss  (%r14,%rbx,4),%xmm0
addss  0x4(%rbp,%rbx,4),%xmm0
movss  0x10(%rsp),%xmm1
subss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movsd  0x18(%rsp),%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
addss  0xc(%rsp),%xmm0
mulss  0x74(%rsp),%xmm0
movss  0x14(%rsp),%xmm2
mulss  0x24(%rsp),%xmm2
mov    0x38(%rsp),%rax
movss  0x10(%rsp),%xmm1
subss  (%rax,%rbx,4),%xmm1
addss  %xmm0,%xmm2
movss  %xmm2,0x14(%rsp)
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
mov    0x48(%rsp),%rax
movss  (%rax,%rbx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1030 <psi@plt>
mov    %r14,%rcx
mov    %r12,%rsi
mov    %r13,%rdi
mov    %r15,%rdx
mov    0x98(%rsp),%rax
mov    0xa0(%rsp),%r11
movss  0x68(%rsp),%xmm4
movss  0x6c(%rsp),%xmm1
mulss  0x70(%rsp),%xmm0
addss  0x14(%rsp),%xmm0
mulss  0x28(%rsp),%xmm0
movss  0x2c(%rsp),%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,0x0(%rbp,%rbx,4)
add    $0x2,%rbx
cmp    0xc8(%rsp),%rbx
jl     1340 <sor_1_iter+0x210>
jmp    1280 <sor_1_iter+0x150>
cmpl   $0x3,0x88(%rsp)
03
mov    0x78(%rsp),%r13
mov    0x100(%rsp),%r12
mov    0xf8(%rsp),%r9
mov    0xf0(%rsp),%rdx
jl     18cc <sor_1_iter+0x79c>
mov    0xe0(%rsp),%rbx
mov    0xd0(%rsp),%rdi
lea    (%rdi,%rbx,4),%r11
lea    (%rdx,%rbx,4),%r10
mov    0xd8(%rsp),%rdx
lea    (%rdx,%rbx,4),%rcx
lea    0x0(%r13,%rbx,8),%r15
mov    0x80(%rsp),%r14
lea    (%r14,%rbx,8),%rax
lea    0x0(%r13,%rbx,4),%r8
lea    (%r14,%rbx,4),%rbp
mov    0xe8(%rsp),%rdx
lea    (%rdx,%rbx,4),%rdi
mov    0xa8(%rsp),%rdx
lea    (%rdx,%rbx,4),%rdx
lea    (%r9,%rbx,4),%rsi
lea    (%r12,%rbx,4),%rbx
mov    $0x1,%r9d
mov    %r14,%r12
jmp    165f <sor_1_iter+0x52f>
nopl   0x0(%rax,%rax,1)
mov    0x90(%rsp),%r9
add    $0x1,%r9
mov    0xb8(%rsp),%rax
mov    0x40(%rsp),%r11
add    %rax,%r11
mov    0x58(%rsp),%r10
add    %rax,%r10
add    %rax,%rcx
mov    0x38(%rsp),%r15
add    %rax,%r15
mov    0x50(%rsp),%r14
add    %rax,%r14
mov    0x78(%rsp),%r13
add    %rax,%r13
mov    0x48(%rsp),%r12
add    %rax,%r12
mov    0x30(%rsp),%r8
add    %rax,%r8
add    %rax,%rbp
add    %rax,%rdi
add    %rax,%rdx
add    %rax,%rsi
add    %rax,%rbx
cmp    0xc0(%rsp),%r9
mov    %r14,%rax
je     18cc <sor_1_iter+0x79c>
mov    %r8,0x30(%rsp)
mov    %r15,0x38(%rsp)
mov    %r11,0x40(%rsp)
mov    %r12,0x48(%rsp)
mov    %rax,0x50(%rsp)
mov    %r10,0x58(%rsp)
mov    %r13,0x78(%rsp)
mov    %r9,0x90(%rsp)
mov    %r9d,%eax
and    $0x1,%eax
or     $0x2,%eax
cmp    0x64(%rsp),%eax
jge    15f0 <sor_1_iter+0x4c0>
mov    %eax,%r13d
mov    %rcx,0xa0(%rsp)
mov    %rdi,0x98(%rsp)
movss  0x0(%rbp,%r13,4),%xmm0
movss  %xmm0,0x2c(%rsp)
mov    0x30(%rsp),%r12
movss  (%r12,%r13,4),%xmm2
addss  %xmm0,%xmm2
movss  %xmm2,0x10(%rsp)
movss  (%rdi,%r13,4),%xmm0
cvtss2sd %xmm0,%xmm0
movaps %xmm0,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtsd2ss %xmm2,%xmm0
movss  %xmm0,0x18(%rsp)
movss  (%rsi,%r13,4),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,0x24(%rsp)
movss  (%rbx,%r13,4),%xmm0
cvtss2sd %xmm0,%xmm0
movss  (%rcx,%r13,4),%xmm1
cvtss2sd %xmm1,%xmm1
movss  (%rdx,%r13,4),%xmm2
mulss  %xmm4,%xmm2
movss  %xmm2,0x28(%rsp)
mov    $0x2,%al
mov    %rdx,%r15
mov    %rsi,%r14
call   1030 <psi@plt>
mov    0x78(%rsp),%rax
movss  (%rax,%r13,4),%xmm1
mov    0x48(%rsp),%rax
addss  (%rax,%r13,4),%xmm1
movss  %xmm0,0x14(%rsp)
movss  0x10(%rsp),%xmm0
subss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movss  0x18(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  %xmm1,0x18(%rsp)
mov    $0x2,%al
call   1050 <psi2@plt>
movss  %xmm0,0xc(%rsp)
movss  -0x4(%r12,%r13,4),%xmm0
addss  -0x4(%rbp,%r13,4),%xmm0
movss  0x10(%rsp),%xmm1
subss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movsd  0x18(%rsp),%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
addss  0xc(%rsp),%xmm0
movss  %xmm0,0xc(%rsp)
mov    0x38(%rsp),%rax
movss  (%rax,%r13,4),%xmm0
mov    0x50(%rsp),%rax
addss  (%rax,%r13,4),%xmm0
movss  0x10(%rsp),%xmm1
subss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movsd  0x18(%rsp),%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
addss  0xc(%rsp),%xmm0
movss  %xmm0,0xc(%rsp)
movss  0x4(%r12,%r13,4),%xmm0
addss  0x4(%rbp,%r13,4),%xmm0
movss  0x10(%rsp),%xmm1
subss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movsd  0x18(%rsp),%xmm1
mov    $0x2,%al
call   1050 <psi2@plt>
addss  0xc(%rsp),%xmm0
mulss  0x74(%rsp),%xmm0
movss  0x14(%rsp),%xmm2
mulss  0x24(%rsp),%xmm2
mov    0x58(%rsp),%rax
movss  0x10(%rsp),%xmm1
subss  (%rax,%r13,4),%xmm1
addss  %xmm0,%xmm2
movss  %xmm2,0x14(%rsp)
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
mov    0x40(%rsp),%rax
movss  (%rax,%r13,4),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1030 <psi@plt>
mov    %r14,%rsi
mov    %r15,%rdx
mov    0x98(%rsp),%rdi
mov    0xa0(%rsp),%rcx
movss  0x68(%rsp),%xmm4
movss  0x6c(%rsp),%xmm1
mulss  0x70(%rsp),%xmm0
addss  0x14(%rsp),%xmm0
mulss  0x28(%rsp),%xmm0
movss  0x2c(%rsp),%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,0x0(%rbp,%r13,4)
add    $0x2,%r13
cmp    0xc8(%rsp),%r13
jl     16b0 <sor_1_iter+0x580>
jmp    15f0 <sor_1_iter+0x4c0>
mov    0x80(%rsp),%rdi
mov    0xb0(%rsp),%rsi
mov    0x88(%rsp),%rdx
xor    %eax,%eax
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <fix_flow@plt>

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

LLVM IR (-O2):
define void @sor_1_iter(float* nocapture %0, float* nocapture %1, float* %2, float* nocapture %3, float* nocapture %4, i32 %5, i32 %6, float %7, float* nocapture %8, float %9, float %10, float %11, float* nocapture %12, float* nocapture %13, float* nocapture %14) {
%16 = add i32 %6, -1
%17 = icmp sgt i32 %6, 2
br i1 %17, label %18, label %221
18: ; preds = %15
%19 = add nsw i32 %5, -1
%20 = sext i32 %5 to i64
%21 = sext i32 %19 to i64
%22 = zext i32 %16 to i64
br label %29
23: ; preds = %123
br i1 %17, label %24, label %221
24: ; preds = %23
%25 = add nsw i32 %5, -1
%26 = sext i32 %5 to i64
%27 = sext i32 %25 to i64
%28 = zext i32 %16 to i64
br label %127
29: ; preds = %18, %123
%30 = phi i64 [ 1, %18 ], [ %124, %123 ]
%31 = phi i32 [ 1, %18 ], [ %125, %123 ]
%32 = trunc i64 %30 to i32
%33 = and i32 %32, 1
%34 = add nuw nsw i32 %33, 1
%35 = icmp slt i32 %34, %19
br i1 %35, label %36, label %123
36: ; preds = %29
%37 = and i32 %31, 1
%38 = add nuw nsw i32 %37, 1
%39 = zext i32 %38 to i64
%40 = mul nsw i64 %30, %20
br label %41
41: ; preds = %36, %41
%42 = phi i64 [ %39, %36 ], [ %121, %41 ]
%43 = add nsw i64 %42, %40
%44 = getelementptr inbounds float, float* %3, i64 %43
%45 = load float, float* %44
%46 = getelementptr inbounds float, float* %2, i64 %43
%47 = load float, float* %46
%48 = fadd float %45, %47
%49 = getelementptr inbounds float, float* %13, i64 %43
%50 = load float, float* %49
%51 = fpext float %50 to double
%52 = fmul double %51, 2.000000e+00
%53 = fmul double %52, %51
%54 = fptrunc double %53 to float
%55 = getelementptr inbounds float, float* %8, i64 %43
%56 = load float, float* %55
%57 = getelementptr inbounds float, float* %0, i64 %43
%58 = load float, float* %57
%59 = fmul float %58, %9
%60 = getelementptr inbounds float, float* %1, i64 %43
%61 = load float, float* %60
%62 = fpext float %61 to double
%63 = getelementptr inbounds float, float* %12, i64 %43
%64 = load float, float* %63
%65 = fpext float %64 to double
%66 = tail call float (double, double, ...) bitcast (float (...)* @psi to float (double, double, ...)*)(double %62, double %65)
%67 = sub nsw i64 %43, %20
%68 = getelementptr inbounds float, float* %3, i64 %67
%69 = load float, float* %68
%70 = getelementptr inbounds float, float* %2, i64 %67
%71 = load float, float* %70
%72 = fadd float %69, %71
%73 = fsub float %48, %72
%74 = fpext float %73 to double
%75 = fpext float %54 to double
%76 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %74, double %75)
%77 = add nsw i64 %43, -1
%78 = getelementptr inbounds float, float* %3, i64 %77
%79 = load float, float* %78
%80 = getelementptr inbounds float, float* %2, i64 %77
%81 = load float, float* %80
%82 = fadd float %79, %81
%83 = fsub float %48, %82
%84 = fpext float %83 to double
%85 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %84, double %75)
%86 = fadd float %76, %85
%87 = add nsw i64 %43, %20
%88 = getelementptr inbounds float, float* %3, i64 %87
%89 = load float, float* %88
%90 = getelementptr inbounds float, float* %2, i64 %87
%91 = load float, float* %90
%92 = fadd float %89, %91
%93 = fsub float %48, %92
%94 = fpext float %93 to double
%95 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %94, double %75)
%96 = fadd float %86, %95
%97 = add nsw i64 %43, 1
%98 = getelementptr inbounds float, float* %3, i64 %97
%99 = load float, float* %98
%100 = getelementptr inbounds float, float* %2, i64 %97
%101 = load float, float* %100
%102 = fadd float %99, %101
%103 = fsub float %48, %102
%104 = fpext float %103 to double
%105 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %104, double %75)
%106 = fadd float %96, %105
%107 = fmul float %106, %10
%108 = tail call float @llvm.fmuladd.f32(float %59, float %66, float %107)
%109 = getelementptr inbounds float, float* %4, i64 %43
%110 = load float, float* %109
%111 = fsub float %48, %110
%112 = fpext float %111 to double
%113 = getelementptr inbounds float, float* %14, i64 %43
%114 = load float, float* %113
%115 = fpext float %114 to double
%116 = tail call float (double, double, ...) bitcast (float (...)* @psi to float (double, double, ...)*)(double %112, double %115)
%117 = tail call float @llvm.fmuladd.f32(float %11, float %116, float %108)
%118 = fneg float %56
%119 = fmul float %118, %7
%120 = tail call float @llvm.fmuladd.f32(float %119, float %117, float %47)
store float %120, float* %46
%121 = add nuw i64 %42, 2
%122 = icmp slt i64 %121, %21
br i1 %122, label %41, label %123
123: ; preds = %41, %29
%124 = add nuw nsw i64 %30, 1
%125 = add nuw nsw i32 %31, 1
%126 = icmp eq i64 %124, %22
br i1 %126, label %23, label %29
127: ; preds = %24, %218
%128 = phi i64 [ 1, %24 ], [ %219, %218 ]
%129 = trunc i64 %128 to i32
%130 = and i32 %129, 1
%131 = or i32 %130, 2
%132 = icmp slt i32 %131, %25
br i1 %132, label %133, label %218
133: ; preds = %127
%134 = mul nsw i64 %128, %26
%135 = zext i32 %131 to i64
br label %136
136: ; preds = %133, %136
%137 = phi i64 [ %135, %133 ], [ %216, %136 ]
%138 = add nsw i64 %137, %134
%139 = getelementptr inbounds float, float* %3, i64 %138
%140 = load float, float* %139
%141 = getelementptr inbounds float, float* %2, i64 %138
%142 = load float, float* %141
%143 = fadd float %140, %142
%144 = getelementptr inbounds float, float* %13, i64 %138
%145 = load float, float* %144
%146 = fpext float %145 to double
%147 = fmul double %146, 2.000000e+00
%148 = fmul double %147, %146
%149 = fptrunc double %148 to float
%150 = getelementptr inbounds float, float* %8, i64 %138
%151 = load float, float* %150
%152 = getelementptr inbounds float, float* %0, i64 %138
%153 = load float, float* %152
%154 = fmul float %153, %9
%155 = getelementptr inbounds float, float* %1, i64 %138
%156 = load float, float* %155
%157 = fpext float %156 to double
%158 = getelementptr inbounds float, float* %12, i64 %138
%159 = load float, float* %158
%160 = fpext float %159 to double
%161 = tail call float (double, double, ...) bitcast (float (...)* @psi to float (double, double, ...)*)(double %157, double %160)
%162 = sub nsw i64 %138, %26
%163 = getelementptr inbounds float, float* %3, i64 %162
%164 = load float, float* %163
%165 = getelementptr inbounds float, float* %2, i64 %162
%166 = load float, float* %165
%167 = fadd float %164, %166
%168 = fsub float %143, %167
%169 = fpext float %168 to double
%170 = fpext float %149 to double
%171 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %169, double %170)
%172 = add nsw i64 %138, -1
%173 = getelementptr inbounds float, float* %3, i64 %172
%174 = load float, float* %173
%175 = getelementptr inbounds float, float* %2, i64 %172
%176 = load float, float* %175
%177 = fadd float %174, %176
%178 = fsub float %143, %177
%179 = fpext float %178 to double
%180 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %179, double %170)
%181 = fadd float %171, %180
%182 = add nsw i64 %138, %26
%183 = getelementptr inbounds float, float* %3, i64 %182
%184 = load float, float* %183
%185 = getelementptr inbounds float, float* %2, i64 %182
%186 = load float, float* %185
%187 = fadd float %184, %186
%188 = fsub float %143, %187
%189 = fpext float %188 to double
%190 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %189, double %170)
%191 = fadd float %181, %190
%192 = add nsw i64 %138, 1
%193 = getelementptr inbounds float, float* %3, i64 %192
%194 = load float, float* %193
%195 = getelementptr inbounds float, float* %2, i64 %192
%196 = load float, float* %195
%197 = fadd float %194, %196
%198 = fsub float %143, %197
%199 = fpext float %198 to double
%200 = tail call float (double, double, ...) bitcast (float (...)* @psi2 to float (double, double, ...)*)(double %199, double %170)
%201 = fadd float %191, %200
%202 = fmul float %201, %10
%203 = tail call float @llvm.fmuladd.f32(float %154, float %161, float %202)
%204 = getelementptr inbounds float, float* %4, i64 %138
%205 = load float, float* %204
%206 = fsub float %143, %205
%207 = fpext float %206 to double
%208 = getelementptr inbounds float, float* %14, i64 %138
%209 = load float, float* %208
%210 = fpext float %209 to double
%211 = tail call float (double, double, ...) bitcast (float (...)* @psi to float (double, double, ...)*)(double %207, double %210)
%212 = tail call float @llvm.fmuladd.f32(float %11, float %211, float %203)
%213 = fneg float %151
%214 = fmul float %213, %7
%215 = tail call float @llvm.fmuladd.f32(float %214, float %212, float %142)
store float %215, float* %141
%216 = add nuw i64 %137, 2
%217 = icmp slt i64 %216, %27
br i1 %217, label %136, label %218
218: ; preds = %136, %127
%219 = add nuw nsw i64 %128, 1
%220 = icmp eq i64 %219, %28
br i1 %220, label %221, label %127
221: ; preds = %218, %15, %23
%222 = tail call i32 (float*, i32, i32, ...) bitcast (i32 (...)* @fix_flow to i32 (float*, i32, i32, ...)*)(float* %2, i32 %5, i32 %6)
ret void
}