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

```asm
<glbChiSpectrumCalib>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  mov    %edi,%r15d
  push   %r14
  mov    %esi,%r14d
  push   %r13
  mov    %r8,%r13
  push   %r12
  push   %rbx
  sub    $0x78,%rsp
  mov    %rcx,-0x80(%rbp)
  mov    %esi,-0x68(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  call   <glbGetNumberOfBins@plt>
  mov    %r14d,%esi
  mov    %r15d,%edi
  mov    %eax,%ebx
  mov    %eax,-0x84(%rbp)
  call   <glbGetRuleRatePtr@plt>
  mov    %rsp,%rcx
  mov    %rax,%r14
  movslq %ebx,%rax
  lea    0xf(,%rax,8),%rax
  mov    %rax,%rsi
  mov    %rax,%rdx
  and    $0xfffffffffffff000,%rsi
  and    $0xfffffffffffffff0,%rdx
  sub    %rsi,%rcx
  cmp    %rcx,%rsp
  je     L46
L91:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L91
L46:
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  jne    L75
L95:
  mov    %rax,%rdx
  mov    %rsp,%rcx
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rbx
  sub    %rax,%rcx
  and    $0xfffffffffffffff0,%rdx
  cmp    %rcx,%rsp
  je     L87
L90:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L90
L87:
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  jne    L0
L105:
  mov    %r15d,%edi
  lea    -0x40(%rbp),%rdx
  lea    -0x48(%rbp),%rsi
  mov    %rsp,%r12
  call   <glbGetEminEmax@plt>
  mov    -0x68(%rbp),%esi
  mov    %r15d,%edi
  lea    -0x5c(%rbp),%rcx
  lea    -0x60(%rbp),%rdx
  call   <glbGetEnergyWindowBins@plt>
  mov    -0x68(%rbp),%esi
  lea    -0x50(%rbp),%rcx
  lea    -0x58(%rbp),%rdx
  mov    %r15d,%edi
  call   <glbGetBGCenters@plt>
  mov    -0x80(%rbp),%rax
  mov    -0x68(%rbp),%esi
  mov    %r15d,%edi
  movsd  D0(%rip),%xmm0
  movsd  -0x40(%rbp),%xmm2
  movsd  (%rax),%xmm6
  movsd  -0x48(%rbp),%xmm1
  movsd  %xmm2,-0xa0(%rbp)
  addsd  %xmm0,%xmm6
  addsd  0x10(%rax),%xmm0
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm1,-0x98(%rbp)
  movsd  %xmm6,-0x78(%rbp)
  movsd  %xmm0,-0x70(%rbp)
  movsd  0x18(%rax),%xmm0
  addsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x90(%rbp)
  call   <glbGetSignalFitRatePtr@plt>
  mov    -0x84(%rbp),%edx
  movsd  -0xa0(%rbp),%xmm2
  mov    %rbx,%rsi
  mov    %eax,%edi
  mov    -0x80(%rbp),%rax
  movsd  -0x98(%rbp),%xmm1
  movsd  0x8(%rax),%xmm0
  call   <glbShiftEnergyScale@plt>
  movsd  -0x40(%rbp),%xmm2
  mov    -0x68(%rbp),%esi
  mov    %r15d,%edi
  movsd  -0x48(%rbp),%xmm1
  movsd  %xmm2,-0xa0(%rbp)
  movsd  %xmm1,-0x98(%rbp)
  call   <glbGetBGFitRatePtr@plt>
  movsd  -0xa0(%rbp),%xmm2
  mov    %r12,%rsi
  movsd  -0x98(%rbp),%xmm1
  mov    -0x84(%rbp),%edx
  movsd  -0x90(%rbp),%xmm0
  mov    %eax,%edi
  call   <glbShiftEnergyScale@plt>
  movslq -0x60(%rbp),%r15
  cmp    -0x5c(%rbp),%r15d
  jg     L45
  pxor   %xmm5,%xmm5
  movsd  %xmm5,-0x68(%rbp)
  nopl   0x0(%rax)
L110:
  movsd  -0x78(%rbp),%xmm2
  movsd  -0x70(%rbp),%xmm1
  mulsd  (%rbx,%r15,8),%xmm2
  movsd  (%r14,%r15,8),%xmm0
  mulsd  (%r12,%r15,8),%xmm1
  add    $0x1,%r15
  addsd  %xmm2,%xmm1
  call   <glb_likelihood@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  addsd  -0x68(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  cmp    %r15d,-0x5c(%rbp)
  jge    L110
L54:
  mov    -0x80(%rbp),%r14
  movsd  0x0(%r13),%xmm2
  pxor   %xmm1,%xmm1
  movsd  (%r14),%xmm0
  call   <glb_prior@plt>
  movsd  0x8(%r14),%xmm0
  movsd  0x8(%r13),%xmm2
  pxor   %xmm1,%xmm1
  mov    %rax,%rbx
  call   <glb_prior@plt>
  movsd  0x10(%r13),%xmm2
  movsd  -0x58(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  add    %rax,%rbx
  call   <glb_prior@plt>
  movsd  0x18(%r13),%xmm2
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x90(%rbp),%xmm0
  add    %rax,%rbx
  call   <glb_prior@plt>
  pxor   %xmm0,%xmm0
  add    %rax,%rbx
  cvtsi2sd %rbx,%xmm0
  addsd  -0x68(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L5
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
  nopl   0x0(%rax)
L75:
  orq    $0x0,-0x8(%rsp,%rdx,1)
  jmp    L95
L45:
  pxor   %xmm7,%xmm7
  movsd  %xmm7,-0x68(%rbp)
  jmp    L54
L0:
  orq    $0x0,-0x8(%rsp,%rdx,1)
  jmp    L105
L5:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: tag: D0
data_type: f64

tool: D0: 1

assistant: ```c
double glbChiSpectrumCalib(int exp, int rule, int n_params, double *x, double *errors, void *user_data) {
  int n_bins = glbGetNumberOfBins(exp);
  double *true_rates = glbGetRuleRatePtr(exp, rule);
  double signal_fit_rates[n_bins];
  double bg_fit_rates[n_bins];
  double signal_norm;
  double bg_norm_center, bg_tilt_center;
  double bg_norm, bg_tilt;
  int ew_low, ew_high;
  double emin, emax;
  double fit_rate, true_rate;
  double chi2 = 0.0;
  int i;
  glbGetEminEmax(exp, &emin, &emax);
  glbGetEnergyWindowBins(exp, rule, &ew_low, &ew_high);
  glbGetBGCenters(exp, rule, &bg_norm_center, &bg_tilt_center);
  signal_norm = 1.0 + x[0];
  bg_norm = bg_norm_center * (1.0 + x[2]);
  bg_tilt = bg_tilt_center + x[3];
  glbShiftEnergyScale(x[1], glbGetSignalFitRatePtr(exp, rule), signal_fit_rates, n_bins, emin, emax);
  glbShiftEnergyScale(bg_tilt, glbGetBGFitRatePtr(exp, rule), bg_fit_rates, n_bins, emin, emax);
  for (i = ew_low; i <= ew_high; i++) {
    fit_rate = signal_norm * signal_fit_rates[i] + bg_norm * bg_fit_rates[i];
    chi2 += glb_likelihood(true_rates[i], fit_rate);
  }
  chi2 += glb_prior(x[0], 0.0, errors[0]) + glb_prior(x[1], 0.0, errors[1]) + glb_prior(bg_norm, bg_norm_center, errors[2]) + glb_prior(bg_tilt, bg_tilt_center, errors[3]);
  return chi2;
}
```