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

```asm
<gsw_enthalpy_sso_0_p>:
  endbr64
  sub    $0x18,%rsp
  movsd  D68(%rip),%xmm1
  movsd  %xmm0,(%rsp)
  mulsd  %xmm0,%xmm1
  movapd %xmm1,%xmm2
  subsd  D19(%rip),%xmm2
  movsd  %xmm1,0x8(%rsp)
  mulsd  %xmm0,%xmm2
  divsd  D133(%rip),%xmm2
  addsd  D20(%rip),%xmm2
  cvttsd2si %xmm2,%edi
  call   <log@plt>
  movsd  (%rsp),%xmm3
  movsd  D63(%rip),%xmm2
  mulsd  D17(%rip),%xmm0
  movsd  0x8(%rsp),%xmm1
  mulsd  %xmm3,%xmm2
  addsd  D90(%rip),%xmm2
  mulsd  %xmm3,%xmm2
  divsd  D68(%rip),%xmm2
  addsd  %xmm0,%xmm2
  movsd  D104(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D116(%rip),%xmm1
  mulsd  D105(%rip),%xmm1
  movsd  %xmm2,(%rsp)
  divsd  %xmm1,%xmm0
  addsd  D20(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <log@plt>
  mulsd  D120(%rip),%xmm0
  addsd  (%rsp),%xmm0
  mulsd  D52(%rip),%xmm0
  add    $0x18,%rsp
  ret
```

tool_calls: name: read_data
arguments: data_tag: D68
type: byte[8]
name: read_data
arguments: data_tag: D19
type: byte[8]
name: read_data
arguments: data_tag: D133
type: byte[8]
name: read_data
arguments: data_tag: D20
type: f64
name: read_data
arguments: data_tag: D63
type: byte[8]
name: read_data
arguments: data_tag: D17
type: byte[8]
name: read_data
arguments: data_tag: D90
type: byte[8]
name: read_data
arguments: data_tag: D104
type: byte[8]
name: read_data
arguments: data_tag: D116
type: byte[8]
name: read_data
arguments: data_tag: D105
type: byte[8]
name: read_data
arguments: data_tag: D120
type: byte[8]
name: read_data
arguments: data_tag: D52
type: f64

tool: {"D68": [0x01,0xf1,0xb3,0xea,0xc8,0x59,0x94,0xbe]}

tool: {"D19": [0x7a,0xaa,0xb4,0xa4,0xef,0x5f,0x9d,0x3f]}

tool: {"D133": [0xc1,0x96,0x4d,0xfa,0xe8,0x9f,0x87,0x40]}

tool: D20: 1

tool: {"D63": [0xeb,0xc5,0x77,0x63,0x5d,0xa0,0xc5,0xbc]}

tool: {"D17": [0xb0,0xff,0xfb,0x75,0x46,0x68,0x4c,0x40]}

tool: {"D90": [0x80,0xe4,0x37,0xf6,0x18,0x3f,0x7d,0x3d]}

tool: {"D104": [0xdb,0xfa,0x1a,0x35,0x50,0x5b,0xa5,0x3f]}

tool: {"D116": [0x78,0x96,0x02,0x8b,0x61,0xad,0x7a,0x3f]}

tool: {"D105": [0x0c,0xa8,0xba,0x03,0xa4,0x05,0xa2,0xbf]}

tool: {"D120": [0x81,0x6f,0xe4,0xfb,0x6c,0x66,0x4c,0x40]}

tool: D52: 10000

assistant: ```c
double gsw_enthalpy_sso_0_p(double p) {
  double v01 = 9.998420897506056e2, v05 = -6.698001071123802e0, v08 = -3.988822378968490e-2, v12 = -2.233269627352527e-2, v15 = -1.806789763745328e-4, v17 = -3.087032500374211e-7, v20 = 1.550932729220080e-10, v21 = 1.0e0, v26 = -7.521448093615448e-3, v31 = -3.303308871386421e-5, v36 = 5.419326551148740e-6, v37 = -2.742185394906099e-5, v41 = -1.105097577149576e-7, v43 = -1.119011592875110e-10, v47 = -1.200507748551599e-15, db2pa = 1e4, sso = 35.16504e0, sqrtsso = 5.930011804372737e0;
  double a0, a1, a2, a3, b0, b1, b2, b1sq, sqrt_disc, n, m, a, b, part;
  a0 = v21 + sso * (v26 + v36 * sso + v31 * sqrtsso);
  a1 = v37 + v41 * sso;
  a2 = v43;
  a3 = v47;
  b0 = v01 + sso * (v05 + v08 * sqrtsso);
  b1 = 0.5 * (v12 + v15 * sso);
  b2 = v17 + v20 * sso;
  b1sq = b1 * b1;
  sqrt_disc = sqrt(b1sq - b0 * b2);
  n = a0 + (2e0 * a3 * b0 * b1 / b2 - a2 * b0) / b2;
  m = a1 + (4e0 * a3 * b1sq / b2 - a3 * b0 - 2 * a2 * b1) / b2;
  a = b1 - sqrt_disc;
  b = b1 + sqrt_disc;
  part = (n * b2 - m * b1) / (b2 * (b - a));
  return (db2pa * (p * (a2 - 2e0 * a3 * b1 / b2 + 0.5e0 * a3 * p) / b2 + (m / (2e0 * b2)) * log(1e0 + p * (2e0 * b1 + b2 * p) / b0) + part * log(1e0 + (b2 * p * (b - a)) / (a * (b + b2 * p)))));
}
```