user: 请从给定的汇编指令推导出C源代码:

```asm
<gsw_pt0_from_t>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  add    $0xffffffffffffff80,%rsp
  movsd  %xmm0,-0x68(%rbp)
  movsd  %xmm1,-0x70(%rbp)
  movsd  %xmm2,-0x78(%rbp)
  movl   $0x0,-0x50(%rbp)
  movl   $0x2,-0x4c(%rbp)
  movsd  D3(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  D143(%rip),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x68(%rbp),%xmm1
  movsd  D125(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  divsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm1
  movsd  D93(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D22(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  -0x78(%rbp),%xmm2
  movsd  D135(%rip),%xmm1
  mulsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x28(%rbp),%xmm2
  movsd  D151(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  movsd  D148(%rip),%xmm3
  movapd %xmm0,%xmm2
  subsd  %xmm3,%xmm2
  movsd  -0x70(%rbp),%xmm3
  movsd  D122(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm2
  movsd  -0x78(%rbp),%xmm3
  movsd  D38(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  mulsd  -0x70(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x78(%rbp),%xmm0
  movsd  -0x70(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  D27(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm3
  divsd  -0x30(%rbp),%xmm3
  movsd  D129(%rip),%xmm0
  movapd %xmm0,%xmm2
  subsd  %xmm3,%xmm2
  movsd  D39(%rip),%xmm0
  mulsd  %xmm0,%xmm2
  movsd  D129(%rip),%xmm0
  subsd  %xmm2,%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x78(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  mov    -0x68(%rbp),%rax
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <gsw_entropy_part@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x20(%rbp)
  movl   $0x1,-0x54(%rbp)
  jmp    L1
L15:
  movsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x18(%rbp),%xmm0
  mov    -0x68(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <gsw_entropy_part_zerop@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  subsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x40(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x18(%rbp),%xmm1
  movsd  D152(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  mov    -0x68(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <gsw_gibbs_pt0_pt0@plt>
  movq   %xmm0,%rax
  movq   D42(%rip),%xmm0
  movq   %rax,%xmm4
  xorpd  %xmm0,%xmm4
  movapd %xmm4,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x40(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  addl   $0x1,-0x54(%rbp)
L1:
  cmpl   $0x2,-0x54(%rbp)
  jle    L15
  movsd  -0x48(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: read_data
arguments: data_label: D3
type_guessed: f64
name: read_data
arguments: data_label: D143
type_guessed: f64
name: read_data
arguments: data_label: D125
type_guessed: f64
name: read_data
arguments: data_label: D93
type_guessed: byte[8]
name: read_data
arguments: data_label: D22
type_guessed: byte[8]
name: read_data
arguments: data_label: D135
type_guessed: byte[8]
name: read_data
arguments: data_label: D151
type_guessed: byte[8]
name: read_data
arguments: data_label: D148
type_guessed: byte[8]
name: read_data
arguments: data_label: D122
type_guessed: byte[8]
name: read_data
arguments: data_label: D38
type_guessed: byte[8]
name: read_data
arguments: data_label: D27
type_guessed: f64
name: read_data
arguments: data_label: D129
type_guessed: f64
name: read_data
arguments: data_label: D39
type_guessed: f64
name: read_data
arguments: data_label: D152
type_guessed: f64
name: read_data
arguments: data_label: D42
type_guessed: byte[8]

tool: D3: 3991.86795711963

tool: D143: 35.16504

tool: D125: 35

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

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

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

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

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

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

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

tool: D27: 273.15

tool: D129: 1

tool: D39: 0.05

tool: D152: 0.5

tool: {"D42": [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);
}
```