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

```asm
<BatchInterpolatedGradY3DTranMirr>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %r9,%r13
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  push   %rbx
  mov    %rdx,%rbx
  sub    $0x88,%rsp
  movsd  %xmm0,0x20(%rsp)
  mulsd  %xmm3,%xmm0
  movsd  %xmm1,0x28(%rsp)
  movapd %xmm3,%xmm1
  movsd  %xmm2,0x30(%rsp)
  movapd %xmm6,%xmm2
  mov    %rsi,0x58(%rsp)
  mov    %rcx,0x48(%rsp)
  mulsd  %xmm6,%xmm0
  mov    %r8,0x70(%rsp)
  movsd  %xmm4,0x60(%rsp)
  movsd  %xmm5,0x50(%rsp)
  movsd  %xmm7,0x68(%rsp)
  movapd %xmm0,%xmm3
  movsd  %xmm0,0x78(%rsp)
  pxor   %xmm0,%xmm0
  comisd %xmm0,%xmm3
  jbe    L84
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L72:
  movl   $0x0,(%r12,%rax,4)
  pxor   %xmm0,%xmm0
  add    $0x1,%rax
  cvtsi2sd %rax,%xmm0
  comisd %xmm0,%xmm3
  ja     L72
  cmpq   $0x0,0x70(%rsp)
  jle    L91
L117:
  movq   $0x0,0x40(%rsp)
  movq   $0x0,0x38(%rsp)
  nopl   0x0(%rax)
L96:
  xor    %r14d,%r14d
  cmpq   $0x0,0x48(%rsp)
  jle    L23
  nop
L50:
  test   %rbx,%rbx
  jle    L102
  mov    0x58(%rsp),%rax
  pxor   %xmm5,%xmm5
  pxor   %xmm6,%xmm6
  xor    %r15d,%r15d
  mov    0x38(%rsp),%rdx
  cvtsi2sd %r14,%xmm5
  subsd  0x60(%rsp),%xmm5
  divsd  0x50(%rsp),%xmm5
  cvtsi2sdq 0x40(%rsp),%xmm6
  subsd  0x68(%rsp),%xmm6
  divsd  0xc0(%rsp),%xmm6
  lea    (%rax,%rdx,4),%rbp
  nopl   0x0(%rax)
L3:
  movsd  0x20(%rsp),%xmm0
  pxor   %xmm4,%xmm4
  mov    %r13,%rsi
  mov    %r12,%rdi
  movss  0x0(%rbp,%r15,4),%xmm3
  cvtsi2sd %r15,%xmm4
  add    $0x1,%r15
  subsd  0x28(%rsp),%xmm4
  movsd  %xmm6,0x18(%rsp)
  divsd  0x30(%rsp),%xmm4
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm1,(%rsp)
  call   <InterpolatedGradY3DTranMirr@plt>
  cmp    %r15,%rbx
  movsd  (%rsp),%xmm1
  movsd  0x8(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm6
  jne    L3
  add    %rbx,0x38(%rsp)
L102:
  add    $0x1,%r14
  cmp    %r14,0x48(%rsp)
  jne    L50
L23:
  addq   $0x1,0x40(%rsp)
  mov    0x40(%rsp),%rax
  cmp    %rax,0x70(%rsp)
  jg     L96
L10:
  movsd  0x78(%rsp),%xmm0
  pxor   %xmm7,%xmm7
  comisd %xmm7,%xmm0
  jbe    L114
L91:
  movsd  0x78(%rsp),%xmm1
  movsd  0x50(%rsp),%xmm2
  xor    %eax,%eax
  nopl   0x0(%rax)
L112:
  pxor   %xmm0,%xmm0
  cvtss2sd (%r12,%rax,4),%xmm0
  divsd  %xmm2,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r12,%rax,4)
  add    $0x1,%rax
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  comisd %xmm0,%xmm1
  ja     L112
L114:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L84:
  cmpq   $0x0,0x70(%rsp)
  jg     L117
  jmp    L10
```

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