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

```asm
<metastable_cooling_SI>:
  endbr64
  sub    $0xd8,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0xc8(%rsp)
  xor    %eax,%eax
  movq   $0x0,0x20(%rsp)
  mov    D63(%rip),%rax
  mov    %rax,0x28(%rsp)
  mov    D91(%rip),%rax
  mov    %rax,0x30(%rsp)
  mov    D29(%rip),%rax
  mov    %rax,0x38(%rsp)
  movl   $0x0,(%rsp)
  movl   $0x2,0x4(%rsp)
  movl   $0x1,0x8(%rsp)
  movl   $0x1,0xc(%rsp)
  movl   $0x0,0x10(%rsp)
  movl   $0x3,0x14(%rsp)
  movl   $0x2,0x18(%rsp)
  movl   $0x3,0x1c(%rsp)
  movq   $0x0,0x40(%rsp)
  mov    D16(%rip),%rax
  mov    %rax,0x48(%rsp)
  mov    D69(%rip),%rax
  mov    %rax,0x50(%rsp)
  mov    D92(%rip),%rax
  mov    %rax,0x58(%rsp)
  movq   $0x0,0x60(%rsp)
  mov    D24(%rip),%rax
  mov    %rax,0x68(%rsp)
  mov    D14(%rip),%rax
  mov    %rax,0x70(%rsp)
  mov    D52(%rip),%rax
  mov    %rax,0x78(%rsp)
  movq   $0x0,0x80(%rsp)
  mov    D30(%rip),%rax
  mov    %rax,0x88(%rsp)
  mov    D110(%rip),%rax
  mov    %rax,0x90(%rsp)
  mov    D32(%rip),%rax
  mov    %rax,0x98(%rsp)
  movsd  D125(%rip),%xmm4
  comisd %xmm3,%xmm4
  jbe    L22
  movsd  D127(%rip),%xmm4
  movsd  %xmm4,0xa8(%rsp)
  movsd  %xmm4,0xb0(%rsp)
  movsd  %xmm4,0xb8(%rsp)
L22:
  comisd D125(%rip),%xmm3
  jb     L21
  mov    D55(%rip),%rax
  mov    %rax,0xa8(%rsp)
  movq   $0x0,0xb0(%rsp)
  movq   $0x0,0xb8(%rsp)
L21:
  mov    %rsp,%rcx
  lea    0x20(%rsp),%rdx
  sub    $0x8,%rsp
  lea    0xa8(%rsp),%rax
  push   %rax
  lea    0x90(%rsp),%rax
  push   %rax
  lea    0x78(%rsp),%rax
  push   %rax
  lea    0x60(%rsp),%r9
  lea    0x30(%rsp),%r8
  mov    $0x3,%esi
  mov    $0x3,%edi
  call   <computeMetastableCoolingRate@plt>
  add    $0x20,%rsp
  mov    0xc8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L23
  add    $0xd8,%rsp
  ret
L23:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_label: D63
datatype: f64
name: unpack_data
arguments: data_label: D91
datatype: f64
name: unpack_data
arguments: data_label: D29
datatype: f64
name: unpack_data
arguments: data_label: D16
datatype: f64
name: unpack_data
arguments: data_label: D69
datatype: f64
name: unpack_data
arguments: data_label: D92
datatype: f64
name: unpack_data
arguments: data_label: D24
datatype: byte[8]
name: unpack_data
arguments: data_label: D14
datatype: f64
name: unpack_data
arguments: data_label: D52
datatype: f64
name: unpack_data
arguments: data_label: D30
datatype: byte[8]
name: unpack_data
arguments: data_label: D110
datatype: byte[8]
name: unpack_data
arguments: data_label: D32
datatype: byte[8]
name: unpack_data
arguments: data_label: D125
datatype: f64[2]
name: unpack_data
arguments: data_label: D127
datatype: f64
name: unpack_data
arguments: data_label: D55
datatype: byte[8]

tool: D63: 1

tool: D91: 5

tool: D29: 9

tool: D16: 13000

tool: D69: 19000

tool: D92: 32000

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

tool: D14: 1.8

tool: D52: 0.36

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

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

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

tool: {"D125": [10000.0,0.57]}

tool: D127: 0.57

tool: {"D55": [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);
}
```