user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<step_forward_k>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,%r12
  mov    %rdi,0x30(%rsp)
  mov    %rsi,%rbx
  mov    %rdx,0x28(%rsp)
  mov    %ecx,%r15d
  mov    %r8d,%r13d
  mov    %r8d,0x74(%rsp)
  movss  %xmm0,0x38(%rsp)
  movss  %xmm1,0x3c(%rsp)
  movss  %xmm3,0x40(%rsp)
  movss  %xmm4,0x44(%rsp)
  movss  %xmm5,0x48(%rsp)
  movss  %xmm6,0x4c(%rsp)
  lea    0x1(%r9),%r14d
  mov    %r14d,0x70(%rsp)
  mov    %r8d,%ebp
  imul   %ecx,%ebp
  shl    $0x2,%ebp
  mov    %ebp,%edi
  call   <malloc@plt>
  mov    %rax,0x68(%rsp)
  mov    %ebp,%edx
  mov    %r12,%rsi
  mov    %rax,%rbp
  mov    %rax,%rdi
  call   <memcpy@plt>
  mov    %r13d,%eax
  sub    %r14d,%eax
  cmp    %eax,%r14d
  jge    L36
  mov    %r14d,%ecx
  mov    %r14d,%eax
  imul   %r15d,%eax
  mov    %eax,0x24(%rsp)
  movslq %r15d,%rax
  mov    %rax,0x78(%rsp)
  shl    $0x2,%rax
  mov    %rax,0x50(%rsp)
  sub    %rbp,%rax
  mov    0x28(%rsp),%rdx
  lea    (%rdx,%rax,1),%r14
  add    %rbx,%rax
  mov    %rax,0x58(%rsp)
  mov    %ecx,0x20(%rsp)
  mov    %rbp,%rax
  neg    %rax
  mov    %rax,0x10(%rsp)
  mov    %rbx,%rax
  sub    %rbp,%rax
  mov    %rax,0x60(%rsp)
  jmp    L16
L64:
  mov    0x18(%rsp),%r13
L25:
  mov    0xc(%rsp),%ebp
  add    0x24(%rsp),%ebp
  movss  (%r14,%r12,1),%xmm0
  pxor   %xmm7,%xmm7
  ucomiss %xmm0,%xmm7
  ja     L75
  sqrtss %xmm0,%xmm0
  movss  %xmm0,0x20(%rsp)
L50:
  movss  0x4c(%rsp),%xmm1
  movss  0x48(%rsp),%xmm0
  call   <fmaxf@plt>
  movaps %xmm0,%xmm1
  movss  0x20(%rsp),%xmm0
  mulss  0x40(%rsp),%xmm0
  mulss  0x44(%rsp),%xmm1
  divss  %xmm1,%xmm0
  call   <floor@plt>
  lea    0x1(%rax),%esi
  mov    %esi,%edx
  imul   %esi,%edx
  pxor   %xmm0,%xmm0
  cvtsi2ss %edx,%xmm0
  movss  0x38(%rsp),%xmm3
  divss  %xmm0,%xmm3
  movss  0x3c(%rsp),%xmm5
  divss  %xmm0,%xmm5
  mov    %r15d,%edi
  imul   %esi,%edi
  add    $0x2,%eax
  mov    %eax,0x20(%rsp)
  lea    (%rsi,%rbp,1),%ecx
  mov    %ebp,%edx
  sub    %esi,%edx
  lea    (%rdi,%rbp,1),%eax
  sub    %edi,%ebp
  lea    0x4(%r13),%rsi
  mov    %rsi,0x18(%rsp)
  mov    0x58(%rsp),%rdi
  movss  (%rdi,%r12,1),%xmm2
  mulss  D25(%rip),%xmm2
  addss  0x4(%rdi,%r12,1),%xmm2
  addss  (%rdi,%r13,1),%xmm2
  addss  -0x4(%rdi,%r12,1),%xmm2
  mov    0x60(%rsp),%rsi
  addss  (%rsi,%r12,1),%xmm2
  movss  D33(%rip),%xmm4
  mulss  %xmm4,%xmm2
  movslq %eax,%rdi
  movslq %ebp,%rsi
  movslq %ecx,%r9
  movslq %edx,%r8
  mov    0x30(%rsp),%r11
  add    0x10(%rsp),%r11
  movss  (%rbx,%r9,4),%xmm0
  mulss  D25(%rip),%xmm0
  addss  0x4(%rbx,%r9,4),%xmm0
  lea    (%r15,%rcx,1),%r10d
  movslq %r10d,%r10
  addss  (%rbx,%r10,4),%xmm0
  addss  -0x4(%rbx,%r9,4),%xmm0
  sub    %r15d,%ecx
  movslq %ecx,%rcx
  addss  (%rbx,%rcx,4),%xmm0
  mulss  %xmm4,%xmm0
  movss  (%rbx,%r8,4),%xmm1
  mulss  D25(%rip),%xmm1
  addss  0x4(%rbx,%r8,4),%xmm1
  lea    (%r15,%rdx,1),%ecx
  movslq %ecx,%rcx
  addss  (%rbx,%rcx,4),%xmm1
  addss  -0x4(%rbx,%r8,4),%xmm1
  sub    %r15d,%edx
  movslq %edx,%rdx
  addss  (%rbx,%rdx,4),%xmm1
  mulss  %xmm4,%xmm1
  addss  %xmm1,%xmm0
  mulss  %xmm3,%xmm0
  movss  (%rbx,%rdi,4),%xmm1
  mulss  D25(%rip),%xmm1
  addss  0x4(%rbx,%rdi,4),%xmm1
  lea    (%r15,%rax,1),%edx
  movslq %edx,%rdx
  addss  (%rbx,%rdx,4),%xmm1
  addss  -0x4(%rbx,%rdi,4),%xmm1
  sub    %r15d,%eax
  cltq
  addss  (%rbx,%rax,4),%xmm1
  mulss  %xmm4,%xmm1
  movss  (%rbx,%rsi,4),%xmm6
  mulss  D25(%rip),%xmm6
  addss  0x4(%rbx,%rsi,4),%xmm6
  lea    (%r15,%rbp,1),%eax
  cltq
  addss  (%rbx,%rax,4),%xmm6
  addss  -0x4(%rbx,%rsi,4),%xmm6
  sub    %r15d,%ebp
  movslq %ebp,%rbp
  addss  (%rbx,%rbp,4),%xmm6
  mulss  %xmm4,%xmm6
  addss  %xmm6,%xmm1
  mulss  %xmm5,%xmm1
  addss  %xmm1,%xmm0
  addss  %xmm5,%xmm3
  addss  %xmm3,%xmm3
  mulss  %xmm2,%xmm3
  subss  %xmm3,%xmm0
  movss  (%r14,%r12,1),%xmm1
  mulss  D25(%rip),%xmm1
  addss  0x4(%r14,%r12,1),%xmm1
  addss  (%r14,%r13,1),%xmm1
  addss  -0x4(%r14,%r12,1),%xmm1
  mov    0x28(%rsp),%rax
  add    0x10(%rsp),%rax
  addss  (%rax,%r12,1),%xmm1
  mulss  %xmm4,%xmm1
  mulss  %xmm1,%xmm0
  addss  %xmm2,%xmm2
  movss  0x0(%r13),%xmm1
  mulss  D25(%rip),%xmm1
  addss  0x4(%r13),%xmm1
  mov    0x50(%rsp),%rcx
  addss  0x0(%r13,%rcx,1),%xmm1
  addss  -0x4(%r13),%xmm1
  addss  (%r12),%xmm1
  mulss  %xmm4,%xmm1
  subss  %xmm1,%xmm2
  addss  %xmm2,%xmm0
  movss  %xmm0,(%r11,%r13,1)
  addl   $0x1,0xc(%rsp)
  mov    0xc(%rsp),%esi
  add    $0x4,%r12
  mov    %r15d,%eax
  sub    0x20(%rsp),%eax
  cmp    %esi,%eax
  jg     L64
L11:
  addl   $0x1,0x70(%rsp)
  mov    0x70(%rsp),%ecx
  add    %r15d,0x24(%rsp)
  mov    0x74(%rsp),%eax
  sub    0x20(%rsp),%eax
  cmp    %ecx,%eax
  jle    L36
L16:
  mov    %r15d,%eax
  mov    0x20(%rsp),%ecx
  sub    %ecx,%eax
  cmp    %ecx,%eax
  jle    L11
  movslq 0x24(%rsp),%rdx
  mov    0x20(%rsp),%esi
  movslq %esi,%rax
  lea    (%rax,%rdx,1),%rcx
  mov    0x68(%rsp),%rdi
  lea    (%rdi,%rcx,4),%r13
  sub    0x78(%rsp),%rax
  add    %rdx,%rax
  lea    (%rdi,%rax,4),%r12
  mov    %esi,0xc(%rsp)
  jmp    L25
L75:
  call   <sqrtf@plt>
  movss  %xmm0,0x20(%rsp)
  jmp    L50
L36:
  mov    0x68(%rsp),%rdi
  call   <free@plt>
  mov    $0x0,%eax
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D25
data_type: byte[4]
name: extract_symbolic_data
arguments: label: D33
data_type: byte[4]

tool: {"D25": [0x00,0x00,0x80,0x40]}

tool: {"D33": [0x00,0x00,0x00,0x3e]}

assistant: ```c
int step_forward_k(float *tgt_p, float *src_p, float *v, int nz, int nx, float rz, float rx, float s, int k, float dt, float CFL, float dx, float dz) {
  int iz;
  int ix;
  int ioff;
  float two = 2.0;
  float k2 = k * k, k_rz = rz / k2, k_rx = rx / k2, k_s = 2.f * (k_rz + k_rx);
  int k_nz = k * nz, k1_nz = (k - 1) * nz, k1 = k - 1, kp1 = k + 1;
  float src_iz_plus_1, src_iz_minus_1, src_ix_iz, src_ix_plus_1, src_ix_minus_1, tgt_ix_iz, vel;
  int ioffkzp, ioffkzm, ioffkxp, ioffkxm, num_cfl = 0;
  float *tgt = (float *)malloc(nx * nz * sizeof(float));
  memcpy(tgt, tgt_p, nx * nz * sizeof(float));
  for (ix = kp1; ix < nx - kp1; ix++) {
    for (iz = kp1; iz < nz - kp1; iz++) {
      ioff = iz + ix * nz;
      vel = sqrtf(v[ioff]);
      k = floor(vel * dt / (CFL * fmaxf(dx, dz))) + 1;
      k2 = k * k;
      k_rz = rz / k2;
      k_rx = rx / k2;
      k_s = 2.f * (k_rz + k_rx);
      k_nz = k * nz;
      k1_nz = (k - 1) * nz;
      k1 = k - 1;
      kp1 = k + 1;
      ioffkzp = ioff + k;
      ioffkzm = ioff - k;
      ioffkxp = ioff + k_nz;
      ioffkxm = ioff - k_nz;
      tgt_ix_iz = 0.125f * (4.f * tgt[ioff] + tgt[ioff + 1] + tgt[ioff + nz] + tgt[ioff - 1] + tgt[ioff - nz]);
      src_ix_iz = 0.125f * (4.f * src_p[ioff] + src_p[ioff + 1] + src_p[ioff + nz] + src_p[ioff - 1] + src_p[ioff - nz]);
      src_ix_plus_1 = 0.125f * (4.f * src_p[ioffkxp] + src_p[ioffkxp + 1] + src_p[ioffkxp + nz] + src_p[ioffkxp - 1] + src_p[ioffkxp - nz]);
      src_ix_minus_1 = 0.125f * (4.f * src_p[ioffkxm] + src_p[ioffkxm + 1] + src_p[ioffkxm + nz] + src_p[ioffkxm - 1] + src_p[ioffkxm - nz]);
      src_iz_plus_1 = 0.125f * (4.f * src_p[ioffkzp] + src_p[ioffkzp + 1] + src_p[ioffkzp + nz] + src_p[ioffkzp - 1] + src_p[ioffkzp - nz]);
      src_iz_minus_1 = 0.125f * (4.f * src_p[ioffkzm] + src_p[ioffkzm + 1] + src_p[ioffkzm + nz] + src_p[ioffkzm - 1] + src_p[ioffkzm - nz]);
      vel = 0.125f * (4.f * v[ioff] + v[ioff + 1] + v[ioff + nz] + v[ioff - 1] + v[ioff - nz]);
      tgt_p[ioff] = two * src_ix_iz - tgt_ix_iz + vel * (k_rz * (src_iz_plus_1 + src_iz_minus_1) + k_rx * (src_ix_plus_1 + src_ix_minus_1) - k_s * src_ix_iz);
    }
  }
  free(tgt);
  return num_cfl;
}
```