user: 请将以下的汇编代码转换成C语言代码:

```asm
<convolve_kernel_inside_fast>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %rsi,-0x60(%rbp)
  mov    %edx,-0x64(%rbp)
  mov    %ecx,-0x68(%rbp)
  mov    %r8d,-0x6c(%rbp)
  mov    %r9,-0x78(%rbp)
  movsd  %xmm0,-0x80(%rbp)
  mov    0x10(%rbp),%eax
  sub    $0x1,%eax
  shr    %eax
  mov    %eax,-0x38(%rbp)
  mov    0x18(%rbp),%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x34(%rbp)
  mov    0x18(%rbp),%rax
  add    $0x4,%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x30(%rbp)
  mov    0x18(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x64(%rbp),%eax
  imul   -0x68(%rbp),%eax
  imul   -0x6c(%rbp),%eax
  mov    %eax,%eax
  mov    %rax,-0x28(%rbp)
  movsd  D1(%rip),%xmm0
  divsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x58(%rbp),%rax
  mov    (%rax),%rax
  mov    -0x28(%rbp),%rdx
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <zero_vect@plt>
  mov    -0x2c(%rbp),%eax
  mov    %eax,-0x48(%rbp)
  jmp    L90
L88:
  mov    -0x30(%rbp),%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L127
L136:
  mov    -0x34(%rbp),%eax
  mov    %eax,-0x50(%rbp)
  jmp    L70
L104:
  mov    -0x64(%rbp),%edi
  mov    -0x68(%rbp),%ecx
  mov    -0x50(%rbp),%edx
  mov    -0x4c(%rbp),%esi
  mov    -0x48(%rbp),%eax
  mov    %edi,%r8d
  mov    %eax,%edi
  call   <gidz_general@plt>
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x20(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x18(%rbp),%xmm0
  jp     L82
  pxor   %xmm0,%xmm0
  ucomisd -0x18(%rbp),%xmm0
  je     L26
L82:
  mov    -0x38(%rbp),%eax
  neg    %eax
  mov    %eax,-0x3c(%rbp)
  jmp    L22
L48:
  mov    -0x38(%rbp),%eax
  neg    %eax
  mov    %eax,-0x40(%rbp)
  jmp    L13
L138:
  mov    -0x38(%rbp),%eax
  neg    %eax
  mov    %eax,-0x44(%rbp)
  jmp    L46
L124:
  mov    -0x44(%rbp),%edx
  mov    -0x38(%rbp),%eax
  add    %eax,%edx
  mov    -0x40(%rbp),%ecx
  mov    -0x38(%rbp),%eax
  lea    (%rcx,%rax,1),%esi
  mov    -0x3c(%rbp),%ecx
  mov    -0x38(%rbp),%eax
  add    %ecx,%eax
  mov    0x10(%rbp),%ecx
  mov    %eax,%edi
  call   <gidz_cube@plt>
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  mov    -0x58(%rbp),%rax
  mov    (%rax),%rbx
  mov    -0x44(%rbp),%edx
  mov    -0x50(%rbp),%eax
  add    %eax,%edx
  mov    -0x40(%rbp),%ecx
  mov    -0x4c(%rbp),%eax
  lea    (%rcx,%rax,1),%esi
  mov    -0x3c(%rbp),%ecx
  mov    -0x48(%rbp),%eax
  lea    (%rcx,%rax,1),%edi
  mov    -0x64(%rbp),%ecx
  mov    -0x68(%rbp),%eax
  mov    %ecx,%r8d
  mov    %eax,%ecx
  call   <gidz_general@plt>
  movslq %eax,%rdx
  shl    $0x3,%rdx
  add    %rbx,%rdx
  movsd  (%rdx),%xmm0
  mov    -0x58(%rbp),%rdx
  mov    (%rdx),%rdx
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  addsd  -0x88(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x44(%rbp)
L46:
  mov    -0x44(%rbp),%eax
  cmp    -0x38(%rbp),%eax
  jle    L124
  addl   $0x1,-0x40(%rbp)
L13:
  mov    -0x40(%rbp),%eax
  cmp    -0x38(%rbp),%eax
  jle    L138
  addl   $0x1,-0x3c(%rbp)
L22:
  mov    -0x3c(%rbp),%eax
  cmp    -0x38(%rbp),%eax
  jle    L48
L26:
  addl   $0x1,-0x50(%rbp)
L70:
  mov    -0x34(%rbp),%edx
  mov    -0x64(%rbp),%eax
  sub    %edx,%eax
  cmp    %eax,-0x50(%rbp)
  jb     L104
  addl   $0x1,-0x4c(%rbp)
L127:
  mov    -0x30(%rbp),%edx
  mov    -0x68(%rbp),%eax
  sub    %edx,%eax
  cmp    %eax,-0x4c(%rbp)
  jb     L136
  addl   $0x1,-0x48(%rbp)
L90:
  mov    -0x2c(%rbp),%edx
  mov    -0x6c(%rbp),%eax
  sub    %edx,%eax
  cmp    %eax,-0x48(%rbp)
  jb     L88
  nop
  nop
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: label: D1
type_guessed: f64

tool: D1: 1

assistant: ```c
void convolve_kernel_inside_fast(double **outmap, double *inmap, unsigned in_extx, unsigned in_exty, unsigned in_extz, double *kernel, unsigned kernel_size, double normfac, unsigned ignored[3]) {
  double dval, inv_normfac;
  unsigned long out_nvox;
  unsigned indx, indy, indz;
  int indx2, indy2, indz2, margin, marginx, marginy, marginz;
  margin = (kernel_size - 1) / 2;
  marginx = margin + ignored[0];
  marginy = margin + ignored[1];
  marginz = margin + ignored[2];
  out_nvox = in_extx * in_exty * in_extz;
  inv_normfac = 1.0 / normfac;
  zero_vect(*outmap, out_nvox);
  for (indz = marginz; indz < in_extz - marginz; indz++)
    for (indy = marginy; indy < in_exty - marginy; indy++)
      for (indx = marginx; indx < in_extx - marginx; indx++) {
        dval = (*(inmap + gidz_general(indz, indy, indx, in_exty, in_extx))) * inv_normfac;
        if (dval != 0.0)
          for (indz2 = -margin; indz2 <= margin; indz2++)
            for (indy2 = -margin; indy2 <= margin; indy2++)
              for (indx2 = -margin; indx2 <= margin; indx2++) {
                *(*outmap + gidz_general(indz + indz2, indy + indy2, indx + indx2, in_exty, in_extx)) += *(kernel + gidz_cube(indz2 + margin, indy2 + margin, indx2 + margin, kernel_size)) * dval;
              }
      }
}
```