user: What is the c source code of the assembly code below:

```asm
<BatchInterpolatedGradX3DMirrWarp>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  movsd  %xmm5,0x40(%rsp)
  pxor   %xmm5,%xmm5
  comisd %xmm5,%xmm6
  mov    %rdi,0x20(%rsp)
  mov    %r9,0x30(%rsp)
  movsd  %xmm0,0x10(%rsp)
  movsd  %xmm1,0x28(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm3,0x60(%rsp)
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm6,0x68(%rsp)
  movsd  %xmm7,0x48(%rsp)
  jbe    L0
  movq   $0x0,0x58(%rsp)
  xor    %r15d,%r15d
  comisd %xmm5,%xmm3
  jbe    L0
  mov    %rsi,%rbx
  mov    %rdx,%rbp
  mov    %rcx,%r12
  mov    %r8,%r13
  mov    %r15,%r14
  nop
L40:
  movq   $0x0,0x50(%rsp)
  mov    %r14,%r15
  nopl   0x0(%rax)
L33:
  movsd  0x10(%rsp),%xmm6
  pxor   %xmm7,%xmm7
  comisd %xmm7,%xmm6
  jbe    L13
  pxor   %xmm5,%xmm5
  pxor   %xmm6,%xmm6
  xor    %r14d,%r14d
  mov    %rbp,%rax
  cvtsi2ssq 0x50(%rsp),%xmm5
  mov    %r14,%rbp
  mov    %rbx,%r14
  mov    %r15,%rbx
  cvtsi2ssq 0x58(%rsp),%xmm6
  mov    %rax,%r15
  movss  %xmm5,0x18(%rsp)
  movss  %xmm6,0x1c(%rsp)
  nopl   0x0(%rax,%rax,1)
L57:
  movss  0x1c(%rsp),%xmm2
  pxor   %xmm0,%xmm0
  mov    %r13,%rcx
  mov    %r12,%rdx
  mov    0xc0(%rsp),%rax
  cvtsi2ss %rbp,%xmm0
  mov    %r15,%rsi
  movss  0x18(%rsp),%xmm1
  mov    0xc8(%rsp),%r8
  mov    %r14,%rdi
  add    $0x1,%rbp
  addss  (%rax,%rbx,4),%xmm2
  mov    0xb0(%rsp),%rax
  addss  (%rax,%rbx,4),%xmm1
  mov    0x30(%rsp),%rax
  cvtss2sd %xmm2,%xmm2
  subsd  0x48(%rsp),%xmm2
  divsd  0xb8(%rsp),%xmm2
  addss  (%rax,%rbx,4),%xmm0
  add    $0x1,%rbx
  cvtss2sd %xmm1,%xmm1
  subsd  0x38(%rsp),%xmm1
  divsd  0x40(%rsp),%xmm1
  cvtss2sd %xmm0,%xmm0
  subsd  0x28(%rsp),%xmm0
  divsd  0x8(%rsp),%xmm0
  call   <InterpolatedGradX3DMirr@plt>
  mov    0x20(%rsp),%rax
  divsd  0x8(%rsp),%xmm0
  cvtsd2ss %xmm0,%xmm0
  movsd  0x10(%rsp),%xmm4
  movss  %xmm0,-0x4(%rax,%rbx,4)
  pxor   %xmm0,%xmm0
  cvtsi2sd %rbp,%xmm0
  comisd %xmm0,%xmm4
  ja     L57
  mov    %r15,%rbp
  mov    %rbx,%r15
  mov    %r14,%rbx
L13:
  addq   $0x1,0x50(%rsp)
  mov    0x50(%rsp),%rax
  pxor   %xmm0,%xmm0
  movsd  0x60(%rsp),%xmm5
  cvtsi2sd %rax,%xmm0
  comisd %xmm0,%xmm5
  ja     L33
  addq   $0x1,0x58(%rsp)
  mov    0x58(%rsp),%rax
  pxor   %xmm0,%xmm0
  mov    %r15,%r14
  movsd  0x68(%rsp),%xmm7
  cvtsi2sd %rax,%xmm0
  comisd %xmm0,%xmm7
  ja     L40
L0:
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
extern void BatchInterpolatedGradX3DMirrWarp(float *Bout, float *Bcoeff, long Width, long Height, long Slice, double nx, double offx, double mx, float *WarpX, double ny, double offy, double my, float *WarpY, double nz, double offz, double mz, float *WarpZ, long SplineDegree) {
  double x, y, z;
  long idx;
  long idcoordx, idcoordy, idcoordz;
  idx = 0;
  for (idcoordz = 0; idcoordz < nz; idcoordz++) {
    for (idcoordy = 0; idcoordy < ny; idcoordy++) {
      for (idcoordx = 0; idcoordx < nx; idcoordx++) {
        x = (idcoordx + WarpX[idx] - offx) / mx;
        y = (idcoordy + WarpY[idx] - offy) / my;
        z = (idcoordz + WarpZ[idx] - offz) / mz;
        Bout[idx++] = InterpolatedGradX3DMirr(Bcoeff, Width, Height, Slice, x, y, z, SplineDegree) / mx;
      }
    }
  }
}
```