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

```asm
<BatchInterpolatedThree3DZero>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x160,%rsp
  mov    %rdi,-0xa8(%rbp)
  mov    %rsi,-0xb0(%rbp)
  mov    %rdx,-0xb8(%rbp)
  mov    %rcx,-0xc0(%rbp)
  mov    %r8,-0xc8(%rbp)
  mov    %r9,-0xd0(%rbp)
  movsd  %xmm0,-0x120(%rbp)
  movsd  %xmm1,-0x128(%rbp)
  movsd  %xmm2,-0x130(%rbp)
  movsd  %xmm3,-0x138(%rbp)
  movsd  %xmm4,-0x140(%rbp)
  movsd  %xmm5,-0x148(%rbp)
  movsd  %xmm6,-0x150(%rbp)
  movsd  %xmm7,-0x158(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xd8(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xe0(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0xe8(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0xf0(%rbp)
  mov    0x30(%rbp),%rax
  mov    %rax,-0xf8(%rbp)
  mov    0x38(%rbp),%rax
  mov    %rax,-0x100(%rbp)
  mov    0x40(%rbp),%rax
  mov    %rax,-0x108(%rbp)
  mov    0x48(%rbp),%rax
  mov    %rax,-0x110(%rbp)
  mov    0x50(%rbp),%rax
  mov    %rax,-0x118(%rbp)
  mov    0x78(%rbp),%rax
  mov    %rax,-0x160(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movq   $0x0,-0x40(%rbp)
  movq   $0x0,-0x28(%rbp)
  jmp    L22
L33:
  movq   $0x0,-0x30(%rbp)
  jmp    L8
L85:
  movq   $0x0,-0x38(%rbp)
  jmp    L3
L83:
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x160(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  ucomiss %xmm1,%xmm0
  jp     L81
  pxor   %xmm1,%xmm1
  ucomiss %xmm1,%xmm0
  jne    L81
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xe8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xf0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x1(%rax),%rdx
  mov    %rdx,-0x40(%rbp)
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  jmp    L89
L81:
  pxor   %xmm0,%xmm0
  cvtsi2sdq -0x38(%rbp),%xmm0
  subsd  -0x128(%rbp),%xmm0
  divsd  -0x130(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdq -0x30(%rbp),%xmm0
  subsd  -0x140(%rbp),%xmm0
  divsd  -0x148(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdq -0x28(%rbp),%xmm0
  subsd  -0x158(%rbp),%xmm0
  divsd  0x70(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mov    -0x20(%rbp),%rdi
  lea    -0x78(%rbp),%r9
  lea    -0x80(%rbp),%r10
  lea    -0x88(%rbp),%rcx
  lea    -0x90(%rbp),%rdx
  lea    -0x98(%rbp),%rsi
  lea    -0xa0(%rbp),%rax
  sub    $0x8,%rsp
  push   0x80(%rbp)
  push   0x68(%rbp)
  push   0x60(%rbp)
  push   0x58(%rbp)
  push   -0x118(%rbp)
  push   -0x110(%rbp)
  push   -0x108(%rbp)
  lea    -0x48(%rbp),%r8
  push   %r8
  lea    -0x50(%rbp),%r8
  push   %r8
  lea    -0x58(%rbp),%r8
  push   %r8
  lea    -0x60(%rbp),%r8
  push   %r8
  lea    -0x68(%rbp),%r8
  push   %r8
  lea    -0x70(%rbp),%r8
  push   %r8
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %rax,%rdi
  call   <InterpolatedThree3DZero@plt>
  add    $0x70,%rsp
  movsd  -0xa0(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x98(%rbp),%xmm0
  divsd  -0x130(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x90(%rbp),%xmm0
  divsd  -0x148(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x88(%rbp),%xmm0
  divsd  0x70(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x80(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x78(%rbp),%xmm0
  divsd  -0x130(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x70(%rbp),%xmm0
  divsd  -0x148(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x68(%rbp),%xmm0
  divsd  0x70(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x60(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xe8(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x58(%rbp),%xmm0
  divsd  -0x130(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xf0(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm0
  divsd  -0x148(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x48(%rbp),%xmm0
  divsd  0x70(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  lea    0x1(%rax),%rdx
  mov    %rdx,-0x40(%rbp)
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
L89:
  addq   $0x1,-0x38(%rbp)
L3:
  pxor   %xmm1,%xmm1
  cvtsi2sdq -0x38(%rbp),%xmm1
  movsd  -0x120(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L83
  addq   $0x1,-0x30(%rbp)
L8:
  pxor   %xmm1,%xmm1
  cvtsi2sdq -0x30(%rbp),%xmm1
  movsd  -0x138(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L85
  addq   $0x1,-0x28(%rbp)
L22:
  pxor   %xmm1,%xmm1
  cvtsi2sdq -0x28(%rbp),%xmm1
  movsd  -0x150(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L33
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L76
  call   <__stack_chk_fail@plt>
L76:
  leave
  ret
```

assistant: ```c
extern void BatchInterpolatedThree3DZero(float *Bo1, float *Bo1gx, float *Bo1gy, float *Bo1gz, float *Bo2, float *Bo2gx, float *Bo2gy, float *Bo2gz, float *Bo3, float *Bo3gx, float *Bo3gy, float *Bo3gz, float *Bcoeff1, float *Bcoeff2, float *Bcoeff3, long Width, long Height, long Slice, double nx, double offx, double mx, double ny, double offy, double my, double nz, double offz, double mz, float *mask, long SplineDegree) {
  double x, y, z;
  double out1, gradx1, grady1, gradz1;
  double out2, gradx2, grady2, gradz2;
  double out3, gradx3, grady3, gradz3;
  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++) {
        if (mask[idx] == 0) {
          Bo1[idx] = 0;
          Bo1gx[idx] = 0;
          Bo1gy[idx] = 0;
          Bo1gz[idx] = 0;
          Bo2[idx] = 0;
          Bo2gx[idx] = 0;
          Bo2gy[idx] = 0;
          Bo2gz[idx] = 0;
          Bo3[idx] = 0;
          Bo3gx[idx] = 0;
          Bo3gy[idx] = 0;
          Bo3gz[idx++] = 0;
        } else {
          x = (idcoordx - offx) / mx;
          y = (idcoordy - offy) / my;
          z = (idcoordz - offz) / mz;
          InterpolatedThree3DZero(&out1, &gradx1, &grady1, &gradz1, &out2, &gradx2, &grady2, &gradz2, &out3, &gradx3, &grady3, &gradz3, Bcoeff1, Bcoeff2, Bcoeff3, Width, Height, Slice, x, y, z, SplineDegree);
          Bo1[idx] = out1;
          Bo1gx[idx] = gradx1 / mx;
          Bo1gy[idx] = grady1 / my;
          Bo1gz[idx] = gradz1 / mz;
          Bo2[idx] = out2;
          Bo2gx[idx] = gradx2 / mx;
          Bo2gy[idx] = grady2 / my;
          Bo2gz[idx] = gradz2 / mz;
          Bo3[idx] = out3;
          Bo3gx[idx] = gradx3 / mx;
          Bo3gy[idx] = grady3 / my;
          Bo3gz[idx++] = gradz3 / mz;
        }
      }
    }
  }
}
```