user: 请从给定的汇编指令推导出C源代码：

```asm
<bwpvtipseudop>:
  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    L4
  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)
L1:
  cmpl   $0x4,0x30(%rsp)
  mov    $0x14,%ebx
  mov    $0x5,%edi
  jle    L20
  nopl   (%rax)
L18:
  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),%xmm8
  mov    0xa8(%rsp),%rax
  mov    %edi,0xc(%rsp)
  mov    -0x28(%rax,%r15,1),%rax
  mulss  %xmm8,%xmm8
  movss  -0x14(%rax,%rbx,1),%xmm9
  mov    0xb0(%rsp),%rax
  mov    -0x28(%rax,%r15,1),%rax
  mulss  %xmm9,%xmm9
  movaps %xmm8,%xmm1
  movss  %xmm8,0x8(%rsp)
  movss  -0x14(%rax,%rbx,1),%xmm0
  mov    0xb8(%rsp),%rax
  addss  %xmm0,%xmm0
  addss  D2(%rip),%xmm0
  mov    -0x28(%rax,%r15,1),%rax
  subss  %xmm9,%xmm1
  movss  %xmm9,0x4(%rsp)
  mulss  %xmm8,%xmm0
  movss  %xmm0,(%rsp)
  movss  -0x14(%rax,%rbx,1),%xmm0
  addss  %xmm0,%xmm0
  addss  D2(%rip),%xmm0
  mulss  %xmm8,%xmm0
  subss  %xmm9,%xmm0
  mulss  %xmm1,%xmm0
  call   <sqrt@plt>
  mov    (%r12),%r8
  mov    0x0(%rbp),%rsi
  xor    %eax,%eax
  pxor   %xmm10,%xmm10
  mov    0xc(%rsp),%edi
  mov    0x10(%rsp),%r10
  movaps %xmm0,%xmm7
  mov    0x18(%rsp),%rdx
  mov    0x20(%rsp),%r11
  movaps %xmm10,%xmm4
  movaps %xmm10,%xmm2
  mov    0x28(%rsp),%r9
  movaps %xmm10,%xmm3
  movaps %xmm10,%xmm5
  add    %rbx,%r8
  movss  0x8(%rsp),%xmm8
  movss  0x4(%rsp),%xmm9
  add    %rbx,%rsi
  nopw   0x0(%rax,%rax,1)
L8:
  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,%xmm5
  movss  (%r14,%rax,1),%xmm0
  addss  %xmm1,%xmm2
  mulss  %xmm0,%xmm6
  mulss  -0x14(%rsi,%rax,1),%xmm0
  add    $0x4,%rax
  addss  %xmm6,%xmm3
  addss  %xmm0,%xmm4
  cmp    $0x2c,%rax
  jne    L8
  mov    (%r10,%r15,1),%rax
  mulss  %xmm4,%xmm8
  mov    (%rdx,%r15,1),%rcx
  add    $0x1,%edi
  movss  (%r8),%xmm0
  addss  %xmm0,%xmm0
  movaps %xmm0,%xmm1
  movss  (%rsp),%xmm0
  subss  (%rax,%rbx,1),%xmm1
  mov    (%r11,%r15,1),%rax
  mulss  %xmm5,%xmm0
  movaps %xmm9,%xmm5
  mulss  %xmm3,%xmm5
  mulss  %xmm2,%xmm9
  addss  %xmm5,%xmm0
  movaps %xmm7,%xmm5
  mulss  %xmm2,%xmm5
  mulss  %xmm3,%xmm7
  addss  %xmm5,%xmm0
  movss  0x34(%rsp),%xmm5
  addss  %xmm9,%xmm7
  mulss  %xmm5,%xmm0
  addss  %xmm8,%xmm7
  mulss  %xmm5,%xmm7
  addss  %xmm1,%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    L18
L20:
  addl   $0x1,0x38(%rsp)
  add    $0x8,%r12
  mov    0x38(%rsp),%eax
  add    $0x8,%rbp
  add    $0x8,%r15
  cmp    %eax,0x3c(%rsp)
  jge    L1
L4:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: label: D2
type_infered: byte[4]

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

assistant: ```c
void bwpvtipseudop(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 * qx);
      q3[i][j] = 2 * q2[i][j] - q1[i][j] + dt2 * (coef * pz + vs2 * qx + vp2 * qz);
    }
  }
}
```