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