user: Converte this assembly snippet into its equivalent C code:

```asm
<fgamma_setf>:
  endbr64
  push   %r14
  pxor   %xmm5,%xmm5
  movaps %xmm0,%xmm2
  push   %r13
  cvtss2sd %xmm0,%xmm5
  push   %r12
  push   %rbp
  movslq %edi,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x10,%rsp
  movsd  D61(%rip),%xmm0
  comisd %xmm5,%xmm0
  jbe    L93
  test   %ebp,%ebp
  js     L80
  mulss  %xmm2,%xmm2
  pxor   %xmm0,%xmm0
  xor    %eax,%eax
  movsd  D1(%rip),%xmm3
  movsd  D32(%rip),%xmm7
  movsd  D11(%rip),%xmm6
  cvtss2sd %xmm2,%xmm0
  nopl   0x0(%rax)
L9:
  lea    (%rax,%rax,1),%edx
  pxor   %xmm2,%xmm2
  movapd %xmm3,%xmm1
  cvtsi2sd %edx,%xmm2
  lea    0x0(,%rax,4),%edx
  movapd %xmm2,%xmm4
  addsd  %xmm7,%xmm2
  addsd  %xmm3,%xmm4
  divsd  %xmm4,%xmm1
  movapd %xmm5,%xmm4
  divsd  %xmm2,%xmm4
  pxor   %xmm2,%xmm2
  cvtsi2sd %edx,%xmm2
  mov    %rax,%rdx
  addsd  %xmm6,%xmm2
  subsd  %xmm4,%xmm1
  movapd %xmm0,%xmm4
  divsd  %xmm2,%xmm4
  addsd  %xmm4,%xmm1
  cvtsd2ss %xmm1,%xmm1
  movss  %xmm1,(%rbx,%rax,4)
  add    $0x1,%rax
  cmp    %rdx,%rbp
  jne    L9
L80:
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax)
L93:
  movaps %xmm2,%xmm0
  movss  %xmm2,0x8(%rsp)
  call   <fgamma_0f@plt>
  movss  %xmm0,(%rbx)
  test   %ebp,%ebp
  jle    L80
  movss  0x8(%rsp),%xmm2
  movss  D17(%rip),%xmm4
  movaps %xmm2,%xmm0
  movaps %xmm4,%xmm5
  addss  %xmm2,%xmm0
  divss  %xmm0,%xmm5
  movaps %xmm2,%xmm0
  xorps  D54(%rip),%xmm0
  movss  %xmm5,0xc(%rsp)
  call   <exp@plt>
  mov    $0x1,%edx
  pxor   %xmm1,%xmm1
  movss  0xc(%rsp),%xmm5
  cvtsi2ss %edx,%xmm1
  movss  D17(%rip),%xmm4
  movaps %xmm0,%xmm3
  mov    $0x1,%eax
  mulss  %xmm5,%xmm3
  movss  0x8(%rsp),%xmm2
  mov    %eax,%r12d
  mulss  %xmm5,%xmm1
  comiss %xmm4,%xmm1
  ja     L86
  nopl   0x0(%rax,%rax,1)
L35:
  mulss  -0x4(%rbx,%rax,4),%xmm1
  add    $0x2,%edx
  subss  %xmm3,%xmm1
  movss  %xmm1,(%rbx,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%ebp
  jl     L80
  pxor   %xmm1,%xmm1
  mov    %eax,%r12d
  cvtsi2ss %edx,%xmm1
  mulss  %xmm5,%xmm1
  comiss %xmm4,%xmm1
  jbe    L35
L86:
  movss  %xmm3,0xc(%rsp)
  movss  %xmm5,0x8(%rsp)
  cmp    %eax,%ebp
  jl     L80
  mov    %ebp,%edi
  movaps %xmm2,%xmm0
  movslq %ebp,%r13
  sub    $0x1,%ebp
  call   <fgammaf@plt>
  movss  0x8(%rsp),%xmm5
  shl    $0x2,%r13
  movss  D17(%rip),%xmm4
  movss  0xc(%rsp),%xmm3
  lea    (%rbx,%r13,1),%r14
  divss  %xmm5,%xmm4
  movss  %xmm0,(%r14)
  mulss  %xmm4,%xmm3
  cmp    %r12d,%ebp
  jl     L80
  lea    -0x4(%rbx,%r13,1),%rdx
  lea    0x1(%rbp,%rbp,1),%eax
  lea    -0x1(%r12,%r12,1),%ecx
  cs nopw 0x0(%rax,%rax,1)
L22:
  movss  0x4(%rdx),%xmm1
  pxor   %xmm0,%xmm0
  sub    $0x4,%rdx
  cvtsi2ss %eax,%xmm0
  sub    $0x2,%eax
  mulss  %xmm4,%xmm1
  addss  %xmm3,%xmm1
  divss  %xmm0,%xmm1
  movss  %xmm1,0x4(%rdx)
  cmp    %eax,%ecx
  jne    L22
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
```

tool_calls: name: read_data
arguments: data_label: D61
type_guessed: f64
name: read_data
arguments: data_label: D1
type_guessed: f64
name: read_data
arguments: data_label: D32
type_guessed: f64
name: read_data
arguments: data_label: D11
type_guessed: f64
name: read_data
arguments: data_label: D17
type_guessed: byte[4]
name: read_data
arguments: data_label: D54
type_guessed: byte[16]

tool: D61: 0.0005

tool: D1: 1

tool: D32: 3

tool: D11: 10

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

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