user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<CE_pHM_P>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x40,%rsp
  movsd  %xmm0,-0x18(%rbp)
  movsd  %xmm1,-0x20(%rbp)
  movsd  %xmm2,-0x28(%rbp)
  movsd  %xmm3,-0x30(%rbp)
  movsd  %xmm4,-0x38(%rbp)
  movsd  %xmm5,-0x40(%rbp)
  movsd  D9(%rip),%xmm0
  subsd  -0x18(%rbp),%xmm0
  cvttsd2si %xmm0,%edx
  mov    D30(%rip),%rax
  mov    %edx,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x40(%rbp),%rax
  mov    $0x2,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movsd  -0x28(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0x40(%rbp),%rax
  mov    $0x4,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movsd  -0x30(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x40(%rbp),%rax
  mov    $0x2,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movq   D17(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mulsd  -0x28(%rbp),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  movsd  -0x28(%rbp),%xmm1
  mulsd  -0x30(%rbp),%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x20(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  -0x38(%rbp),%xmm0
  subsd  %xmm0,%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x38(%rbp),%xmm1
  mulsd  -0x28(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  mulsd  -0x28(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movq   D17(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  -0x30(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movapd %xmm0,%xmm2
  addsd  -0x10(%rbp),%xmm2
  divsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  divsd  -0x10(%rbp),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_label: D9
datatype: byte[8]
name: parse_data
arguments: data_label: D30
datatype: byte[8]
name: parse_data
arguments: data_label: D17
datatype: byte[8]

tool: {"D9": [0x00,0x00,0x00,0x00,0x00,0x00,0x08,0x40]}

tool: {"D30": [0x00,0x00,0x00,0x00,0x00,0x00,0x24,0x40]}

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

assistant: ```c
double CE_pHM_P(double pH, double M, double K1, double K2, double Kw, double f) {
  double H3O, t0;
  H3O = pow(10, (3 - pH)) / f;
  K1 = K1 / pow(f, 2);
  K2 = K2 / pow(f, 4);
  Kw = Kw / pow(f, 2);
  t0 = -(-Kw * K1 * K2 - K1 * K2 * H3O * M + K1 * K2 * H3O * H3O - H3O * H3O * Kw + H3O * H3O * H3O * M + H3O * H3O * H3O * H3O - Kw * K1 * H3O + H3O * H3O * H3O * K1) / (2.0 * K2 + H3O) / H3O / K1;
  return t0;
}
```