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

```asm
<IncompleteGamma>:
  endbr64
  pxor   %xmm9,%xmm9
  movapd %xmm0,%xmm3
  ucomisd %xmm9,%xmm0
  jp     L83
  movapd %xmm9,%xmm0
  je     L137
L83:
  sub    $0x28,%rsp
  comisd %xmm3,%xmm9
  movsd  D23(%rip),%xmm0
  movsd  %xmm2,(%rsp)
  ja     L126
  comisd %xmm1,%xmm9
  movsd  %xmm1,0x8(%rsp)
  jae    L126
  movapd %xmm3,%xmm0
  movsd  %xmm3,0x18(%rsp)
  call   <log@plt>
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm1
  movsd  0x18(%rsp),%xmm3
  cvtsi2sd %rax,%xmm0
  movsd  (%rsp),%xmm2
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm3,0x8(%rsp)
  mulsd  %xmm1,%xmm0
  subsd  %xmm3,%xmm0
  subsd  %xmm2,%xmm0
  call   <exp@plt>
  movsd  0x8(%rsp),%xmm3
  pxor   %xmm7,%xmm7
  movsd  D63(%rip),%xmm0
  cvtsi2sd %rax,%xmm7
  movsd  0x10(%rsp),%xmm1
  comisd %xmm0,%xmm3
  movsd  %xmm7,(%rsp)
  jbe    L41
  comisd %xmm1,%xmm3
  jb     L41
  movapd %xmm0,%xmm7
  movapd %xmm3,%xmm4
  movapd %xmm3,%xmm5
  movsd  D41(%rip),%xmm11
  subsd  %xmm1,%xmm7
  movapd %xmm3,%xmm6
  pxor   %xmm9,%xmm9
  movq   D34(%rip),%xmm12
  addsd  %xmm0,%xmm6
  movapd %xmm0,%xmm1
  movapd %xmm9,%xmm8
  movsd  D35(%rip),%xmm10
  addsd  %xmm7,%xmm4
  movapd %xmm6,%xmm13
  addsd  %xmm0,%xmm4
  mulsd  %xmm4,%xmm5
  divsd  %xmm5,%xmm13
  jmp    L24
  nopl   0x0(%rax)
L101:
  movapd %xmm5,%xmm3
  movapd %xmm1,%xmm5
  movapd %xmm6,%xmm1
  movapd %xmm15,%xmm6
L24:
  addsd  %xmm0,%xmm7
  addsd  %xmm0,%xmm8
  addsd  D37(%rip),%xmm4
  movapd %xmm7,%xmm2
  mulsd  %xmm8,%xmm2
  movapd %xmm4,%xmm15
  mulsd  %xmm6,%xmm15
  mulsd  %xmm2,%xmm1
  mulsd  %xmm3,%xmm2
  subsd  %xmm1,%xmm15
  movapd %xmm4,%xmm1
  mulsd  %xmm5,%xmm1
  subsd  %xmm2,%xmm1
  ucomisd %xmm9,%xmm1
  jp     L91
  je     L77
L91:
  movapd %xmm15,%xmm3
  movapd %xmm13,%xmm2
  divsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  andpd  %xmm12,%xmm2
  comisd %xmm11,%xmm2
  ja     L48
  movapd %xmm3,%xmm14
  mulsd  %xmm11,%xmm14
  comisd %xmm2,%xmm14
  jae    L90
L48:
  movapd %xmm3,%xmm13
L77:
  movapd %xmm15,%xmm2
  andpd  %xmm12,%xmm2
  comisd %xmm2,%xmm10
  ja     L101
  divsd  %xmm10,%xmm6
  divsd  %xmm10,%xmm5
  divsd  %xmm10,%xmm15
  divsd  %xmm10,%xmm1
  jmp    L101
  nopl   (%rax)
L137:
  ret
  nopl   0x0(%rax)
L41:
  movapd %xmm1,%xmm4
  movapd %xmm0,%xmm2
  movapd %xmm0,%xmm5
  movsd  D41(%rip),%xmm11
  nopl   (%rax)
L37:
  addsd  %xmm0,%xmm4
  movapd %xmm3,%xmm6
  divsd  %xmm4,%xmm6
  mulsd  %xmm6,%xmm2
  comisd %xmm11,%xmm2
  addsd  %xmm2,%xmm5
  ja     L37
  movsd  (%rsp),%xmm0
  divsd  %xmm1,%xmm0
  mulsd  %xmm5,%xmm0
L126:
  add    $0x28,%rsp
  ret
  nopl   (%rax)
L90:
  mulsd  (%rsp),%xmm13
  add    $0x28,%rsp
  subsd  %xmm13,%xmm0
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D23
type_guessed: f32[2]
name: retrieve_data
arguments: tag: D63
type_guessed: f32[2]
name: retrieve_data
arguments: tag: D41
type_guessed: byte[8]
name: retrieve_data
arguments: tag: D34
type_guessed: qword
name: retrieve_data
arguments: tag: D35
type_guessed: f64
name: retrieve_data
arguments: tag: D37
type_guessed: f32[2]

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

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

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

tool: D34: 9223372036854775807

tool: D35: 1e+30

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

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