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

```asm
<dcs_photonuclear_f2a_BM>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xe0,%rsp
  movsd  %xmm0,-0xb8(%rbp)
  movsd  %xmm1,-0xc0(%rbp)
  movsd  %xmm2,-0xc8(%rbp)
  movsd  %xmm3,-0xd0(%rbp)
  movsd  %xmm4,-0xd8(%rbp)
  movsd  -0xd8(%rbp),%xmm1
  movsd  D43(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd8(%rbp)
  movsd  D176(%rip),%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  D193(%rip),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  D213(%rip),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  D116(%rip),%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  D185(%rip),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  D67(%rip),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  D42(%rip),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  D100(%rip),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  D12(%rip),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  D10(%rip),%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0xd8(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  -0x98(%rbp),%xmm2
  movsd  -0xd8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D16(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  D41(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0xd8(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  -0xd8(%rbp),%xmm1
  movapd %xmm1,%xmm2
  addsd  -0x90(%rbp),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D16(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x88(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movq   D5(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0xd0(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm5
  mulsd  -0x70(%rbp),%xmm5
  movsd  %xmm5,-0xe0(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  D16(%rip),%xmm0
  movapd %xmm1,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0xd8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0xa8(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movapd %xmm2,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  mulsd  -0xe0(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x58(%rbp),%xmm1
  movsd  D109(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  D16(%rip),%xmm0
  subsd  -0xd0(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <pow@plt>
  movsd  -0x48(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x78(%rbp),%xmm1
  movsd  D16(%rip),%xmm0
  subsd  -0xd0(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <pow@plt>
  movsd  -0x48(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  D16(%rip),%xmm0
  subsd  -0x80(%rbp),%xmm0
  mov    -0xd0(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movsd  %xmm0,-0xe0(%rbp)
  movsd  D16(%rip),%xmm0
  movapd %xmm0,%xmm6
  subsd  -0xd0(%rbp),%xmm6
  movq   %xmm6,%rax
  movsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movapd %xmm0,%xmm7
  mulsd  -0xe0(%rbp),%xmm7
  movsd  %xmm7,-0xe0(%rbp)
  movsd  -0xd8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0xa0(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movsd  -0x80(%rbp),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  mulsd  -0xe0(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x68(%rbp),%xmm0
  mulsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x60(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x48(%rbp),%xmm1
  movsd  D16(%rip),%xmm0
  subsd  -0xd0(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x30(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movsd  D109(%rip),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  -0x28(%rbp),%xmm2
  movsd  D109(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x40(%rbp),%xmm0
  addsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x38(%rbp),%xmm0
  addsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0xb8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  -0xc0(%rbp),%xmm0
  subsd  -0xb8(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0xc8(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: parse_data
arguments: label: D43
data_type: byte[8]
name: parse_data
arguments: label: D176
data_type: byte[8]
name: parse_data
arguments: label: D193
data_type: byte[8]
name: parse_data
arguments: label: D213
data_type: byte[8]
name: parse_data
arguments: label: D116
data_type: byte[8]
name: parse_data
arguments: label: D185
data_type: f64
name: parse_data
arguments: label: D67
data_type: f64
name: parse_data
arguments: label: D42
data_type: f64
name: parse_data
arguments: label: D100
data_type: f64
name: parse_data
arguments: label: D12
data_type: f64
name: parse_data
arguments: label: D10
data_type: f64
name: parse_data
arguments: label: D16
data_type: f64
name: parse_data
arguments: label: D41
data_type: f64
name: parse_data
arguments: label: D5
data_type: byte[8]
name: parse_data
arguments: label: D109
data_type: f64

tool: {"D43": [0x00,0x00,0x00,0x00,0x80,0x84,0x2e,0x41]}

tool: {"D176": [0x00,0x00,0x00,0x00,0x20,0xad,0x0e,0x41]}

tool: {"D193": [0x00,0x00,0x00,0x00,0xd0,0xe0,0x22,0x41]}

tool: {"D213": [0x00,0x00,0x00,0x00,0x68,0x1c,0x47,0x41]}

tool: {"D116": [0x00,0x00,0x00,0x00,0xe4,0x9b,0x36,0x41]}

tool: D185: 0.0988

tool: D67: 0.4056

tool: D42: 1.8152

tool: D100: 0.12

tool: D12: 1.2437

tool: D10: 0.1853

tool: D16: 1

tool: D41: 1.5

tool: {"D5": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

tool: D109: 4

assistant: ```c
double dcs_photonuclear_f2a_BM(double Z, double A, double shadowing, double x, double Q2) {
  Q2 *= 1E+06;
  const double a = 0.2513E+06;
  const double b = 0.6186E+06;
  const double c = 3.0292E+06;
  const double d = 1.4817E+06;
  const double intercept_pomeron = 0.0988;
  const double intercept_reggeon = 0.4056;
  const double tau = 1.8152;
  const double A_s = 0.12;
  const double B_up = 1.2437;
  const double B_down = 0.1853;
  const double n = 1.5 * (1. + Q2 / (Q2 + c));
  const double dl = intercept_pomeron * (1. + 2. * Q2 / (Q2 + d));
  double aux = A_s * pow(x, -dl) * pow(Q2 / (Q2 + a), 1. + dl);
  const double F_proton_singlet = aux * pow(1. - x, n + 4.);
  const double F_neutron_singlet = aux * pow(1. - x, n + tau);
  aux = pow(x, 1. - intercept_reggeon) * pow(1. - x, n) * pow(Q2 / (Q2 + b), intercept_reggeon);
  const double Up_valence = B_up * aux;
  const double Down_valence = B_down * aux * (1. - x);
  const double F_proton_non_singlet = Up_valence + Down_valence;
  const double F_neutron_non_singlet = Up_valence / 4. + Down_valence * 4.;
  const double structure_function_proton = F_proton_singlet + F_proton_non_singlet;
  const double structure_function_neutron = F_neutron_singlet + F_neutron_non_singlet;
  return shadowing * (Z * structure_function_proton + (A - Z) * structure_function_neutron);
}
```