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

```asm
<fwpvtipseudosv>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    0x90(%rsp),%eax
  mov    0x80(%rsp),%r13
  movss  %xmm0,0x34(%rsp)
  mov    0x88(%rsp),%r14
  add    $0x4,%eax
  mov    %eax,0x3c(%rsp)
  cmp    $0x4,%eax
  jle    L25
  mov    0x98(%rsp),%eax
  mov    %rdi,%r15
  mov    %rsi,%r12
  mov    %r8,%rbp
  mov    $0x28,%edi
  mov    %r15,%r10
  mov    %rcx,%r11
  add    $0x28,%r12
  add    $0x4,%eax
  movl   $0x5,0x38(%rsp)
  add    $0x28,%rbp
  mov    %rdi,%r15
  mov    %eax,0x30(%rsp)
  nopl   0x0(%rax)
L45:
  cmpl   $0x4,0x30(%rsp)
  mov    $0x14,%ebx
  mov    $0x5,%edi
  jle    L21
  nopl   (%rax)
L40:
  mov    0xa0(%rsp),%rax
  mov    %r9,0x28(%rsp)
  mov    %r11,0x20(%rsp)
  mov    -0x28(%rax,%r15,1),%rax
  mov    %rdx,0x18(%rsp)
  mov    %r10,0x10(%rsp)
  movss  -0x14(%rax,%rbx,1),%xmm10
  mov    0xa8(%rsp),%rax
  mov    %edi,0xc(%rsp)
  mov    -0x28(%rax,%r15,1),%rax
  mulss  %xmm10,%xmm10
  movss  -0x14(%rax,%rbx,1),%xmm7
  mov    0xb0(%rsp),%rax
  mov    -0x28(%rax,%r15,1),%rax
  mulss  %xmm7,%xmm7
  movaps %xmm10,%xmm1
  movss  %xmm10,0x4(%rsp)
  movss  -0x14(%rax,%rbx,1),%xmm8
  mov    0xb8(%rsp),%rax
  mov    -0x28(%rax,%r15,1),%rax
  addss  %xmm8,%xmm8
  addss  D5(%rip),%xmm8
  subss  %xmm7,%xmm1
  movss  %xmm7,(%rsp)
  movss  -0x14(%rax,%rbx,1),%xmm0
  mulss  %xmm10,%xmm8
  addss  %xmm0,%xmm0
  addss  D5(%rip),%xmm0
  mulss  %xmm10,%xmm0
  movss  %xmm8,0x8(%rsp)
  subss  %xmm7,%xmm0
  mulss  %xmm1,%xmm0
  call   <sqrt@plt>
  mov    (%r12),%r8
  mov    0x0(%rbp),%rsi
  xor    %eax,%eax
  pxor   %xmm9,%xmm9
  mov    0xc(%rsp),%edi
  mov    0x10(%rsp),%r10
  movaps %xmm0,%xmm11
  movss  0x8(%rsp),%xmm8
  mov    0x18(%rsp),%rdx
  movaps %xmm9,%xmm3
  movaps %xmm9,%xmm5
  mov    0x20(%rsp),%r11
  mov    0x28(%rsp),%r9
  add    %rbx,%r8
  add    %rbx,%rsi
  movss  0x4(%rsp),%xmm10
  movss  (%rsp),%xmm7
  movaps %xmm9,%xmm2
  movaps %xmm9,%xmm4
  nop
L1:
  mov    -0x28(%r12,%rax,2),%rcx
  movss  0x0(%r13,%rax,1),%xmm1
  movss  -0x14(%r8,%rax,1),%xmm6
  movss  (%rcx,%rbx,1),%xmm0
  mov    -0x28(%rbp,%rax,2),%rcx
  mulss  %xmm1,%xmm0
  mulss  (%rcx,%rbx,1),%xmm1
  addss  %xmm0,%xmm4
  movss  (%r14,%rax,1),%xmm0
  addss  %xmm1,%xmm5
  mulss  %xmm0,%xmm6
  mulss  -0x14(%rsi,%rax,1),%xmm0
  add    $0x4,%rax
  addss  %xmm6,%xmm2
  addss  %xmm0,%xmm3
  cmp    $0x2c,%rax
  jne    L1
  mulss  %xmm7,%xmm2
  movaps %xmm11,%xmm1
  mov    (%r10,%r15,1),%rax
  mov    (%rdx,%r15,1),%rcx
  mulss  %xmm4,%xmm8
  movss  (%r8),%xmm0
  add    $0x1,%edi
  mulss  %xmm3,%xmm1
  mulss  %xmm4,%xmm11
  addss  %xmm0,%xmm0
  subss  (%rax,%rbx,1),%xmm0
  mov    (%r11,%r15,1),%rax
  mulss  %xmm5,%xmm7
  mulss  %xmm3,%xmm10
  addss  %xmm2,%xmm8
  movss  0x34(%rsp),%xmm2
  subss  %xmm1,%xmm8
  subss  %xmm11,%xmm7
  mulss  %xmm2,%xmm8
  addss  %xmm10,%xmm7
  mulss  %xmm2,%xmm7
  addss  %xmm8,%xmm0
  movss  %xmm0,(%rcx,%rbx,1)
  movss  (%rsi),%xmm0
  mov    (%r9,%r15,1),%rcx
  addss  %xmm0,%xmm0
  subss  (%rax,%rbx,1),%xmm0
  addss  %xmm7,%xmm0
  movss  %xmm0,(%rcx,%rbx,1)
  add    $0x4,%rbx
  cmp    0x30(%rsp),%edi
  jle    L40
L21:
  addl   $0x1,0x38(%rsp)
  add    $0x8,%r12
  mov    0x38(%rsp),%eax
  add    $0x8,%rbp
  add    $0x8,%r15
  cmp    %eax,0x3c(%rsp)
  jge    L45
L25:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: tag: D5
type_infered: byte[4]

tool: {"D5": [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);
    }
  }
}
```