user: 下面的汇编指令对应的C源代码是什么:

```asm
<fgamma_set>:
  endbr64
  push   %r14
  movapd %xmm0,%xmm3
  pxor   %xmm0,%xmm0
  push   %r13
  push   %r12
  push   %rbp
  movslq %edi,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x20,%rsp
  comisd %xmm3,%xmm0
  ja     L72
L79:
  movsd  D17(%rip),%xmm0
  comisd %xmm3,%xmm0
  jbe    L22
  test   %ebp,%ebp
  js     L36
  movapd %xmm3,%xmm7
  movsd  D41(%rip),%xmm2
  xor    %eax,%eax
  movsd  D69(%rip),%xmm6
  mulsd  %xmm3,%xmm7
  movsd  D25(%rip),%xmm5
L8:
  lea    (%rax,%rax,1),%edx
  pxor   %xmm1,%xmm1
  movapd %xmm2,%xmm0
  cvtsi2sd %edx,%xmm1
  lea    0x0(,%rax,4),%edx
  movapd %xmm1,%xmm4
  addsd  %xmm6,%xmm1
  addsd  %xmm2,%xmm4
  divsd  %xmm4,%xmm0
  movapd %xmm3,%xmm4
  divsd  %xmm1,%xmm4
  pxor   %xmm1,%xmm1
  cvtsi2sd %edx,%xmm1
  mov    %rax,%rdx
  addsd  %xmm5,%xmm1
  subsd  %xmm4,%xmm0
  movapd %xmm7,%xmm4
  divsd  %xmm1,%xmm4
  addsd  %xmm4,%xmm0
  movsd  %xmm0,(%rbx,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rbp
  jne    L8
L36:
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax,%rax,1)
L22:
  movapd %xmm3,%xmm0
  movsd  %xmm3,0x8(%rsp)
  call   <fgamma_0@plt>
  movsd  %xmm0,(%rbx)
  test   %ebp,%ebp
  jle    L36
  movsd  0x8(%rsp),%xmm3
  movsd  D41(%rip),%xmm2
  movapd %xmm3,%xmm0
  movapd %xmm2,%xmm4
  movsd  %xmm2,0x18(%rsp)
  addsd  %xmm3,%xmm0
  divsd  %xmm0,%xmm4
  movapd %xmm3,%xmm0
  xorpd  D46(%rip),%xmm0
  movsd  %xmm4,0x10(%rsp)
  call   <exp@plt>
  movsd  0x10(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm2
  mov    $0x1,%edx
  movapd %xmm0,%xmm1
  movsd  0x8(%rsp),%xmm3
  mov    $0x1,%eax
  mulsd  %xmm4,%xmm1
  cs nopw 0x0(%rax,%rax,1)
L51:
  pxor   %xmm0,%xmm0
  mov    %eax,%r12d
  cvtsi2sd %edx,%xmm0
  mulsd  %xmm4,%xmm0
  comisd %xmm2,%xmm0
  ja     L4
  mulsd  -0x8(%rbx,%rax,8),%xmm0
  add    $0x2,%edx
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rbx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%ebp
  jge    L51
  jmp    L36
  nopl   0x0(%rax)
L72:
  lea    D32(%rip),%rdi
  movsd  %xmm3,0x8(%rsp)
  call   <printf@plt>
  movsd  0x8(%rsp),%xmm3
  jmp    L79
  nopl   (%rax)
L4:
  movsd  %xmm2,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm4,0x8(%rsp)
  cmp    %eax,%ebp
  jl     L36
  mov    %ebp,%edi
  movapd %xmm3,%xmm0
  movslq %ebp,%r13
  sub    $0x1,%ebp
  call   <fgamma@plt>
  movsd  0x8(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm2
  shl    $0x3,%r13
  movsd  0x10(%rsp),%xmm1
  lea    (%rbx,%r13,1),%r14
  divsd  %xmm4,%xmm2
  movsd  %xmm0,(%r14)
  mulsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  cmp    %r12d,%ebp
  jl     L36
  lea    -0x8(%rbx,%r13,1),%rdx
  lea    0x1(%rbp,%rbp,1),%eax
  lea    -0x1(%r12,%r12,1),%ecx
  nopl   (%rax)
L46:
  movsd  0x8(%rdx),%xmm1
  pxor   %xmm3,%xmm3
  sub    $0x8,%rdx
  cvtsi2sd %eax,%xmm3
  sub    $0x2,%eax
  mulsd  %xmm2,%xmm1
  addsd  %xmm0,%xmm1
  divsd  %xmm3,%xmm1
  movsd  %xmm1,0x8(%rdx)
  cmp    %eax,%ecx
  jne    L46
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
```

tool_calls: name: read_data
arguments: data_label: D17
data_type: f64
name: read_data
arguments: data_label: D41
data_type: f64
name: read_data
arguments: data_label: D69
data_type: f64
name: read_data
arguments: data_label: D25
data_type: f64
name: read_data
arguments: data_label: D46
data_type: f32[4]
name: read_data
arguments: data_label: D32
data_type: string

tool: D17: 0.0005

tool: D41: 1

tool: D69: 3

tool: D25: 10

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

tool: D32: fgamma_set - Warning : approximate values for x<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);
}
```