user: Please show me the C implementation that matches the assembly code listed below:

```asm
<kernel_CPU_06_mod_3DRhoCte>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %r8d,%r12d
  push   %rbp
  imul   %ecx,%r12d
  push   %rbx
  mov    0x40(%rsp),%rbx
  mov    %rdx,-0x20(%rsp)
  mov    0x38(%rsp),%edx
  mov    %rsi,-0x28(%rsp)
  sub    %r9d,%edx
  mov    %r8d,-0x18(%rsp)
  mov    %edx,-0x4(%rsp)
  test   %edx,%edx
  jle    L29
  mov    %ecx,%r14d
  lea    0x0(,%rcx,4),%edx
  mov    %ecx,%esi
  mov    %rdi,%r13
  lea    (%r12,%r12,1),%ebp
  sub    %edx,%esi
  lea    (%r14,%r14,1),%r11d
  add    %r14d,%edx
  lea    0x0(%rbp,%r12,1),%r15d
  mov    %r11d,-0x38(%rsp)
  lea    (%rsi,%rsi,1),%ecx
  lea    (%r15,%r15,1),%r11d
  lea    0x0(,%r12,4),%r10d
  mov    %edx,-0x54(%rsp)
  neg    %edx
  mov    %r11d,-0x10(%rsp)
  lea    (%r10,%r12,1),%r11d
  mov    %esi,-0x6c(%rsp)
  mov    %ecx,-0x14(%rsp)
  test   %r8d,%r8d
  jle    L29
  test   %r14d,%r14d
  jle    L29
  mov    %r9d,%eax
  movaps %xmm0,%xmm5
  pxor   %xmm4,%xmm4
  mov    %r15d,%edi
  imul   %r12d,%eax
  movl   $0x0,-0xc(%rsp)
  mov    %r12d,%r15d
  mov    %edi,%r12d
  mov    %eax,-0x8(%rsp)
  movslq %r14d,%rax
  shl    $0x2,%rax
  mov    %rax,-0x30(%rsp)
  mov    %r14d,%eax
  sub    %ecx,%eax
  mov    %eax,-0x70(%rsp)
  lea    (%rdx,%r14,1),%eax
  mov    %eax,-0x50(%rsp)
  mov    %r14d,%eax
  sub    %edx,%eax
  mov    %eax,-0x4c(%rsp)
  lea    (%rsi,%r14,1),%eax
  mov    %eax,-0x48(%rsp)
  mov    %r11d,%eax
  mov    %r10d,%r11d
  mov    %r14d,%r10d
  mov    %ebp,%r14d
  mov    %eax,%ebp
  nopl   0x0(%rax)
L37:
  movslq -0x8(%rsp),%rax
  movl   $0x0,-0x44(%rsp)
  mov    %rax,%rdi
  shl    $0x2,%rax
  mov    %edi,-0x74(%rsp)
  mov    %rax,-0x40(%rsp)
  mov    %ebp,%eax
  mov    %r12d,%ebp
  mov    %eax,%r12d
  nopl   (%rax)
L1:
  mov    -0x38(%rsp),%ecx
  mov    %r10d,%esi
  sub    -0x6c(%rsp),%esi
  mov    -0x40(%rsp),%rax
  mov    -0x20(%rsp),%r9
  mov    %esi,-0x68(%rsp)
  mov    %r10d,%esi
  sub    %ecx,%esi
  mov    -0x28(%rsp),%r8
  mov    -0x14(%rsp),%edi
  mov    %esi,-0x64(%rsp)
  lea    (%r10,%rcx,1),%esi
  add    %rax,%r9
  lea    0x0(%r13,%rax,1),%rdx
  mov    %esi,-0x60(%rsp)
  mov    -0x10(%rsp),%esi
  add    %rax,%r8
  mov    -0x74(%rsp),%eax
  lea    (%r10,%rsi,1),%ecx
  mov    %ecx,-0x5c(%rsp)
  mov    %r10d,%ecx
  add    %eax,%edi
  sub    %r10d,%eax
  sub    %esi,%ecx
  mov    %eax,-0x34(%rsp)
  mov    %ecx,-0x58(%rsp)
  nopl   0x0(%rax,%rax,1)
L36:
  movss  (%r9),%xmm6
  comiss %xmm4,%xmm6
  jbe    L0
  mov    -0x50(%rsp),%esi
  mov    -0x4c(%rsp),%ecx
  movss  0x18(%rbx),%xmm13
  movss  0x14(%rbx),%xmm12
  add    %eax,%esi
  add    %eax,%ecx
  movss  0x10(%rbx),%xmm11
  movss  0xc(%rbx),%xmm10
  movslq %esi,%rsi
  movslq %ecx,%rcx
  movss  0x8(%rbx),%xmm9
  movss  (%rbx),%xmm7
  movss  0x0(%r13,%rsi,4),%xmm0
  mov    -0x70(%rsp),%esi
  addss  0x0(%r13,%rcx,4),%xmm0
  movslq %edi,%rcx
  movss  0x4(%rbx),%xmm8
  add    %eax,%esi
  movss  (%rdx),%xmm14
  movslq %esi,%rsi
  movss  0x0(%r13,%rsi,4),%xmm3
  addss  0x0(%r13,%rcx,4),%xmm3
  mulss  %xmm12,%xmm0
  mov    -0x6c(%rsp),%ecx
  mulss  %xmm14,%xmm7
  addss  %xmm14,%xmm14
  mulss  %xmm13,%xmm3
  lea    (%rcx,%rax,1),%esi
  mov    -0x54(%rsp),%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  movslq %ecx,%rcx
  addss  %xmm3,%xmm0
  movss  0x0(%r13,%rsi,4),%xmm3
  addss  0x0(%r13,%rcx,4),%xmm3
  mov    -0x48(%rsp),%ecx
  mulss  %xmm11,%xmm3
  lea    (%rcx,%rax,1),%esi
  mov    -0x68(%rsp),%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  movslq %ecx,%rcx
  addss  %xmm3,%xmm0
  movss  0x0(%r13,%rsi,4),%xmm3
  addss  0x0(%r13,%rcx,4),%xmm3
  mov    -0x64(%rsp),%ecx
  mulss  %xmm10,%xmm3
  lea    (%rcx,%rax,1),%esi
  mov    -0x60(%rsp),%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  movslq %ecx,%rcx
  addss  %xmm3,%xmm0
  movss  0x0(%r13,%rsi,4),%xmm3
  mov    -0x70(%rsp),%esi
  addss  0x0(%r13,%rcx,4),%xmm3
  movslq %eax,%rcx
  add    %edi,%esi
  movslq %esi,%rsi
  mulss  %xmm9,%xmm3
  addss  %xmm3,%xmm0
  movss  0x0(%r13,%rsi,4),%xmm3
  addss  0x0(%r13,%rcx,4),%xmm3
  mulss  %xmm8,%xmm3
  addss  %xmm3,%xmm0
  movaps %xmm6,%xmm3
  mulss  %xmm1,%xmm3
  addss  %xmm7,%xmm0
  mulss  %xmm6,%xmm3
  mulss  %xmm3,%xmm0
  movss  -0x18(%rdx),%xmm3
  addss  0x18(%rdx),%xmm3
  movss  -0x14(%rdx),%xmm15
  addss  0x14(%rdx),%xmm15
  subss  (%r8),%xmm14
  mov    -0x5c(%rsp),%esi
  mulss  %xmm13,%xmm3
  mov    -0x58(%rsp),%ecx
  mulss  %xmm12,%xmm15
  add    %eax,%esi
  add    %eax,%ecx
  movslq %esi,%rsi
  movslq %ecx,%rcx
  addss  %xmm15,%xmm3
  movss  -0x10(%rdx),%xmm15
  addss  0x10(%rdx),%xmm15
  mulss  %xmm11,%xmm15
  addss  %xmm15,%xmm3
  movss  -0xc(%rdx),%xmm15
  addss  0xc(%rdx),%xmm15
  mulss  %xmm10,%xmm15
  addss  %xmm15,%xmm3
  movss  -0x8(%rdx),%xmm15
  addss  0x8(%rdx),%xmm15
  mulss  %xmm9,%xmm15
  addss  %xmm15,%xmm3
  movss  -0x4(%rdx),%xmm15
  addss  0x4(%rdx),%xmm15
  mulss  %xmm8,%xmm15
  addss  %xmm15,%xmm3
  movaps %xmm6,%xmm15
  mulss  %xmm5,%xmm15
  addss  %xmm7,%xmm3
  mulss  %xmm6,%xmm15
  mulss  %xmm15,%xmm3
  addss  %xmm14,%xmm3
  addss  %xmm3,%xmm0
  movss  0x0(%r13,%rsi,4),%xmm3
  lea    (%r10,%r12,1),%esi
  addss  0x0(%r13,%rcx,4),%xmm3
  mov    %r10d,%ecx
  add    %eax,%esi
  sub    %r12d,%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  mulss  %xmm13,%xmm3
  movslq %ecx,%rcx
  movss  0x0(%r13,%rsi,4),%xmm13
  lea    (%r10,%r11,1),%esi
  addss  0x0(%r13,%rcx,4),%xmm13
  mov    %r10d,%ecx
  add    %eax,%esi
  sub    %r11d,%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  mulss  %xmm13,%xmm12
  movslq %ecx,%rcx
  addss  %xmm12,%xmm3
  movss  0x0(%r13,%rsi,4),%xmm12
  lea    (%r10,%rbp,1),%esi
  addss  0x0(%r13,%rcx,4),%xmm12
  mov    %r10d,%ecx
  add    %eax,%esi
  sub    %ebp,%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  mulss  %xmm12,%xmm11
  movslq %ecx,%rcx
  addss  %xmm11,%xmm3
  movss  0x0(%r13,%rsi,4),%xmm11
  lea    (%r10,%r14,1),%esi
  addss  0x0(%r13,%rcx,4),%xmm11
  mov    %r10d,%ecx
  add    %eax,%esi
  sub    %r14d,%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  mulss  %xmm11,%xmm10
  movslq %ecx,%rcx
  addss  %xmm10,%xmm3
  movss  0x0(%r13,%rsi,4),%xmm10
  lea    (%r10,%r15,1),%esi
  addss  0x0(%r13,%rcx,4),%xmm10
  mov    %r10d,%ecx
  add    %eax,%esi
  sub    %r15d,%ecx
  movslq %esi,%rsi
  add    %eax,%ecx
  mulss  %xmm10,%xmm9
  movslq %ecx,%rcx
  addss  %xmm9,%xmm3
  movss  0x0(%r13,%rsi,4),%xmm9
  addss  0x0(%r13,%rcx,4),%xmm9
  mulss  %xmm9,%xmm8
  addss  %xmm8,%xmm3
  addss  %xmm7,%xmm3
  movaps %xmm6,%xmm7
  mulss  %xmm2,%xmm7
  mulss  %xmm7,%xmm6
  mulss  %xmm6,%xmm3
  addss  %xmm3,%xmm0
  movss  %xmm0,(%r8)
L0:
  add    $0x4,%r9
  add    $0x4,%rdx
  add    $0x4,%r8
  add    $0x1,%edi
  add    $0x1,%eax
  cmp    %eax,-0x74(%rsp)
  jne    L36
  mov    -0x34(%rsp),%edi
  add    -0x38(%rsp),%edi
  addl   $0x1,-0x44(%rsp)
  mov    -0x30(%rsp),%rdx
  mov    %edi,-0x74(%rsp)
  mov    -0x44(%rsp),%eax
  add    %rdx,-0x40(%rsp)
  cmp    %eax,-0x18(%rsp)
  jne    L1
  mov    %r12d,%eax
  addl   $0x1,-0xc(%rsp)
  mov    %ebp,%r12d
  add    %r15d,-0x8(%rsp)
  mov    %eax,%ebp
  mov    -0xc(%rsp),%eax
  cmp    %eax,-0x4(%rsp)
  jne    L37
L29:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
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) collapse(2) private(index_X, index_Y, index, k)
  for (k = 0; k < k1 - k0; k++)
    for (index_Y = 0; index_Y < nnoj; index_Y++) {
#pragma GCC ivdep
      for (index_X = 0; index_X < nnoi; index_X++) {
        index = (index_Y * nnoi + index_X) + (k0 + k) * stride;
        if (gVorg[index] > 0.0f) {
          gU1[index] = 2.0f * gU0[index] - gU1[index] + FATMDFX * gVorg[index] * gVorg[index] * (+W[6] * (gU0[index - 6] + gU0[index + 6]) + W[5] * (gU0[index - 5] + gU0[index + 5]) + W[4] * (gU0[index - 4] + gU0[index + 4]) + W[3] * (gU0[index - 3] + gU0[index + 3]) + W[2] * (gU0[index - 2] + gU0[index + 2]) + W[1] * (gU0[index - 1] + gU0[index + 1]) + W[0] * gU0[index]) + FATMDFY * gVorg[index] * gVorg[index] * (+W[6] * (gU0[index - 6 * nnoi] + gU0[index + 6 * nnoi]) + W[5] * (gU0[index - 5 * nnoi] + gU0[index + 5 * nnoi]) + W[4] * (gU0[index - 4 * nnoi] + gU0[index + 4 * nnoi]) + W[3] * (gU0[index - 3 * nnoi] + gU0[index + 3 * nnoi]) + W[2] * (gU0[index - 2 * nnoi] + gU0[index + 2 * nnoi]) + W[1] * (gU0[index - nnoi] + gU0[index + nnoi]) + W[0] * gU0[index]) + FATMDFZ * gVorg[index] * gVorg[index] * (+W[6] * (gU0[index + 6 * stride] + gU0[index - 6 * stride]) + W[5] * (gU0[index + 5 * stride] + gU0[index - 5 * stride]) + W[4] * (gU0[index + 4 * stride] + gU0[index - 4 * stride]) + W[3] * (gU0[index + 3 * stride] + gU0[index - 3 * stride]) + W[2] * (gU0[index + 2 * stride] + gU0[index - 2 * stride]) + W[1] * (gU0[index + stride] + gU0[index - stride]) + W[0] * gU0[index]);
        }
      }
    }
}
```