user: 请将以下的汇编代码转换成C语言代码：

```asm
<fgamma_lenz>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %edi,-0x44(%rbp)
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x44(%rbp),%xmm1
  movsd  D166(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  D117(%rip),%xmm0
  addsd  %xmm1,%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  D72(%rip),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  D117(%rip),%xmm0
  divsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movl   $0x1,-0x3c(%rbp)
  jmp    L25
L19:
  pxor   %xmm2,%xmm2
  cvtsi2sdl -0x3c(%rbp),%xmm2
  movsd  -0x18(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x3c(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x38(%rbp),%xmm1
  movsd  D110(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  -0x28(%rbp),%xmm0
  movsd  -0x38(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm0
  movq   D125(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D112(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L66
  movsd  D112(%rip),%xmm0
  movsd  %xmm0,-0x28(%rbp)
L66:
  movsd  -0x10(%rbp),%xmm0
  divsd  -0x30(%rbp),%xmm0
  movsd  -0x38(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movq   D125(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D112(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L37
  movsd  D112(%rip),%xmm0
  movsd  %xmm0,-0x30(%rbp)
L37:
  movsd  D117(%rip),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x20(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movsd  D117(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movq   D125(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D197(%rip),%xmm0
  comisd %xmm1,%xmm0
  ja     L16
  addl   $0x1,-0x3c(%rbp)
L25:
  cmpl   $0x64,-0x3c(%rbp)
  jle    L19
  jmp    L61
L16:
  nop
L61:
  movsd  D117(%rip),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  %xmm0,%xmm1
  movsd  D117(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movl   $0x1,-0x3c(%rbp)
  jmp    L22
L21:
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x3c(%rbp),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  movsd  -0x28(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  addl   $0x2,-0x3c(%rbp)
L22:
  mov    -0x44(%rbp),%eax
  add    %eax,%eax
  cmp    %eax,-0x3c(%rbp)
  jl     L21
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movsd  D153(%rip),%xmm1
  divsd  %xmm0,%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movq   D91(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  movsd  D166(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  movsd  -0x28(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: label: D166
datatype: f64
name: disassemble_data
arguments: label: D117
datatype: f64
name: disassemble_data
arguments: label: D72
datatype: byte[8]
name: disassemble_data
arguments: label: D110
datatype: byte[8]
name: disassemble_data
arguments: label: D125
datatype: qword
name: disassemble_data
arguments: label: D112
datatype: byte[8]
name: disassemble_data
arguments: label: D197
datatype: byte[8]
name: disassemble_data
arguments: label: D153
datatype: f64
name: disassemble_data
arguments: label: D91
datatype: byte[8]

tool: D166: 0.5

tool: D117: 1

tool: {"D72": [0xe9,0x8c,0xa0,0x39,0x59,0x3e,0x29,0x46]}

tool: {"D110": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x40]}

tool: D125: 9223372036854775807

tool: {"D112": [0xa0,0xc2,0xeb,0xfe,0x4b,0x48,0xb4,0x39]}

tool: {"D197": [0x16,0x56,0xe7,0x9e,0xaf,0x03,0xd2,0x3c]}

tool: D153: 0.8862269254528

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

assistant: ```c
double fgamma_lenz(int m, double x) {
  int i;
  double a, an, b, c, d, del, h;
  a = m + 0.5;
  b = x + 1.0 - a;
  c = 1.0 / 1.0E-30;
  d = 1.0 / b;
  h = d;
  for (i = 1; i <= 100; i++) {
    an = (a - i) * i;
    b += 2;
    d = an * d + b;
    if (fabs(d) < 1.0E-30)
      d = 1.0E-30;
    c = b + an / c;
    if (fabs(c) < 1.0E-30)
      c = 1.0E-30;
    d = 1.0 / d;
    del = d * c;
    h *= del;
    if (fabs(del - 1.0) < 1.0e-15)
      break;
  }
  d = 1.0;
  c = 1.0 / (x + x);
  for (i = 1; i <= (m + m - 1); i += 2)
    d *= i * c;
  d *= (double)0.8862269254528 / sqrt(x);
  return d - 0.5 * exp(-x) * h;
}
```