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

```asm
<fwpvtipseudosv>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  movss  %xmm0,0x24(%rsp)
  mov    %rdi,0x28(%rsp)
  mov    %rdx,0x30(%rsp)
  mov    %rcx,0x38(%rsp)
  mov    %r9,0x40(%rsp)
  mov    0x90(%rsp),%rbp
  mov    0x98(%rsp),%r12
  mov    0xa0(%rsp),%ecx
  add    $0x4,%ecx
  mov    %ecx,0x4c(%rsp)
  cmp    $0x4,%ecx
  jle    L23
  mov    %rsi,%rax
  mov    %r8,%r15
  add    $0x28,%r15
  add    $0x28,%rax
  mov    %rax,0x18(%rsp)
  mov    $0x28,%r13d
  movl   $0x5,0x48(%rsp)
  mov    0xa8(%rsp),%eax
  add    $0x4,%eax
  mov    %eax,0x20(%rsp)
L14:
  cmpl   $0x4,0x20(%rsp)
  jle    L17
  mov    $0x14,%ebx
  mov    $0x5,%r14d
L29:
  mov    0xb0(%rsp),%rax
  mov    -0x28(%rax,%r13,1),%rax
  movss  -0x14(%rax,%rbx,1),%xmm2
  movaps %xmm2,%xmm7
  mulss  %xmm2,%xmm7
  movss  %xmm7,0xc(%rsp)
  mov    0xb8(%rsp),%rax
  mov    -0x28(%rax,%r13,1),%rax
  movss  -0x14(%rax,%rbx,1),%xmm0
  movaps %xmm0,%xmm5
  mulss  %xmm0,%xmm5
  movss  %xmm5,0x10(%rsp)
  mov    0xc0(%rsp),%rax
  mov    -0x28(%rax,%r13,1),%rax
  movss  -0x14(%rax,%rbx,1),%xmm0
  addss  %xmm0,%xmm0
  addss  D0(%rip),%xmm0
  mulss  %xmm7,%xmm0
  movss  %xmm0,0x14(%rsp)
  mov    0xc8(%rsp),%rax
  mov    -0x28(%rax,%r13,1),%rax
  movss  -0x14(%rax,%rbx,1),%xmm0
  addss  %xmm0,%xmm0
  addss  D0(%rip),%xmm0
  mulss  %xmm7,%xmm0
  subss  %xmm5,%xmm0
  subss  %xmm5,%xmm7
  mulss  %xmm7,%xmm0
  call   <sqrt@plt>
  movaps %xmm0,%xmm7
  mov    0x18(%rsp),%rdi
  mov    %rdi,%r8
  mov    %r15,%rcx
  mov    $0x0,%eax
  pxor   %xmm2,%xmm2
  movaps %xmm2,%xmm5
  movaps %xmm2,%xmm3
  movaps %xmm2,%xmm4
  mov    %rbx,%rsi
  add    (%rdi),%rsi
  mov    %rbx,%rdi
  add    (%r15),%rdi
L40:
  movss  0x0(%rbp,%rax,1),%xmm1
  mov    -0x28(%r8,%rax,2),%rdx
  movaps %xmm1,%xmm0
  mulss  (%rdx,%rbx,1),%xmm0
  addss  %xmm0,%xmm4
  movss  (%r12,%rax,1),%xmm0
  movaps %xmm0,%xmm6
  mulss  -0x14(%rsi,%rax,1),%xmm6
  addss  %xmm6,%xmm3
  mov    -0x28(%rcx,%rax,2),%rdx
  mulss  (%rdx,%rbx,1),%xmm1
  addss  %xmm1,%xmm5
  mulss  -0x14(%rdi,%rax,1),%xmm0
  addss  %xmm0,%xmm2
  add    $0x4,%rax
  cmp    $0x2c,%rax
  jne    L40
  mov    0x30(%rsp),%rax
  mov    (%rax,%r13,1),%rax
  movss  (%rsi),%xmm1
  addss  %xmm1,%xmm1
  mov    0x28(%rsp),%rsi
  mov    (%rsi,%r13,1),%rdx
  subss  (%rdx,%rbx,1),%xmm1
  movss  0x14(%rsp),%xmm0
  mulss  %xmm4,%xmm0
  movss  0x10(%rsp),%xmm6
  mulss  %xmm6,%xmm3
  addss  %xmm3,%xmm0
  movaps %xmm7,%xmm3
  mulss  %xmm2,%xmm3
  subss  %xmm3,%xmm0
  movss  0x24(%rsp),%xmm3
  mulss  %xmm3,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax,%rbx,1)
  mov    0x40(%rsp),%rax
  mov    (%rax,%r13,1),%rax
  mov    (%rcx),%rdx
  movss  (%rdx,%rbx,1),%xmm1
  addss  %xmm1,%xmm1
  mov    0x38(%rsp),%rcx
  mov    (%rcx,%r13,1),%rdx
  subss  (%rdx,%rbx,1),%xmm1
  mulss  %xmm4,%xmm7
  mulss  %xmm5,%xmm6
  movaps %xmm6,%xmm0
  subss  %xmm7,%xmm0
  mulss  0xc(%rsp),%xmm2
  addss  %xmm2,%xmm0
  mulss  %xmm3,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax,%rbx,1)
  add    $0x1,%r14d
  add    $0x4,%rbx
  cmp    0x20(%rsp),%r14d
  jle    L29
L17:
  addl   $0x1,0x48(%rsp)
  mov    0x48(%rsp),%eax
  add    $0x8,%r15
  addq   $0x8,0x18(%rsp)
  add    $0x8,%r13
  cmp    0x4c(%rsp),%eax
  jle    L14
L23:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D0
type: byte[4]

tool: {"D0": [0x00,0x00,0x80,0x3f]}

assistant: ```c
void fwpvtipseudosv(float dt2, float **p1, float **p2, float **p3, float **q1, float **q2, float **q3, float *coeff_x, float *coeff_z, int nx, int nz, float **vp0, float **vs0, float **epsilon, float **delta) {
  int i, j, k, im, jm, km;
  float px, pz, qx, qz, vp2, vs2, vpx2, vpn2, ep, de, coef;
  for (i = 5; i < nx + 5; i++) {
    im = i - 5;
    for (j = 5; j < nz + 5; j++) {
      jm = j - 5;
      vp2 = vp0[im][jm] * vp0[im][jm];
      vs2 = vs0[im][jm] * vs0[im][jm];
      ep = 1 + 2 * epsilon[im][jm];
      de = 1 + 2 * delta[im][jm];
      vpx2 = vp2 * ep;
      vpn2 = vp2 * de;
      coef = sqrt((vpn2 - vs2) * (vp2 - vs2));
      px = 0;
      qx = 0;
      pz = 0;
      qz = 0;
      for (k = -5; k <= 5; k++) {
        km = k + 5;
        px += coeff_x[km] * p2[i + k][j];
        pz += coeff_z[km] * p2[i][j + k];
        qx += coeff_x[km] * q2[i + k][j];
        qz += coeff_z[km] * q2[i][j + k];
      }
      p3[i][j] = 2 * p2[i][j] - p1[i][j] + dt2 * (vpx2 * px + vs2 * pz - coef * qz);
      q3[i][j] = 2 * q2[i][j] - q1[i][j] + dt2 * (-coef * px + vs2 * qx + vp2 * qz);
    }
  }
}
```