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

```asm
<metastable_cooling_SI>:
  endbr64
  sub    $0xd8,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0xc8(%rsp)
  xor    %eax,%eax
  mov    D106(%rip),%rax
  movq   $0x0,0x20(%rsp)
  movq   $0x0,0x40(%rsp)
  mov    %rax,0x28(%rsp)
  mov    D34(%rip),%rax
  movq   $0x0,0x60(%rsp)
  mov    %rax,0x30(%rsp)
  mov    D36(%rip),%rax
  movq   $0x0,0x80(%rsp)
  mov    %rax,0x38(%rsp)
  movabs $0x200000000,%rax
  mov    %rax,(%rsp)
  movabs $0x100000001,%rax
  mov    %rax,0x8(%rsp)
  movabs $0x300000000,%rax
  mov    %rax,0x10(%rsp)
  add    $0x2,%rax
  mov    %rax,0x18(%rsp)
  mov    D142(%rip),%rax
  mov    %rax,0x48(%rsp)
  mov    D66(%rip),%rax
  mov    %rax,0x50(%rsp)
  mov    D149(%rip),%rax
  mov    %rax,0x58(%rsp)
  mov    D127(%rip),%rax
  mov    %rax,0x68(%rsp)
  mov    D17(%rip),%rax
  mov    %rax,0x70(%rsp)
  mov    D47(%rip),%rax
  mov    %rax,0x78(%rsp)
  mov    D67(%rip),%rax
  mov    %rax,0x88(%rsp)
  mov    D7(%rip),%rax
  mov    %rax,0x90(%rsp)
  movsd  D12(%rip),%xmm4
  mov    D91(%rip),%rax
  comisd %xmm3,%xmm4
  mov    %rax,0x98(%rsp)
  jbe    L4
  movsd  D135(%rip),%xmm5
  movsd  %xmm5,0xa8(%rsp)
  movsd  %xmm5,0xb0(%rsp)
  movsd  %xmm5,0xb8(%rsp)
L4:
  comisd %xmm4,%xmm3
  jb     L28
  mov    D6(%rip),%rax
  movq   $0x0,0xb0(%rsp)
  movq   $0x0,0xb8(%rsp)
  mov    %rax,0xa8(%rsp)
L28:
  mov    %rsp,%rcx
  lea    0x20(%rsp),%rdx
  sub    $0x8,%rsp
  mov    $0x3,%esi
  lea    0xa8(%rsp),%rax
  mov    $0x3,%edi
  push   %rax
  lea    0x90(%rsp),%rax
  push   %rax
  lea    0x78(%rsp),%rax
  push   %rax
  lea    0x60(%rsp),%r9
  lea    0x30(%rsp),%r8
  call   <computeMetastableCoolingRate@plt>
  add    $0x20,%rsp
  mov    0xc8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L5
  add    $0xd8,%rsp
  ret
L5:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_label: D106
type: f64
name: retrieve_data
arguments: data_label: D34
type: f64
name: retrieve_data
arguments: data_label: D36
type: f64
name: retrieve_data
arguments: data_label: D142
type: f64
name: retrieve_data
arguments: data_label: D66
type: f64
name: retrieve_data
arguments: data_label: D149
type: f64
name: retrieve_data
arguments: data_label: D127
type: byte[8]
name: retrieve_data
arguments: data_label: D17
type: f64
name: retrieve_data
arguments: data_label: D47
type: f64
name: retrieve_data
arguments: data_label: D67
type: byte[8]
name: retrieve_data
arguments: data_label: D7
type: byte[8]
name: retrieve_data
arguments: data_label: D12
type: f64
name: retrieve_data
arguments: data_label: D91
type: byte[8]
name: retrieve_data
arguments: data_label: D135
type: f64
name: retrieve_data
arguments: data_label: D6
type: byte[8]

tool: D106: 1

tool: D34: 5

tool: D36: 9

tool: D142: 13000

tool: D66: 19000

tool: D149: 32000

tool: {"D127": [0x3b,0xdf,0x4f,0x8d,0x97,0x6e,0xa2,0x3f]}

tool: D17: 1.8

tool: D47: 0.36

tool: {"D67": [0xfd,0xe5,0xfd,0xf9,0x62,0x61,0x6b,0x3e]}

tool: {"D7": [0x76,0xd4,0xb7,0x1d,0x31,0xb2,0x58,0x3e]}

tool: D12: 10000

tool: {"D91": [0x48,0xaf,0xbc,0x9a,0xf2,0xd7,0x5a,0x3e]}

tool: D135: 0.57

tool: {"D6": [0xa4,0x70,0x3d,0x0a,0xd7,0xa3,0xc0,0xbf]}

assistant: ```c
double metastable_cooling_SI(const double densThis, const double Hdens, const double elecDens, const double T) {
  const double weights[3 + 1] = {0.0, 1.0, 5.0, 9.0};
  const int lower[3 + 1] = {0, 2, 1, 1};
  const int upper[3 + 1] = {0, 3, 2, 3};
  const double energy[3 + 1] = {0.0, 1.3e4, 1.9e4, 3.2e4};
  const double EinsA[3 + 1] = {0.0, 3.6e-2, 1.8, 0.36};
  const double elec_a[3 + 1] = {0.0, 5.1e-8, 2.3e-8, 2.5e-8};
  double elec_b[3 + 1];
  if (T < 1e4) {
    elec_b[1] = 0.57;
    elec_b[2] = 0.57;
    elec_b[3] = 0.57;
  }
  if (T >= 1e4) {
    elec_b[1] = -0.13;
    elec_b[2] = 0.0;
    elec_b[3] = 0.0;
  }
  return computeMetastableCoolingRate(3, 3, weights, lower, upper, energy, EinsA, elec_a, elec_b, densThis, Hdens, elecDens, T);
}
```