user: Please show me the C implementation that matches the assembly code listed below:

```asm
<acd_2d_2_d>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x38(%rsp),%r15
  mov    0x40(%rsp),%r13
  mov    0x48(%rsp),%r11
  mov    0x4(%r15),%r14d
  cmp    0x4(%r13),%r14d
  jg     L1
  mov    %rsi,%r10
  mov    %rdx,%rbx
  mov    %rcx,%r12
  movaps %xmm0,%xmm2
  movslq %r14d,%rax
  shl    $0x3,%rax
  lea    (%rdi,%rax,1),%rcx
  add    %rax,%r10
  add    %rax,%r12
  lea    (%r9,%rax,1),%rbp
  add    %rax,%r8
  add    %rax,%rbx
L24:
  mov    (%r15),%edi
  cmp    0x0(%r13),%edi
  jg     L25
  movslq %edi,%rdx
  shl    $0x2,%rdx
  mov    %r14d,-0x4(%rsp)
L22:
  mov    (%rcx),%rsi
  movaps %xmm2,%xmm3
  mulss  (%rsi,%rdx,1),%xmm3
  movss  (%r11),%xmm4
  mov    %rdx,%rax
  add    $0x4,%rdx
  movss  0x4(%r11),%xmm5
  mov    -0x8(%rcx),%r9
  mov    0x8(%rcx),%r14
  movss  (%r14,%rax,1),%xmm1
  addss  (%r9,%rax,1),%xmm1
  mulss  %xmm5,%xmm1
  movss  (%rsi,%rdx,1),%xmm0
  addss  -0x4(%rsi,%rax,1),%xmm0
  mulss  %xmm4,%xmm0
  addss  %xmm3,%xmm0
  addss  %xmm0,%xmm1
  mov    (%r10),%r9
  movss  (%r9,%rax,1),%xmm3
  mov    %rax,%rsi
  add    (%r12),%rsi
  movss  (%r9,%rdx,1),%xmm0
  addss  -0x4(%r9,%rax,1),%xmm0
  mulss  %xmm4,%xmm0
  movaps %xmm3,%xmm4
  mulss  %xmm2,%xmm4
  addss  %xmm4,%xmm0
  mov    -0x8(%r10),%r9
  mov    0x8(%r10),%r14
  movss  (%r14,%rax,1),%xmm4
  addss  (%r9,%rax,1),%xmm4
  mulss  %xmm5,%xmm4
  addss  %xmm4,%xmm0
  mov    (%r8),%r9
  mulss  (%r9,%rax,1),%xmm0
  cvtss2sd %xmm0,%xmm0
  mov    0x0(%rbp),%r9
  movaps %xmm1,%xmm4
  mulss  (%r9,%rax,1),%xmm4
  cvtss2sd %xmm4,%xmm4
  cvtss2sd %xmm3,%xmm3
  addsd  %xmm3,%xmm3
  pxor   %xmm5,%xmm5
  cvtss2sd (%rsi),%xmm5
  subsd  %xmm5,%xmm3
  addsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rsi)
  mov    %rax,%rsi
  add    (%rbx),%rsi
  mov    (%rcx),%r9
  pxor   %xmm0,%xmm0
  cvtss2sd (%r9,%rax,1),%xmm0
  addsd  %xmm0,%xmm0
  pxor   %xmm3,%xmm3
  cvtss2sd (%rsi),%xmm3
  subsd  %xmm3,%xmm0
  mov    (%r8),%r9
  mulss  (%r9,%rax,1),%xmm1
  cvtss2sd %xmm1,%xmm1
  addsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rsi)
  add    $0x1,%edi
  cmp    %edi,0x0(%r13)
  jge    L22
  mov    -0x4(%rsp),%r14d
L25:
  add    $0x1,%r14d
  add    $0x8,%rcx
  add    $0x8,%r10
  add    $0x8,%r12
  add    $0x8,%rbp
  add    $0x8,%r8
  add    $0x8,%rbx
  cmp    %r14d,0x4(%r13)
  jge    L24
L1:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void acd_2d_2_d(float **uc, float **ucd, float **up, float **upd, float **csq, float **csqd, int *s, int *e, float c0, float *c1) {
  int i0, i1;
  for (i1 = s[1]; i1 < e[1] + 1; ++i1)
    for (i0 = s[0]; i0 < e[0] + 1; ++i0) {
      float lap = (c0 * uc[i1][i0] + c1[0] * (uc[i1][i0 + 1] + uc[i1][i0 - 1]) + c1[1] * (uc[i1 + 1][i0] + uc[i1 - 1][i0]));
      float lapd = c0 * ucd[i1][i0] + c1[0] * (ucd[i1][i0 + 1] + ucd[i1][i0 - 1]) + c1[1] * (ucd[i1 + 1][i0] + ucd[i1 - 1][i0]);
      upd[i1][i0] = 2.0 * ucd[i1][i0] - upd[i1][i0] + csqd[i1][i0] * lap + csq[i1][i0] * lapd;
      up[i1][i0] = 2.0 * uc[i1][i0] - up[i1][i0] + csq[i1][i0] * lap;
    }
}
```