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

```asm
<kEpsilon_epsilon_3D_Evaluate_sd>:
  endbr64
  push   %r15
  movslq %edi,%rax
  push   %r14
  push   %r13
  push   %r12
  mov    %r9,%r12
  push   %rbp
  mov    %r8,%rbp
  push   %rbx
  mov    0x38(%rsp),%r14
  mov    0x68(%rsp),%r10
  mov    0x70(%rsp),%r11
  mov    0x78(%rsp),%r8
  mov    0x80(%rsp),%r9
  mov    %esi,-0xc(%rsp)
  test   %eax,%eax
  jle    L31
  movapd %xmm2,%xmm10
  mov    %esi,%edi
  mov    %rax,-0x8(%rsp)
  mov    %rdx,%rsi
  addsd  %xmm2,%xmm10
  movslq %edi,%r13
  lea    -0x1(%rdi),%edx
  mov    %rcx,%rbx
  mov    0x40(%rsp),%rcx
  movapd %xmm0,%xmm6
  shl    $0x3,%r13
  xor    %edi,%edi
  movq   D18(%rip),%xmm9
  movsd  D3(%rip),%xmm4
  xor    %r15d,%r15d
  lea    0x8(,%rdx,8),%rdx
  movsd  D20(%rip),%xmm7
  xorpd  %xmm9,%xmm1
  movapd %xmm1,%xmm8
  pxor   %xmm1,%xmm1
  nopl   (%rax)
L70:
  movsd  (%r14,%r15,8),%xmm12
  movsd  (%rcx),%xmm0
  movapd %xmm12,%xmm11
  movapd %xmm0,%xmm13
  mulsd  %xmm3,%xmm11
  addsd  %xmm4,%xmm13
  mulsd  %xmm12,%xmm11
  movapd %xmm11,%xmm12
  divsd  %xmm13,%xmm12
  comisd %xmm12,%xmm1
  ja     L37
  movapd %xmm0,%xmm13
  xorpd  %xmm9,%xmm11
  mulsd  %xmm0,%xmm13
  addsd  %xmm4,%xmm13
  divsd  %xmm13,%xmm11
L50:
  mov    0x48(%rsp),%rax
  movsd  %xmm0,(%rax,%r15,8)
  mov    0x50(%rsp),%rax
  movsd  %xmm7,(%rax,%r15,8)
  mov    -0xc(%rsp),%eax
  test   %eax,%eax
  jle    L98
  mov    %rdi,%rax
  cs nopw 0x0(%rax,%rax,1)
L49:
  movsd  (%rcx),%xmm0
  mulsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,(%r10,%rax,1)
  movsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,(%r11,%rax,1)
  add    $0x8,%rax
  cmp    %rax,%rdx
  jne    L49
  movsd  (%rcx),%xmm0
  mov    0x58(%rsp),%rax
  divsd  %xmm6,%xmm12
  movsd  %xmm0,(%rax,%r15,8)
  mov    0x60(%rsp),%rax
  movsd  %xmm7,(%rax,%r15,8)
  mov    %rdi,%rax
  divsd  %xmm6,%xmm11
  addsd  %xmm5,%xmm12
  nop
L96:
  movsd  %xmm12,(%r8,%rax,1)
  movsd  %xmm11,(%r9,%rax,1)
  add    $0x8,%rax
  cmp    %rax,%rdx
  jne    L96
L33:
  movsd  (%rbx,%rdi,1),%xmm0
  movsd  0x8(%rbp,%rdi,1),%xmm15
  movsd  0x10(%r12,%rdi,1),%xmm14
  movsd  0x8(%rbx,%rdi,1),%xmm13
  mulsd  %xmm0,%xmm0
  addsd  0x0(%rbp,%rdi,1),%xmm13
  movsd  0x10(%rbx,%rdi,1),%xmm12
  mulsd  %xmm15,%xmm15
  addsd  (%r12,%rdi,1),%xmm12
  movsd  0x10(%rbp,%rdi,1),%xmm11
  mulsd  %xmm14,%xmm14
  addsd  0x8(%r12,%rdi,1),%xmm11
  mulsd  %xmm13,%xmm13
  mulsd  %xmm12,%xmm12
  mulsd  %xmm11,%xmm11
  addsd  %xmm15,%xmm0
  addsd  %xmm14,%xmm0
  addsd  %xmm0,%xmm0
  addsd  %xmm13,%xmm0
  addsd  %xmm12,%xmm0
  movsd  (%rcx),%xmm12
  movapd %xmm12,%xmm13
  mulsd  %xmm2,%xmm13
  addsd  %xmm11,%xmm0
  movsd  (%r14,%r15,8),%xmm11
  movapd %xmm11,%xmm14
  addsd  %xmm4,%xmm14
  mulsd  %xmm12,%xmm13
  divsd  %xmm14,%xmm13
  comisd %xmm13,%xmm1
  ja     L29
  mulsd  %xmm10,%xmm12
  divsd  %xmm14,%xmm12
L95:
  mulsd  %xmm8,%xmm11
  add    $0x8,%rcx
  add    %r13,%rdi
  add    %r13,%rdx
  mov    0x88(%rsp),%rax
  mulsd  %xmm0,%xmm11
  addsd  %xmm13,%xmm11
  movsd  %xmm11,(%rax,%r15,8)
  mov    0x90(%rsp),%rax
  movsd  %xmm12,(%rax,%r15,8)
  add    $0x1,%r15
  cmp    %r15,-0x8(%rsp)
  jne    L70
L31:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L29:
  movapd %xmm1,%xmm12
  movapd %xmm1,%xmm13
  jmp    L95
  nopl   0x0(%rax)
L37:
  movapd %xmm1,%xmm11
  movapd %xmm1,%xmm12
  jmp    L50
  nopw   0x0(%rax,%rax,1)
L98:
  movsd  (%rcx),%xmm0
  mov    0x58(%rsp),%rax
  movsd  %xmm0,(%rax,%r15,8)
  mov    0x60(%rsp),%rax
  movsd  %xmm7,(%rax,%r15,8)
  jmp    L33
```

tool_calls: name: retrieve_data
arguments: label: D18
type: f32[2]
name: retrieve_data
arguments: label: D3
type: byte[8]
name: retrieve_data
arguments: label: D20
type: f64

tool: {"D18": [0.0,-0.0]}

tool: {"D3": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xb0,0x3e]}

tool: D20: 1

assistant: ```c
void kEpsilon_epsilon_3D_Evaluate_sd(int nPoints, int nSpace, double sigma_e, double c_1, double c_2, double c_mu, double c_e, double nu, double *velocity, double *gradu, double *gradv, double *gradw, double *k, double *epsilon, double *m_e, double *dm_e, double *phi_e, double *dphi_e, double *f_e, double *df_e, double *a_e, double *da_e_de, double *r_e, double *dr_e_de) {
  int i, I;
  double nu_t, dnu_t_de, PiD4, eval, kval, disp, ddisp_de;
  const double div_eps = 1.0e-6;
  for (i = 0; i < nPoints; i++) {
    eval = epsilon[i] >= 0.0 ? epsilon[i] : 0.0;
    kval = k[i] >= 0.0 ? k[i] : 0.0;
    nu_t = c_mu * k[i] * k[i] / (epsilon[i] + div_eps);
    dnu_t_de = -c_mu * k[i] * k[i] / (epsilon[i] * epsilon[i] + div_eps);
    if (nu_t < 0.0) {
      nu_t = 0.0;
      dnu_t_de = 0.0;
    }
    m_e[i] = epsilon[i];
    dm_e[i] = 1.0;
    for (I = 0; I < nSpace; I++) {
      f_e[i * nSpace + I] = epsilon[i] * velocity[i * nSpace + I];
      df_e[i * nSpace + I] = velocity[i * nSpace + I];
    }
    phi_e[i] = epsilon[i];
    dphi_e[i] = 1.0;
    for (I = 0; I < nSpace; I++) {
      a_e[i * nSpace + I] = nu_t / sigma_e + nu;
      da_e_de[i * nSpace + I] = dnu_t_de / sigma_e;
    }
    PiD4 = 2.0 * (gradu[i * nSpace + 0] * gradu[i * nSpace + 0] + gradv[i * nSpace + 1] * gradv[i * nSpace + 1] + gradw[i * nSpace + 2] * gradw[i * nSpace + 2]) + (gradu[i * nSpace + 1] + gradv[i * nSpace + 0]) * (gradu[i * nSpace + 1] + gradv[i * nSpace + 0]) + (gradu[i * nSpace + 2] + gradw[i * nSpace + 0]) * (gradu[i * nSpace + 2] + gradw[i * nSpace + 0]) + (gradv[i * nSpace + 2] + gradw[i * nSpace + 1]) * (gradv[i * nSpace + 2] + gradw[i * nSpace + 1]);
    disp = c_2 * epsilon[i] * epsilon[i] / (k[i] + div_eps);
    ddisp_de = 2.0 * c_2 * epsilon[i] / (k[i] + div_eps);
    if (disp < 0.0) {
      disp = 0.0;
      ddisp_de = 0.0;
    }
    r_e[i] = -c_1 * k[i] * PiD4 + disp;
    dr_e_de[i] = ddisp_de;
  }
}
```