void kernel_CPU_06_mod_3DRhoCte(float *gU0, float *gU1, float *gVorg, int nnoi,
                                int nnoj, int k0, int k1, float FATMDFX,
                                float FATMDFY, float FATMDFZ, float *W) {
  int index_X, index_Y;
  int stride = nnoi * nnoj;
  int index, k;
#pragma omp parallel for default(shared) private(index_X, index_Y, index, k)
  for (index_Y = 0; index_Y < nnoj; index_Y++)
    for (index_X = 0; index_X < nnoi; index_X++) {
#pragma simd
      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,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    16e6 <kernel_CPU_06_mod_3DRhoCte+0x5e6>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%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    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1155 <kernel_CPU_06_mod_3DRhoCte+0x55>
jmp    16d8 <kernel_CPU_06_mod_3DRhoCte+0x5d8>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1142 <kernel_CPU_06_mod_3DRhoCte+0x42>
pop    %rbp
ret

x86 assembly (clang, -O2):
<kernel_CPU_06_mod_3DRhoCte>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
movss  %xmm2,0x24(%rsp)
mov    %rcx,0x28(%rsp)
test   %r8d,%r8d
jle    16b8 <kernel_CPU_06_mod_3DRhoCte+0x5b8>
mov    0x90(%rsp),%ecx
mov    %r8d,%eax
mov    %r8d,-0x68(%rsp)
mov    0x28(%rsp),%r8
imul   %r8d,%eax
sub    %r9d,%ecx
mov    %rcx,0x30(%rsp)
lea    (%r8,%r8,1),%ecx
lea    (%rcx,%rcx,2),%r14d
lea    (%r8,%r8,2),%r10d
lea    (%rax,%rax,1),%ebx
lea    (%rbx,%rbx,2),%r12d
lea    (%rax,%rax,4),%r13d
mov    %rdi,%r11
lea    (%rax,%rax,2),%edi
movslq %r9d,%r15
movslq %eax,%rbp
mov    %rbp,-0x78(%rsp)
shl    $0x2,%eax
movslq %ebx,%rbx
mov    %rbx,0x18(%rsp)
movslq %edi,%rdi
mov    %rdi,0x10(%rsp)
cltq
mov    %rax,0x8(%rsp)
movslq %r13d,%rax
mov    %rax,(%rsp)
movslq %r12d,%r9
movslq %ecx,%rax
movslq %r10d,%rbx
lea    0x0(,%r8,4),%ecx
movslq %ecx,%r10
lea    (%r8,%r8,4),%ecx
movslq %ecx,%rcx
movslq %r14d,%r14
mov    %r15,%r13
imul   %rbp,%r13
lea    (%rdx,%r13,4),%rdx
lea    (%rsi,%r13,4),%rsi
mov    %r13,%rdi
sub    %r14,%rdi
mov    %rdi,-0x8(%rsp)
add    %r13,%r14
mov    %r14,-0x10(%rsp)
mov    %r13,%rdi
sub    %rcx,%rdi
mov    %rdi,-0x18(%rsp)
add    %r13,%rcx
mov    %rcx,-0x20(%rsp)
mov    %r13,%rdi
sub    %r10,%rdi
mov    %rdi,-0x28(%rsp)
add    %r13,%r10
mov    %r10,-0x30(%rsp)
mov    %r13,%rcx
sub    %rbx,%rcx
mov    %rcx,-0x50(%rsp)
add    %r13,%rbx
mov    %rbx,-0x38(%rsp)
mov    %r13,%r12
sub    %rax,%r12
add    %r13,%rax
mov    %rax,-0x40(%rsp)
movslq %r8d,%rcx
mov    %r13,%rbx
sub    %rcx,%rbx
mov    %ecx,%eax
mov    %rax,0x48(%rsp)
lea    0x0(,%rcx,4),%rax
mov    %rax,0x40(%rsp)
add    %r13,%rcx
mov    %r13,%r14
sub    %r9,%r14
lea    (%r9,%r13,1),%rax
mov    %rax,-0x48(%rsp)
mov    %r13,%r10
mov    (%rsp),%rax
sub    %rax,%r10
add    %r13,%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,-0x80(%rsp)
lea    (%r11,%r13,4),%rax
mov    0x18(%rsp),%rbp
mov    %r13,%rdi
add    %rbp,%rdi
mov    %rdi,-0x70(%rsp)
sub    %rbp,%r13
mov    -0x8(%rsp),%rbp
lea    (%r11,%rbp,4),%rbp
mov    %rbp,0x18(%rsp)
mov    -0x10(%rsp),%rbp
lea    (%r11,%rbp,4),%rbp
mov    %rbp,0x10(%rsp)
mov    -0x18(%rsp),%rbp
lea    (%r11,%rbp,4),%rbp
mov    %rbp,0x8(%rsp)
mov    -0x20(%rsp),%rbp
lea    (%r11,%rbp,4),%rbp
mov    %rbp,(%rsp)
mov    -0x28(%rsp),%rbp
lea    (%r11,%rbp,4),%rbp
mov    %rbp,-0x8(%rsp)
mov    -0x30(%rsp),%rbp
lea    (%r11,%rbp,4),%rbp
mov    %rbp,-0x10(%rsp)
mov    -0x50(%rsp),%rdi
lea    (%r11,%rdi,4),%rbp
mov    %rbp,-0x18(%rsp)
mov    -0x38(%rsp),%rdi
lea    (%r11,%rdi,4),%rbp
mov    %rbp,-0x20(%rsp)
lea    (%r11,%r12,4),%rbp
mov    %rbp,-0x28(%rsp)
mov    -0x40(%rsp),%rdi
lea    (%r11,%rdi,4),%rbp
mov    %rbp,-0x30(%rsp)
lea    (%r11,%rbx,4),%rbx
mov    %rbx,-0x38(%rsp)
lea    (%r11,%rcx,4),%rcx
mov    %rcx,-0x40(%rsp)
mov    -0x48(%rsp),%rcx
lea    (%r11,%rcx,4),%rcx
mov    %rcx,-0x48(%rsp)
lea    (%r11,%r14,4),%rcx
mov    %rcx,-0x50(%rsp)
mov    -0x58(%rsp),%rcx
lea    (%r11,%rcx,4),%rcx
mov    %rcx,-0x58(%rsp)
lea    (%r11,%r10,4),%rcx
mov    %rcx,0x50(%rsp)
mov    -0x60(%rsp),%rcx
lea    (%r11,%rcx,4),%rcx
mov    %rcx,-0x60(%rsp)
lea    (%r11,%r9,4),%rbp
mov    -0x80(%rsp),%rcx
lea    (%r11,%rcx,4),%r10
lea    (%r11,%r8,4),%r8
mov    -0x70(%rsp),%rcx
lea    (%r11,%rcx,4),%rcx
lea    (%r11,%r13,4),%r9
lea    0x1(%r15),%rdi
mov    -0x78(%rsp),%r14
imul   %r14,%rdi
lea    (%r11,%rdi,4),%rdi
add    $0xffffffffffffffff,%r15
imul   %r14,%r15
lea    (%r11,%r15,4),%r11
mov    -0x68(%rsp),%ebx
mov    %rbx,0x38(%rsp)
lea    0x0(,%r14,4),%r14
xorps  %xmm2,%xmm2
mov    0x98(%rsp),%r15
xor    %ebx,%ebx
mov    %rbx,-0x80(%rsp)
xor    %r12d,%r12d
jmp    13f3 <kernel_CPU_06_mod_3DRhoCte+0x2f3>
nopw   0x0(%rax,%rax,1)
mov    -0x70(%rsp),%r12
add    $0x1,%r12
mov    -0x80(%rsp),%rbx
add    0x40(%rsp),%rbx
mov    %rbx,-0x80(%rsp)
cmp    0x38(%rsp),%r12
je     16b8 <kernel_CPU_06_mod_3DRhoCte+0x5b8>
mov    %r12,-0x70(%rsp)
cmpl   $0x0,0x28(%rsp)
jle    13d0 <kernel_CPU_06_mod_3DRhoCte+0x2d0>
mov    -0x80(%rsp),%rbx
mov    %rbx,-0x78(%rsp)
xor    %ebx,%ebx
jmp    1426 <kernel_CPU_06_mod_3DRhoCte+0x326>
nopl   (%rax)
mov    -0x68(%rsp),%rbx
add    $0x1,%rbx
addq   $0x4,-0x78(%rsp)
cmp    0x48(%rsp),%rbx
je     13d0 <kernel_CPU_06_mod_3DRhoCte+0x2d0>
mov    %rbx,-0x68(%rsp)
cmpl   $0x0,0x30(%rsp)
jle    1410 <kernel_CPU_06_mod_3DRhoCte+0x310>
mov    0x30(%rsp),%r12
mov    -0x78(%rsp),%r13
jmp    1449 <kernel_CPU_06_mod_3DRhoCte+0x349>
xchg   %ax,%ax
add    %r14,%r13
add    $0xffffffffffffffff,%r12
je     1410 <kernel_CPU_06_mod_3DRhoCte+0x310>
movss  (%rdx,%r13,1),%xmm3
ucomiss %xmm2,%xmm3
jbe    1440 <kernel_CPU_06_mod_3DRhoCte+0x340>
movss  (%rax,%r13,1),%xmm15
movaps %xmm15,%xmm7
addss  %xmm15,%xmm7
subss  (%rsi,%r13,1),%xmm7
movaps %xmm3,%xmm5
mulss  %xmm0,%xmm5
mulss  %xmm3,%xmm5
movss  0x18(%r15),%xmm9
movss  -0x18(%rax,%r13,1),%xmm6
movss  -0x14(%rax,%r13,1),%xmm4
addss  0x18(%rax,%r13,1),%xmm6
movss  0x14(%r15),%xmm11
addss  0x14(%rax,%r13,1),%xmm4
mulss  %xmm11,%xmm4
mulss  %xmm9,%xmm6
addss  %xmm4,%xmm6
movss  0x10(%r15),%xmm12
movss  -0x10(%rax,%r13,1),%xmm4
addss  0x10(%rax,%r13,1),%xmm4
mulss  %xmm12,%xmm4
addss  %xmm6,%xmm4
movss  0xc(%r15),%xmm13
movss  -0xc(%rax,%r13,1),%xmm6
addss  0xc(%rax,%r13,1),%xmm6
mulss  %xmm13,%xmm6
addss  %xmm4,%xmm6
movss  0x4(%r15),%xmm10
movss  0x8(%r15),%xmm14
movss  -0x8(%rax,%r13,1),%xmm4
addss  0x8(%rax,%r13,1),%xmm4
mulss  %xmm14,%xmm4
addss  %xmm6,%xmm4
movss  -0x4(%rax,%r13,1),%xmm6
addss  0x4(%rax,%r13,1),%xmm6
mulss  %xmm10,%xmm6
addss  %xmm4,%xmm6
mulss  (%r15),%xmm15
addss  %xmm15,%xmm6
mulss  %xmm5,%xmm6
addss  %xmm7,%xmm6
movaps %xmm3,%xmm5
mulss  %xmm1,%xmm5
mov    0x18(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm4
mov    0x10(%rsp),%rbx
addss  (%rbx,%r13,1),%xmm4
mov    0x8(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm7
mov    (%rsp),%rbx
addss  (%rbx,%r13,1),%xmm7
mulss  %xmm11,%xmm7
mulss  %xmm9,%xmm4
addss  %xmm7,%xmm4
mov    -0x8(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm7
mov    -0x10(%rsp),%rbx
addss  (%rbx,%r13,1),%xmm7
mulss  %xmm3,%xmm5
mulss  %xmm12,%xmm7
addss  %xmm4,%xmm7
mov    -0x18(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm4
mov    -0x20(%rsp),%rbx
addss  (%rbx,%r13,1),%xmm4
mulss  %xmm13,%xmm4
addss  %xmm7,%xmm4
mov    -0x28(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm8
mov    -0x30(%rsp),%rbx
addss  (%rbx,%r13,1),%xmm8
mulss  %xmm14,%xmm8
addss  %xmm4,%xmm8
mov    -0x38(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm7
mov    -0x40(%rsp),%rbx
addss  (%rbx,%r13,1),%xmm7
mulss  %xmm10,%xmm7
addss  %xmm8,%xmm7
addss  %xmm15,%xmm7
mulss  %xmm5,%xmm7
addss  %xmm6,%xmm7
movaps %xmm3,%xmm5
mov    -0x48(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm4
mov    -0x50(%rsp),%rbx
addss  (%rbx,%r13,1),%xmm4
mulss  0x24(%rsp),%xmm5
mov    -0x58(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm6
mov    0x50(%rsp),%rbx
addss  (%rbx,%r13,1),%xmm6
mulss  %xmm3,%xmm5
mulss  %xmm11,%xmm6
mulss  %xmm9,%xmm4
addss  %xmm6,%xmm4
mov    -0x60(%rsp),%rbx
movss  (%rbx,%r13,1),%xmm3
addss  0x0(%rbp,%r13,1),%xmm3
mulss  %xmm12,%xmm3
addss  %xmm4,%xmm3
movss  (%r10,%r13,1),%xmm4
addss  (%r8,%r13,1),%xmm4
mulss  %xmm13,%xmm4
movss  (%rcx,%r13,1),%xmm6
addss  (%r9,%r13,1),%xmm6
addss  %xmm3,%xmm4
mulss  %xmm14,%xmm6
addss  %xmm4,%xmm6
movss  (%rdi,%r13,1),%xmm3
addss  (%r11,%r13,1),%xmm3
mulss  %xmm10,%xmm3
addss  %xmm6,%xmm3
addss  %xmm15,%xmm3
mulss  %xmm5,%xmm3
addss  %xmm7,%xmm3
movss  %xmm3,(%rsi,%r13,1)
jmp    1440 <kernel_CPU_06_mod_3DRhoCte+0x340>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @kernel_CPU_06_mod_3DRhoCte(float* %0, float* %1, float* %2, i32 %3, i32 %4, i32 %5, i32 %6, float %7, float %8, float %9, float* %10) {
%12 = alloca float*
%13 = alloca float*
%14 = alloca float*
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca float
%20 = alloca float
%21 = alloca float
%22 = alloca float*
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
store float* %0, float** %12
store float* %1, float** %13
store float* %2, float** %14
store i32 %3, i32* %15
store i32 %4, i32* %16
store i32 %5, i32* %17
store i32 %6, i32* %18
store float %7, float* %19
store float %8, float* %20
store float %9, float* %21
store float* %10, float** %22
%28 = load i32, i32* %15
%29 = load i32, i32* %16
%30 = mul nsw i32 %28, %29
store i32 %30, i32* %25
store i32 0, i32* %24
br label %31
31: ; preds = %510, %11
%32 = load i32, i32* %24
%33 = load i32, i32* %16
%34 = icmp slt i32 %32, %33
br i1 %34, label %35, label %513
35: ; preds = %31
store i32 0, i32* %23
br label %36
36: ; preds = %506, %35
%37 = load i32, i32* %23
%38 = load i32, i32* %15
%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* %23
%508 = add nsw i32 %507, 1
store i32 %508, i32* %23
br label %36
509: ; preds = %36
br label %510
510: ; preds = %509
%511 = load i32, i32* %24
%512 = add nsw i32 %511, 1
store i32 %512, i32* %24
br label %31
513: ; preds = %31
ret void
}

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