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

```asm
<Stokes_2D_Evaluate>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %edi,%r11d
  mov    %rsi,%rdx
  mov    0x58(%rsp),%rsi
  mov    0x70(%rsp),%rdi
  mov    0x78(%rsp),%r10
  mov    0x80(%rsp),%rbx
  mov    0x88(%rsp),%rbp
  mov    0x90(%rsp),%r12
  mov    0x98(%rsp),%r13
  mov    0xa0(%rsp),%r14
  mov    0xa8(%rsp),%r15
  test   %r11d,%r11d
  jle    L8
  movapd %xmm0,%xmm3
  movsd  D10(%rip),%xmm4
  divsd  %xmm0,%xmm4
  mov    %r11d,%r11d
  shl    $0x3,%r11
  mov    $0x0,%eax
  movsd  D10(%rip),%xmm0
  movq   D5(%rip),%xmm5
  mov    %r11,-0x8(%rsp)
L14:
  movsd  (%r8,%rax,1),%xmm2
  mov    0x38(%rsp),%r11
  movsd  %xmm2,(%r11,%rax,1)
  mov    0x40(%rsp),%r11
  movsd  %xmm0,(%r11,%rax,1)
  movsd  (%r9,%rax,1),%xmm2
  mov    0x48(%rsp),%r11
  movsd  %xmm2,(%r11,%rax,1)
  mov    0x50(%rsp),%r11
  movsd  %xmm0,(%r11,%rax,1)
  movsd  (%r8,%rax,1),%xmm2
  movsd  %xmm2,(%rsi,%rax,2)
  movsd  (%r9,%rax,1),%xmm2
  movsd  %xmm2,0x8(%rsi,%rax,2)
  mov    0x60(%rsp),%r11
  movsd  %xmm0,(%r11,%rax,2)
  mov    0x68(%rsp),%r11
  movsd  %xmm0,0x8(%r11,%rax,2)
  movsd  %xmm1,(%rdi,%rax,4)
  movsd  %xmm1,0x18(%rdi,%rax,4)
  movsd  %xmm1,(%r10,%rax,4)
  movsd  %xmm1,0x18(%r10,%rax,4)
  movsd  (%rdx),%xmm2
  xorpd  %xmm5,%xmm2
  movsd  %xmm2,(%rbx,%rax,1)
  movsd  0x8(%rdx),%xmm2
  xorpd  %xmm5,%xmm2
  movsd  %xmm2,0x0(%rbp,%rax,1)
  movsd  (%rcx,%rax,2),%xmm2
  divsd  %xmm3,%xmm2
  movsd  %xmm2,(%r12,%rax,1)
  movsd  %xmm4,0x0(%r13,%rax,2)
  movsd  0x8(%rcx,%rax,2),%xmm2
  divsd  %xmm3,%xmm2
  movsd  %xmm2,(%r14,%rax,1)
  movsd  %xmm4,0x8(%r15,%rax,2)
  add    $0x8,%rax
  cmp    -0x8(%rsp),%rax
  jne    L14
L8:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_label: D10
type: f64
name: read_data
arguments: data_label: D5
type: byte[8]

tool: D10: 1

tool: {"D5": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

assistant: ```c
void Stokes_2D_Evaluate(const int nPoints, const double rho, const double nu, const double *g, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p) {
  int k;
  for (k = 0; k < nPoints; k++) {
    mom_u_acc[k] = u[k];
    dmom_u_acc_u[k] = 1.0;
    mom_v_acc[k] = v[k];
    dmom_v_acc_v[k] = 1.0;
    mass_adv[k * 2 + 0] = u[k];
    mass_adv[k * 2 + 1] = v[k];
    dmass_adv_u[k * 2 + 0] = 1.0;
    dmass_adv_v[k * 2 + 1] = 1.0;
    mom_u_diff_ten[k * 4 + 0] = nu;
    mom_u_diff_ten[k * 4 + 3] = nu;
    mom_v_diff_ten[k * 4 + 0] = nu;
    mom_v_diff_ten[k * 4 + 3] = nu;
    mom_u_source[k] = -g[0];
    mom_v_source[k] = -g[1];
    mom_u_ham[k] = grad_p[k * 2 + 0] / rho;
    dmom_u_ham_grad_p[k * 2 + 0] = 1.0 / rho;
    mom_v_ham[k] = grad_p[k * 2 + 1] / rho;
    dmom_v_ham_grad_p[k * 2 + 1] = 1.0 / rho;
  }
}
```