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

```asm
<fgamma_set>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x20,%rsp
  mov    %edi,%ebp
  movapd %xmm0,%xmm7
  movsd  %xmm0,0x8(%rsp)
  mov    %rsi,%rbx
  pxor   %xmm0,%xmm0
  comisd %xmm7,%xmm0
  ja     L101
L29:
  movsd  D50(%rip),%xmm0
  movsd  0x8(%rsp),%xmm5
  comisd %xmm5,%xmm0
  jbe    L128
  test   %ebp,%ebp
  js     L73
  movapd %xmm5,%xmm6
  mulsd  %xmm5,%xmm6
  lea    0x1(%rbp),%ecx
  mov    $0x0,%eax
  movsd  D91(%rip),%xmm2
  movsd  D75(%rip),%xmm5
  movsd  D71(%rip),%xmm4
L2:
  lea    (%rax,%rax,1),%edx
  pxor   %xmm1,%xmm1
  cvtsi2sd %edx,%xmm1
  movapd %xmm1,%xmm3
  addsd  %xmm2,%xmm3
  movapd %xmm2,%xmm0
  divsd  %xmm3,%xmm0
  addsd  %xmm5,%xmm1
  movsd  0x8(%rsp),%xmm3
  divsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm0
  lea    0x0(,%rax,4),%edx
  pxor   %xmm1,%xmm1
  cvtsi2sd %edx,%xmm1
  addsd  %xmm4,%xmm1
  movapd %xmm6,%xmm3
  divsd  %xmm1,%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%rbx,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rcx
  jne    L2
L73:
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L101:
  lea    D14(%rip),%rdi
  call   <printf@plt>
  jmp    L29
L128:
  movsd  0x8(%rsp),%xmm0
  call   <fgamma_0@plt>
  movsd  %xmm0,(%rbx)
  test   %ebp,%ebp
  jle    L73
  movsd  0x8(%rsp),%xmm6
  movapd %xmm6,%xmm0
  addsd  %xmm6,%xmm0
  movsd  D91(%rip),%xmm4
  divsd  %xmm0,%xmm4
  movsd  %xmm4,0x10(%rsp)
  xorpd  D57(%rip),%xmm6
  movapd %xmm6,%xmm0
  call   <exp@plt>
  movsd  0x10(%rsp),%xmm6
  mulsd  %xmm6,%xmm0
  movsd  %xmm0,0x18(%rsp)
  comisd D91(%rip),%xmm6
  ja     L109
  movapd %xmm6,%xmm0
  mov    $0x3,%edx
  mov    $0x2,%eax
  movsd  D91(%rip),%xmm1
L44:
  mulsd  -0x10(%rbx,%rax,8),%xmm0
  subsd  0x18(%rsp),%xmm0
  movsd  %xmm0,-0x8(%rbx,%rax,8)
  mov    %eax,%r12d
  cmp    %eax,%ebp
  jl     L73
  pxor   %xmm0,%xmm0
  cvtsi2sd %edx,%xmm0
  mulsd  0x10(%rsp),%xmm0
  add    $0x1,%rax
  add    $0x2,%edx
  comisd %xmm1,%xmm0
  jbe    L44
L74:
  movslq %ebp,%r13
  shl    $0x3,%r13
  lea    (%rbx,%r13,1),%r14
  movsd  0x8(%rsp),%xmm0
  mov    %ebp,%edi
  call   <fgamma@plt>
  movsd  %xmm0,(%r14)
  movsd  D91(%rip),%xmm3
  divsd  0x10(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm1
  mulsd  %xmm3,%xmm1
  sub    $0x1,%ebp
  cmp    %r12d,%ebp
  jl     L73
  lea    -0x8(%rbx,%r13,1),%rdx
  lea    0x1(%rbp,%rbp,1),%eax
  lea    -0x1(%r12,%r12,1),%ecx
L21:
  movapd %xmm3,%xmm0
  mulsd  0x8(%rdx),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%rdx)
  sub    $0x8,%rdx
  sub    $0x2,%eax
  cmp    %ecx,%eax
  jne    L21
  jmp    L73
L109:
  mov    $0x1,%r12d
  jmp    L74
```

tool_calls: name: read_data
arguments: data_label: D50
type_infered: f64
name: read_data
arguments: data_label: D91
type_infered: f64[2]
name: read_data
arguments: data_label: D75
type_infered: f64
name: read_data
arguments: data_label: D71
type_infered: f64
name: read_data
arguments: data_label: D14
type_infered: string
name: read_data
arguments: data_label: D57
type_infered: f32[4]

tool: D50: 0.0005

tool: {"D91": [1.0,3.0]}

tool: D75: 3

tool: D71: 10

tool: D14: fgamma_set - Warning : approximate values for x<0


tool: {"D57": [0.0,-0.0,0.0,0.0]}

assistant: ```c
void fgamma_set(int max, double x, double *fset) {
  double a, b, c;
  int i, j;
  if (x < 0.0) {
    printf("fgamma_set - Warning : approximate values for x<0\n");
  }
  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_0(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] = fgamma(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);
}
```