user: Please show me the C implementation that matches the assembly code listed below:

```asm
<IncompleteGamma>:
  endbr64
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  movq   %xmm0,%rbx
  movsd  %xmm1,(%rsp)
  movq   %xmm2,%rbp
  mov    %fs:0x28,%rax
  mov    %rax,0x48(%rsp)
  xor    %eax,%eax
  ucomisd D43(%rip),%xmm0
  jp     L68
  pxor   %xmm0,%xmm0
  je     L197
L68:
  pxor   %xmm0,%xmm0
  movq   %rbx,%xmm7
  comisd %xmm7,%xmm0
  ja     L22
  movsd  D107(%rip),%xmm0
  pxor   %xmm1,%xmm1
  comisd (%rsp),%xmm1
  jae    L197
  movq   %rbx,%xmm0
  call   <log@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  mulsd  (%rsp),%xmm0
  movq   %rbx,%xmm7
  subsd  %xmm7,%xmm0
  movq   %rbp,%xmm7
  subsd  %xmm7,%xmm0
  call   <exp@plt>
  pxor   %xmm7,%xmm7
  cvtsi2sd %rax,%xmm7
  movsd  %xmm7,0x8(%rsp)
  movq   %rbx,%xmm7
  comisd D92(%rip),%xmm7
  jbe    L13
  movsd  (%rsp),%xmm6
  comisd %xmm6,%xmm7
  jb     L20
  movsd  D92(%rip),%xmm3
  movapd %xmm3,%xmm2
  subsd  %xmm6,%xmm2
  movq   %rbx,%xmm1
  addsd  %xmm2,%xmm1
  addsd  %xmm3,%xmm1
  movsd  %xmm3,0x10(%rsp)
  mov    %rbx,0x18(%rsp)
  addsd  %xmm7,%xmm3
  movsd  %xmm3,0x20(%rsp)
  movq   %rbx,%xmm4
  mulsd  %xmm1,%xmm4
  movsd  %xmm4,0x28(%rsp)
  movapd %xmm3,%xmm0
  divsd  %xmm4,%xmm0
  pxor   %xmm4,%xmm4
  movsd  D92(%rip),%xmm6
  movsd  D143(%rip),%xmm9
  movapd %xmm4,%xmm8
  movsd  D53(%rip),%xmm11
  movapd %xmm11,%xmm7
  jmp    L96
L13:
  movsd  (%rsp),%xmm2
  movsd  D92(%rip),%xmm1
  movapd %xmm1,%xmm3
  jmp    L164
L20:
  movsd  (%rsp),%xmm2
  movsd  D92(%rip),%xmm1
  movapd %xmm1,%xmm3
L164:
  movsd  D92(%rip),%xmm6
  movsd  D65(%rip),%xmm5
L162:
  addsd  %xmm6,%xmm2
  movq   %rbx,%xmm4
  divsd  %xmm2,%xmm4
  mulsd  %xmm4,%xmm1
  addsd  %xmm1,%xmm3
  comisd %xmm5,%xmm1
  ja     L162
  movsd  0x8(%rsp),%xmm0
  divsd  (%rsp),%xmm0
  mulsd  %xmm3,%xmm0
L197:
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L5
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  ret
L121:
  movsd  0x8(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  D92(%rip),%xmm0
  subsd  %xmm1,%xmm0
  jmp    L197
L176:
  movapd %xmm5,%xmm0
L48:
  movsd  %xmm12,0x10(%rsp)
  movsd  %xmm13,0x18(%rsp)
  movsd  %xmm15,0x20(%rsp)
  movsd  %xmm14,0x28(%rsp)
  movapd %xmm15,%xmm3
  andpd  D6(%rip),%xmm3
  comisd %xmm3,%xmm11
  jbe    L42
L96:
  addsd  %xmm6,%xmm2
  addsd  %xmm9,%xmm1
  addsd  %xmm6,%xmm4
  movapd %xmm2,%xmm5
  mulsd  %xmm4,%xmm5
  movsd  0x20(%rsp),%xmm12
  movapd %xmm1,%xmm15
  mulsd  %xmm12,%xmm15
  movapd %xmm5,%xmm3
  mulsd  0x10(%rsp),%xmm3
  subsd  %xmm3,%xmm15
  movsd  %xmm15,0x30(%rsp)
  movsd  0x28(%rsp),%xmm13
  movapd %xmm13,%xmm14
  mulsd  %xmm1,%xmm14
  mulsd  0x18(%rsp),%xmm5
  subsd  %xmm5,%xmm14
  movsd  %xmm14,0x38(%rsp)
  ucomisd %xmm8,%xmm14
  jp     L4
  ucomisd D43(%rip),%xmm14
  je     L48
L4:
  movapd %xmm15,%xmm5
  divsd  %xmm14,%xmm5
  movapd %xmm0,%xmm3
  subsd  %xmm5,%xmm3
  andpd  D6(%rip),%xmm3
  comisd D65(%rip),%xmm3
  ja     L176
  movapd %xmm5,%xmm10
  mulsd  D65(%rip),%xmm10
  comisd %xmm3,%xmm10
  jae    L121
  movapd %xmm5,%xmm0
  jmp    L48
L42:
  lea    0x10(%rsp),%rax
  lea    0x30(%rsp),%rdx
L84:
  movsd  (%rax),%xmm3
  divsd  %xmm7,%xmm3
  movsd  %xmm3,(%rax)
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L84
  jmp    L96
L22:
  movsd  D107(%rip),%xmm0
  jmp    L197
L5:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_tag: D43
type_infered: f64
name: retrieve_data
arguments: data_tag: D107
type_infered: byte[8]
name: retrieve_data
arguments: data_tag: D92
type_infered: byte[16]
name: retrieve_data
arguments: data_tag: D143
type_infered: byte[8]
name: retrieve_data
arguments: data_tag: D53
type_infered: f64
name: retrieve_data
arguments: data_tag: D65
type_infered: byte[16]
name: retrieve_data
arguments: data_tag: D6
type_infered: qword[2]

tool: D43: 0

tool: {"D107": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: {"D92": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f,0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

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

tool: D53: 1e+30

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

tool: {"D6": [0x7fffffffffffffff,0x0000000000000000]}

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);
}
```