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

```asm
<wave13pt>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %edi,-0x58(%rsp)
  mov    %esi,-0x10(%rsp)
  mov    %r8,-0x60(%rsp)
  cmp    $0x4,%edx
  jle    L16
  mov    %edi,%ebx
  mov    %edx,%eax
  movapd %xmm0,%xmm4
  mov    %rcx,%r13
  mov    %r9,%r14
  imul   %esi,%edi
  mov    %edi,-0x18(%rsp)
  lea    (%rdi,%rdi,1),%ecx
  mov    %ecx,-0x2c(%rsp)
  add    %edi,%ecx
  mov    %ecx,-0x28(%rsp)
  lea    0x0(,%rdi,4),%ecx
  mov    %ecx,-0x20(%rsp)
  sub    $0x2,%eax
  mov    %eax,-0x14(%rsp)
  mov    %ebx,%eax
  lea    (%rbx,%rbx,2),%ebx
  mov    %ebx,-0xc(%rsp)
  mov    %edi,-0x24(%rsp)
  movl   $0x0,-0x1c(%rsp)
  movl   $0x2,-0x30(%rsp)
  lea    -0x2(%rsi),%ebx
  mov    %ebx,-0x48(%rsp)
  lea    -0x5(%rax),%eax
  mov    %rax,-0x8(%rsp)
  lea    0x8(%r8),%rax
  mov    %rax,-0x38(%rsp)
  jmp    L64
L55:
  mov    -0x58(%rsp),%eax
  lea    (%rax,%rbp,1),%ecx
  movslq %ecx,%rcx
  mov    -0x50(%rsp),%rdx
  lea    (%rcx,%rdx,1),%rax
  mov    -0x60(%rsp),%r15
  lea    (%r15,%rax,8),%rax
  lea    0x2(%rcx,%rdx,1),%rdx
  shl    $0x3,%rdx
  mov    -0x40(%rsp),%r9
  lea    (%r9,%rcx,1),%rsi
  mov    -0x38(%rsp),%r8
  lea    (%r8,%rsi,8),%r9
  neg    %rcx
  shl    $0x3,%rcx
  movslq %r12d,%rsi
  lea    (%rcx,%rsi,8),%r8
  add    %r15,%r8
  movslq %ebp,%rsi
  lea    (%rcx,%rsi,8),%rdi
  add    %r15,%rdi
  mov    -0x58(%rsp),%esi
  add    %r12d,%esi
  movslq %esi,%rsi
  lea    (%rcx,%rsi,8),%rsi
  add    %r15,%rsi
  mov    %ebp,%r15d
  sub    -0x58(%rsp),%r15d
  movslq %r15d,%r15
  lea    (%rcx,%r15,8),%r15
  add    -0x60(%rsp),%r15
  movslq -0x44(%rsp),%rcx
  sub    -0x50(%rsp),%rcx
L19:
  movsd  0x18(%rax),%xmm0
  addsd  0x8(%rax),%xmm0
  addsd  (%r8,%rdx,1),%xmm0
  addsd  (%rdi,%rdx,1),%xmm0
  addsd  0x10(%rax,%rbx,8),%xmm0
  addsd  0x10(%rax,%r11,8),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm4,%xmm3
  mulsd  0x10(%rax),%xmm3
  subsd  0x0(%r13,%rdx,1),%xmm3
  addsd  %xmm3,%xmm0
  movsd  0x20(%rax),%xmm3
  addsd  (%rax),%xmm3
  addsd  (%rsi,%rdx,1),%xmm3
  addsd  (%r15,%rdx,1),%xmm3
  addsd  0x10(%rax,%r10,8),%xmm3
  addsd  0x10(%rax,%rcx,8),%xmm3
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%r14,%rdx,1)
  add    $0x8,%rax
  add    $0x8,%rdx
  cmp    %r9,%rax
  jne    L19
L70:
  addl   $0x1,-0x54(%rsp)
  mov    -0x54(%rsp),%eax
  mov    -0x58(%rsp),%edi
  add    %edi,%r12d
  add    %edi,%ebp
  cmp    -0x48(%rsp),%eax
  je     L82
L33:
  cmpl   $0x4,-0x58(%rsp)
  jg     L55
  jmp    L70
L82:
  addl   $0x1,-0x30(%rsp)
  mov    -0x30(%rsp),%ebx
  mov    -0x18(%rsp),%eax
  add    %eax,-0x2c(%rsp)
  add    %eax,-0x28(%rsp)
  add    %eax,-0x24(%rsp)
  add    %eax,-0x20(%rsp)
  add    %eax,-0x1c(%rsp)
  cmp    -0x14(%rsp),%ebx
  je     L16
L64:
  cmpl   $0x4,-0x10(%rsp)
  jle    L82
  mov    -0x1c(%rsp),%eax
  mov    %eax,-0x44(%rsp)
  mov    -0x58(%rsp),%ebp
  mov    -0xc(%rsp),%r12d
  movl   $0x2,-0x54(%rsp)
  movslq -0x2c(%rsp),%rax
  mov    %rax,-0x50(%rsp)
  mov    -0x8(%rsp),%rbx
  add    %rax,%rbx
  mov    %rbx,-0x40(%rsp)
  movslq -0x28(%rsp),%rbx
  sub    %rax,%rbx
  movslq -0x24(%rsp),%r11
  sub    %rax,%r11
  movslq -0x20(%rsp),%r10
  sub    %rax,%r10
  jmp    L33
L16:
  mov    $0x0,%eax
  pop    %rbx
  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;
}
```