user: 这段汇编代码用C语言表示是什么样的？

```asm
<C_advect_velocity>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %rdi,0x8(%rsp)
  mov    %rdx,0x30(%rsp)
  mov    %rcx,0x68(%rsp)
  mov    %r9d,0x54(%rsp)
  movsd  %xmm0,0x38(%rsp)
  movsd  %xmm1,0x40(%rsp)
  movsd  %xmm2,0x48(%rsp)
  mov    %r9d,%eax
  imul   0xb0(%rsp),%eax
  mov    %eax,0x64(%rsp)
  test   %r9d,%r9d
  jle    L28
  mov    %rsi,%r15
  mov    %r8,%r14
  movl   $0x0,0x60(%rsp)
  movl   $0x0,0x50(%rsp)
  jmp    L80
L58:
  pxor   %xmm1,%xmm1
L101:
  pxor   %xmm7,%xmm7
  comiss %xmm0,%xmm7
  ja     L32
  pxor   %xmm2,%xmm2
  cvtsi2sdl 0xb0(%rsp),%xmm2
  subsd  D18(%rip),%xmm2
  pxor   %xmm3,%xmm3
  cvtss2sd %xmm0,%xmm3
  comisd %xmm2,%xmm3
  jb     L85
  pxor   %xmm0,%xmm0
  cvtsd2ss %xmm2,%xmm0
  jmp    L85
L32:
  pxor   %xmm0,%xmm0
L85:
  cvttss2si %xmm1,%ebp
  cvttss2si %xmm0,%eax
  cvtss2sd %xmm1,%xmm1
  pxor   %xmm2,%xmm2
  cvtsi2sd %ebp,%xmm2
  subsd  %xmm2,%xmm1
  pxor   %xmm5,%xmm5
  cvtsd2ss %xmm1,%xmm5
  movss  %xmm5,(%rsp)
  cvtss2sd %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm6,%xmm6
  cvtsd2ss %xmm0,%xmm6
  movss  %xmm6,0x4(%rsp)
  mov    0x20(%rsp),%rcx
  mov    %ebp,(%rcx,%rbx,4)
  mov    0x28(%rsp),%rsi
  mov    %eax,(%rsi,%rbx,4)
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm5,%xmm0
  movsd  %xmm0,(%r14,%r13,1)
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm6,%xmm0
  movsd  %xmm0,(%r14,%r12,1)
  mov    0x30(%rsp),%rdx
  add    0x10(%rsp),%rdx
  cmpb   $0x0,(%rdx,%rbx,1)
  je     L44
  movsd  (%r15,%r13,1),%xmm0
  mov    0x8(%rsp),%rax
  movsd  %xmm0,(%rax,%r13,1)
  movsd  (%r15,%r12,1),%xmm0
  movsd  %xmm0,(%rax,%r12,1)
L106:
  add    $0x1,%rbx
  add    $0x8,%r13
  add    $0x8,%r12
  cmp    0x18(%rsp),%rbx
  je     L57
L1:
  pxor   %xmm1,%xmm1
  cvtsi2sdl 0x50(%rsp),%xmm1
  movsd  0x48(%rsp),%xmm4
  movapd %xmm4,%xmm0
  mulsd  (%r15,%r13,1),%xmm0
  divsd  0x38(%rsp),%xmm0
  subsd  %xmm0,%xmm1
  cvtsd2ss %xmm1,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  mulsd  (%r15,%r12,1),%xmm4
  movapd %xmm4,%xmm2
  divsd  0x40(%rsp),%xmm2
  subsd  %xmm2,%xmm0
  cvtsd2ss %xmm0,%xmm0
  pxor   %xmm3,%xmm3
  comiss %xmm1,%xmm3
  ja     L58
  pxor   %xmm2,%xmm2
  cvtsi2sdl 0x54(%rsp),%xmm2
  subsd  D18(%rip),%xmm2
  pxor   %xmm3,%xmm3
  cvtss2sd %xmm1,%xmm3
  comisd %xmm2,%xmm3
  jb     L101
  pxor   %xmm1,%xmm1
  cvtsd2ss %xmm2,%xmm1
  jmp    L101
L44:
  imul   0xb0(%rsp),%ebp
  add    %eax,%ebp
  movslq %ebp,%rbp
  lea    (%r15,%rbp,8),%rdi
  movaps %xmm6,%xmm1
  movaps %xmm5,%xmm0
  mov    0xb0(%rsp),%esi
  call   <advect_sample@plt>
  mov    0x8(%rsp),%rax
  movsd  %xmm0,(%rax,%r13,1)
  add    0x58(%rsp),%rbp
  lea    (%r15,%rbp,8),%rdi
  movss  0x4(%rsp),%xmm1
  movss  (%rsp),%xmm0
  mov    0xb0(%rsp),%esi
  call   <advect_sample@plt>
  mov    0x8(%rsp),%rax
  movsd  %xmm0,(%rax,%r12,1)
  jmp    L106
L57:
  addl   $0x1,0x50(%rsp)
  mov    0x50(%rsp),%eax
  mov    0xb0(%rsp),%edi
  add    %edi,0x60(%rsp)
  cmp    %eax,0x54(%rsp)
  je     L28
L80:
  cmpl   $0x0,0xb0(%rsp)
  jle    L57
  movslq 0x60(%rsp),%rax
  mov    %rax,%rdi
  mov    %rax,0x10(%rsp)
  lea    0x0(,%rax,8),%r13
  movslq 0x64(%rsp),%rcx
  mov    %rcx,0x58(%rsp)
  add    %rcx,%rax
  lea    0x0(,%rax,8),%r12
  mov    0xb0(%rsp),%ecx
  mov    %rcx,0x18(%rsp)
  mov    0x68(%rsp),%rsi
  lea    (%rsi,%rdi,4),%rdi
  mov    %rdi,0x20(%rsp)
  lea    (%rsi,%rax,4),%rax
  mov    %rax,0x28(%rsp)
  mov    $0x0,%ebx
  jmp    L1
L28:
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: disassemble_data
arguments: label: D18
datatype: f64

tool: D18: 1.01

assistant: ```c
void C_advect_velocity(double *v, const double *const v0, const unsigned char *bound, int *advect_indexes, double *advect_lerps, const int Nx, const int Ny, const double dx, const double dy, const double dt) {
  int x, y, idx, xi, yi, iidx;
  int vyidx = Nx * Ny;
  float xa, ya, s, t;
#pragma omp for schedule(auto)
  for (x = 0; x < Nx; ++x) {
    for (y = 0; y < Ny; ++y) {
      idx = y + x * Ny;
      xa = (double)x - dt * v0[0 + idx] / dx;
      ya = (double)y - dt * v0[vyidx + idx] / dy;
      xa = xa < 0.0 ? 0.0 : (xa >= Nx - 1.01) ? (Nx - 1.01) : xa;
      ya = ya < 0.0 ? 0.0 : (ya >= Ny - 1.01) ? (Ny - 1.01) : ya;
      xi = (int)xa;
      yi = (int)ya;
      s = xa - (double)xi;
      t = ya - (double)yi;
      advect_indexes[0 + idx] = xi;
      advect_indexes[vyidx + idx] = yi;
      advect_lerps[0 + idx] = s;
      advect_lerps[vyidx + idx] = t;
      if (!bound[idx]) {
        iidx = yi + xi * Ny;
        v[0 + idx] = advect_sample(v0 + iidx, Ny, s, t);
        v[vyidx + idx] = advect_sample(v0 + vyidx + iidx, Ny, s, t);
      } else {
        v[0 + idx] = v0[0 + idx];
        v[vyidx + idx] = v0[vyidx + idx];
      }
    }
  }
}
```