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

```asm
<gsw_pt_from_ct>:
  endbr64
  sub    $0x38,%rsp
  movapd %xmm1,%xmm4
  movsd  D146(%rip),%xmm3
  movsd  D144(%rip),%xmm8
  movsd  D42(%rip),%xmm7
  movapd %xmm0,%xmm5
  movsd  D114(%rip),%xmm1
  movsd  %xmm4,0x10(%rsp)
  mulsd  %xmm4,%xmm8
  movsd  D31(%rip),%xmm6
  movsd  %xmm5,0x18(%rsp)
  mulsd  %xmm0,%xmm7
  movsd  D148(%rip),%xmm0
  divsd  D106(%rip),%xmm7
  mulsd  %xmm4,%xmm1
  mulsd  %xmm7,%xmm0
  addsd  D69(%rip),%xmm0
  mulsd  %xmm7,%xmm6
  addsd  %xmm8,%xmm3
  addsd  D53(%rip),%xmm6
  movapd %xmm3,%xmm2
  addsd  %xmm1,%xmm6
  addsd  %xmm8,%xmm3
  mulsd  %xmm4,%xmm2
  addsd  %xmm6,%xmm1
  addsd  %xmm2,%xmm0
  movsd  D75(%rip),%xmm2
  mulsd  %xmm7,%xmm2
  subsd  D10(%rip),%xmm2
  mulsd  %xmm7,%xmm2
  movapd %xmm4,%xmm7
  mulsd  %xmm6,%xmm7
  subsd  D139(%rip),%xmm2
  addsd  %xmm7,%xmm2
  divsd  %xmm0,%xmm2
  mulsd  %xmm2,%xmm3
  movsd  %xmm2,0x20(%rsp)
  subsd  %xmm3,%xmm1
  divsd  %xmm1,%xmm0
  movapd %xmm2,%xmm1
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm5,%xmm0
  call   <gsw_ct_from_pt@plt>
  movsd  0x10(%rsp),%xmm4
  movsd  0x20(%rsp),%xmm2
  movapd %xmm0,%xmm6
  movsd  0x18(%rsp),%xmm5
  subsd  %xmm4,%xmm6
  movapd %xmm2,%xmm1
  movsd  %xmm4,0x28(%rsp)
  movapd %xmm5,%xmm0
  movapd %xmm6,%xmm4
  divsd  0x8(%rsp),%xmm4
  subsd  %xmm4,%xmm1
  movsd  %xmm6,0x10(%rsp)
  addsd  %xmm2,%xmm1
  mulsd  D9(%rip),%xmm1
  movsd  %xmm1,0x8(%rsp)
  call   <gsw_gibbs_pt0_pt0@plt>
  movsd  0x8(%rsp),%xmm1
  addsd  D104(%rip),%xmm1
  movsd  0x10(%rsp),%xmm6
  movsd  0x20(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm5
  movapd %xmm1,%xmm3
  xorpd  D97(%rip),%xmm3
  mulsd  %xmm0,%xmm3
  movapd %xmm5,%xmm0
  divsd  D151(%rip),%xmm3
  divsd  %xmm3,%xmm6
  movsd  %xmm3,0x10(%rsp)
  subsd  %xmm6,%xmm2
  movapd %xmm2,%xmm1
  movsd  %xmm2,0x8(%rsp)
  call   <gsw_ct_from_pt@plt>
  movsd  0x28(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm1
  add    $0x38,%rsp
  subsd  %xmm4,%xmm0
  divsd  %xmm3,%xmm0
  subsd  %xmm0,%xmm1
  movapd %xmm1,%xmm0
  ret
```

tool_calls: name: disassemble_data
arguments: data_tag: D146
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D144
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D42
type_infered: f64
name: disassemble_data
arguments: data_tag: D114
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D31
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D148
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D106
type_infered: f64
name: disassemble_data
arguments: data_tag: D69
type_infered: f64
name: disassemble_data
arguments: data_tag: D53
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D75
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D10
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D139
type_infered: byte[8]
name: disassemble_data
arguments: data_tag: D9
type_infered: f64
name: disassemble_data
arguments: data_tag: D104
type_infered: f64
name: disassemble_data
arguments: data_tag: D97
type_infered: f32[4]
name: disassemble_data
arguments: data_tag: D151
type_infered: f64

tool: {"D146": [0xaa,0x4d,0x91,0x02,0xb3,0x60,0x6f,0x3f]}

tool: {"D144": [0x1d,0x5a,0x30,0x8a,0x4f,0xef,0xb4,0x3e]}

tool: D42: 35

tool: {"D114": [0xb4,0xd3,0xbf,0x68,0xaa,0x5d,0x6f,0x3f]}

tool: {"D31": [0x0e,0xd5,0xe5,0x31,0xac,0xbf,0x61,0x3f]}

tool: {"D148": [0xeb,0x0a,0x81,0xb7,0xb5,0x51,0x45,0x3f]}

tool: D106: 35.16504

tool: D69: 1

tool: {"D53": [0xe9,0xbd,0x5e,0xca,0x05,0x54,0xee,0x3f]}

tool: {"D75": [0xae,0x71,0x0b,0xab,0xc0,0xd9,0x1b,0x3f]}

tool: {"D10": [0x27,0xa7,0x4a,0x41,0xbc,0x13,0x6b,0x3f]}

tool: {"D139": [0x32,0x4f,0xa6,0xa9,0x46,0x9d,0x8d,0x3f]}

tool: D9: 0.5

tool: D104: 273.15

tool: {"D97": [0.0,-0.0,0.0,0.0]}

tool: D151: 3991.86795711963

assistant: ```c
double gsw_pt_from_ct(double sa, double ct) {
  int n0, n2;
  double s1, p0, cp0;
  double a0, a1, a2, a3, a4, a5, b0, b1, b2, b3;
  double a5ct, b3ct, ct_factor, pt_num, pt_den, ct_diff;
  double pt, pt_old, ptm, dct_dpt;
  cp0 = 3991.86795711963e0;
  n0 = 0;
  n2 = 2;
  s1 = sa * 35.e0 / 35.16504e0;
  p0 = 0.e0;
  a0 = -1.446013646344788e-2;
  a1 = -3.305308995852924e-3;
  a2 = 1.062415929128982e-4;
  a3 = 9.477566673794488e-1;
  a4 = 2.166591947736613e-3;
  a5 = 3.828842955039902e-3;
  b0 = 1.000000000000000e0;
  b1 = 6.506097115635800e-4;
  b2 = 3.830289486850898e-3;
  b3 = 1.247811760368034e-6;
  a5ct = a5 * ct;
  b3ct = b3 * ct;
  ;
  ct_factor = (a3 + a4 * s1 + a5ct);
  pt_num = a0 + s1 * (a1 + a2 * s1) + ct * ct_factor;
  pt_den = b0 + b1 * s1 + ct * (b2 + b3ct);
  pt = (pt_num) / (pt_den);
  dct_dpt = (pt_den) / (ct_factor + a5ct - (b2 + b3ct + b3ct) * pt);
  ct_diff = gsw_ct_from_pt(sa, pt) - ct;
  pt_old = pt;
  pt = pt_old - (ct_diff) / dct_dpt;
  ptm = 0.5e0 * (pt + pt_old);
  dct_dpt = -(ptm + 273.15e0) * gsw_gibbs_pt0_pt0(sa, ptm) / cp0;
  pt = pt_old - (ct_diff) / dct_dpt;
  ct_diff = gsw_ct_from_pt(sa, pt) - ct;
  pt_old = pt;
  return (pt_old - (ct_diff) / dct_dpt);
}
```