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

```asm
<fgamma_setf>:
  endbr64
  push   %r14
  pxor   %xmm10,%xmm10
  movaps %xmm0,%xmm2
  push   %r13
  cvtss2sd %xmm0,%xmm10
  push   %r12
  push   %rbp
  mov    %edi,%ebp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x10,%rsp
  movsd  D116(%rip),%xmm0
  comisd %xmm10,%xmm0
  jbe    L120
  test   %edi,%edi
  js     L77
  mulss  %xmm2,%xmm2
  pxor   %xmm0,%xmm0
  lea    0x1(%rdi),%ecx
  cvtss2sd %xmm2,%xmm0
  cmp    $0x2,%edi
  jle    L115
  mov    %ecx,%edx
  movapd %xmm10,%xmm9
  movapd %xmm0,%xmm8
  mov    %rsi,%rax
  shr    $0x2,%edx
  movdqa D47(%rip),%xmm5
  movdqa D89(%rip),%xmm11
  unpcklpd %xmm9,%xmm9
  sub    $0x1,%edx
  movapd D105(%rip),%xmm4
  movapd D63(%rip),%xmm7
  unpcklpd %xmm8,%xmm8
  shl    $0x4,%rdx
  movapd D124(%rip),%xmm6
  lea    0x10(%rsi,%rdx,1),%rdx
  cs nopw 0x0(%rax,%rax,1)
L118:
  movdqa %xmm5,%xmm1
  movapd %xmm4,%xmm3
  movapd %xmm4,%xmm14
  add    $0x10,%rax
  movdqa %xmm1,%xmm2
  pslld  $0x2,%xmm1
  paddd  %xmm11,%xmm5
  pslld  $0x1,%xmm2
  cvtdq2pd %xmm2,%xmm12
  movapd %xmm12,%xmm13
  addpd  %xmm7,%xmm12
  pshufd $0xee,%xmm2,%xmm2
  addpd  %xmm4,%xmm13
  cvtdq2pd %xmm2,%xmm2
  divpd  %xmm13,%xmm3
  movapd %xmm9,%xmm13
  divpd  %xmm12,%xmm13
  cvtdq2pd %xmm1,%xmm12
  addpd  %xmm6,%xmm12
  pshufd $0xee,%xmm1,%xmm1
  cvtdq2pd %xmm1,%xmm1
  addpd  %xmm6,%xmm1
  subpd  %xmm13,%xmm3
  movapd %xmm8,%xmm13
  divpd  %xmm12,%xmm13
  movapd %xmm2,%xmm12
  addpd  %xmm7,%xmm2
  addpd  %xmm4,%xmm12
  divpd  %xmm12,%xmm14
  movapd %xmm8,%xmm12
  addpd  %xmm13,%xmm3
  movapd %xmm9,%xmm13
  divpd  %xmm2,%xmm13
  movapd %xmm14,%xmm2
  divpd  %xmm1,%xmm12
  cvtpd2ps %xmm3,%xmm1
  subpd  %xmm13,%xmm2
  addpd  %xmm12,%xmm2
  cvtpd2ps %xmm2,%xmm2
  movlhps %xmm2,%xmm1
  movups %xmm1,-0x10(%rax)
  cmp    %rdx,%rax
  jne    L118
  mov    %ecx,%eax
  and    $0xfffffffc,%eax
  and    $0x3,%ecx
  je     L77
L23:
  lea    (%rax,%rax,1),%esi
  pxor   %xmm3,%xmm3
  movapd %xmm10,%xmm5
  movslq %eax,%rdi
  cvtsi2sd %esi,%xmm3
  movsd  D105(%rip),%xmm1
  movapd %xmm0,%xmm6
  lea    0x0(,%rax,4),%ecx
  lea    0x0(,%rdi,4),%rdx
  movapd %xmm1,%xmm2
  movapd %xmm3,%xmm4
  addsd  %xmm1,%xmm4
  divsd  %xmm4,%xmm2
  movsd  D63(%rip),%xmm4
  addsd  %xmm4,%xmm3
  divsd  %xmm3,%xmm5
  pxor   %xmm3,%xmm3
  cvtsi2sd %ecx,%xmm3
  subsd  %xmm5,%xmm2
  movsd  D124(%rip),%xmm5
  addsd  %xmm5,%xmm3
  divsd  %xmm3,%xmm6
  addsd  %xmm6,%xmm2
  cvtsd2ss %xmm2,%xmm2
  movss  %xmm2,(%rbx,%rdi,4)
  lea    0x1(%rax),%edi
  cmp    %edi,%ebp
  jl     L77
  lea    0x2(%rsi),%edi
  pxor   %xmm3,%xmm3
  movapd %xmm1,%xmm2
  add    $0x2,%eax
  cvtsi2sd %edi,%xmm3
  lea    0x4(%rcx),%edi
  movapd %xmm3,%xmm6
  addsd  %xmm4,%xmm3
  addsd  %xmm1,%xmm6
  divsd  %xmm6,%xmm2
  movapd %xmm10,%xmm6
  divsd  %xmm3,%xmm6
  pxor   %xmm3,%xmm3
  cvtsi2sd %edi,%xmm3
  addsd  %xmm5,%xmm3
  subsd  %xmm6,%xmm2
  movapd %xmm0,%xmm6
  divsd  %xmm3,%xmm6
  addsd  %xmm6,%xmm2
  cvtsd2ss %xmm2,%xmm2
  movss  %xmm2,0x4(%rbx,%rdx,1)
  cmp    %eax,%ebp
  jl     L77
  add    $0x4,%esi
  pxor   %xmm2,%xmm2
  add    $0x8,%ecx
  cvtsi2sd %esi,%xmm2
  movapd %xmm2,%xmm3
  addsd  %xmm4,%xmm2
  addsd  %xmm1,%xmm3
  divsd  %xmm2,%xmm10
  pxor   %xmm2,%xmm2
  cvtsi2sd %ecx,%xmm2
  addsd  %xmm5,%xmm2
  divsd  %xmm3,%xmm1
  divsd  %xmm2,%xmm0
  subsd  %xmm10,%xmm1
  addsd  %xmm0,%xmm1
  cvtsd2ss %xmm1,%xmm1
  movss  %xmm1,0x8(%rbx,%rdx,1)
L77:
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax)
L120:
  movaps %xmm2,%xmm0
  movss  %xmm2,0x8(%rsp)
  call   <fgamma_0f@plt>
  movss  %xmm0,(%rbx)
  test   %ebp,%ebp
  jle    L77
  movss  0x8(%rsp),%xmm2
  movss  D57(%rip),%xmm4
  movaps %xmm2,%xmm0
  movaps %xmm4,%xmm5
  addss  %xmm2,%xmm0
  divss  %xmm0,%xmm5
  movaps %xmm2,%xmm0
  xorps  D98(%rip),%xmm0
  movss  %xmm5,0xc(%rsp)
  call   <exp@plt>
  mov    $0x1,%edx
  pxor   %xmm1,%xmm1
  movss  0xc(%rsp),%xmm5
  cvtsi2ss %edx,%xmm1
  movss  D57(%rip),%xmm4
  mov    $0x1,%eax
  movss  0x8(%rsp),%xmm2
  mulss  %xmm5,%xmm0
  mov    %eax,%r12d
  mulss  %xmm5,%xmm1
  movaps %xmm0,%xmm3
  comiss %xmm4,%xmm1
  ja     L83
  nopl   0x0(%rax,%rax,1)
L66:
  mulss  -0x4(%rbx,%rax,4),%xmm1
  add    $0x2,%edx
  subss  %xmm3,%xmm1
  movss  %xmm1,(%rbx,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%ebp
  jl     L77
  pxor   %xmm1,%xmm1
  mov    %eax,%r12d
  cvtsi2ss %edx,%xmm1
  mulss  %xmm5,%xmm1
  comiss %xmm4,%xmm1
  jbe    L66
L83:
  movss  %xmm3,0xc(%rsp)
  movss  %xmm5,0x8(%rsp)
  cmp    %eax,%ebp
  jl     L77
  mov    %ebp,%edi
  movaps %xmm2,%xmm0
  movslq %ebp,%r13
  sub    $0x1,%ebp
  call   <fgammaf@plt>
  movss  0x8(%rsp),%xmm5
  shl    $0x2,%r13
  movss  D57(%rip),%xmm4
  movss  0xc(%rsp),%xmm3
  lea    (%rbx,%r13,1),%r14
  divss  %xmm5,%xmm4
  movss  %xmm0,(%r14)
  mulss  %xmm4,%xmm3
  cmp    %r12d,%ebp
  jl     L77
  lea    0x1(%rbp,%rbp,1),%eax
  lea    -0x4(%rbx,%r13,1),%rdx
  lea    -0x1(%r12,%r12,1),%ecx
  cs nopw 0x0(%rax,%rax,1)
L91:
  mulss  %xmm4,%xmm0
  pxor   %xmm1,%xmm1
  sub    $0x4,%rdx
  cvtsi2ss %eax,%xmm1
  sub    $0x2,%eax
  addss  %xmm3,%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,0x4(%rdx)
  cmp    %ecx,%eax
  jne    L91
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L115:
  xor    %eax,%eax
  jmp    L23
```

tool_calls: name: unpack_data
arguments: data_label: D116
type_infered: f64
name: unpack_data
arguments: data_label: D47
type_infered: byte[16]
name: unpack_data
arguments: data_label: D89
type_infered: byte[16]
name: unpack_data
arguments: data_label: D105
type_infered: f64[2]
name: unpack_data
arguments: data_label: D63
type_infered: f64[2]
name: unpack_data
arguments: data_label: D124
type_infered: f64[2]
name: unpack_data
arguments: data_label: D57
type_infered: byte[4]
name: unpack_data
arguments: data_label: D98
type_infered: byte[16]

tool: D116: 0.0005

tool: {"D47": [0x00,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x03,0x00,0x00,0x00]}

tool: {"D89": [0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00]}

tool: {"D105": [1.0,1.0]}

tool: {"D63": [3.0,3.0]}

tool: {"D124": [10.0,10.0]}

tool: {"D57": [0x00,0x00,0x80,0x3f]}

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

assistant: ```c
void fgamma_setf(int max, float x, float *fset) {
  float a, b, c;
  int i, j;
  if (x < 0.0005) {
    for (i = 0; i <= max; i++)
      fset[i] = 1.0 / (i + i + 1.0) - x / (i + i + 3.0) + x * x / (i + i + i + i + 10.0);
    return;
  }
  fset[0] = fgamma_0f(x);
  if (max < 1)
    return;
  a = 1.0 / (x + x);
  c = a * exp(-x);
  for (i = 1; i <= max; i++) {
    if ((b = (i + i - 1) * a) > 1.0)
      break;
    fset[i] = b * fset[i - 1] - c;
  }
  if (max < i)
    return;
  fset[max] = fgammaf(max, x);
  a = 1.0 / a;
  c = c * a;
  for (j = max - 1; j >= i; j--)
    fset[j] = (a * fset[j + 1] + c) / (j + j + 1);
}
```