user: 请从给定的汇编指令推导出C源代码:

```asm
<radiativeTransitions>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xa8,%rsp
  mov    %edi,0x74(%rsp)
  mov    0xf8(%rsp),%r14
  mov    %rsi,0x78(%rsp)
  mov    %rdx,0x80(%rsp)
  mov    %rcx,0x88(%rsp)
  mov    %r8,0x90(%rsp)
  mov    %r9,0x98(%rsp)
  test   %edi,%edi
  jle    L13
  mov    0xe0(%rsp),%rax
  mov    $0x1,%ebp
  lea    0x8(%rax),%rbx
  mov    0xf0(%rsp),%eax
  sub    $0x1,%eax
  lea    0x8(,%rax,8),%rax
  mov    %rax,0x30(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L37:
  mov    0x80(%rsp),%rax
  mov    0x88(%rsp),%rcx
  mov    0x78(%rsp),%rsi
  movsd  D33(%rip),%xmm1
  movslq (%rcx,%rbp,4),%rcx
  movslq (%rax,%rbp,4),%rax
  movsd  (%rsi,%rcx,8),%xmm0
  movsd  (%rsi,%rax,8),%xmm5
  lea    0x0(,%rcx,8),%r13
  lea    0x0(,%rax,8),%r15
  mov    0x90(%rsp),%rsi
  mov    %rcx,0x68(%rsp)
  movsd  %xmm0,0x28(%rsp)
  movsd  (%rsi,%rbp,8),%xmm4
  mov    0x98(%rsp),%rsi
  movsd  %xmm5,0x20(%rsp)
  mov    %rax,0x60(%rsp)
  movsd  (%rsi,%rcx,8),%xmm0
  subsd  (%rsi,%rax,8),%xmm0
  movsd  %xmm4,0x18(%rsp)
  mulsd  D91(%rip),%xmm0
  mov    %r15,0x38(%rsp)
  divsd  D28(%rip),%xmm0
  divsd  %xmm0,%xmm1
  movsd  %xmm0,(%rbx)
  mov    %r13,0x40(%rsp)
  movapd %xmm1,%xmm0
  call   <calc_background_radiation@plt>
  mov    D33(%rip),%rax
  mov    $0x2,%edi
  movsd  %xmm0,0x48(%rsp)
  movq   %rax,%xmm0
  call   <pow@plt>
  movsd  (%rbx),%xmm0
  mov    $0x3,%edi
  mov    %eax,%r12d
  call   <pow@plt>
  pxor   %xmm1,%xmm1
  pxor   %xmm0,%xmm0
  movsd  0x18(%rsp),%xmm7
  cvtsi2sd %eax,%xmm1
  mulsd  D32(%rip),%xmm1
  mov    0xf0(%rsp),%eax
  cvtsi2sd %r12d,%xmm0
  divsd  %xmm1,%xmm0
  mulsd  %xmm0,%xmm7
  movsd  0x28(%rsp),%xmm0
  divsd  0x20(%rsp),%xmm0
  mulsd  %xmm7,%xmm0
  movsd  %xmm7,0x50(%rsp)
  movsd  %xmm0,0x58(%rsp)
  test   %eax,%eax
  jle    L67
  mov    0x100(%rsp),%r12
  add    0x100(%rsp),%r15
  movq   $0x0,0x8(%rsp)
  add    %r13,%r12
  mov    %r15,%r13
  xor    %r15d,%r15d
  cs nopw 0x0(%rax,%rax,1)
L94:
  mov    D33(%rip),%rax
  mov    $0x3,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  pxor   %xmm1,%xmm1
  movsd  (%rbx),%xmm0
  mov    $0x3,%edi
  cvtsi2sd %eax,%xmm1
  mulsd  0x18(%rsp),%xmm1
  movsd  %xmm1,0x10(%rsp)
  call   <pow@plt>
  mov    (%r12),%rdx
  pxor   %xmm5,%xmm5
  movsd  0x20(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm1
  movsd  (%rdx,%r15,1),%xmm3
  mov    0x0(%r13),%rdx
  mulsd  %xmm3,%xmm2
  mulsd  (%rdx,%r15,1),%xmm0
  divsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  mulsd  D58(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  D127(%rip),%xmm0
  mulsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm0
  mulsd  (%r14,%r15,1),%xmm0
  divsd  D87(%rip),%xmm0
  comisd %xmm0,%xmm5
  jbe    L66
  movsd  0x8(%rsp),%xmm6
  add    $0x8,%r15
  addsd  %xmm5,%xmm6
  movsd  %xmm6,0x8(%rsp)
  cmp    %r15,0x30(%rsp)
  jne    L94
  movapd %xmm6,%xmm0
L71:
  call   <beta_tau_finestructure@plt>
  mov    D95(%rip),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  comisd %xmm0,%xmm1
  jbe    L30
L73:
  movsd  0x48(%rsp),%xmm6
  movsd  0x58(%rsp),%xmm1
  add    $0x8,%rbx
  mov    0xe8(%rsp),%rax
  mov    0x60(%rsp),%rsi
  mulsd  %xmm6,%xmm1
  mov    0x38(%rsp),%rcx
  mov    0xe8(%rsp),%rdi
  movsd  %xmm0,(%rax,%rbp,8)
  mov    0x108(%rsp),%rax
  movsd  %xmm6,(%rax,%rbp,8)
  mov    0x110(%rsp),%rax
  mulsd  0x50(%rsp),%xmm6
  mov    (%rax,%rsi,8),%rax
  mov    0xe8(%rsp),%rsi
  movsd  (%rsi,%rbp,8),%xmm2
  lea    (%rax,%rcx,1),%rdx
  mov    0x40(%rsp),%rsi
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm2
  add    %rsi,%rax
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rdx)
  movsd  0x18(%rsp),%xmm0
  movsd  (%rdi,%rbp,8),%xmm2
  mov    0x68(%rsp),%rdi
  addsd  %xmm6,%xmm0
  mulsd  %xmm0,%xmm2
  addsd  (%rax),%xmm2
  movsd  %xmm2,(%rax)
  mov    0x110(%rsp),%rax
  mov    (%rax,%rdi,8),%rdx
  lea    (%rdx,%rcx,1),%rax
  mov    0xe8(%rsp),%rcx
  mulsd  (%rcx,%rbp,8),%xmm1
  addsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  mulsd  (%rcx,%rbp,8),%xmm0
  mov    %rsi,%rax
  add    $0x1,%rbp
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax)
  cmp    %ebp,0x74(%rsp)
  jge    L37
L13:
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L66:
  mov    D33(%rip),%rax
  mov    $0x3,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  pxor   %xmm1,%xmm1
  movsd  (%rbx),%xmm0
  mov    $0x3,%edi
  cvtsi2sd %eax,%xmm1
  mulsd  0x18(%rsp),%xmm1
  movsd  %xmm1,0x10(%rsp)
  call   <pow@plt>
  mov    (%r12),%rdx
  movsd  0x20(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm1
  movsd  (%rdx,%r15,1),%xmm3
  mov    0x0(%r13),%rdx
  mulsd  %xmm3,%xmm2
  mulsd  (%rdx,%r15,1),%xmm0
  divsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  mulsd  D58(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  D127(%rip),%xmm0
  mulsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm0
  mulsd  (%r14,%r15,1),%xmm0
  add    $0x8,%r15
  divsd  D87(%rip),%xmm0
  addsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x8(%rsp)
  cmp    0x30(%rsp),%r15
  jne    L94
  jmp    L71
  nopl   0x0(%rax,%rax,1)
L30:
  movsd  0x8(%rsp),%xmm0
  call   <beta_tau_finestructure@plt>
  jmp    L73
L67:
  movq   $0x0,0x8(%rsp)
  movsd  0x8(%rsp),%xmm0
  jmp    L71
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D33
datatype: f64
name: extract_symbolic_data
arguments: data_label: D91
datatype: byte[8]
name: extract_symbolic_data
arguments: data_label: D28
datatype: byte[8]
name: extract_symbolic_data
arguments: data_label: D32
datatype: byte[8]
name: extract_symbolic_data
arguments: data_label: D58
datatype: byte[8]
name: extract_symbolic_data
arguments: data_label: D127
datatype: f64
name: extract_symbolic_data
arguments: data_label: D87
datatype: f64
name: extract_symbolic_data
arguments: data_label: D95
datatype: f64

tool: D33: 29979250000

tool: {"D91": [0xcf,0xa7,0x60,0x9f,0xfd,0xe6,0xa3,0x3c]}

tool: {"D28": [0x79,0xe1,0xc8,0x2c,0xd7,0x67,0x80,0x3a]}

tool: {"D32": [0x79,0xe1,0xc8,0x2c,0xd7,0x67,0x90,0x3a]}

tool: {"D58": [0xe0,0x86,0x44,0x54,0xfb,0x21,0x39,0x40]}

tool: D127: 1

tool: D87: 100000

tool: D95: 0.5

assistant: ```c
void radiativeTransitions(const int transitions, const double *weights, const int *lower, const int *upper, const double *EinsA, const double *energy, double *nu, double *betatau, const int slabInd, const double *dx, double **popDensThis, double *Pnu, double **DFDX) {
  int i, k;
  double lw, uw;
  int li, ui;
  double Aul, lambda, background, Bul, Blu, tau;
  for (i = 1; i <= transitions; i++) {
    li = lower[i];
    lw = weights[li];
    ui = upper[i];
    uw = weights[ui];
    Aul = EinsA[i];
    nu[i] = 1.381e-16 * (energy[ui] - energy[li]) / 6.62618e-27;
    lambda = 2.997925e10 / nu[i];
    background = calc_background_radiation(lambda);
    Bul = pow(2.997925e10, 2e0) / (2e0 * 6.62618e-27 * pow(nu[i], 3e0)) * Aul;
    Blu = uw / lw * Bul;
    tau = 0e0;
    for (k = 0; k < slabInd; k++) {
      tau += ((0e0) > ((Aul * pow(2.997925e10, 3e0)) / (8e0 * 3.1415926536 * pow(nu[i], 3e0)) * popDensThis[ui][k] * ((popDensThis[li][k] * uw) / (popDensThis[ui][k] * lw) - 1e0) * dx[k] / (1.0e0 * 1e5)) ? (0e0) : ((Aul * pow(2.997925e10, 3e0)) / (8e0 * 3.1415926536 * pow(nu[i], 3e0)) * popDensThis[ui][k] * ((popDensThis[li][k] * uw) / (popDensThis[ui][k] * lw) - 1e0) * dx[k] / (1.0e0 * 1e5)));
    }
    betatau[i] = ((0.5) < (beta_tau_finestructure(tau)) ? (0.5) : (beta_tau_finestructure(tau)));
    Pnu[i] = background;
    DFDX[li][li] -= Blu * background * betatau[i];
    DFDX[li][ui] += (Aul + Bul * background) * betatau[i];
    DFDX[ui][li] += Blu * background * betatau[i];
    DFDX[ui][ui] -= (Aul + Bul * background) * betatau[i];
  }
}
```