user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<coreGTRCATPROTMULT>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x180,%rsp
  mov    %rdi,-0x128(%rbp)
  movsd  %xmm0,-0x130(%rbp)
  mov    %esi,-0x134(%rbp)
  mov    %rdx,-0x140(%rbp)
  mov    %rcx,-0x148(%rbp)
  mov    %r8d,-0x138(%rbp)
  mov    %r9d,-0x14c(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x158(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0x160(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0x168(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0x170(%rbp)
  mov    0x30(%rbp),%rax
  mov    %rax,-0x178(%rbp)
  mov    0x40(%rbp),%rax
  mov    %rax,-0x180(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x108(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x100(%rbp)
  mov    0x38(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0xd8(%rbp)
  mov    -0xd8(%rbp),%rax
  mov    %rax,-0xd0(%rbp)
  mov    0x38(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0xc8(%rbp)
  mov    -0xc8(%rbp),%rax
  mov    %rax,-0xc0(%rbp)
  mov    0x38(%rbp),%eax
  imul   -0x134(%rbp),%eax
  mov    %eax,%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0xb8(%rbp)
  mov    -0xb8(%rbp),%rax
  mov    %rax,-0xf8(%rbp)
  movl   $0x0,-0x10c(%rbp)
  jmp    L91
L123:
  movl   $0x0,-0x110(%rbp)
  jmp    L21
L45:
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x110(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x128(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x110(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x128(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x110(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x110(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x128(%rbp),%rax
  lea    (%rdx,%rax,1),%rcx
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x110(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rcx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %eax,%edx
  mov    -0x110(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x130(%rbp),%xmm0
  mov    -0x110(%rbp),%eax
  cltq
  movsd  %xmm0,-0xa0(%rbp,%rax,8)
  addl   $0x1,-0x110(%rbp)
L21:
  cmpl   $0x12,-0x110(%rbp)
  jle    L45
  movl   $0x0,-0x114(%rbp)
  jmp    L109
L4:
  movl   $0x0,-0x110(%rbp)
  jmp    L108
L126:
  mov    -0x110(%rbp),%eax
  cltq
  movsd  -0xa0(%rbp,%rax,8),%xmm1
  mov    -0x114(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x140(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  mov    %eax,%ecx
  mov    -0xf8(%rbp),%rax
  lea    0x8(%rax),%rdx
  mov    %rdx,-0xf8(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sd %ecx,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x110(%rbp)
L108:
  cmpl   $0x12,-0x110(%rbp)
  jle    L126
  addl   $0x1,-0x114(%rbp)
L109:
  mov    -0x114(%rbp),%eax
  cmp    -0x134(%rbp),%eax
  jl     L4
  addl   $0x1,-0x10c(%rbp)
L91:
  mov    -0x10c(%rbp),%eax
  cmp    0x38(%rbp),%eax
  jl     L123
  mov    -0x138(%rbp),%eax
  mov    %eax,-0x114(%rbp)
  jmp    L129
L128:
  mov    -0x114(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x180(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x10c(%rbp)
  mov    -0x114(%rbp),%edx
  mov    %edx,%eax
  shl    $0x2,%eax
  add    %edx,%eax
  shl    $0x2,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x178(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0xb0(%rbp)
  mov    -0x10c(%rbp),%eax
  imul   -0x134(%rbp),%eax
  mov    %eax,%edx
  mov    -0x114(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x148(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  add    %eax,%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0xf8(%rbp)
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0xd0(%rbp)
  mov    -0x10c(%rbp),%edx
  mov    %edx,%eax
  shl    $0x3,%eax
  add    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0xc0(%rbp)
  mov    -0xb0(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  mov    -0xf8(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0xb0(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  movsd  -0xf0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  mov    -0xc0(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  -0xa8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  mov    -0xd0(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  -0xa8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  movl   $0x1,-0x110(%rbp)
  jmp    L116
L62:
  mov    -0x110(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x110(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  movsd  -0xf0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  mov    -0x110(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xa8(%rbp),%xmm0
  movsd  -0xe8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  mov    -0x110(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xa8(%rbp),%xmm0
  movsd  -0xe0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  addl   $0x1,-0x110(%rbp)
L116:
  cmpl   $0x12,-0x110(%rbp)
  jle    L62
  movsd  D7(%rip),%xmm0
  divsd  -0xf0(%rbp),%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  movsd  -0xe8(%rbp),%xmm0
  mulsd  -0xf0(%rbp),%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  movsd  -0xe0(%rbp),%xmm0
  mulsd  -0xf0(%rbp),%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  mov    -0x114(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x158(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xe8(%rbp),%xmm0
  movsd  -0x108(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x108(%rbp)
  mov    -0x114(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x160(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0xe8(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0xe0(%rbp),%xmm0
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x100(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x100(%rbp)
  addl   $0x1,-0x114(%rbp)
L129:
  mov    -0x114(%rbp),%eax
  cmp    -0x14c(%rbp),%eax
  jl     L128
  mov    -0x168(%rbp),%rax
  movsd  -0x108(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x170(%rbp),%rax
  movsd  -0x100(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xb8(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    -0xc8(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    -0xd8(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L11
  call   <__stack_chk_fail@plt>
L11:
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D7
type_guessed: f64

tool: D7: 1

assistant: ```c
void coreGTRCATPROTMULT(double *EIGN, double lz, int numberOfCategories, double *rptr, int *cptr, int lower, int upper, double *wrptr, double *wr2ptr, double *ext_dlnLdlz, double *ext_d2lnLdlz2, double *sumtable, int numberOfModels, int *modelptr) {
  double *sum;
  int i, l;
  double dlnLdlz = 0.0;
  double d2lnLdlz2 = 0.0;
  double *d1, *d_start;
  double *e, *s, *e_start, *s_start;
  int model;
  double dd[19];
  double tmp;
  double inv_Li, dlnLidlz, d2lnLidlz2;
  e = e_start = (double *)malloc(19 * numberOfModels * sizeof(double));
  s = s_start = (double *)malloc(19 * numberOfModels * sizeof(double));
  d1 = d_start = (double *)malloc(numberOfModels * numberOfCategories * 19 * sizeof(double));
  for (model = 0; model < numberOfModels; model++) {
    for (l = 0; l < 19; l++) {
      e[model * 19 + l] = EIGN[model * 19 + l] * EIGN[model * 19 + l];
      s[model * 19 + l] = EIGN[model * 19 + l];
      dd[l] = s[model * 19 + l] * lz;
    }
    for (i = 0; i < numberOfCategories; i++) {
      for (l = 0; l < 19; l++)
        *d1++ = exp(dd[l] * rptr[i]);
    }
  }
  for (i = lower; i < upper; i++) {
    model = modelptr[i];
    sum = &sumtable[20 * i];
    d1 = &(d_start[model * numberOfCategories * 19 + 19 * cptr[i]]);
    e = &(e_start[model * 19]);
    s = &(s_start[model * 19]);
    inv_Li = sum[0];
    inv_Li += (tmp = d1[0] * sum[1]);
    dlnLidlz = tmp * s[0];
    d2lnLidlz2 = tmp * e[0];
    for (l = 1; l < 19; l++) {
      inv_Li += (tmp = d1[l] * sum[l + 1]);
      dlnLidlz += tmp * s[l];
      d2lnLidlz2 += tmp * e[l];
    }
    inv_Li = 1.0 / inv_Li;
    dlnLidlz *= inv_Li;
    d2lnLidlz2 *= inv_Li;
    dlnLdlz += wrptr[i] * dlnLidlz;
    d2lnLdlz2 += wr2ptr[i] * (d2lnLidlz2 - dlnLidlz * dlnLidlz);
  }
  *ext_dlnLdlz = dlnLdlz;
  *ext_d2lnLdlz2 = d2lnLdlz2;
  free(d_start);
  free(s_start);
  free(e_start);
}
```