user: Please show me the C implementation that matches the assembly code listed below:

```asm
<gsw_enthalpy>:
  endbr64
  pxor   %xmm5,%xmm5
  sub    $0x88,%rsp
  movapd %xmm0,%xmm4
  ucomisd %xmm0,%xmm5
  ja     L32
  movapd %xmm0,%xmm6
  sqrtsd %xmm6,%xmm6
L6:
  movsd  D117(%rip),%xmm3
  movsd  D24(%rip),%xmm0
  movsd  D141(%rip),%xmm11
  movsd  D256(%rip),%xmm7
  movsd  D70(%rip),%xmm8
  mulsd  %xmm1,%xmm3
  addsd  D69(%rip),%xmm3
  mulsd  %xmm4,%xmm0
  mulsd  %xmm1,%xmm11
  subsd  D548(%rip),%xmm11
  mulsd  %xmm1,%xmm7
  mulsd  %xmm1,%xmm3
  subsd  D172(%rip),%xmm7
  addsd  D52(%rip),%xmm3
  mulsd  %xmm1,%xmm8
  subsd  D222(%rip),%xmm8
  mulsd  %xmm1,%xmm11
  addsd  D340(%rip),%xmm11
  mulsd  %xmm1,%xmm3
  subsd  D506(%rip),%xmm3
  mulsd  %xmm1,%xmm11
  mulsd  %xmm1,%xmm3
  subsd  D484(%rip),%xmm3
  addsd  %xmm0,%xmm3
  movsd  D37(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D201(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D350(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D365(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D410(%rip),%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm0,%xmm3
  movsd  D235(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D296(%rip),%xmm0
  mulsd  %xmm4,%xmm3
  mulsd  %xmm1,%xmm0
  subsd  D382(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D111(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D405(%rip),%xmm0
  addsd  %xmm0,%xmm3
  movsd  D472(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm3,0x18(%rsp)
  movsd  D544(%rip),%xmm3
  subsd  D312(%rip),%xmm0
  mulsd  %xmm1,%xmm3
  addsd  D455(%rip),%xmm3
  mulsd  %xmm4,%xmm0
  mulsd  %xmm1,%xmm3
  subsd  D77(%rip),%xmm3
  mulsd  %xmm1,%xmm3
  subsd  D528(%rip),%xmm3
  addsd  %xmm0,%xmm3
  movsd  D246(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  %xmm0,%xmm7
  mulsd  %xmm1,%xmm7
  subsd  D417(%rip),%xmm7
  subsd  D212(%rip),%xmm11
  movsd  D156(%rip),%xmm0
  movsd  D445(%rip),%xmm9
  mulsd  %xmm1,%xmm0
  subsd  D55(%rip),%xmm0
  mulsd  %xmm6,%xmm11
  movsd  D278(%rip),%xmm6
  mulsd  %xmm1,%xmm9
  subsd  D29(%rip),%xmm9
  mulsd  %xmm1,%xmm6
  mulsd  %xmm1,%xmm0
  subsd  D318(%rip),%xmm6
  subsd  D479(%rip),%xmm0
  mulsd  %xmm1,%xmm9
  subsd  D16(%rip),%xmm9
  mulsd  %xmm1,%xmm6
  addsd  %xmm0,%xmm11
  movsd  D478(%rip),%xmm0
  subsd  D19(%rip),%xmm6
  mulsd  %xmm1,%xmm0
  subsd  D198(%rip),%xmm0
  mulsd  %xmm4,%xmm11
  mulsd  %xmm1,%xmm0
  addsd  D160(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D522(%rip),%xmm0
  addsd  %xmm0,%xmm11
  movsd  D531(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D183(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  mulsd  D22(%rip),%xmm4
  addsd  %xmm0,%xmm6
  mulsd  D314(%rip),%xmm6
  movapd %xmm4,%xmm0
  movapd %xmm9,%xmm4
  movapd %xmm11,%xmm9
  addsd  %xmm0,%xmm4
  mulsd  %xmm4,%xmm9
  movapd %xmm6,%xmm10
  mulsd  %xmm6,%xmm10
  movapd %xmm10,%xmm0
  subsd  %xmm9,%xmm0
  ucomisd %xmm0,%xmm5
  ja     L20
  sqrtsd %xmm0,%xmm0
L33:
  movsd  D400(%rip),%xmm5
  movsd  %xmm7,0x60(%rsp)
  movapd %xmm8,%xmm9
  addsd  %xmm8,%xmm9
  movsd  %xmm1,0x78(%rsp)
  mulsd  %xmm8,%xmm5
  movsd  %xmm8,0x68(%rsp)
  movsd  %xmm4,0x40(%rsp)
  movsd  %xmm9,0x70(%rsp)
  movsd  %xmm6,0x38(%rsp)
  movsd  %xmm11,0x20(%rsp)
  mulsd  %xmm10,%xmm5
  movapd %xmm8,%xmm10
  movsd  %xmm2,0x28(%rsp)
  mulsd  %xmm11,%xmm10
  divsd  %xmm4,%xmm5
  subsd  %xmm10,%xmm5
  movapd %xmm7,%xmm10
  addsd  %xmm7,%xmm10
  mulsd  %xmm6,%xmm10
  subsd  %xmm10,%xmm5
  movapd %xmm6,%xmm10
  subsd  %xmm0,%xmm10
  addsd  %xmm6,%xmm0
  divsd  %xmm4,%xmm5
  movsd  %xmm0,0x50(%rsp)
  movapd %xmm0,%xmm7
  movapd %xmm6,%xmm0
  addsd  %xmm6,%xmm0
  subsd  %xmm10,%xmm7
  movsd  %xmm10,0x48(%rsp)
  movsd  %xmm7,0x8(%rsp)
  addsd  %xmm3,%xmm5
  movapd %xmm4,%xmm3
  mulsd  %xmm2,%xmm3
  movsd  %xmm5,0x58(%rsp)
  addsd  %xmm3,%xmm0
  movsd  %xmm3,0x30(%rsp)
  mulsd  %xmm2,%xmm0
  divsd  %xmm11,%xmm0
  addsd  D405(%rip),%xmm0
  call   <log@plt>
  movsd  0x30(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movsd  0x50(%rsp),%xmm12
  movsd  0x48(%rsp),%xmm10
  movsd  %xmm0,0x10(%rsp)
  mulsd  %xmm3,%xmm2
  addsd  %xmm12,%xmm3
  mulsd  %xmm10,%xmm3
  movapd %xmm2,%xmm0
  divsd  %xmm3,%xmm0
  addsd  D405(%rip),%xmm0
  call   <log@plt>
  movsd  0x70(%rsp),%xmm9
  movsd  0x20(%rsp),%xmm11
  movsd  0x38(%rsp),%xmm6
  movsd  0x40(%rsp),%xmm4
  movapd %xmm0,%xmm10
  movapd %xmm9,%xmm5
  movsd  0x60(%rsp),%xmm7
  movsd  0x68(%rsp),%xmm8
  mulsd  %xmm11,%xmm5
  movsd  0x28(%rsp),%xmm2
  movapd %xmm6,%xmm1
  mulsd  %xmm7,%xmm11
  mulsd  %xmm6,%xmm9
  mulsd  D314(%rip),%xmm8
  movapd %xmm5,%xmm3
  movsd  0x18(%rsp),%xmm5
  mulsd  %xmm6,%xmm3
  divsd  %xmm4,%xmm9
  mulsd  %xmm2,%xmm8
  divsd  %xmm4,%xmm3
  subsd  %xmm9,%xmm7
  addsd  %xmm8,%xmm7
  mulsd  %xmm2,%xmm7
  divsd  %xmm4,%xmm7
  subsd  %xmm11,%xmm3
  divsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm5
  movsd  0x8(%rsp),%xmm3
  mulsd  %xmm4,%xmm3
  movapd %xmm5,%xmm0
  movsd  0x58(%rsp),%xmm5
  mulsd  %xmm4,%xmm0
  addsd  %xmm4,%xmm4
  mulsd  %xmm5,%xmm1
  divsd  %xmm4,%xmm5
  subsd  %xmm1,%xmm0
  movsd  0x78(%rsp),%xmm1
  mulsd  D126(%rip),%xmm1
  divsd  %xmm3,%xmm0
  mulsd  0x10(%rsp),%xmm5
  addsd  %xmm7,%xmm5
  mulsd  %xmm10,%xmm0
  addsd  %xmm5,%xmm0
  mulsd  D27(%rip),%xmm0
  add    $0x88,%rsp
  addsd  %xmm1,%xmm0
  ret
L32:
  movsd  %xmm2,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm0,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  pxor   %xmm5,%xmm5
  movsd  0x8(%rsp),%xmm4
  movapd %xmm0,%xmm6
  jmp    L6
L20:
  movsd  %xmm2,0x50(%rsp)
  movsd  %xmm1,0x48(%rsp)
  movsd  %xmm10,0x40(%rsp)
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm6,0x30(%rsp)
  movsd  %xmm11,0x28(%rsp)
  movsd  %xmm8,0x20(%rsp)
  movsd  %xmm7,0x10(%rsp)
  movsd  %xmm3,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x50(%rsp),%xmm2
  movsd  0x48(%rsp),%xmm1
  movsd  0x40(%rsp),%xmm10
  movsd  0x38(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm6
  movsd  0x28(%rsp),%xmm11
  movsd  0x20(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm7
  movsd  0x8(%rsp),%xmm3
  jmp    L33
```

tool_calls: name: unpack_data
arguments: data_tag: D117
data_type: byte[8]
name: unpack_data
arguments: data_tag: D24
data_type: byte[8]
name: unpack_data
arguments: data_tag: D141
data_type: byte[8]
name: unpack_data
arguments: data_tag: D256
data_type: byte[8]
name: unpack_data
arguments: data_tag: D70
data_type: byte[8]
name: unpack_data
arguments: data_tag: D69
data_type: byte[8]
name: unpack_data
arguments: data_tag: D548
data_type: byte[8]
name: unpack_data
arguments: data_tag: D172
data_type: byte[8]
name: unpack_data
arguments: data_tag: D52
data_type: byte[8]
name: unpack_data
arguments: data_tag: D222
data_type: byte[8]
name: unpack_data
arguments: data_tag: D340
data_type: byte[8]
name: unpack_data
arguments: data_tag: D506
data_type: byte[8]
name: unpack_data
arguments: data_tag: D484
data_type: byte[8]
name: unpack_data
arguments: data_tag: D37
data_type: byte[8]
name: unpack_data
arguments: data_tag: D201
data_type: byte[8]
name: unpack_data
arguments: data_tag: D350
data_type: byte[8]
name: unpack_data
arguments: data_tag: D365
data_type: byte[8]
name: unpack_data
arguments: data_tag: D410
data_type: byte[8]
name: unpack_data
arguments: data_tag: D235
data_type: byte[8]
name: unpack_data
arguments: data_tag: D296
data_type: byte[8]
name: unpack_data
arguments: data_tag: D382
data_type: byte[8]
name: unpack_data
arguments: data_tag: D111
data_type: byte[8]
name: unpack_data
arguments: data_tag: D405
data_type: f64
name: unpack_data
arguments: data_tag: D472
data_type: byte[8]
name: unpack_data
arguments: data_tag: D544
data_type: byte[8]
name: unpack_data
arguments: data_tag: D312
data_type: byte[8]
name: unpack_data
arguments: data_tag: D455
data_type: byte[8]
name: unpack_data
arguments: data_tag: D77
data_type: byte[8]
name: unpack_data
arguments: data_tag: D528
data_type: byte[8]
name: unpack_data
arguments: data_tag: D246
data_type: byte[8]
name: unpack_data
arguments: data_tag: D417
data_type: byte[8]
name: unpack_data
arguments: data_tag: D212
data_type: byte[8]
name: unpack_data
arguments: data_tag: D156
data_type: byte[8]
name: unpack_data
arguments: data_tag: D445
data_type: byte[8]
name: unpack_data
arguments: data_tag: D55
data_type: byte[8]
name: unpack_data
arguments: data_tag: D278
data_type: byte[8]
name: unpack_data
arguments: data_tag: D29
data_type: byte[8]
name: unpack_data
arguments: data_tag: D318
data_type: byte[8]
name: unpack_data
arguments: data_tag: D479
data_type: f64
name: unpack_data
arguments: data_tag: D16
data_type: byte[8]
name: unpack_data
arguments: data_tag: D478
data_type: byte[8]
name: unpack_data
arguments: data_tag: D19
data_type: byte[8]
name: unpack_data
arguments: data_tag: D198
data_type: byte[8]
name: unpack_data
arguments: data_tag: D160
data_type: f64
name: unpack_data
arguments: data_tag: D522
data_type: byte[8]
name: unpack_data
arguments: data_tag: D531
data_type: byte[8]
name: unpack_data
arguments: data_tag: D183
data_type: byte[8]
name: unpack_data
arguments: data_tag: D22
data_type: byte[8]
name: unpack_data
arguments: data_tag: D314
data_type: f64
name: unpack_data
arguments: data_tag: D400
data_type: f64
name: unpack_data
arguments: data_tag: D126
data_type: f64
name: unpack_data
arguments: data_tag: D27
data_type: f64

tool: {"D117": [0x25,0x50,0x24,0x19,0xed,0xe9,0xb3,0xbd]}

tool: {"D24": [0xf6,0x68,0xb1,0x23,0xf5,0xba,0xd6,0x3e]}

tool: {"D141": [0x73,0x33,0x84,0xef,0x4f,0x14,0x86,0x3e]}

tool: {"D256": [0x10,0xf7,0x4c,0x08,0xff,0xfe,0x14,0xbd]}

tool: {"D70": [0x2c,0x35,0x32,0x6d,0xf1,0x75,0x91,0x3c]}

tool: {"D69": [0xda,0x54,0xdf,0x5d,0x67,0x5e,0x10,0x3e]}

tool: {"D548": [0xd6,0xa6,0x84,0xcc,0x10,0xed,0xed,0x3e]}

tool: {"D172": [0x61,0xf3,0x4c,0xa2,0x49,0x59,0xb5,0x3d]}

tool: {"D52": [0x9b,0x60,0x15,0x69,0x6b,0xf3,0x80,0x3e]}

tool: {"D222": [0xeb,0xc5,0x77,0x63,0x5d,0xa0,0xd5,0x3c]}

tool: {"D340": [0xbf,0xb2,0x6b,0xaa,0x59,0xb2,0x40,0x3f]}

tool: {"D506": [0xbc,0xe6,0xbd,0x4a,0x61,0xfc,0xfc,0x3e]}

tool: {"D484": [0x0f,0x4d,0x4a,0x59,0xcf,0xce,0x7e,0x3f]}

tool: {"D37": [0x8b,0x87,0x54,0x54,0x4b,0x95,0x87,0x3d]}

tool: {"D201": [0x61,0xf7,0x2f,0x26,0x08,0x7a,0xc9,0x3d]}

tool: {"D350": [0x57,0x27,0x46,0x6b,0x35,0x7a,0x40,0x3e]}

tool: {"D365": [0x2d,0xaf,0xa5,0xc9,0x09,0x83,0x99,0x3e]}

tool: {"D410": [0xe4,0x53,0x23,0x49,0xa0,0x51,0x01,0x3f]}

tool: {"D235": [0x40,0x56,0x29,0xd6,0xcf,0x2b,0x07,0x3e]}

tool: {"D296": [0x42,0xd1,0x39,0xf4,0x45,0xc8,0xb2,0x3e]}

tool: {"D382": [0x63,0xc6,0xe2,0xf0,0x2d,0xa3,0xf8,0x3e]}

tool: {"D111": [0xb6,0xcd,0x91,0xdc,0x8a,0xbd,0x66,0x3f]}

tool: D405: 1

tool: {"D472": [0xe8,0x94,0x78,0xf4,0xa0,0x57,0x05,0x3e]}

tool: {"D544": [0xb7,0x34,0xcf,0x2a,0x41,0xfc,0x80,0xbd]}

tool: {"D312": [0x8f,0x77,0xaa,0x37,0x2c,0xaa,0x7d,0x3e]}

tool: {"D455": [0x58,0xd0,0xf2,0x4f,0xfe,0x69,0x2b,0x3e]}

tool: {"D77": [0xa4,0xa8,0x06,0xd8,0x74,0x8f,0x95,0x3e]}

tool: {"D528": [0x9b,0xcc,0x5a,0x74,0xff,0xc0,0xfc,0x3e]}

tool: {"D246": [0x3e,0xd3,0xdc,0x5f,0x9b,0x35,0x09,0x3d]}

tool: {"D417": [0x5a,0x66,0xbb,0x13,0x58,0xc2,0xde,0x3d]}

tool: {"D212": [0xfc,0x2d,0x58,0xb1,0x3a,0x6c,0xa4,0x3f]}

tool: {"D156": [0xfc,0x2c,0x45,0x01,0x00,0x83,0x2e,0x3f]}

tool: {"D445": [0x18,0xe6,0x61,0x34,0xd2,0x57,0xa7,0xbd]}

tool: {"D55": [0xf2,0xa0,0x4f,0x84,0xed,0x94,0x9e,0x3f]}

tool: {"D278": [0x6e,0xb3,0xdd,0x0b,0xa2,0x41,0xcf,0x3e]}

tool: {"D29": [0x40,0x72,0xfa,0x98,0x94,0x59,0x55,0x3e]}

tool: {"D318": [0xd1,0x33,0xb0,0x73,0xcd,0x84,0x36,0x3f]}

tool: D479: 6.698001071123802

tool: {"D16": [0xb4,0x8a,0x75,0x40,0x7b,0xb7,0x94,0x3e]}

tool: {"D478": [0x1f,0x14,0xa5,0x0a,0xda,0x5c,0x53,0x3f]}

tool: {"D19": [0xb8,0x52,0x81,0xe0,0x61,0xde,0x96,0x3f]}

tool: {"D198": [0x8e,0x15,0x07,0xbe,0xd4,0x1d,0xa0,0x3f]}

tool: D160: 2.839940833161907

tool: {"D522": [0x2f,0x19,0x8d,0x99,0xbc,0x3e,0x8f,0x40]}

tool: {"D531": [0x0d,0x0d,0x3e,0x7f,0x28,0x13,0xa7,0x3e]}

tool: {"D183": [0xea,0x40,0xd6,0x96,0x94,0xae,0x27,0x3f]}

tool: {"D22": [0x3d,0xae,0x24,0x03,0xdc,0x50,0xe5,0x3d]}

tool: D314: 0.5

tool: D400: 4

tool: D126: 3991.86795711963

tool: D27: 10000

assistant: ```c
double gsw_enthalpy(double sa, double ct, double p) {
  double v01 = 9.998420897506056e+2, v02 = 2.839940833161907e0, v03 = -3.147759265588511e-2, v04 = 1.181805545074306e-3, v05 = -6.698001071123802e0, v06 = -2.986498947203215e-2, v07 = 2.327859407479162e-4, v08 = -3.988822378968490e-2, v09 = 5.095422573880500e-4, v10 = -1.426984671633621e-5, v11 = 1.645039373682922e-7, v12 = -2.233269627352527e-2, v13 = -3.436090079851880e-4, v14 = 3.726050720345733e-6, v15 = -1.806789763745328e-4, v16 = 6.876837219536232e-7, v17 = -3.087032500374211e-7, v18 = -1.988366587925593e-8, v19 = -1.061519070296458e-11, v20 = 1.550932729220080e-10, v21 = 1.0e0, v22 = 2.775927747785646e-3, v23 = -2.349607444135925e-5, v24 = 1.119513357486743e-6, v25 = 6.743689325042773e-10, v26 = -7.521448093615448e-3, v27 = -2.764306979894411e-5, v28 = 1.262937315098546e-7, v29 = 9.527875081696435e-10, v30 = -1.811147201949891e-11, v31 = -3.303308871386421e-5, v32 = 3.801564588876298e-7, v33 = -7.672876869259043e-9, v34 = -4.634182341116144e-11, v35 = 2.681097235569143e-12, v36 = 5.419326551148740e-6, v37 = -2.742185394906099e-5, v38 = -3.212746477974189e-7, v39 = 3.191413910561627e-9, v40 = -1.931012931541776e-12, v41 = -1.105097577149576e-7, v42 = 6.211426728363857e-10, v43 = -1.119011592875110e-10, v44 = -1.941660213148725e-11, v45 = -1.864826425365600e-14, v46 = 1.119522344879478e-14, v47 = -1.200507748551599e-15, v48 = 6.057902487546866e-17;
  double db2pa, cp0, sqrtsa, a0, a1, a2, a3, b0, b1, b2, b1sq;
  double sqrt_disc, ca, cb, cn, cm, part;
  db2pa = 1e4;
  cp0 = 3991.86795711963e0;
  sqrtsa = sqrt(sa);
  a0 = v21 + ct * (v22 + ct * (v23 + ct * (v24 + v25 * ct))) + sa * (v26 + ct * (v27 + ct * (v28 + ct * (v29 + v30 * ct))) + v36 * sa + sqrtsa * (v31 + ct * (v32 + ct * (v33 + ct * (v34 + v35 * ct)))));
  a1 = v37 + ct * (v38 + ct * (v39 + v40 * ct)) + sa * (v41 + v42 * ct);
  a2 = v43 + ct * (v44 + v45 * ct + v46 * sa);
  a3 = v47 + v48 * ct;
  b0 = v01 + ct * (v02 + ct * (v03 + v04 * ct)) + sa * (v05 + ct * (v06 + v07 * ct) + sqrtsa * (v08 + ct * (v09 + ct * (v10 + v11 * ct))));
  b1 = 0.5e0 * (v12 + ct * (v13 + v14 * ct) + sa * (v15 + v16 * ct));
  b2 = v17 + ct * (v18 + v19 * ct) + v20 * sa;
  b1sq = b1 * b1;
  sqrt_disc = sqrt(b1sq - b0 * b2);
  cn = a0 + (2.0 * a3 * b0 * b1 / b2 - a2 * b0) / b2;
  cm = a1 + (4.0 * a3 * b1sq / b2 - a3 * b0 - 2 * a2 * b1) / b2;
  ca = b1 - sqrt_disc;
  cb = b1 + sqrt_disc;
  part = (cn * b2 - cm * b1) / (b2 * (cb - ca));
  return (cp0 * ct + db2pa * (p * (a2 - 2.0 * a3 * b1 / b2 + 0.5 * a3 * p) / b2 + (cm / (2.0 * b2)) * log(1 + p * (2.0 * b1 + b2 * p) / b0) + part * log(1.0 + (b2 * p * (cb - ca)) / (ca * (cb + b2 * p)))));
}
```