user: 请显示与下面列出的汇编代码相匹配的C语言实现:

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

tool_calls: name: disassemble_data
arguments: label: D13
data_type: byte[8]
name: disassemble_data
arguments: label: D14
data_type: f64
name: disassemble_data
arguments: label: D2
data_type: f32[2]
name: disassemble_data
arguments: label: D11
data_type: byte[8]
name: disassemble_data
arguments: label: D37
data_type: byte[8]

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

tool: D14: 1

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

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

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