user: What is the c source code of the assembly code below?

```asm
<gss3tar>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  movsd  %xmm0,0x8(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm2,0x18(%rsp)
  mov    %rdi,%rbx
  mov    %rsi,%r15
  mov    %rdx,%r14
  mov    %rcx,%r13
  mov    %r8,%r12
  mov    %r9,%rbp
  movsd  D23(%rip),%xmm1
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movsd  0x8(%rsp),%xmm13
  movapd %xmm13,%xmm9
  mulsd  %xmm0,%xmm9
  movsd  0x10(%rsp),%xmm12
  movapd %xmm12,%xmm0
  divsd  %xmm1,%xmm0
  movsd  0x18(%rsp),%xmm6
  movapd %xmm6,%xmm2
  mulsd  %xmm6,%xmm0
  movsd  D38(%rip),%xmm5
  movapd %xmm0,%xmm4
  addsd  %xmm5,%xmm4
  divsd  %xmm4,%xmm0
  movapd %xmm4,%xmm3
  mulsd  %xmm4,%xmm3
  movapd %xmm6,%xmm8
  divsd  %xmm3,%xmm8
  mulsd  D43(%rip),%xmm2
  mulsd  %xmm6,%xmm2
  mulsd  %xmm4,%xmm3
  divsd  %xmm3,%xmm2
  movsd  D22(%rip),%xmm3
  movapd %xmm12,%xmm6
  mulsd  %xmm3,%xmm6
  divsd  %xmm9,%xmm6
  divsd  %xmm1,%xmm5
  movapd %xmm8,%xmm4
  mulsd  %xmm6,%xmm4
  movapd %xmm8,%xmm7
  mulsd  %xmm5,%xmm7
  movapd %xmm2,%xmm10
  mulsd  %xmm6,%xmm10
  movapd %xmm12,%xmm1
  mulsd  D5(%rip),%xmm1
  movapd %xmm13,%xmm11
  mulsd  %xmm9,%xmm11
  divsd  %xmm11,%xmm1
  mulsd  %xmm8,%xmm1
  mulsd  %xmm10,%xmm6
  addsd  %xmm6,%xmm1
  mulsd  %xmm5,%xmm2
  mulsd  %xmm5,%xmm2
  divsd  %xmm9,%xmm3
  mulsd  %xmm8,%xmm3
  mulsd  %xmm5,%xmm10
  addsd  %xmm10,%xmm3
  movapd %xmm0,%xmm5
  mulsd  0x8(%rbx),%xmm5
  addsd  (%rbx),%xmm5
  movapd %xmm0,%xmm6
  mulsd  0x10(%rbx),%xmm6
  mulsd  %xmm0,%xmm6
  addsd  %xmm6,%xmm5
  movsd  %xmm5,(%r15)
  movsd  0x10(%rbx),%xmm5
  addsd  %xmm5,%xmm5
  mulsd  %xmm0,%xmm5
  mulsd  %xmm4,%xmm5
  movapd %xmm4,%xmm8
  mulsd  0x8(%rbx),%xmm8
  pxor   %xmm6,%xmm6
  addsd  %xmm6,%xmm8
  addsd  %xmm8,%xmm5
  movsd  %xmm5,(%r14)
  movsd  0x10(%rbx),%xmm5
  addsd  %xmm5,%xmm5
  mulsd  %xmm0,%xmm5
  mulsd  %xmm7,%xmm5
  movapd %xmm7,%xmm8
  mulsd  0x8(%rbx),%xmm8
  addsd  %xmm6,%xmm8
  addsd  %xmm8,%xmm5
  movsd  %xmm5,0x0(%r13)
  movapd %xmm4,%xmm5
  addsd  %xmm4,%xmm5
  addsd  %xmm0,%xmm0
  mulsd  %xmm5,%xmm4
  movapd %xmm0,%xmm8
  mulsd  %xmm1,%xmm8
  addsd  %xmm8,%xmm4
  mulsd  0x10(%rbx),%xmm4
  mulsd  0x8(%rbx),%xmm1
  addsd  %xmm6,%xmm1
  addsd  %xmm1,%xmm4
  movsd  %xmm4,(%r12)
  movapd %xmm7,%xmm1
  addsd  %xmm7,%xmm1
  mulsd  %xmm7,%xmm1
  movapd %xmm0,%xmm4
  mulsd  %xmm2,%xmm4
  addsd  %xmm4,%xmm1
  mulsd  0x10(%rbx),%xmm1
  mulsd  0x8(%rbx),%xmm2
  addsd  %xmm6,%xmm2
  addsd  %xmm2,%xmm1
  mov    0x60(%rsp),%rax
  movsd  %xmm1,(%rax)
  mulsd  %xmm7,%xmm5
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm5
  mulsd  0x10(%rbx),%xmm5
  mulsd  0x8(%rbx),%xmm3
  addsd  %xmm6,%xmm3
  addsd  %xmm3,%xmm5
  movsd  %xmm5,0x0(%rbp)
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: data_label: D23
type_guessed: byte[8]
name: parse_data
arguments: data_label: D38
type_guessed: f64
name: parse_data
arguments: data_label: D43
type_guessed: f32[2]
name: parse_data
arguments: data_label: D22
type_guessed: byte[8]
name: parse_data
arguments: data_label: D5
type_guessed: byte[8]

tool: {"D23": [0x55,0x55,0x55,0x55,0x55,0x55,0x05,0x40]}

tool: D38: 1

tool: {"D43": [0.0,-2.0]}

tool: {"D22": [0x55,0x55,0x55,0x55,0x55,0x55,0x05,0xc0]}

tool: {"D5": [0xe4,0x38,0x8e,0xe3,0x38,0x8e,0x23,0x40]}

assistant: ```c
void gss3tar(double r, double g, double alpha, double *a, double *f, double *fr, double *fg, double *frr, double *frg, double *fgg) {
  double r83, r113, r143, s2, s2p1, u;
  r83 = pow(r, 8.0 / 3.0);
  r113 = r83 * r;
  r143 = r113 * r;
  s2 = g / r83;
  s2p1 = 1.0 + alpha * s2;
  u = alpha * s2 / s2p1;
  double D1uDs21 = alpha / (s2p1 * s2p1);
  double D2uDs22 = -2.0 * alpha * alpha / (s2p1 * s2p1 * s2p1);
  double D1s2Dr1 = (-8.0 / 3.0) * g / r113;
  double D2s2Dr2 = (+88.0 / 9.0) * g / r143;
  double D1s2Dg1 = 1.0 / r83;
  double D2s2Dr1g1 = (-8.0 / 3.0) / r113;
  double D1uDr1 = D1uDs21 * D1s2Dr1;
  double D1uDg1 = D1uDs21 * D1s2Dg1;
  double D2uDr2 = D2uDs22 * D1s2Dr1 * D1s2Dr1 + D1uDs21 * D2s2Dr2;
  double D2uDg2 = D2uDs22 * D1s2Dg1 * D1s2Dg1;
  double D2uDr1g1 = D2uDs22 * D1s2Dr1 * D1s2Dg1 + D1uDs21 * D2s2Dr1g1;
  *f = a[0] + a[1] * u + a[2] * u * u;
  *fr = 0.0 + a[1] * D1uDr1 + a[2] * 2.0 * u * D1uDr1;
  *fg = 0.0 + a[1] * D1uDg1 + a[2] * 2.0 * u * D1uDg1;
  *frr = 0.0 + a[1] * D2uDr2 + a[2] * (2.0 * D1uDr1 * D1uDr1 + 2.0 * u * D2uDr2);
  *fgg = 0.0 + a[1] * D2uDg2 + a[2] * (2.0 * D1uDg1 * D1uDg1 + 2.0 * u * D2uDg2);
  *frg = 0.0 + a[1] * D2uDr1g1 + a[2] * (2.0 * D1uDr1 * D1uDg1 + 2.0 * u * D2uDr1g1);
  return;
}
```