user: What would the C source look like for the assembly code provided here:

```asm
<wave13pt>:
  endbr64
  push   %r15
  mov    %esi,%eax
  lea    -0x2(%rdx),%esi
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %edi,-0x54(%rsp)
  mov    %esi,-0x10(%rsp)
  cmp    $0x2,%esi
  jle    L53
  mov    %edi,%esi
  mov    %edi,%ebx
  movl   $0x0,-0x3c(%rsp)
  movapd %xmm0,%xmm4
  imul   %eax,%esi
  sub    $0x2,%eax
  mov    %rcx,%r13
  mov    %r8,%r15
  mov    %eax,-0x40(%rsp)
  mov    %r9,%r14
  movl   $0x2,-0x28(%rsp)
  lea    (%rsi,%rsi,1),%eax
  mov    %esi,-0x14(%rsp)
  lea    (%rax,%rsi,1),%edi
  mov    %esi,-0x20(%rsp)
  mov    %edi,-0x18(%rsp)
  mov    %edi,-0x1c(%rsp)
  lea    0x0(,%rsi,4),%edi
  mov    %edi,-0x24(%rsp)
  lea    (%rbx,%rbx,2),%edi
  sub    $0x5,%ebx
  mov    %rbx,-0x8(%rsp)
  lea    0x8(%r8),%rbx
  mov    %edi,-0xc(%rsp)
  mov    %rbx,-0x30(%rsp)
  nopl   0x0(%rax,%rax,1)
L11:
  cmpl   $0x2,-0x40(%rsp)
  jle    L58
  mov    -0x8(%rsp),%rbx
  cltq
  movslq -0x1c(%rsp),%rbp
  movl   $0x2,-0x48(%rsp)
  movslq -0x24(%rsp),%r11
  mov    -0x54(%rsp),%esi
  mov    %rax,-0x50(%rsp)
  add    %rax,%rbx
  mov    -0xc(%rsp),%r12d
  sub    %rax,%rbp
  mov    %rbx,-0x38(%rsp)
  movslq -0x20(%rsp),%rbx
  sub    %rax,%r11
  sub    %rax,%rbx
  nopl   0x0(%rax,%rax,1)
L2:
  addl   $0x1,-0x48(%rsp)
  cmpl   $0x4,-0x54(%rsp)
  mov    -0x54(%rsp),%eax
  jle    L74
  lea    (%rax,%rsi,1),%ecx
  mov    -0x50(%rsp),%rdi
  mov    -0x38(%rsp),%r10
  mov    %ecx,-0x44(%rsp)
  movslq %ecx,%rcx
  lea    (%rcx,%rdi,1),%rax
  lea    0x2(%rcx,%rdi,1),%rdx
  lea    (%r10,%rcx,1),%rdi
  mov    -0x30(%rsp),%r10
  neg    %rcx
  lea    (%r15,%rax,8),%rax
  shl    $0x3,%rcx
  shl    $0x3,%rdx
  lea    (%r10,%rdi,8),%r10
  movslq %r12d,%rdi
  add    -0x54(%rsp),%r12d
  lea    (%rcx,%rdi,8),%r9
  movslq %esi,%rdi
  sub    -0x54(%rsp),%esi
  lea    (%rcx,%rdi,8),%r8
  movslq %esi,%rsi
  movslq %r12d,%rdi
  add    %r15,%r9
  lea    (%rcx,%rdi,8),%rdi
  lea    (%rcx,%rsi,8),%rsi
  add    %r15,%r8
  movslq -0x3c(%rsp),%rcx
  add    %r15,%rdi
  add    %r15,%rsi
  sub    -0x50(%rsp),%rcx
  nopl   0x0(%rax,%rax,1)
L61:
  movsd  0x18(%rax),%xmm0
  addsd  0x8(%rax),%xmm0
  addsd  (%r9,%rdx,1),%xmm0
  addsd  (%r8,%rdx,1),%xmm0
  addsd  0x10(%rax,%rbp,8),%xmm0
  addsd  0x10(%rax,%rbx,8),%xmm0
  movsd  0x10(%rax),%xmm3
  mulsd  %xmm1,%xmm0
  mulsd  %xmm4,%xmm3
  subsd  0x0(%r13,%rdx,1),%xmm3
  addsd  %xmm3,%xmm0
  movsd  0x20(%rax),%xmm3
  addsd  (%rax),%xmm3
  addsd  (%rdi,%rdx,1),%xmm3
  addsd  (%rsi,%rdx,1),%xmm3
  addsd  0x10(%rax,%r11,8),%xmm3
  addsd  0x10(%rax,%rcx,8),%xmm3
  add    $0x8,%rax
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%r14,%rdx,1)
  add    $0x8,%rdx
  cmp    %r10,%rax
  jne    L61
  mov    -0x44(%rsp),%esi
  mov    -0x48(%rsp),%edi
  cmp    %edi,-0x40(%rsp)
  jne    L2
L58:
  mov    -0x18(%rsp),%ecx
  mov    -0x14(%rsp),%ebx
  addl   $0x1,-0x28(%rsp)
  mov    -0x28(%rsp),%esi
  add    %ebx,-0x1c(%rsp)
  mov    %ecx,%eax
  add    %ebx,-0x20(%rsp)
  add    %ebx,-0x24(%rsp)
  add    %ebx,-0x3c(%rsp)
  cmp    -0x10(%rsp),%esi
  je     L53
  add    %ecx,%ebx
  mov    %ebx,-0x18(%rsp)
  jmp    L11
  nopl   0x0(%rax,%rax,1)
L74:
  add    %eax,%esi
  mov    -0x48(%rsp),%edi
  add    %eax,%r12d
  mov    %esi,-0x44(%rsp)
  mov    -0x44(%rsp),%esi
  cmp    %edi,-0x40(%rsp)
  jne    L2
  jmp    L58
L53:
  pop    %rbx
  xor    %eax,%eax
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
int wave13pt(int nx, int ny, int ns, const double c0, const double c1, const double c2, double *w0, double *w1, double *w2) {
  int i, j, k;
  unsigned int szarray = (unsigned int)nx * ny * ns;
#pragma acc data present(w0 [0:szarray], w1 [0:szarray], w2 [0:szarray])
  {
#pragma acc kernels
    {
#pragma acc loop independent vector(2)
      for (k = 2; k < ns - 2; k++) {
#pragma acc loop independent gang(31) vector(4)
        for (j = 2; j < ny - 2; j++) {
#pragma acc loop independent gang(2) vector(64)
          for (i = 2; i < nx - 2; i++) {
            (w2[(i) + nx * (j) + nx * ny * (k)]) = c0 * (w1[(i) + nx * (j) + nx * ny * (k)]) - (w0[(i) + nx * (j) + nx * ny * (k)]) + c1 * ((w1[(i + 1) + nx * (j) + nx * ny * (k)]) + (w1[(i - 1) + nx * (j) + nx * ny * (k)]) + (w1[(i) + nx * (j + 1) + nx * ny * (k)]) + (w1[(i) + nx * (j - 1) + nx * ny * (k)]) + (w1[(i) + nx * (j) + nx * ny * (k + 1)]) + (w1[(i) + nx * (j) + nx * ny * (k - 1)])) + c2 * ((w1[(i + 2) + nx * (j) + nx * ny * (k)]) + (w1[(i - 2) + nx * (j) + nx * ny * (k)]) + (w1[(i) + nx * (j + 2) + nx * ny * (k)]) + (w1[(i) + nx * (j - 2) + nx * ny * (k)]) + (w1[(i) + nx * (j) + nx * ny * (k + 2)]) + (w1[(i) + nx * (j) + nx * ny * (k - 2)]));
          }
        }
      }
    }
  }
  return 0;
}
```