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

```asm
<fwpvtipseudosv>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  movss  %xmm0,-0x54(%rbp)
  mov    %rdi,-0x60(%rbp)
  mov    %rsi,-0x68(%rbp)
  mov    %rdx,-0x70(%rbp)
  mov    %rcx,-0x78(%rbp)
  mov    %r8,-0x80(%rbp)
  mov    %r9,-0x88(%rbp)
  movl   $0x5,-0x44(%rbp)
  jmp    L5
L32:
  mov    -0x44(%rbp),%eax
  sub    $0x5,%eax
  mov    %eax,-0x28(%rbp)
  movl   $0x5,-0x40(%rbp)
  jmp    L12
L42:
  mov    -0x40(%rbp),%eax
  sub    $0x5,%eax
  mov    %eax,-0x24(%rbp)
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x20(%rbp)
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x1c(%rbp)
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x40(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  D17(%rip),%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x18(%rbp)
  mov    -0x28(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  D17(%rip),%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x14(%rbp)
  movss  -0x20(%rbp),%xmm0
  mulss  -0x18(%rbp),%xmm0
  movss  %xmm0,-0x10(%rbp)
  movss  -0x20(%rbp),%xmm0
  mulss  -0x14(%rbp),%xmm0
  movss  %xmm0,-0xc(%rbp)
  movss  -0xc(%rbp),%xmm0
  movaps %xmm0,%xmm1
  subss  -0x1c(%rbp),%xmm1
  movss  -0x20(%rbp),%xmm0
  subss  -0x1c(%rbp),%xmm0
  mulss  %xmm0,%xmm1
  movd   %xmm1,%eax
  movd   %eax,%xmm0
  call   <sqrt@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x8(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x34(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x2c(%rbp)
  movl   $0xfffffffb,-0x3c(%rbp)
  jmp    L1
L50:
  mov    -0x3c(%rbp),%eax
  add    $0x5,%eax
  mov    %eax,-0x4(%rbp)
  mov    -0x4(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x44(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  -0x38(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x38(%rbp)
  mov    -0x4(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%ecx
  mov    -0x3c(%rbp),%eax
  add    %ecx,%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  -0x34(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x34(%rbp)
  mov    -0x4(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x44(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  -0x30(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x30(%rbp)
  mov    -0x4(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%ecx
  mov    -0x3c(%rbp),%eax
  add    %ecx,%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  -0x2c(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x2c(%rbp)
  addl   $0x1,-0x3c(%rbp)
L1:
  cmpl   $0x5,-0x3c(%rbp)
  jle    L50
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  addss  %xmm0,%xmm0
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm2
  movss  -0x10(%rbp),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x38(%rbp),%xmm1
  movss  -0x1c(%rbp),%xmm0
  mulss  -0x34(%rbp),%xmm0
  addss  %xmm1,%xmm0
  movss  -0x8(%rbp),%xmm1
  mulss  -0x2c(%rbp),%xmm1
  subss  %xmm1,%xmm0
  mulss  -0x54(%rbp),%xmm0
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  addss  %xmm2,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  addss  %xmm0,%xmm0
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm2
  movaps %xmm0,%xmm1
  subss  %xmm2,%xmm1
  movss  -0x8(%rbp),%xmm0
  movss  D9(%rip),%xmm2
  xorps  %xmm2,%xmm0
  movaps %xmm0,%xmm2
  mulss  -0x38(%rbp),%xmm2
  movss  -0x1c(%rbp),%xmm0
  mulss  -0x30(%rbp),%xmm0
  addss  %xmm0,%xmm2
  movss  -0x20(%rbp),%xmm0
  mulss  -0x2c(%rbp),%xmm0
  addss  %xmm2,%xmm0
  mulss  -0x54(%rbp),%xmm0
  mov    -0x44(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x40(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x40(%rbp)
L12:
  mov    0x28(%rbp),%eax
  add    $0x4,%eax
  cmp    %eax,-0x40(%rbp)
  jle    L42
  addl   $0x1,-0x44(%rbp)
L5:
  mov    0x20(%rbp),%eax
  add    $0x4,%eax
  cmp    %eax,-0x44(%rbp)
  jle    L32
  nop
  nop
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D17
type_infered: byte[4]
name: extract_symbolic_data
arguments: tag: D9
type_infered: byte[4]

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

tool: {"D9": [0x00,0x00,0x00,0x80]}

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);
    }
  }
}
```