user: 请将这段汇编代码片段转换成等效的C代码:

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