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

```asm
<gsw_alpha>:
  endbr64
  movapd %xmm0,%xmm5
  pxor   %xmm0,%xmm0
  sub    $0x28,%rsp
  ucomisd %xmm5,%xmm0
  ja     L1
  movapd %xmm5,%xmm9
  sqrtsd %xmm9,%xmm9
L15:
  movsd  D89(%rip),%xmm3
  movsd  D614(%rip),%xmm0
  movapd %xmm1,%xmm8
  movsd  D15(%rip),%xmm4
  movsd  D171(%rip),%xmm6
  movsd  D363(%rip),%xmm14
  mulsd  %xmm1,%xmm3
  movsd  D248(%rip),%xmm15
  subsd  D671(%rip),%xmm3
  mulsd  %xmm1,%xmm0
  movsd  D97(%rip),%xmm11
  subsd  D677(%rip),%xmm0
  mulsd  %xmm1,%xmm4
  subsd  D347(%rip),%xmm4
  mulsd  %xmm5,%xmm6
  mulsd  %xmm1,%xmm3
  addsd  D528(%rip),%xmm3
  mulsd  %xmm1,%xmm0
  subsd  D587(%rip),%xmm0
  mulsd  %xmm5,%xmm4
  mulsd  %xmm5,%xmm11
  mulsd  %xmm1,%xmm3
  subsd  D722(%rip),%xmm3
  mulsd  %xmm9,%xmm3
  addsd  %xmm0,%xmm3
  movsd  D551(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D304(%rip),%xmm0
  mulsd  %xmm5,%xmm3
  mulsd  %xmm1,%xmm0
  addsd  D261(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D728(%rip),%xmm0
  addsd  %xmm0,%xmm3
  movsd  D265(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D372(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D250(%rip),%xmm0
  addsd  %xmm4,%xmm0
  movsd  D2(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  subsd  D518(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  subsd  D342(%rip),%xmm4
  addsd  %xmm6,%xmm4
  movapd %xmm1,%xmm6
  mulsd  %xmm2,%xmm4
  addsd  %xmm4,%xmm0
  movsd  D313(%rip),%xmm4
  mulsd  %xmm5,%xmm4
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm3
  movsd  D293(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D329(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D125(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  %xmm14,%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D154(%rip),%xmm0
  addsd  %xmm4,%xmm0
  movsd  D291(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  subsd  D284(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  subsd  D678(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  %xmm15,%xmm4
  mulsd  %xmm1,%xmm4
  subsd  D13(%rip),%xmm4
  mulsd  %xmm9,%xmm4
  addsd  %xmm4,%xmm0
  movsd  D19(%rip),%xmm4
  movsd  D674(%rip),%xmm13
  movsd  D665(%rip),%xmm7
  movsd  D108(%rip),%xmm12
  movsd  D732(%rip),%xmm10
  mulsd  %xmm1,%xmm4
  addsd  D25(%rip),%xmm4
  mulsd  %xmm5,%xmm0
  mulsd  %xmm7,%xmm6
  subsd  D107(%rip),%xmm6
  mulsd  %xmm10,%xmm8
  subsd  D729(%rip),%xmm8
  mulsd  %xmm1,%xmm4
  subsd  D663(%rip),%xmm4
  mulsd  %xmm5,%xmm7
  mulsd  %xmm5,%xmm6
  mulsd  %xmm2,%xmm8
  mulsd  %xmm1,%xmm4
  mulsd  %xmm2,%xmm10
  addsd  %xmm13,%xmm4
  mulsd  %xmm1,%xmm4
  addsd  D6(%rip),%xmm4
  addsd  %xmm4,%xmm0
  movsd  D226(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  D64(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  subsd  %xmm12,%xmm4
  mulsd  %xmm1,%xmm4
  subsd  D29(%rip),%xmm4
  addsd  %xmm6,%xmm4
  movsd  D498(%rip),%xmm6
  mulsd  %xmm1,%xmm6
  subsd  D713(%rip),%xmm6
  addsd  %xmm11,%xmm6
  mulsd  %xmm1,%xmm6
  subsd  D3(%rip),%xmm6
  addsd  %xmm8,%xmm6
  movsd  D690(%rip),%xmm8
  mulsd  %xmm1,%xmm8
  subsd  D297(%rip),%xmm8
  mulsd  %xmm2,%xmm6
  mulsd  %xmm1,%xmm8
  addsd  D528(%rip),%xmm8
  addsd  %xmm6,%xmm4
  movsd  D8(%rip),%xmm6
  mulsd  %xmm1,%xmm6
  mulsd  %xmm2,%xmm4
  addsd  D148(%rip),%xmm6
  mulsd  %xmm9,%xmm8
  mulsd  %xmm1,%xmm6
  addsd  D680(%rip),%xmm6
  addsd  %xmm0,%xmm4
  movsd  D478(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D493(%rip),%xmm0
  mulsd  %xmm1,%xmm6
  mulsd  %xmm1,%xmm0
  subsd  D396(%rip),%xmm0
  subsd  %xmm14,%xmm6
  mulsd  %xmm1,%xmm0
  addsd  %xmm15,%xmm0
  mulsd  %xmm9,%xmm0
  addsd  %xmm6,%xmm0
  movsd  D701(%rip),%xmm6
  mulsd  %xmm1,%xmm6
  addsd  D206(%rip),%xmm6
  mulsd  %xmm5,%xmm0
  mulsd  %xmm1,%xmm6
  subsd  D300(%rip),%xmm6
  mulsd  %xmm1,%xmm6
  addsd  %xmm13,%xmm6
  addsd  %xmm0,%xmm6
  movsd  D364(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D104(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  %xmm12,%xmm0
  addsd  %xmm0,%xmm7
  movsd  D354(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  subsd  D713(%rip),%xmm0
  addsd  %xmm11,%xmm0
  addsd  %xmm10,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm7
  mulsd  %xmm2,%xmm7
  addsd  %xmm7,%xmm6
  movapd %xmm6,%xmm0
  movsd  D548(%rip),%xmm6
  mulsd  %xmm3,%xmm0
  mulsd  %xmm1,%xmm6
  subsd  D677(%rip),%xmm6
  mulsd  %xmm4,%xmm3
  addsd  %xmm6,%xmm8
  movsd  D708(%rip),%xmm6
  mulsd  %xmm1,%xmm6
  subsd  D472(%rip),%xmm6
  mulsd  %xmm5,%xmm8
  mulsd  D15(%rip),%xmm5
  mulsd  %xmm1,%xmm6
  addsd  D261(%rip),%xmm6
  addsd  %xmm6,%xmm8
  movsd  D609(%rip),%xmm6
  mulsd  %xmm1,%xmm6
  mulsd  D152(%rip),%xmm1
  subsd  D372(%rip),%xmm6
  subsd  D518(%rip),%xmm1
  add    $0x28,%rsp
  addsd  %xmm6,%xmm5
  mulsd  %xmm2,%xmm1
  addsd  %xmm5,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm8
  mulsd  %xmm4,%xmm8
  subsd  %xmm8,%xmm0
  divsd  %xmm3,%xmm0
  ret
L1:
  movapd %xmm5,%xmm0
  movsd  %xmm2,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm5,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm5
  movapd %xmm0,%xmm9
  jmp    L15
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D89
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D614
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D15
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D171
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D363
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D248
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D671
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D97
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D677
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D347
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D528
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D587
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D722
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D551
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D304
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D261
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D728
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D265
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D372
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D250
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D2
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D518
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D342
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D313
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D293
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D329
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D125
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D154
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D291
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D284
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D678
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D13
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D19
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D674
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D665
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D108
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D732
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D25
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D107
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D729
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D663
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D6
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D226
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D64
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D29
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D498
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D713
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D3
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D690
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D297
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D8
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D148
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D680
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D478
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D493
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D396
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D701
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D206
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D300
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D364
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D104
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D354
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D548
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D708
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D472
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D609
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D152
type_guessed: byte[8]

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

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

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

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

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

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

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

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

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

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

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

tool: D587: 6.698001071123802

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

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

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

tool: D261: 2.839940833161907

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

tool: D6: 1

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

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

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

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

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

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

tool: {"D690": [0x98,0x26,0xa3,0xf3,0x3b,0x8f,0xa0,0x3e]}

tool: {"D297": [0xd3,0xa6,0x84,0xcc,0x10,0xed,0xfd,0x3e]}

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

tool: {"D148": [0x48,0xff,0xce,0x0c,0x9b,0x8d,0x28,0x3e]}

tool: {"D680": [0x99,0x60,0x15,0x69,0x6b,0xf3,0x90,0x3e]}

tool: {"D478": [0x89,0x87,0x54,0x54,0x4b,0x95,0xa7,0x3d]}

tool: {"D493": [0x88,0xf9,0xa3,0x1c,0x86,0x1b,0xe3,0x3d]}

tool: {"D396": [0x58,0x27,0x46,0x6b,0x35,0x7a,0x50,0x3e]}

tool: {"D701": [0x3f,0x56,0x29,0xd6,0xcf,0x2b,0x27,0x3e]}

tool: {"D206": [0xe6,0xb9,0x56,0xee,0x68,0x2c,0xcc,0x3e]}

tool: {"D300": [0x63,0xc6,0xe2,0xf0,0x2d,0xa3,0x08,0x3f]}

tool: {"D364": [0x14,0xcf,0x36,0xc0,0x61,0x7a,0x99,0xbd]}

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

tool: {"D354": [0x11,0xf7,0x4c,0x08,0xff,0xfe,0x24,0xbd]}

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

tool: {"D708": [0x2f,0x9e,0xf7,0x0f,0x47,0x0b,0x6d,0x3f]}

tool: {"D472": [0x8f,0x15,0x07,0xbe,0xd4,0x1d,0xb0,0x3f]}

tool: {"D609": [0x6f,0xb3,0xdd,0x0b,0xa2,0x41,0xdf,0x3e]}

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

assistant: ```c
double gsw_alpha(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, a01 = 2.839940833161907e0, a02 = -6.295518531177023e-2, a03 = 3.545416635222918e-3, a04 = -2.986498947203215e-2, a05 = 4.655718814958324e-4, a06 = 5.095422573880500e-4, a07 = -2.853969343267241e-5, a08 = 4.935118121048767e-7, a09 = -3.436090079851880e-4, a10 = 7.452101440691467e-6, a11 = 6.876837219536232e-7, a12 = -1.988366587925593e-8, a13 = -2.123038140592916e-11, a14 = 2.775927747785646e-3, a15 = -4.699214888271850e-5, a16 = 3.358540072460230e-6, a17 = 2.697475730017109e-9, a18 = -2.764306979894411e-5, a19 = 2.525874630197091e-7, a20 = 2.858362524508931e-9, a21 = -7.244588807799565e-11, a22 = 3.801564588876298e-7, a23 = -1.534575373851809e-8, a24 = -1.390254702334843e-10, a25 = 1.072438894227657e-11, a26 = -3.212746477974189e-7, a27 = 6.382827821123254e-9, a28 = -5.793038794625329e-12, a29 = 6.211426728363857e-10, a30 = -1.941660213148725e-11, a31 = -3.729652850731201e-14, a32 = 1.119522344879478e-14, a33 = 6.057902487546866e-17;
  double sqrtsa, v_hat_denominator, v_hat_numerator;
  double dvhatden_dct, dvhatnum_dct;
  sqrtsa = sqrt(sa);
  v_hat_denominator = v01 + ct * (v02 + ct * (v03 + v04 * ct)) + sa * (v05 + ct * (v06 + v07 * ct) + sqrtsa * (v08 + ct * (v09 + ct * (v10 + v11 * ct)))) + p * (v12 + ct * (v13 + v14 * ct) + sa * (v15 + v16 * ct) + p * (v17 + ct * (v18 + v19 * ct) + v20 * sa));
  v_hat_numerator = 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))))) + p * (v37 + ct * (v38 + ct * (v39 + v40 * ct)) + sa * (v41 + v42 * ct) + p * (v43 + ct * (v44 + v45 * ct + v46 * sa) + p * (v47 + v48 * ct)));
  dvhatden_dct = a01 + ct * (a02 + a03 * ct) + sa * (a04 + a05 * ct + sqrtsa * (a06 + ct * (a07 + a08 * ct))) + p * (a09 + a10 * ct + a11 * sa + p * (a12 + a13 * ct));
  dvhatnum_dct = a14 + ct * (a15 + ct * (a16 + a17 * ct)) + sa * (a18 + ct * (a19 + ct * (a20 + a21 * ct)) + sqrtsa * (a22 + ct * (a23 + ct * (a24 + a25 * ct)))) + p * (a26 + ct * (a27 + a28 * ct) + a29 * sa + p * (a30 + a31 * ct + a32 * sa + a33 * p));
  return ((v_hat_denominator * dvhatnum_dct - v_hat_numerator * dvhatden_dct) / (v_hat_numerator * v_hat_denominator));
}
```