user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<VolumeAveragedVOFCoefficientsEvaluate>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x48(%rsp),%r12
  mov    0x50(%rsp),%r13
  mov    %esi,-0x4c(%rsp)
  test   %edi,%edi
  jle    L77
  mov    %esi,%ebp
  movslq %edi,%rbx
  mov    %rdx,%rsi
  mov    %r9,%rcx
  mov    %ebp,%edi
  mov    %rbx,-0x28(%rsp)
  mov    %r8,%rax
  and    $0xfffffffe,%ebp
  mov    %edi,%edx
  mov    %rdi,-0x20(%rsp)
  xor    %r14d,%r14d
  movq   $0x0,-0x60(%rsp)
  shr    %edx
  sub    $0x1,%edx
  add    $0x1,%rdx
  shl    $0x4,%rdx
  mov    %rdx,-0x18(%rsp)
  lea    -0x1(%rdi),%edx
  lea    0x1(%rdx),%rbx
  mov    %rbx,-0x10(%rsp)
  nopl   (%rax)
L29:
  movsd  (%rcx),%xmm0
  mulsd  (%rax),%xmm0
  mov    0x38(%rsp),%rbx
  mov    -0x60(%rsp),%rdi
  mov    -0x4c(%rsp),%edx
  movsd  %xmm0,(%rbx,%rdi,8)
  mov    0x40(%rsp),%rbx
  movsd  (%rax),%xmm1
  movsd  %xmm1,(%rbx,%rdi,8)
  test   %edx,%edx
  jle    L8
  movslq %r14d,%r10
  lea    0x8(%rcx),%r11
  lea    0x0(,%r10,8),%rdx
  mov    %r11,-0x58(%rsp)
  lea    0x10(,%r10,8),%r15
  lea    (%rsi,%rdx,1),%rdi
  lea    -0x8(%r15),%rbx
  mov    %rdi,-0x38(%rsp)
  mov    -0x20(%rsp),%rdi
  lea    0x0(%r13,%rdx,1),%r8
  lea    (%r12,%rdx,1),%r9
  lea    (%rdi,%r10,1),%r11
  shl    $0x3,%r11
  lea    (%r12,%r11,1),%rdi
  add    %r13,%r11
  mov    %r11,-0x40(%rsp)
  lea    (%rsi,%rbx,1),%r11
  cmp    %r11,%r8
  mov    %rdi,-0x30(%rsp)
  lea    0x8(%rax),%rdi
  setne  %r11b
  add    %r12,%rbx
  mov    %rdi,-0x48(%rsp)
  mov    -0x30(%rsp),%rdi
  cmp    %rbx,%r8
  setne  %bl
  and    %ebx,%r11d
  cmpl   $0x1,-0x4c(%rsp)
  setne  %bl
  and    %ebx,%r11d
  lea    (%r12,%r15,1),%rbx
  cmp    %rbx,-0x38(%rsp)
  setae  %bl
  add    %rsi,%r15
  cmp    %r15,%r9
  setae  %r15b
  or     %r15d,%ebx
  and    %r11d,%ebx
  cmp    -0x58(%rsp),%r9
  setae  %r11b
  cmp    %rcx,%rdi
  setbe  %r15b
  or     %r15d,%r11d
  mov    %ebx,%r15d
  and    %r11d,%r15d
  cmp    %rax,%rdi
  mov    -0x48(%rsp),%rdi
  setbe  %bl
  cmp    %rdi,%r9
  setae  %r11b
  or     %ebx,%r11d
  and    %r15d,%r11d
  cmp    -0x58(%rsp),%r8
  setae  %bl
  cmp    -0x40(%rsp),%rcx
  setae  %r15b
  or     %r15d,%ebx
  test   %bl,%r11b
  je     L14
  cmp    %rdi,%r8
  setae  %r11b
  cmp    -0x40(%rsp),%rax
  setae  %bl
  or     %r11b,%bl
  je     L14
  movsd  (%rcx),%xmm2
  mov    -0x38(%rsp),%rdi
  unpcklpd %xmm1,%xmm1
  xor    %edx,%edx
  mov    -0x18(%rsp),%r10
  unpcklpd %xmm2,%xmm2
  nopl   0x0(%rax,%rax,1)
L83:
  movupd (%rdi,%rdx,1),%xmm0
  mulpd  %xmm1,%xmm0
  mulpd  %xmm2,%xmm0
  movups %xmm0,(%r9,%rdx,1)
  movupd (%rdi,%rdx,1),%xmm0
  mulpd  %xmm1,%xmm0
  movups %xmm0,(%r8,%rdx,1)
  add    $0x10,%rdx
  cmp    %r10,%rdx
  jne    L83
  cmp    %ebp,-0x4c(%rsp)
  je     L15
  lea    (%r14,%rbp,1),%edx
  movsd  (%rax),%xmm0
  movslq %edx,%rdx
  lea    (%rsi,%rdx,8),%rdi
  mulsd  (%rdi),%xmm0
  mulsd  (%rcx),%xmm0
  movsd  %xmm0,(%r12,%rdx,8)
  movsd  (%rax),%xmm0
  mulsd  (%rdi),%xmm0
  movsd  %xmm0,0x0(%r13,%rdx,8)
L15:
  addq   $0x1,-0x60(%rsp)
  mov    -0x58(%rsp),%rcx
  mov    -0x60(%rsp),%rbx
  mov    -0x48(%rsp),%rax
  add    -0x4c(%rsp),%r14d
  cmp    -0x28(%rsp),%rbx
  jne    L29
L77:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L14:
  add    -0x10(%rsp),%r10
  shl    $0x3,%r10
  nopl   0x0(%rax)
L43:
  movsd  (%rax),%xmm0
  mulsd  (%rsi,%rdx,1),%xmm0
  mulsd  (%rcx),%xmm0
  movsd  %xmm0,(%r12,%rdx,1)
  movsd  (%rax),%xmm0
  mulsd  (%rsi,%rdx,1),%xmm0
  movsd  %xmm0,0x0(%r13,%rdx,1)
  add    $0x8,%rdx
  cmp    %rdx,%r10
  jne    L43
  jmp    L15
  xchg   %ax,%ax
L8:
  add    $0x8,%rcx
  add    $0x8,%rax
  mov    %rcx,-0x58(%rsp)
  mov    %rax,-0x48(%rsp)
  jmp    L15
```

assistant: ```c
void VolumeAveragedVOFCoefficientsEvaluate(int nPoints, int nSpace, double eps, double *v, double *phi, double *porosity, double *u, double *m, double *dm, double *f, double *df) {
  int i, I;
  for (i = 0; i < nPoints; i++) {
    m[i] = u[i] * porosity[i];
    dm[i] = porosity[i];
    for (I = 0; I < nSpace; I++) {
      f[i * nSpace + I] = porosity[i] * v[i * nSpace + I] * u[i];
      df[i * nSpace + I] = porosity[i] * v[i * nSpace + I];
    }
  }
}
```