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_X = 0; index_X < nnoi; index_X++)
    for (k = 0; k < k1 - k0; k++) {
#pragma simd
      for (index_Y = 0; index_Y < nnoj; index_Y++) {
        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,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%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,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%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    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%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    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%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    $0x80,%rsp
movss  %xmm2,-0x5c(%rsp)
movss  %xmm1,-0x60(%rsp)
movss  %xmm0,-0x64(%rsp)
mov    %r8d,-0x74(%rsp)
test   %ecx,%ecx
jle    1692 <kernel_CPU_06_mod_3DRhoCte+0x592>
mov    0xb8(%rsp),%ebp
mov    -0x74(%rsp),%r8d
mov    %r8d,%r10d
imul   %ecx,%r10d
sub    %r9d,%ebp
mov    %ebp,-0x70(%rsp)
mov    %r9d,%ebp
imul   %r8d,%ebp
mov    %ebp,%ebx
imul   %ecx,%ebx
lea    -0x6(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x70(%rsp)
lea    0x6(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x68(%rsp)
lea    -0x5(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x60(%rsp)
lea    0x5(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x58(%rsp)
lea    -0x4(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x50(%rsp)
lea    0x4(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x48(%rsp)
lea    -0x3(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x40(%rsp)
lea    0x3(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x38(%rsp)
lea    -0x2(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x30(%rsp)
lea    0x2(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x28(%rsp)
lea    -0x1(%rbp),%eax
imul   %ecx,%eax
mov    %rax,0x20(%rsp)
add    $0x1,%ebp
imul   %ecx,%ebp
mov    %rbp,0x78(%rsp)
lea    0x6(%r9),%eax
imul   %r10d,%eax
mov    %rax,0x18(%rsp)
lea    -0x6(%r9),%eax
imul   %r10d,%eax
mov    %rax,0x10(%rsp)
lea    0x5(%r9),%eax
imul   %r10d,%eax
mov    %rax,0x8(%rsp)
lea    -0x5(%r9),%eax
imul   %r10d,%eax
mov    %rax,(%rsp)
lea    0x4(%r9),%eax
imul   %r10d,%eax
mov    %rax,-0x8(%rsp)
lea    -0x4(%r9),%eax
imul   %r10d,%eax
mov    %rax,-0x10(%rsp)
lea    0x3(%r9),%eax
imul   %r10d,%eax
mov    %rax,-0x18(%rsp)
lea    -0x3(%r9),%eax
imul   %r10d,%eax
mov    %rax,-0x20(%rsp)
lea    0x2(%r9),%eax
imul   %r10d,%eax
mov    %rax,-0x28(%rsp)
lea    -0x2(%r9),%eax
imul   %r10d,%eax
mov    %rax,-0x30(%rsp)
lea    0x1(%r9),%eax
imul   %r10d,%eax
mov    %rax,-0x38(%rsp)
add    $0xffffffff,%r9d
mov    %r10d,-0x6c(%rsp)
imul   %r10d,%r9d
mov    %r9,-0x50(%rsp)
xorps  %xmm0,%xmm0
mov    0xc0(%rsp),%r13
mov    %r8d,%eax
mov    %rax,-0x58(%rsp)
lea    -0x1(%rbx),%eax
mov    %rax,-0x40(%rsp)
lea    -0x2(%rbx),%eax
mov    %rax,-0x48(%rsp)
lea    -0x3(%rbx),%r10d
lea    -0x4(%rbx),%r8d
lea    -0x5(%rbx),%r11d
lea    -0x6(%rbx),%r14d
movl   $0x0,-0x7c(%rsp)
jmp    12d3 <kernel_CPU_06_mod_3DRhoCte+0x1d3>
nopw   0x0(%rax,%rax,1)
mov    -0x7c(%rsp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rsp)
cmp    %ecx,%eax
je     1692 <kernel_CPU_06_mod_3DRhoCte+0x592>
cmpl   $0x0,-0x70(%rsp)
jle    12c0 <kernel_CPU_06_mod_3DRhoCte+0x1c0>
xor    %eax,%eax
mov    -0x7c(%rsp),%ebp
mov    %ebp,-0x78(%rsp)
jmp    1309 <kernel_CPU_06_mod_3DRhoCte+0x209>
cs nopw 0x0(%rax,%rax,1)
mov    -0x68(%rsp),%eax
add    $0x1,%eax
mov    -0x78(%rsp),%ebp
add    -0x6c(%rsp),%ebp
mov    %ebp,-0x78(%rsp)
cmp    -0x70(%rsp),%eax
je     12c0 <kernel_CPU_06_mod_3DRhoCte+0x1c0>
mov    %eax,-0x68(%rsp)
cmpl   $0x0,-0x74(%rsp)
jle    12f0 <kernel_CPU_06_mod_3DRhoCte+0x1f0>
mov    -0x58(%rsp),%r9
mov    -0x78(%rsp),%eax
mov    %eax,%r15d
jmp    1339 <kernel_CPU_06_mod_3DRhoCte+0x239>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    %ecx,%r15d
add    $0xffffffffffffffff,%r9
je     12f0 <kernel_CPU_06_mod_3DRhoCte+0x1f0>
lea    (%rbx,%r15,1),%eax
movslq %eax,%r12
movss  (%rdx,%r12,4),%xmm7
ucomiss %xmm0,%xmm7
jbe    1330 <kernel_CPU_06_mod_3DRhoCte+0x230>
movss  (%rdi,%r12,4),%xmm15
movaps %xmm15,%xmm3
addss  %xmm15,%xmm3
subss  (%rsi,%r12,4),%xmm3
movaps %xmm7,%xmm5
mulss  -0x64(%rsp),%xmm5
mulss  %xmm7,%xmm5
movss  0x18(%r13),%xmm9
lea    (%r14,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm2
lea    0x6(%rbx,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm2
movss  0x14(%r13),%xmm10
lea    (%r11,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm4
lea    0x5(%rbx,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm4
mulss  %xmm10,%xmm4
mulss  %xmm9,%xmm2
addss  %xmm4,%xmm2
movss  0x10(%r13),%xmm11
lea    (%r8,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm4
lea    0x4(%rbx,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm4
mulss  %xmm11,%xmm4
addss  %xmm2,%xmm4
movss  0xc(%r13),%xmm12
lea    (%r10,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm2
lea    0x3(%rbx,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm2
mulss  %xmm12,%xmm2
addss  %xmm4,%xmm2
movss  0x4(%r13),%xmm13
movss  0x8(%r13),%xmm14
mov    -0x48(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm4
lea    0x2(%rbx,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm4
mulss  %xmm14,%xmm4
addss  %xmm2,%xmm4
mov    -0x40(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm6
lea    0x1(%rbx,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm6
mulss  %xmm13,%xmm6
addss  %xmm4,%xmm6
mulss  0x0(%r13),%xmm15
addss  %xmm15,%xmm6
mulss  %xmm5,%xmm6
addss  %xmm3,%xmm6
movaps %xmm7,%xmm5
mulss  -0x60(%rsp),%xmm5
mulss  %xmm7,%xmm5
mov    0x70(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm2
mov    0x68(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm2
mov    0x60(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm3
mov    0x58(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm3
mulss  %xmm10,%xmm3
mulss  %xmm9,%xmm2
mov    0x50(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm4
mov    0x48(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm4
addss  %xmm3,%xmm2
mulss  %xmm11,%xmm4
addss  %xmm2,%xmm4
mov    0x40(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm3
mov    0x38(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm3
mulss  %xmm12,%xmm3
addss  %xmm4,%xmm3
mov    0x30(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm2
mov    0x28(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm2
mulss  %xmm14,%xmm2
addss  %xmm3,%xmm2
mov    0x20(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm3
mov    0x78(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm3
mulss  %xmm13,%xmm3
mov    0x18(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm4
mov    0x10(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm4
addss  %xmm2,%xmm3
mov    0x8(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm8
mov    (%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm8
addss  %xmm15,%xmm3
mov    -0x8(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm2
mov    -0x10(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm2
mulss  %xmm5,%xmm3
mov    -0x18(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm5
mov    -0x20(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm5
mov    -0x28(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm1
mov    -0x30(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm1
mov    -0x38(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm0
mov    -0x50(%rsp),%rax
lea    (%rax,%r15,1),%eax
cltq
addss  (%rdi,%rax,4),%xmm0
addss  %xmm6,%xmm3
movaps %xmm7,%xmm6
mulss  -0x5c(%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    1330 <kernel_CPU_06_mod_3DRhoCte+0x230>
add    $0x80,%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* %23
br label %31
31: ; preds = %510, %11
%32 = load i32, i32* %23
%33 = load i32, i32* %15
%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* %24
br label %43
43: ; preds = %502, %42
%44 = load i32, i32* %24
%45 = load i32, i32* %16
%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* %24
%504 = add nsw i32 %503, 1
store i32 %504, i32* %24
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* %23
%512 = add nsw i32 %511, 1
store i32 %512, i32* %23
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 %3, 0
br i1 %13, label %14, label %267
14: ; preds = %11
%15 = sub i32 %6, %5
%16 = icmp sgt i32 %15, 0
%17 = icmp sgt i32 %4, 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 nuw 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 %4 to i64
br label %35
35: ; preds = %14, %264
%36 = phi i32 [ 0, %14 ], [ %265, %264 ]
br i1 %16, label %37, label %264
37: ; preds = %35, %261
%38 = phi i32 [ %262, %261 ], [ 0, %35 ]
br i1 %17, label %39, label %261
39: ; preds = %37
%40 = add nsw i32 %38, %5
%41 = mul nsw i32 %40, %12
%42 = add i32 %41, %36
br label %43
43: ; preds = %39, %258
%44 = phi i64 [ 0, %39 ], [ %259, %258 ]
%45 = trunc i64 %44 to i32
%46 = mul i32 %45, %3
%47 = add i32 %42, %46
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %2, i64 %48
%50 = load float, float* %49
%51 = fcmp ogt float %50, 0.000000e+00
br i1 %51, label %52, label %258
52: ; preds = %43
%53 = getelementptr inbounds float, float* %0, i64 %48
%54 = load float, float* %53
%55 = getelementptr inbounds float, float* %1, i64 %48
%56 = load float, float* %55
%57 = fneg float %56
%58 = tail call float @llvm.fmuladd.f32(float %54, float 2.000000e+00, float %57)
%59 = fmul float %50, %7
%60 = fmul float %50, %59
%61 = load float, float* %18
%62 = add nsw i32 %47, -6
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds float, float* %0, i64 %63
%65 = load float, float* %64
%66 = add nsw i32 %47, 6
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds float, float* %0, i64 %67
%69 = load float, float* %68
%70 = fadd float %65, %69
%71 = load float, float* %19
%72 = add nsw i32 %47, -5
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float, float* %0, i64 %73
%75 = load float, float* %74
%76 = add nsw i32 %47, 5
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float, float* %0, i64 %77
%79 = load float, float* %78
%80 = fadd float %75, %79
%81 = fmul float %71, %80
%82 = tail call float @llvm.fmuladd.f32(float %61, float %70, float %81)
%83 = load float, float* %20
%84 = add nsw i32 %47, -4
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds float, float* %0, i64 %85
%87 = load float, float* %86
%88 = add nsw i32 %47, 4
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds float, float* %0, i64 %89
%91 = load float, float* %90
%92 = fadd float %87, %91
%93 = tail call float @llvm.fmuladd.f32(float %83, float %92, float %82)
%94 = load float, float* %21
%95 = add nsw i32 %47, -3
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds float, float* %0, i64 %96
%98 = load float, float* %97
%99 = add nsw i32 %47, 3
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %0, i64 %100
%102 = load float, float* %101
%103 = fadd float %98, %102
%104 = tail call float @llvm.fmuladd.f32(float %94, float %103, float %93)
%105 = load float, float* %22
%106 = add nsw i32 %47, -2
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds float, float* %0, i64 %107
%109 = load float, float* %108
%110 = add nsw i32 %47, 2
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds float, float* %0, i64 %111
%113 = load float, float* %112
%114 = fadd float %109, %113
%115 = tail call float @llvm.fmuladd.f32(float %105, float %114, float %104)
%116 = load float, float* %23
%117 = add nsw i32 %47, -1
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds float, float* %0, i64 %118
%120 = load float, float* %119
%121 = add nsw i32 %47, 1
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds float, float* %0, i64 %122
%124 = load float, float* %123
%125 = fadd float %120, %124
%126 = tail call float @llvm.fmuladd.f32(float %116, float %125, float %115)
%127 = load float, float* %10
%128 = tail call float @llvm.fmuladd.f32(float %127, float %54, float %126)
%129 = tail call float @llvm.fmuladd.f32(float %60, float %128, float %58)
%130 = fmul float %50, %8
%131 = fmul float %50, %130
%132 = sub nsw i32 %47, %24
%133 = sext i32 %132 to i64
%134 = getelementptr inbounds float, float* %0, i64 %133
%135 = load float, float* %134
%136 = add nsw i32 %47, %24
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds float, float* %0, i64 %137
%139 = load float, float* %138
%140 = fadd float %135, %139
%141 = sub nsw i32 %47, %25
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds float, float* %0, i64 %142
%144 = load float, float* %143
%145 = add nsw i32 %47, %25
%146 = sext i32 %145 to i64
%147 = getelementptr inbounds float, float* %0, i64 %146
%148 = load float, float* %147
%149 = fadd float %144, %148
%150 = fmul float %71, %149
%151 = tail call float @llvm.fmuladd.f32(float %61, float %140, float %150)
%152 = sub nsw i32 %47, %26
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds float, float* %0, i64 %153
%155 = load float, float* %154
%156 = add nsw i32 %47, %26
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds float, float* %0, i64 %157
%159 = load float, float* %158
%160 = fadd float %155, %159
%161 = tail call float @llvm.fmuladd.f32(float %83, float %160, float %151)
%162 = sub nsw i32 %47, %27
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds float, float* %0, i64 %163
%165 = load float, float* %164
%166 = add nsw i32 %47, %27
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds float, float* %0, i64 %167
%169 = load float, float* %168
%170 = fadd float %165, %169
%171 = tail call float @llvm.fmuladd.f32(float %94, float %170, float %161)
%172 = sub nsw i32 %47, %28
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds float, float* %0, i64 %173
%175 = load float, float* %174
%176 = add nsw i32 %47, %28
%177 = sext i32 %176 to i64
%178 = getelementptr inbounds float, float* %0, i64 %177
%179 = load float, float* %178
%180 = fadd float %175, %179
%181 = tail call float @llvm.fmuladd.f32(float %105, float %180, float %171)
%182 = sub nsw i32 %47, %3
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds float, float* %0, i64 %183
%185 = load float, float* %184
%186 = add nsw i32 %47, %3
%187 = sext i32 %186 to i64
%188 = getelementptr inbounds float, float* %0, i64 %187
%189 = load float, float* %188
%190 = fadd float %185, %189
%191 = tail call float @llvm.fmuladd.f32(float %116, float %190, float %181)
%192 = tail call float @llvm.fmuladd.f32(float %127, float %54, float %191)
%193 = tail call float @llvm.fmuladd.f32(float %131, float %192, float %129)
%194 = fmul float %50, %9
%195 = fmul float %50, %194
%196 = add nsw i32 %47, %29
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds float, float* %0, i64 %197
%199 = load float, float* %198
%200 = sub nsw i32 %47, %29
%201 = sext i32 %200 to i64
%202 = getelementptr inbounds float, float* %0, i64 %201
%203 = load float, float* %202
%204 = fadd float %199, %203
%205 = add nsw i32 %47, %30
%206 = sext i32 %205 to i64
%207 = getelementptr inbounds float, float* %0, i64 %206
%208 = load float, float* %207
%209 = sub nsw i32 %47, %30
%210 = sext i32 %209 to i64
%211 = getelementptr inbounds float, float* %0, i64 %210
%212 = load float, float* %211
%213 = fadd float %208, %212
%214 = fmul float %71, %213
%215 = tail call float @llvm.fmuladd.f32(float %61, float %204, float %214)
%216 = add nsw i32 %47, %31
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds float, float* %0, i64 %217
%219 = load float, float* %218
%220 = sub nsw i32 %47, %31
%221 = sext i32 %220 to i64
%222 = getelementptr inbounds float, float* %0, i64 %221
%223 = load float, float* %222
%224 = fadd float %219, %223
%225 = tail call float @llvm.fmuladd.f32(float %83, float %224, float %215)
%226 = add nsw i32 %47, %32
%227 = sext i32 %226 to i64
%228 = getelementptr inbounds float, float* %0, i64 %227
%229 = load float, float* %228
%230 = sub nsw i32 %47, %32
%231 = sext i32 %230 to i64
%232 = getelementptr inbounds float, float* %0, i64 %231
%233 = load float, float* %232
%234 = fadd float %229, %233
%235 = tail call float @llvm.fmuladd.f32(float %94, float %234, float %225)
%236 = add nsw i32 %47, %33
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds float, float* %0, i64 %237
%239 = load float, float* %238
%240 = sub nsw i32 %47, %33
%241 = sext i32 %240 to i64
%242 = getelementptr inbounds float, float* %0, i64 %241
%243 = load float, float* %242
%244 = fadd float %239, %243
%245 = tail call float @llvm.fmuladd.f32(float %105, float %244, float %235)
%246 = add nsw i32 %47, %12
%247 = sext i32 %246 to i64
%248 = getelementptr inbounds float, float* %0, i64 %247
%249 = load float, float* %248
%250 = sub nsw i32 %47, %12
%251 = sext i32 %250 to i64
%252 = getelementptr inbounds float, float* %0, i64 %251
%253 = load float, float* %252
%254 = fadd float %249, %253
%255 = tail call float @llvm.fmuladd.f32(float %116, float %254, float %245)
%256 = tail call float @llvm.fmuladd.f32(float %127, float %54, float %255)
%257 = tail call float @llvm.fmuladd.f32(float %195, float %256, float %193)
store float %257, float* %55
br label %258
258: ; preds = %43, %52
%259 = add nuw nsw i64 %44, 1
%260 = icmp eq i64 %259, %34
br i1 %260, label %261, label %43
261: ; preds = %258, %37
%262 = add nuw nsw i32 %38, 1
%263 = icmp eq i32 %262, %15
br i1 %263, label %264, label %37
264: ; preds = %261, %35
%265 = add nuw nsw i32 %36, 1
%266 = icmp eq i32 %265, %3
br i1 %266, label %267, label %35
267: ; preds = %264, %11
ret void
}