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

```asm
<gsw_pt0_from_t>:
  endbr64
  push   %rbx
  movapd %xmm0,%xmm4
  mov    $0x2,%ebx
  sub    $0x30,%rsp
  movsd  D50(%rip),%xmm5
  movsd  D31(%rip),%xmm6
  mulsd  D8(%rip),%xmm0
  movsd  D4(%rip),%xmm3
  movsd  %xmm4,0x10(%rsp)
  mulsd  %xmm1,%xmm6
  movsd  D111(%rip),%xmm7
  mulsd  %xmm2,%xmm7
  divsd  %xmm5,%xmm0
  mulsd  %xmm0,%xmm3
  subsd  D120(%rip),%xmm3
  mulsd  D119(%rip),%xmm0
  addsd  %xmm6,%xmm3
  movsd  D113(%rip),%xmm6
  mulsd  %xmm2,%xmm6
  addsd  %xmm6,%xmm3
  movsd  D139(%rip),%xmm6
  subsd  %xmm0,%xmm6
  movsd  D23(%rip),%xmm0
  mulsd  %xmm1,%xmm3
  subsd  %xmm7,%xmm6
  addsd  %xmm6,%xmm3
  movapd %xmm4,%xmm6
  divsd  %xmm5,%xmm6
  movapd %xmm0,%xmm5
  mulsd  %xmm2,%xmm3
  addsd  %xmm1,%xmm3
  movsd  %xmm3,0x18(%rsp)
  subsd  %xmm6,%xmm5
  mulsd  D108(%rip),%xmm5
  subsd  %xmm5,%xmm0
  movapd %xmm3,%xmm5
  addsd  D18(%rip),%xmm5
  mulsd  %xmm5,%xmm0
  movsd  D138(%rip),%xmm5
  divsd  %xmm0,%xmm5
  movapd %xmm4,%xmm0
  movsd  %xmm5,0x8(%rsp)
  call   <gsw_entropy_part@plt>
  movsd  0x10(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm3
  movsd  %xmm0,0x28(%rsp)
L10:
  movapd %xmm3,%xmm1
  movapd %xmm4,%xmm0
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm4,0x10(%rsp)
  call   <gsw_entropy_part_zerop@plt>
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm4
  movapd %xmm0,%xmm2
  subsd  0x28(%rsp),%xmm2
  movapd %xmm3,%xmm1
  movapd %xmm2,%xmm0
  divsd  0x8(%rsp),%xmm0
  subsd  %xmm0,%xmm1
  movapd %xmm4,%xmm0
  movsd  %xmm2,0x20(%rsp)
  addsd  %xmm3,%xmm1
  mulsd  D135(%rip),%xmm1
  call   <gsw_gibbs_pt0_pt0@plt>
  movq   D17(%rip),%xmm5
  movsd  0x20(%rsp),%xmm2
  cmp    $0x1,%ebx
  movsd  0x18(%rsp),%xmm3
  movapd %xmm0,%xmm7
  movsd  0x10(%rsp),%xmm4
  xorpd  %xmm5,%xmm2
  xorpd  %xmm5,%xmm7
  divsd  %xmm0,%xmm2
  movsd  %xmm7,0x8(%rsp)
  subsd  %xmm2,%xmm3
  jne    L5
  add    $0x30,%rsp
  movapd %xmm3,%xmm0
  pop    %rbx
  ret
  nopw   0x0(%rax,%rax,1)
L5:
  mov    $0x1,%ebx
  jmp    L10
```

tool_calls: name: retrieve_data
arguments: tag: D50
type: f64
name: retrieve_data
arguments: tag: D31
type: byte[8]
name: retrieve_data
arguments: tag: D8
type: f64
name: retrieve_data
arguments: tag: D4
type: byte[8]
name: retrieve_data
arguments: tag: D111
type: byte[8]
name: retrieve_data
arguments: tag: D120
type: byte[8]
name: retrieve_data
arguments: tag: D119
type: byte[8]
name: retrieve_data
arguments: tag: D113
type: byte[8]
name: retrieve_data
arguments: tag: D139
type: byte[8]
name: retrieve_data
arguments: tag: D23
type: f64
name: retrieve_data
arguments: tag: D108
type: f64
name: retrieve_data
arguments: tag: D18
type: f64
name: retrieve_data
arguments: tag: D138
type: f64
name: retrieve_data
arguments: tag: D135
type: f64
name: retrieve_data
arguments: tag: D17
type: byte[8]

tool: D50: 35.16504

tool: {"D31": [0x88,0xfb,0x99,0x2c,0x70,0x17,0x53,0x3e]}

tool: D8: 35

tool: {"D4": [0xf8,0x16,0x2f,0x75,0xb3,0x7c,0x5e,0x3e]}

tool: {"D111": [0x76,0x5c,0x53,0xa9,0x8b,0xb5,0x3f,0x3e]}

tool: {"D120": [0xe0,0x35,0x10,0x26,0x46,0x94,0xe1,0x3e]}

tool: {"D119": [0x55,0x82,0xb8,0x95,0x3b,0xc3,0xb7,0x3e]}

tool: {"D113": [0x4e,0xf8,0x09,0xc4,0x00,0x86,0xe7,0x3d]}

tool: {"D139": [0xe7,0x41,0x75,0x08,0x88,0x26,0xe2,0x3e]}

tool: D23: 1

tool: D108: 0.05

tool: D18: 273.15

tool: D138: 3991.86795711963

tool: D135: 0.5

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

assistant: ```c
double gsw_pt0_from_t(double sa, double t, double p) {
  int n0, n2, no_iter;
  double pt0, pt0_old, dentropy, dentropy_dt, sso, cp0;
  double s1, true_entropy_part, pt0m;
  n0 = 0;
  n2 = 2;
  cp0 = 3991.86795711963e0;
  sso = 35.16504e0;
  s1 = sa * 35e0 / sso;
  pt0 = t + p * (8.65483913395442e-6 - s1 * 1.41636299744881e-6 - p * 7.38286467135737e-9 + t * (-8.38241357039698e-6 + s1 * 2.83933368585534e-8 + t * 1.77803965218656e-8 + p * 1.71155619208233e-10));
  dentropy_dt = cp0 / ((273.15e0 + pt0) * (1e0 - 0.05e0 * (1e0 - sa / sso)));
  true_entropy_part = gsw_entropy_part(sa, t, p);
  for (no_iter = 1; no_iter <= 2; no_iter++) {
    pt0_old = pt0;
    dentropy = gsw_entropy_part_zerop(sa, pt0_old) - true_entropy_part;
    pt0 = pt0_old - dentropy / dentropy_dt;
    pt0m = 0.5e0 * (pt0 + pt0_old);
    dentropy_dt = -gsw_gibbs_pt0_pt0(sa, pt0m);
    pt0 = pt0_old - dentropy / dentropy_dt;
  }
  return (pt0);
}
```