void kernel_CPU_06_mod_3DRhoCte(float *gU0, float *gU1, float *gVorg, int nnoi,
                                int nnoj, int k0, int k1, float FATMDFX,
                                float FATMDFY, float FATMDFZ, float *W) {
  int index_X, index_Y;
  int stride = nnoi * nnoj;
  int index, k;
#pragma omp parallel for default(shared) private(index_X, index_Y, index, k)
  for (index_Y = 0; index_Y < nnoj; index_Y++)
    for (k = 0; k < k1 - k0; k++) {
      for (index_X = 0; index_X < nnoi; index_X++) {
        index = (index_Y * nnoi + index_X) + (k0 + k) * stride;
        if (gVorg[index] > 0.0f) {
          gU1[index] =
              2.0f * gU0[index] - gU1[index] +
              FATMDFX * gVorg[index] * gVorg[index] *
                  (+W[6] * (gU0[index - 6] + gU0[index + 6]) +
                   W[5] * (gU0[index - 5] + gU0[index + 5]) +
                   W[4] * (gU0[index - 4] + gU0[index + 4]) +
                   W[3] * (gU0[index - 3] + gU0[index + 3]) +
                   W[2] * (gU0[index - 2] + gU0[index + 2]) +
                   W[1] * (gU0[index - 1] + gU0[index + 1]) +
                   W[0] * gU0[index]) +
              FATMDFY * gVorg[index] * gVorg[index] *
                  (+W[6] * (gU0[index - 6 * nnoi] + gU0[index + 6 * nnoi]) +
                   W[5] * (gU0[index - 5 * nnoi] + gU0[index + 5 * nnoi]) +
                   W[4] * (gU0[index - 4 * nnoi] + gU0[index + 4 * nnoi]) +
                   W[3] * (gU0[index - 3 * nnoi] + gU0[index + 3 * nnoi]) +
                   W[2] * (gU0[index - 2 * nnoi] + gU0[index + 2 * nnoi]) +
                   W[1] * (gU0[index - nnoi] + gU0[index + nnoi]) +
                   W[0] * gU0[index]) +
              FATMDFZ * gVorg[index] * gVorg[index] *
                  (+W[6] * (gU0[index + 6 * stride] + gU0[index - 6 * stride]) +
                   W[5] * (gU0[index + 5 * stride] + gU0[index - 5 * stride]) +
                   W[4] * (gU0[index + 4 * stride] + gU0[index - 4 * stride]) +
                   W[3] * (gU0[index + 3 * stride] + gU0[index - 3 * stride]) +
                   W[2] * (gU0[index + 2 * stride] + gU0[index - 2 * stride]) +
                   W[1] * (gU0[index + stride] + gU0[index - stride]) +
                   W[0] * gU0[index]);
        }
      }
    }
}


x86 assembly (clang, -O0):
<kernel_CPU_06_mod_3DRhoCte>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movss  %xmm0,-0x28(%rbp)
movss  %xmm1,-0x2c(%rbp)
movss  %xmm2,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    16e6 <kernel_CPU_06_mod_3DRhoCte+0x5e6>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    0x10(%rbp),%ecx
sub    -0x24(%rbp),%ecx
cmp    %ecx,%eax
jge    16d3 <kernel_CPU_06_mod_3DRhoCte+0x5d3>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    16c0 <kernel_CPU_06_mod_3DRhoCte+0x5c0>
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x34(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    -0x44(%rbp),%ecx
imul   -0x3c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    16ad <kernel_CPU_06_mod_3DRhoCte+0x5ad>
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
addss  %xmm1,%xmm1
addss  %xmm0,%xmm1
movss  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    0x18(%rbp),%rax
movss  0x18(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x6,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x6,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mov    0x18(%rbp),%rax
movss  0x14(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x5,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm5
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x5,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm5
mulss  %xmm5,%xmm2
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  0x10(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x4,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mov    0x18(%rbp),%rax
movss  0xc(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x3,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  0x8(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mov    0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x2c(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mov    0x18(%rbp),%rax
movss  0x18(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x6,-0x1c(%rbp),%edx
sub    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x6,-0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mov    0x18(%rbp),%rax
movss  0x14(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x5,-0x1c(%rbp),%edx
sub    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm5
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x5,-0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm5
mulss  %xmm5,%xmm2
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  0x10(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
shl    $0x2,%edx
sub    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mov    0x18(%rbp),%rax
movss  0xc(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x3,-0x1c(%rbp),%edx
sub    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x3,-0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  0x8(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
shl    $0x1,%edx
sub    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mov    0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    0x18(%rbp),%rax
movss  0x18(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x6,-0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x6,-0x3c(%rbp),%edx
sub    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mov    0x18(%rbp),%rax
movss  0x14(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x5,-0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm5
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x5,-0x3c(%rbp),%edx
sub    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm5
mulss  %xmm5,%xmm2
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  0x10(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x3c(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x3c(%rbp),%edx
shl    $0x2,%edx
sub    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mov    0x18(%rbp),%rax
movss  0xc(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x3,-0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   $0x3,-0x3c(%rbp),%edx
sub    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  0x8(%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x3c(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x3c(%rbp),%edx
shl    $0x1,%edx
sub    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mov    0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm4
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp),%rax
movss  (%rax),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm4
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    16b2 <kernel_CPU_06_mod_3DRhoCte+0x5b2>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    116d <kernel_CPU_06_mod_3DRhoCte+0x6d>
jmp    16c5 <kernel_CPU_06_mod_3DRhoCte+0x5c5>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1155 <kernel_CPU_06_mod_3DRhoCte+0x55>
jmp    16d8 <kernel_CPU_06_mod_3DRhoCte+0x5d8>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1142 <kernel_CPU_06_mod_3DRhoCte+0x42>
pop    %rbp
ret

x86 assembly (clang, -O2):
<kernel_CPU_06_mod_3DRhoCte>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x170,%rsp
movss  %xmm2,-0x64(%rsp)
movss  %xmm1,-0x68(%rsp)
movss  %xmm0,-0x6c(%rsp)
mov    %ecx,-0x78(%rsp)
mov    %rdx,0x168(%rsp)
mov    %rsi,0x160(%rsp)
mov    %r8d,-0x7c(%rsp)
test   %r8d,%r8d
jle    19f2 <kernel_CPU_06_mod_3DRhoCte+0x8f2>
mov    0x1a8(%rsp),%esi
mov    -0x7c(%rsp),%ecx
mov    %ecx,%eax
mov    -0x78(%rsp),%edx
imul   %edx,%eax
sub    %r9d,%esi
mov    %esi,-0x74(%rsp)
mov    %r9d,%esi
imul   %ecx,%esi
mov    %esi,%ecx
imul   %edx,%ecx
lea    0x6(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x78(%rsp)
lea    -0x6(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x70(%rsp)
lea    0x5(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x68(%rsp)
lea    -0x5(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x60(%rsp)
lea    0x4(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x58(%rsp)
lea    -0x4(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x50(%rsp)
lea    0x3(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x48(%rsp)
lea    -0x3(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x40(%rsp)
lea    0x2(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x38(%rsp)
lea    -0x2(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x30(%rsp)
lea    0x1(%r9),%ebp
imul   %eax,%ebp
mov    %rbp,0x28(%rsp)
add    $0xffffffff,%r9d
imul   %eax,%r9d
mov    %r9,0x90(%rsp)
mov    %eax,%r15d
lea    -0x6(%rsi),%eax
imul   %edx,%eax
mov    %rax,0x20(%rsp)
lea    0x6(%rsi),%eax
imul   %edx,%eax
mov    %rax,0x18(%rsp)
lea    -0x5(%rsi),%eax
imul   %edx,%eax
mov    %rax,0x10(%rsp)
lea    0x5(%rsi),%eax
imul   %edx,%eax
mov    %rax,0x8(%rsp)
lea    -0x4(%rsi),%eax
imul   %edx,%eax
mov    %rax,(%rsp)
lea    0x4(%rsi),%eax
imul   %edx,%eax
mov    %rax,-0x8(%rsp)
lea    -0x3(%rsi),%eax
imul   %edx,%eax
mov    %rax,-0x10(%rsp)
lea    0x3(%rsi),%eax
imul   %edx,%eax
mov    %rax,-0x18(%rsp)
lea    -0x2(%rsi),%eax
imul   %edx,%eax
mov    %rax,-0x20(%rsp)
lea    0x2(%rsi),%eax
imul   %edx,%eax
mov    %rax,-0x28(%rsp)
lea    -0x1(%rsi),%eax
imul   %edx,%eax
mov    %rax,-0x30(%rsp)
add    $0x1,%esi
imul   %edx,%esi
mov    %rsi,0x88(%rsp)
xorps  %xmm0,%xmm0
mov    0x1b0(%rsp),%rbx
mov    %edx,%r9d
lea    -0x1(%rcx),%eax
mov    %rax,-0x38(%rsp)
lea    -0x2(%rcx),%eax
mov    %rax,-0x40(%rsp)
lea    -0x3(%rcx),%eax
mov    %rax,-0x48(%rsp)
lea    -0x4(%rcx),%eax
mov    %rax,-0x50(%rsp)
lea    -0x5(%rcx),%eax
mov    %rax,-0x58(%rsp)
mov    %rcx,0x80(%rsp)
lea    -0x6(%rcx),%eax
mov    %rax,-0x60(%rsp)
xor    %eax,%eax
mov    0x160(%rsp),%rsi
jmp    1398 <kernel_CPU_06_mod_3DRhoCte+0x298>
nopl   (%rax)
mov    0x158(%rsp),%rax
add    $0x1,%eax
add    %r9,-0x38(%rsp)
add    %r9,-0x40(%rsp)
add    %r9,-0x48(%rsp)
add    %r9,-0x50(%rsp)
add    %r9,-0x58(%rsp)
add    %r9,-0x60(%rsp)
add    %r9,0x80(%rsp)
add    %r9,0x20(%rsp)
add    %r9,0x18(%rsp)
add    %r9,0x10(%rsp)
add    %r9,0x8(%rsp)
add    %r9,(%rsp)
add    %r9,-0x8(%rsp)
add    %r9,-0x10(%rsp)
add    %r9,-0x18(%rsp)
add    %r9,-0x20(%rsp)
add    %r9,-0x28(%rsp)
add    %r9,-0x30(%rsp)
add    %r9,0x88(%rsp)
add    %r9,0x78(%rsp)
add    %r9,0x70(%rsp)
add    %r9,0x68(%rsp)
add    %r9,0x60(%rsp)
add    %r9,0x58(%rsp)
add    %r9,0x50(%rsp)
add    %r9,0x48(%rsp)
add    %r9,0x40(%rsp)
add    %r9,0x38(%rsp)
add    %r9,0x30(%rsp)
add    %r9,0x28(%rsp)
add    %r9,0x90(%rsp)
cmp    -0x7c(%rsp),%eax
je     19f2 <kernel_CPU_06_mod_3DRhoCte+0x8f2>
mov    %rax,0x158(%rsp)
cmpl   $0x0,-0x74(%rsp)
jle    12e0 <kernel_CPU_06_mod_3DRhoCte+0x1e0>
xor    %edx,%edx
mov    0x90(%rsp),%rax
mov    %rax,0x150(%rsp)
mov    0x28(%rsp),%rax
mov    %rax,0x148(%rsp)
mov    0x30(%rsp),%rax
mov    %rax,0x140(%rsp)
mov    0x38(%rsp),%rax
mov    %rax,0x138(%rsp)
mov    0x40(%rsp),%rax
mov    %rax,0x130(%rsp)
mov    0x48(%rsp),%rax
mov    %rax,0x128(%rsp)
mov    0x50(%rsp),%rax
mov    %rax,0x120(%rsp)
mov    0x58(%rsp),%rax
mov    %rax,0x118(%rsp)
mov    0x60(%rsp),%rax
mov    %rax,0x110(%rsp)
mov    0x68(%rsp),%rax
mov    %rax,0x108(%rsp)
mov    0x70(%rsp),%rax
mov    %rax,0x100(%rsp)
mov    0x78(%rsp),%rax
mov    %rax,0xf8(%rsp)
mov    0x88(%rsp),%rax
mov    %rax,0xf0(%rsp)
mov    -0x30(%rsp),%rax
mov    %rax,0xe8(%rsp)
mov    -0x28(%rsp),%rax
mov    %rax,0xe0(%rsp)
mov    -0x20(%rsp),%rax
mov    %rax,0xd8(%rsp)
mov    -0x18(%rsp),%rax
mov    %rax,0xd0(%rsp)
mov    -0x10(%rsp),%rax
mov    %rax,0xc8(%rsp)
mov    -0x8(%rsp),%rax
mov    %rax,0xc0(%rsp)
mov    (%rsp),%rax
mov    %rax,0xb8(%rsp)
mov    0x8(%rsp),%rax
mov    %rax,0xb0(%rsp)
mov    0x10(%rsp),%rax
mov    %rax,0xa8(%rsp)
mov    0x18(%rsp),%rax
mov    %rax,0xa0(%rsp)
mov    0x20(%rsp),%rax
mov    %rax,0x98(%rsp)
mov    0x80(%rsp),%r11
mov    -0x60(%rsp),%r13
mov    -0x58(%rsp),%r14
mov    -0x50(%rsp),%r8
mov    -0x48(%rsp),%rbp
mov    -0x40(%rsp),%rax
mov    -0x38(%rsp),%rcx
jmp    1606 <kernel_CPU_06_mod_3DRhoCte+0x506>
cs nopw 0x0(%rax,%rax,1)
nop
mov    -0x70(%rsp),%edx
add    $0x1,%edx
add    %r15,%rcx
add    %r15,%rax
add    %r15,%rbp
add    %r15,%r8
add    %r15,%r14
add    %r15,%r13
add    %r15,%r11
add    %r15,0x98(%rsp)
add    %r15,0xa0(%rsp)
add    %r15,0xa8(%rsp)
add    %r15,0xb0(%rsp)
add    %r15,0xb8(%rsp)
add    %r15,0xc0(%rsp)
add    %r15,0xc8(%rsp)
add    %r15,0xd0(%rsp)
add    %r15,0xd8(%rsp)
add    %r15,0xe0(%rsp)
add    %r15,0xe8(%rsp)
add    %r15,0xf0(%rsp)
add    %r15,0xf8(%rsp)
add    %r15,0x100(%rsp)
add    %r15,0x108(%rsp)
add    %r15,0x110(%rsp)
add    %r15,0x118(%rsp)
add    %r15,0x120(%rsp)
add    %r15,0x128(%rsp)
add    %r15,0x130(%rsp)
add    %r15,0x138(%rsp)
add    %r15,0x140(%rsp)
add    %r15,0x148(%rsp)
add    %r15,0x150(%rsp)
cmp    -0x74(%rsp),%edx
je     12e0 <kernel_CPU_06_mod_3DRhoCte+0x1e0>
mov    %edx,-0x70(%rsp)
cmpl   $0x0,-0x78(%rsp)
jle    1520 <kernel_CPU_06_mod_3DRhoCte+0x420>
xor    %r10d,%r10d
jmp    162d <kernel_CPU_06_mod_3DRhoCte+0x52d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r10
cmp    %r10,%r9
je     1520 <kernel_CPU_06_mod_3DRhoCte+0x420>
lea    (%r11,%r10,1),%edx
movslq %edx,%r12
mov    0x168(%rsp),%rdx
movss  (%rdx,%r12,4),%xmm7
ucomiss %xmm0,%xmm7
jbe    1620 <kernel_CPU_06_mod_3DRhoCte+0x520>
movss  (%rdi,%r12,4),%xmm15
movaps %xmm15,%xmm3
addss  %xmm15,%xmm3
subss  (%rsi,%r12,4),%xmm3
movaps %xmm7,%xmm5
mulss  -0x6c(%rsp),%xmm5
mulss  %xmm7,%xmm5
movss  0x18(%rbx),%xmm9
lea    0x0(%r13,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm2
lea    0x6(%r11,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm2
movss  0x14(%rbx),%xmm10
lea    (%r14,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm4
lea    0x5(%r11,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm4
mulss  %xmm10,%xmm4
mulss  %xmm9,%xmm2
addss  %xmm4,%xmm2
movss  0x10(%rbx),%xmm11
lea    (%r8,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm4
lea    0x4(%r11,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm4
mulss  %xmm11,%xmm4
addss  %xmm2,%xmm4
movss  0xc(%rbx),%xmm12
lea    0x0(%rbp,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm2
lea    0x3(%r11,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm2
mulss  %xmm12,%xmm2
addss  %xmm4,%xmm2
movss  0x4(%rbx),%xmm13
movss  0x8(%rbx),%xmm14
lea    (%rax,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm4
lea    0x2(%r11,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm4
mulss  %xmm14,%xmm4
addss  %xmm2,%xmm4
lea    (%rcx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm6
lea    0x1(%r11,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm6
mulss  %xmm13,%xmm6
addss  %xmm4,%xmm6
mulss  (%rbx),%xmm15
addss  %xmm15,%xmm6
mulss  %xmm5,%xmm6
addss  %xmm3,%xmm6
movaps %xmm7,%xmm5
mulss  -0x68(%rsp),%xmm5
mulss  %xmm7,%xmm5
mov    0x98(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm2
mov    0xa0(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm2
mov    0xa8(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm3
mov    0xb0(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm3
mulss  %xmm10,%xmm3
mulss  %xmm9,%xmm2
mov    0xb8(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm4
mov    0xc0(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm4
addss  %xmm3,%xmm2
mulss  %xmm11,%xmm4
addss  %xmm2,%xmm4
mov    0xc8(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm3
mov    0xd0(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm3
mulss  %xmm12,%xmm3
addss  %xmm4,%xmm3
mov    0xd8(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm2
mov    0xe0(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm2
mulss  %xmm14,%xmm2
addss  %xmm3,%xmm2
mov    0xe8(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm3
mov    0xf0(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm3
mulss  %xmm13,%xmm3
mov    0xf8(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm4
mov    0x100(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm4
addss  %xmm2,%xmm3
mov    0x108(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm8
mov    0x110(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm8
addss  %xmm15,%xmm3
mov    0x118(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm2
mov    0x120(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm2
mulss  %xmm5,%xmm3
mov    0x128(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm5
mov    0x130(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm5
mov    0x138(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm1
mov    0x140(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm1
mov    0x148(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm0
mov    0x150(%rsp),%rdx
lea    (%rdx,%r10,1),%edx
movslq %edx,%rdx
addss  (%rdi,%rdx,4),%xmm0
addss  %xmm6,%xmm3
movaps %xmm7,%xmm6
mulss  -0x64(%rsp),%xmm6
mulss  %xmm7,%xmm6
mulss  %xmm10,%xmm8
mulss  %xmm9,%xmm4
addss  %xmm8,%xmm4
mulss  %xmm11,%xmm2
addss  %xmm4,%xmm2
mulss  %xmm12,%xmm5
addss  %xmm2,%xmm5
mulss  %xmm14,%xmm1
addss  %xmm5,%xmm1
mulss  %xmm13,%xmm0
addss  %xmm1,%xmm0
addss  %xmm15,%xmm0
mulss  %xmm6,%xmm0
addss  %xmm3,%xmm0
movss  %xmm0,(%rsi,%r12,4)
xorps  %xmm0,%xmm0
jmp    1620 <kernel_CPU_06_mod_3DRhoCte+0x520>
add    $0x170,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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