user: 请从给定的汇编指令推导出C源代码:

```asm
<kernel_CPU_06_mod_3DRhoCte>:
  endbr64
  mov    %r8d,%r11d
  push   %r15
  mov    %ecx,%r15d
  push   %r14
  imul   %ecx,%r11d
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x38(%rsp),%ecx
  mov    0x40(%rsp),%r12
  sub    %r9d,%ecx
  mov    %rsi,-0x28(%rsp)
  mov    %rdx,-0x20(%rsp)
  mov    %r8d,-0x38(%rsp)
  mov    %ecx,-0x8(%rsp)
  test   %ecx,%ecx
  jle    L38
  mov    %r9d,%eax
  mov    %r15d,%edx
  lea    (%r15,%r15,1),%ebx
  mov    %rdi,%r14
  imul   %r11d,%eax
  lea    (%r11,%r11,1),%ebp
  mov    %ebx,-0x18(%rsp)
  mov    %r11d,%r10d
  lea    0x0(,%r15,4),%ecx
  lea    0x0(%rbp,%r11,1),%ebx
  movaps %xmm0,%xmm5
  movl   $0x0,-0xc(%rsp)
  sub    %ecx,%edx
  lea    (%rcx,%r15,1),%edi
  lea    (%rbx,%rbx,1),%esi
  mov    %ecx,-0x4(%rsp)
  mov    %eax,-0x34(%rsp)
  neg    %eax
  lea    0x0(,%r11,4),%r13d
  neg    %ecx
  mov    %edx,-0x2c(%rsp)
  add    %edx,%edx
  pxor   %xmm4,%xmm4
  mov    %eax,-0x30(%rsp)
  mov    %r15d,%eax
  sub    %edx,%eax
  mov    %esi,-0x10(%rsp)
  lea    0x0(%r13,%r11,1),%esi
  mov    %r13d,%r11d
  mov    %eax,-0x6c(%rsp)
  lea    (%r15,%rdi,1),%eax
  mov    %ebx,%r13d
  mov    %edx,-0x14(%rsp)
  mov    %edi,-0x48(%rsp)
  mov    %esi,-0x4c(%rsp)
  mov    %ecx,-0x44(%rsp)
  mov    %eax,-0x40(%rsp)
  xchg   %ax,%ax
L62:
  mov    -0x38(%rsp),%eax
  test   %eax,%eax
  jle    L39
  mov    %r15d,%eax
  sub    -0x4(%rsp),%eax
  movl   $0x0,-0x3c(%rsp)
  mov    %eax,-0x54(%rsp)
  mov    -0x2c(%rsp),%eax
  mov    -0x34(%rsp),%ebx
  add    %r15d,%eax
  mov    %eax,-0x50(%rsp)
  xchg   %ax,%ax
L4:
  test   %r15d,%r15d
  jle    L36
  mov    -0x18(%rsp),%esi
  mov    %r15d,%edi
  sub    -0x2c(%rsp),%edi
  mov    %ebx,-0x70(%rsp)
  mov    %edi,-0x68(%rsp)
  mov    %r15d,%edi
  mov    -0x30(%rsp),%eax
  sub    %esi,%edi
  mov    %edi,-0x64(%rsp)
  lea    (%r15,%rsi,1),%edi
  lea    (%rbx,%rax,1),%ecx
  movslq -0x34(%rsp),%rax
  mov    %edi,-0x60(%rsp)
  mov    -0x10(%rsp),%edi
  movslq %ecx,%rcx
  add    %rax,%rcx
  mov    -0x20(%rsp),%rax
  lea    (%r15,%rdi,1),%esi
  shl    $0x2,%rcx
  mov    %esi,-0x5c(%rsp)
  mov    %r15d,%esi
  lea    (%rax,%rcx,1),%r9
  mov    -0x14(%rsp),%eax
  sub    %edi,%esi
  lea    (%r14,%rcx,1),%rdx
  add    -0x28(%rsp),%rcx
  mov    %esi,-0x58(%rsp)
  lea    (%rax,%rbx,1),%r8d
  mov    %ebx,%eax
  sub    %r15d,%eax
  nopw   0x0(%rax,%rax,1)
L16:
  movss  (%r9),%xmm6
  comiss %xmm4,%xmm6
  jbe    L3
  mov    -0x44(%rsp),%ebx
  movss  0x18(%r12),%xmm13
  movss  0x14(%r12),%xmm12
  movss  0x10(%r12),%xmm11
  lea    (%rbx,%rax,1),%edi
  mov    -0x40(%rsp),%ebx
  movss  0xc(%r12),%xmm10
  movslq %edi,%rdi
  movss  (%rdx),%xmm14
  movss  0x8(%r12),%xmm9
  lea    (%rbx,%rax,1),%esi
  movss  (%r14,%rdi,4),%xmm0
  mov    -0x54(%rsp),%ebx
  movslq %esi,%rsi
  movss  0x4(%r12),%xmm8
  movss  (%r12),%xmm7
  addss  (%r14,%rsi,4),%xmm0
  mov    -0x6c(%rsp),%esi
  mulss  %xmm14,%xmm7
  addss  %xmm14,%xmm14
  lea    (%rax,%rsi,1),%edi
  movslq %r8d,%rsi
  mulss  %xmm12,%xmm0
  movslq %edi,%rdi
  movss  (%r14,%rdi,4),%xmm3
  addss  (%r14,%rsi,4),%xmm3
  lea    (%rbx,%rax,1),%edi
  mov    -0x48(%rsp),%ebx
  movslq %edi,%rdi
  mulss  %xmm13,%xmm3
  lea    (%rbx,%rax,1),%esi
  mov    -0x50(%rsp),%ebx
  movslq %esi,%rsi
  addss  %xmm3,%xmm0
  movss  (%r14,%rdi,4),%xmm3
  addss  (%r14,%rsi,4),%xmm3
  lea    (%rbx,%rax,1),%edi
  mov    -0x68(%rsp),%ebx
  movslq %edi,%rdi
  mulss  %xmm11,%xmm3
  lea    (%rbx,%rax,1),%esi
  mov    -0x64(%rsp),%ebx
  movslq %esi,%rsi
  addss  %xmm3,%xmm0
  movss  (%r14,%rdi,4),%xmm3
  addss  (%r14,%rsi,4),%xmm3
  lea    (%rbx,%rax,1),%edi
  mov    -0x60(%rsp),%ebx
  movslq %edi,%rdi
  mulss  %xmm10,%xmm3
  lea    (%rbx,%rax,1),%esi
  movslq %esi,%rsi
  addss  %xmm3,%xmm0
  movss  (%r14,%rdi,4),%xmm3
  addss  (%r14,%rsi,4),%xmm3
  mov    -0x6c(%rsp),%esi
  mulss  %xmm9,%xmm3
  lea    (%r8,%rsi,1),%edi
  movslq %eax,%rsi
  movslq %edi,%rdi
  addss  %xmm3,%xmm0
  movss  (%r14,%rdi,4),%xmm3
  addss  (%r14,%rsi,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  (%rcx),%xmm14
  mov    -0x5c(%rsp),%ebx
  mulss  %xmm13,%xmm3
  mulss  %xmm12,%xmm15
  lea    (%rbx,%rax,1),%edi
  mov    -0x58(%rsp),%ebx
  movslq %edi,%rdi
  lea    (%rbx,%rax,1),%esi
  mov    %r15d,%ebx
  movslq %esi,%rsi
  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  (%r14,%rdi,4),%xmm3
  addss  (%r14,%rsi,4),%xmm3
  mov    -0x4c(%rsp),%esi
  lea    (%r15,%rsi,1),%edi
  mulss  %xmm13,%xmm3
  sub    %esi,%ebx
  mov    %ebx,%esi
  add    %eax,%edi
  add    %eax,%esi
  movslq %edi,%rdi
  movslq %esi,%rsi
  movss  (%r14,%rdi,4),%xmm13
  lea    (%r15,%r11,1),%edi
  addss  (%r14,%rsi,4),%xmm13
  mov    %r15d,%esi
  add    %eax,%edi
  sub    %r11d,%esi
  movslq %edi,%rdi
  add    %eax,%esi
  mulss  %xmm13,%xmm12
  movslq %esi,%rsi
  addss  %xmm12,%xmm3
  movss  (%r14,%rdi,4),%xmm12
  lea    (%r15,%r13,1),%edi
  addss  (%r14,%rsi,4),%xmm12
  mov    %r15d,%esi
  add    %eax,%edi
  sub    %r13d,%esi
  movslq %edi,%rdi
  add    %eax,%esi
  mulss  %xmm12,%xmm11
  movslq %esi,%rsi
  addss  %xmm11,%xmm3
  movss  (%r14,%rdi,4),%xmm11
  lea    (%r15,%rbp,1),%edi
  addss  (%r14,%rsi,4),%xmm11
  mov    %r15d,%esi
  add    %eax,%edi
  sub    %ebp,%esi
  movslq %edi,%rdi
  add    %eax,%esi
  mulss  %xmm11,%xmm10
  movslq %esi,%rsi
  addss  %xmm10,%xmm3
  movss  (%r14,%rdi,4),%xmm10
  lea    (%r15,%r10,1),%edi
  addss  (%r14,%rsi,4),%xmm10
  mov    %r15d,%esi
  add    %eax,%edi
  sub    %r10d,%esi
  movslq %edi,%rdi
  add    %eax,%esi
  mulss  %xmm10,%xmm9
  movslq %esi,%rsi
  addss  %xmm9,%xmm3
  movss  (%r14,%rdi,4),%xmm9
  addss  (%r14,%rsi,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,(%rcx)
L3:
  add    $0x4,%r9
  add    $0x4,%rdx
  add    $0x4,%rcx
  add    $0x1,%r8d
  add    $0x1,%eax
  cmp    -0x70(%rsp),%eax
  jne    L16
  mov    -0x70(%rsp),%ebx
L36:
  addl   $0x1,-0x3c(%rsp)
  add    %r15d,%ebx
  mov    -0x3c(%rsp),%eax
  cmp    %eax,-0x38(%rsp)
  jne    L4
L39:
  addl   $0x1,-0xc(%rsp)
  mov    -0xc(%rsp),%eax
  add    %r10d,-0x34(%rsp)
  sub    %r10d,-0x30(%rsp)
  cmp    %eax,-0x8(%rsp)
  jne    L62
L38:
  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]);
        }
      }
    }
}
```