user: Translate the below assembly code into C programming language:

```asm
<kEpsilon_3D_Evaluate_sd>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %esi,%eax
  mov    %esi,-0x30(%rsp)
  mov    %rdx,%rsi
  mov    %rcx,%rdx
  mov    %r9,%rcx
  mov    0x88(%rsp),%r9
  mov    0x90(%rsp),%r11
  mov    0xa8(%rsp),%rbx
  mov    0xb0(%rsp),%rbp
  mov    0xb8(%rsp),%r12
  mov    0xc0(%rsp),%r13
  mov    0xc8(%rsp),%r14
  mov    0xd0(%rsp),%r15
  test   %edi,%edi
  jle    L73
  movapd %xmm0,%xmm12
  movapd %xmm4,%xmm5
  addsd  %xmm4,%xmm5
  movsd  %xmm5,-0x8(%rsp)
  mov    0x40(%rsp),%r10
  mov    %r10,-0x58(%rsp)
  mov    0x38(%rsp),%r10
  mov    %r10,-0x50(%rsp)
  cltq
  shl    $0x3,%rax
  mov    %rax,-0x18(%rsp)
  mov    %rdx,-0x48(%rsp)
  mov    %r8,-0x40(%rsp)
  mov    %rcx,-0x38(%rsp)
  mov    %edi,%eax
  mov    %rax,-0x10(%rsp)
  movl   $0x0,-0x2c(%rsp)
  mov    $0x0,%edi
  movq   D70(%rip),%xmm13
  xorpd  %xmm13,%xmm2
  jmp    L2
L109:
  pxor   %xmm5,%xmm5
  movapd %xmm5,%xmm8
  movapd %xmm5,%xmm7
  jmp    L14
L86:
  pxor   %xmm7,%xmm7
  movapd %xmm7,%xmm5
  movapd %xmm7,%xmm8
L51:
  mulsd  %xmm2,%xmm0
  mulsd  %xmm9,%xmm0
  addsd  %xmm8,%xmm0
  mov    0xf0(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  mulsd  %xmm2,%xmm9
  addsd  %xmm5,%xmm9
  mov    0xf8(%rsp),%rax
  movsd  %xmm9,(%rax,%rdi,8)
  mov    0x100(%rsp),%rax
  movsd  %xmm7,(%rax,%rdi,8)
  add    $0x1,%rdi
  addq   $0x8,-0x58(%rsp)
  addq   $0x8,-0x50(%rsp)
  mov    -0x18(%rsp),%rax
  add    %rax,-0x48(%rsp)
  add    %rax,-0x40(%rsp)
  add    %rax,-0x38(%rsp)
  mov    -0x30(%rsp),%ecx
  add    %ecx,-0x2c(%rsp)
  cmp    -0x10(%rsp),%rdi
  je     L73
L2:
  mov    -0x58(%rsp),%rax
  mov    %rax,%r8
  movsd  (%rax),%xmm0
  mov    -0x50(%rsp),%rax
  mov    %rax,%r10
  movsd  (%rax),%xmm9
  movapd %xmm9,%xmm5
  mulsd  %xmm4,%xmm5
  mulsd  %xmm9,%xmm5
  movsd  D80(%rip),%xmm10
  addsd  %xmm0,%xmm10
  movapd %xmm5,%xmm7
  divsd  %xmm10,%xmm7
  pxor   %xmm14,%xmm14
  comisd %xmm7,%xmm14
  ja     L109
  movsd  -0x8(%rsp),%xmm8
  mulsd  %xmm9,%xmm8
  divsd  %xmm10,%xmm8
  xorpd  D70(%rip),%xmm5
  mulsd  %xmm0,%xmm0
  addsd  D80(%rip),%xmm0
  divsd  %xmm0,%xmm5
L14:
  mov    0x48(%rsp),%rax
  movsd  %xmm9,(%rax,%rdi,8)
  movsd  D32(%rip),%xmm0
  mov    0x50(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  movsd  (%r8),%xmm9
  mov    0x58(%rsp),%rax
  movsd  %xmm9,(%rax,%rdi,8)
  mov    0x60(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  mov    -0x30(%rsp),%ecx
  test   %ecx,%ecx
  jle    L0
  movslq -0x2c(%rsp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    %ecx,%ecx
  add    %rax,%rcx
  shl    $0x3,%rcx
  mov    %rdx,%rax
  mov    %rdx,-0x28(%rsp)
  mov    %rdi,-0x20(%rsp)
  mov    0x98(%rsp),%rdx
  mov    0xa0(%rsp),%rdi
L33:
  movsd  (%r10),%xmm0
  mulsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,(%r9,%rax,1)
  movsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,(%r11,%rax,1)
  mulsd  (%r8),%xmm0
  movsd  %xmm0,(%rdx,%rax,1)
  movsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,(%rdi,%rax,1)
  add    $0x8,%rax
  cmp    %rcx,%rax
  jne    L33
  mov    -0x28(%rsp),%rdx
  mov    -0x20(%rsp),%rdi
  movsd  (%r10),%xmm0
  mov    0x68(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  movsd  D32(%rip),%xmm0
  mov    0x70(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  movsd  (%r8),%xmm9
  mov    0x78(%rsp),%rax
  movsd  %xmm9,(%rax,%rdi,8)
  mov    0x80(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  movapd %xmm7,%xmm9
  divsd  %xmm12,%xmm9
  addsd  %xmm6,%xmm9
  movapd %xmm8,%xmm15
  divsd  %xmm12,%xmm15
  movapd %xmm5,%xmm14
  divsd  %xmm12,%xmm14
  movapd %xmm7,%xmm0
  divsd  %xmm1,%xmm0
  addsd  %xmm6,%xmm0
  movapd %xmm8,%xmm11
  divsd  %xmm1,%xmm11
  movapd %xmm5,%xmm10
  divsd  %xmm1,%xmm10
L54:
  movsd  %xmm9,(%rbx,%rdx,1)
  movsd  %xmm15,0x0(%rbp,%rdx,1)
  movsd  %xmm14,(%r12,%rdx,1)
  movsd  %xmm0,0x0(%r13,%rdx,1)
  movsd  %xmm11,(%r14,%rdx,1)
  movsd  %xmm10,(%r15,%rdx,1)
  add    $0x8,%rdx
  cmp    %rcx,%rdx
  jne    L54
L8:
  mov    -0x48(%rsp),%rax
  movsd  (%rax),%xmm15
  mov    -0x40(%rsp),%rcx
  movsd  0x8(%rcx),%xmm0
  mov    -0x38(%rsp),%rdx
  movsd  0x10(%rdx),%xmm14
  movsd  0x8(%rax),%xmm11
  addsd  (%rcx),%xmm11
  movsd  0x10(%rax),%xmm10
  addsd  (%rdx),%xmm10
  movsd  0x10(%rcx),%xmm9
  addsd  0x8(%rdx),%xmm9
  mulsd  %xmm15,%xmm15
  mulsd  %xmm0,%xmm0
  addsd  %xmm15,%xmm0
  mulsd  %xmm14,%xmm14
  addsd  %xmm14,%xmm0
  addsd  %xmm0,%xmm0
  mulsd  %xmm11,%xmm11
  addsd  %xmm11,%xmm0
  mulsd  %xmm10,%xmm10
  addsd  %xmm0,%xmm10
  mulsd  %xmm9,%xmm9
  addsd  %xmm10,%xmm9
  xorpd  %xmm13,%xmm7
  mulsd  %xmm9,%xmm7
  addsd  (%r8),%xmm7
  mov    0xd8(%rsp),%rax
  movsd  %xmm7,(%rax,%rdi,8)
  xorpd  %xmm13,%xmm8
  mulsd  %xmm9,%xmm8
  mov    0xe0(%rsp),%rax
  movsd  %xmm8,(%rax,%rdi,8)
  mulsd  %xmm9,%xmm5
  movsd  D32(%rip),%xmm0
  subsd  %xmm5,%xmm0
  mov    0xe8(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  movsd  (%r8),%xmm10
  movapd %xmm10,%xmm5
  mulsd  %xmm3,%xmm5
  mulsd  %xmm10,%xmm5
  movsd  (%r10),%xmm0
  movapd %xmm0,%xmm11
  addsd  D80(%rip),%xmm11
  movapd %xmm5,%xmm8
  divsd  %xmm11,%xmm8
  pxor   %xmm7,%xmm7
  comisd %xmm8,%xmm7
  ja     L86
  xorpd  D70(%rip),%xmm5
  movapd %xmm0,%xmm7
  mulsd  %xmm0,%xmm7
  addsd  D80(%rip),%xmm7
  divsd  %xmm7,%xmm5
  movapd %xmm3,%xmm7
  addsd  %xmm3,%xmm7
  mulsd  %xmm10,%xmm7
  divsd  %xmm11,%xmm7
  jmp    L51
L0:
  movsd  (%r10),%xmm0
  mov    0x68(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  movsd  D32(%rip),%xmm0
  mov    0x70(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  movsd  (%r8),%xmm9
  mov    0x78(%rsp),%rax
  movsd  %xmm9,(%rax,%rdi,8)
  mov    0x80(%rsp),%rax
  movsd  %xmm0,(%rax,%rdi,8)
  jmp    L8
L73:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: retrieve_data
arguments: label: D70
type: f32[4]
name: retrieve_data
arguments: label: D80
type: byte[8]
name: retrieve_data
arguments: label: D32
type: f64

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

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

tool: D32: 1

assistant: ```c
void kEpsilon_3D_Evaluate_sd(int nPoints, int nSpace, double sigma_k, 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_k, double *dm_k, double *m_e, double *dm_e, double *phi_k, double *dphi_k, double *phi_e, double *dphi_e, double *f_k, double *df_k, double *f_e, double *df_e, double *a_k, double *da_k_dk, double *da_k_de, double *a_e, double *da_e_dk, double *da_e_de, double *r_k, double *dr_k_dk, double *dr_k_de, double *r_e, double *dr_e_dk, double *dr_e_de) {
  int i, I;
  double nu_t, dnu_t_dk, dnu_t_de, PiD4, eval, kval, disp, ddisp_dk, 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_dk = 2.0 * c_mu * 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_dk = 0.0;
      dnu_t_de = 0.0;
    }
    m_k[i] = k[i];
    dm_k[i] = 1.0;
    m_e[i] = epsilon[i];
    dm_e[i] = 1.0;
    for (I = 0; I < nSpace; I++) {
      f_k[i * nSpace + I] = k[i] * velocity[i * nSpace + I];
      df_k[i * nSpace + I] = velocity[i * nSpace + I];
      f_e[i * nSpace + I] = epsilon[i] * velocity[i * nSpace + I];
      df_e[i * nSpace + I] = velocity[i * nSpace + I];
    }
    phi_k[i] = k[i];
    dphi_k[i] = 1.0;
    phi_e[i] = epsilon[i];
    dphi_e[i] = 1.0;
    for (I = 0; I < nSpace; I++) {
      a_k[i * nSpace + I] = nu_t / sigma_k + nu;
      da_k_dk[i * nSpace + I] = dnu_t_dk / sigma_k;
      da_k_de[i * nSpace + I] = dnu_t_de / sigma_k;
      a_e[i * nSpace + I] = nu_t / sigma_e + nu;
      da_e_dk[i * nSpace + I] = dnu_t_dk / sigma_e;
      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]);
    r_k[i] = -nu_t * PiD4 + epsilon[i];
    dr_k_dk[i] = -dnu_t_dk * PiD4;
    dr_k_de[i] = -dnu_t_de * PiD4 + 1.0;
    disp = c_2 * epsilon[i] * epsilon[i] / (k[i] + div_eps);
    ddisp_dk = -c_2 * epsilon[i] * epsilon[i] / (k[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_dk = 0.0;
      ddisp_de = 0.0;
    }
    r_e[i] = -c_1 * k[i] * PiD4 + disp;
    dr_e_dk[i] = -c_1 * PiD4 + ddisp_dk;
    dr_e_de[i] = ddisp_de;
  }
}
```