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

```asm
<IncompleteGamma>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xd0,%rsp
  movsd  %xmm0,-0xb8(%rbp)
  movsd  %xmm1,-0xc0(%rbp)
  movsd  %xmm2,-0xc8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  -0xc0(%rbp),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0xc8(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  D24(%rip),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  D101(%rip),%xmm0
  movsd  %xmm0,-0x60(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0xb8(%rbp),%xmm0
  jp     L143
  pxor   %xmm0,%xmm0
  ucomisd -0xb8(%rbp),%xmm0
  jne    L143
  pxor   %xmm0,%xmm0
  jmp    L195
L143:
  pxor   %xmm0,%xmm0
  comisd -0xb8(%rbp),%xmm0
  ja     L145
  pxor   %xmm0,%xmm0
  comisd -0x78(%rbp),%xmm0
  jb     L21
L145:
  movsd  D106(%rip),%xmm0
  jmp    L195
L21:
  mov    -0xb8(%rbp),%rax
  movq   %rax,%xmm0
  call   <log@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  mulsd  -0x78(%rbp),%xmm0
  subsd  -0xb8(%rbp),%xmm0
  subsd  -0x70(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0xb8(%rbp),%xmm0
  movsd  D42(%rip),%xmm1
  comisd %xmm1,%xmm0
  jbe    L161
  movsd  -0xb8(%rbp),%xmm0
  comisd -0x78(%rbp),%xmm0
  jae    L94
L161:
  movsd  D42(%rip),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  D42(%rip),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x78(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
L225:
  movsd  -0x98(%rbp),%xmm1
  movsd  D42(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0xb8(%rbp),%xmm0
  divsd  -0x98(%rbp),%xmm0
  movsd  -0x80(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  addsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  -0x80(%rbp),%xmm0
  comisd -0x68(%rbp),%xmm0
  jbe    L56
  jmp    L225
L56:
  movsd  -0x48(%rbp),%xmm0
  divsd  -0x78(%rbp),%xmm0
  movsd  -0xa0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  jmp    L95
L94:
  nop
  movsd  D42(%rip),%xmm0
  subsd  -0x78(%rbp),%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0xb8(%rbp),%xmm1
  movsd  D42(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  D42(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0xb8(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0xb8(%rbp),%xmm1
  movsd  D42(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0xb8(%rbp),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movsd  -0x28(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
L121:
  movsd  -0x90(%rbp),%xmm1
  movsd  D42(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x88(%rbp),%xmm1
  movsd  D0(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x80(%rbp),%xmm1
  movsd  D42(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x90(%rbp),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movl   $0x0,-0xa4(%rbp)
  jmp    L226
L23:
  mov    -0xa4(%rbp),%eax
  add    $0x2,%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  mov    -0xa4(%rbp),%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm1
  mulsd  -0x58(%rbp),%xmm1
  mov    -0xa4(%rbp),%eax
  add    $0x4,%eax
  subsd  %xmm1,%xmm0
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  addl   $0x1,-0xa4(%rbp)
L226:
  cmpl   $0x1,-0xa4(%rbp)
  jle    L23
  movsd  -0x18(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  jp     L250
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  je     L251
L250:
  movsd  -0x20(%rbp),%xmm0
  movsd  -0x18(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  subsd  -0x98(%rbp),%xmm0
  movq   D38(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x50(%rbp),%xmm0
  comisd -0x68(%rbp),%xmm0
  ja     L144
  movsd  -0x68(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  comisd -0x50(%rbp),%xmm0
  jae    L75
  jmp    L210
L144:
  nop
L210:
  movsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  jmp    L6
L251:
  nop
L6:
  movl   $0x0,-0xa4(%rbp)
  jmp    L10
L260:
  mov    -0xa4(%rbp),%eax
  add    $0x2,%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  mov    -0xa4(%rbp),%eax
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  addl   $0x1,-0xa4(%rbp)
L10:
  cmpl   $0x3,-0xa4(%rbp)
  jle    L260
  movsd  -0x20(%rbp),%xmm0
  movq   D38(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0x60(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L192
  jmp    L121
L192:
  movl   $0x0,-0xa4(%rbp)
  jmp    L59
L26:
  mov    -0xa4(%rbp),%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  divsd  -0x60(%rbp),%xmm0
  mov    -0xa4(%rbp),%eax
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  addl   $0x1,-0xa4(%rbp)
L59:
  cmpl   $0x3,-0xa4(%rbp)
  jle    L26
  jmp    L121
L75:
  nop
  movsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xa0(%rbp),%xmm1
  movsd  D42(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
L95:
  movsd  -0xa0(%rbp),%xmm0
L195:
  movq   %xmm0,%rax
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L53
  call   <__stack_chk_fail@plt>
L53:
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D24
type_infered: byte[8]
name: unpack_data
arguments: data_label: D101
type_infered: f64
name: unpack_data
arguments: data_label: D106
type_infered: f32[2]
name: unpack_data
arguments: data_label: D42
type_infered: f32[2]
name: unpack_data
arguments: data_label: D0
type_infered: f32[2]
name: unpack_data
arguments: data_label: D38
type_infered: qword

tool: {"D24": [0x3a,0x8c,0x30,0xe2,0x8e,0x79,0x45,0x3e]}

tool: D101: 1e+30

tool: {"D106": [0.0,-1.875]}

tool: {"D42": [0.0,1.875]}

tool: {"D0": [0.0,2.0]}

tool: D38: 9223372036854775807

assistant: ```c
double IncompleteGamma(double x, double alpha, double ln_gamma_alpha) {
  int i;
  double p = alpha, g = ln_gamma_alpha;
  double accurate = 1e-8, overflow = 1e30;
  double factor, gin = 0, rn = 0, a = 0, b = 0, an = 0, dif = 0, term = 0, pn[6];
  if (x == 0)
    return (0);
  if (x < 0 || p <= 0)
    return (-1);
  factor = (double)exp(p * (double)log(x) - x - g);
  if (x > 1 && x >= p)
    goto l30;
  gin = 1;
  term = 1;
  rn = p;
l20:
  rn++;
  term *= x / rn;
  gin += term;
  if (term > accurate)
    goto l20;
  gin *= factor / p;
  goto l50;
l30:
  a = 1 - p;
  b = a + x + 1;
  term = 0;
  pn[0] = 1;
  pn[1] = x;
  pn[2] = x + 1;
  pn[3] = x * b;
  gin = pn[2] / pn[3];
l32:
  a++;
  b += 2;
  term++;
  an = a * term;
  for (i = 0; i < 2; i++)
    pn[i + 4] = b * pn[i + 2] - an * pn[i];
  if (pn[5] == 0)
    goto l35;
  rn = pn[4] / pn[5];
  dif = fabs(gin - rn);
  if (dif > accurate)
    goto l34;
  if (dif <= accurate * rn)
    goto l42;
l34:
  gin = rn;
l35:
  for (i = 0; i < 4; i++)
    pn[i] = pn[i + 2];
  if (fabs(pn[4]) < overflow)
    goto l32;
  for (i = 0; i < 4; i++)
    pn[i] /= overflow;
  goto l32;
l42:
  gin = 1 - factor * gin;
l50:
  return (gin);
}
```