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

```asm
<cephes_igamcf>:
  endbr64
  sub    $0x38,%rsp
  movss  %xmm0,0x4(%rsp)
  movaps %xmm1,%xmm0
  movss  %xmm1,0x8(%rsp)
  call   <log@plt>
  movss  0x4(%rsp),%xmm7
  movss  0x8(%rsp),%xmm1
  mulss  %xmm7,%xmm0
  movss  %xmm1,0xc(%rsp)
  movaps %xmm0,%xmm2
  movaps %xmm7,%xmm0
  subss  %xmm1,%xmm2
  movss  %xmm2,0x8(%rsp)
  call   <cephes_lgamf@plt>
  movss  0x8(%rsp),%xmm2
  subss  %xmm0,%xmm2
  movss  D13(%rip),%xmm0
  comiss %xmm2,%xmm0
  jbe    L37
  pxor   %xmm0,%xmm0
  add    $0x38,%rsp
  ret
  nopw   0x0(%rax,%rax,1)
L37:
  movaps %xmm2,%xmm0
  call   <expf@plt>
  movss  0xc(%rsp),%xmm1
  pxor   %xmm8,%xmm8
  movss  D54(%rip),%xmm9
  movss  %xmm0,0x2c(%rsp)
  movaps %xmm9,%xmm5
  subss  0x4(%rsp),%xmm5
  movaps %xmm1,%xmm4
  movaps %xmm1,%xmm6
  addss  %xmm9,%xmm6
  addss  %xmm5,%xmm4
  movaps %xmm5,%xmm7
  movaps %xmm1,%xmm5
  movaps %xmm6,%xmm3
  addss  %xmm9,%xmm4
  mulss  %xmm4,%xmm5
  divss  %xmm5,%xmm3
  movss  %xmm3,0x20(%rsp)
  movaps %xmm9,%xmm3
  nopl   0x0(%rax,%rax,1)
L62:
  addss  %xmm9,%xmm8
  addss  %xmm9,%xmm7
  movaps %xmm5,%xmm10
  addss  D113(%rip),%xmm4
  movaps %xmm6,%xmm2
  movaps %xmm8,%xmm0
  mulss  %xmm7,%xmm0
  mulss  %xmm4,%xmm10
  mulss  %xmm4,%xmm2
  mulss  %xmm0,%xmm3
  mulss  %xmm1,%xmm0
  subss  %xmm3,%xmm2
  pxor   %xmm3,%xmm3
  subss  %xmm0,%xmm10
  ucomiss %xmm3,%xmm10
  jp     L38
  jne    L38
  movaps %xmm2,%xmm0
  movss  %xmm10,0x1c(%rsp)
  movss  %xmm5,0x18(%rsp)
  movss  %xmm6,0x14(%rsp)
  movss  %xmm4,0x10(%rsp)
  movss  %xmm7,0xc(%rsp)
  movss  %xmm8,0x8(%rsp)
  movss  %xmm2,0x4(%rsp)
  call   <fabs@plt>
  comiss D79(%rip),%xmm0
  movss  0x4(%rsp),%xmm2
  movss  0x8(%rsp),%xmm8
  movss  0xc(%rsp),%xmm7
  movss  0x10(%rsp),%xmm4
  movss  0x14(%rsp),%xmm6
  movss  0x18(%rsp),%xmm5
  movss  0x1c(%rsp),%xmm10
  movss  D54(%rip),%xmm9
  ja     L10
  movaps %xmm5,%xmm1
  movaps %xmm6,%xmm3
  movaps %xmm10,%xmm5
  movaps %xmm2,%xmm6
  jmp    L62
  nopl   0x0(%rax,%rax,1)
L10:
  movaps %xmm9,%xmm11
L17:
  movss  D29(%rip),%xmm3
  movss  D29(%rip),%xmm1
  mulss  %xmm6,%xmm3
  movss  D29(%rip),%xmm6
  mulss  %xmm5,%xmm1
  movss  D29(%rip),%xmm5
  mulss  %xmm2,%xmm6
  mulss  %xmm10,%xmm5
L66:
  comiss D29(%rip),%xmm11
  ja     L62
  movss  0x2c(%rsp),%xmm0
  mulss  0x20(%rsp),%xmm0
  add    $0x38,%rsp
  ret
  nopl   0x0(%rax)
L38:
  movaps %xmm2,%xmm1
  movss  0x20(%rsp),%xmm0
  movss  %xmm5,0x28(%rsp)
  divss  %xmm10,%xmm1
  movss  %xmm6,0x24(%rsp)
  movss  %xmm4,0x1c(%rsp)
  movss  %xmm7,0x18(%rsp)
  movss  %xmm8,0x14(%rsp)
  movss  %xmm10,0x10(%rsp)
  movss  %xmm2,0x4(%rsp)
  subss  %xmm1,%xmm0
  movss  %xmm1,0xc(%rsp)
  divss  %xmm1,%xmm0
  call   <fabs@plt>
  movss  0x4(%rsp),%xmm2
  movss  %xmm0,0x8(%rsp)
  movaps %xmm2,%xmm0
  call   <fabs@plt>
  comiss D79(%rip),%xmm0
  movss  0xc(%rsp),%xmm1
  movss  0x4(%rsp),%xmm2
  movss  0x8(%rsp),%xmm11
  movss  0x10(%rsp),%xmm10
  movss  0x14(%rsp),%xmm8
  movss  %xmm1,0x20(%rsp)
  movss  0x18(%rsp),%xmm7
  movss  0x1c(%rsp),%xmm4
  movss  0x24(%rsp),%xmm6
  movss  0x28(%rsp),%xmm5
  movss  D54(%rip),%xmm9
  ja     L17
  movaps %xmm5,%xmm1
  movaps %xmm6,%xmm3
  movaps %xmm10,%xmm5
  movaps %xmm2,%xmm6
  jmp    L66
```

tool_calls: name: read_data
arguments: data_label: D13
type_guessed: byte[4]
name: read_data
arguments: data_label: D54
type_guessed: f32
name: read_data
arguments: data_label: D113
type_guessed: f32
name: read_data
arguments: data_label: D79
type_guessed: f32
name: read_data
arguments: data_label: D29
type_guessed: byte[4]

tool: {"D13": [0x18,0x72,0xb1,0xc2]}

tool: D54: 1

tool: D113: 2

tool: D79: 16777216

tool: {"D29": [0x00,0x00,0x80,0x33]}

assistant: ```c
float cephes_igamcf(float a, float x) {
  const float MAXLOGF = 88.72283905206835;
  const float MACHEPF = 5.9604644775390625E-8;
  const float BIG = 16777216.;
  float ans, c, yc, ax, y, z;
  float pk, pkm1, pkm2, qk, qkm1, qkm2;
  float r, t;
  ax = a * log(x) - x - cephes_lgamf(a);
  if (ax < -MAXLOGF)
    return 0.0;
  ax = expf(ax);
  y = 1.0 - a;
  z = x + y + 1.0;
  c = 0.0;
  pkm2 = 1.0;
  qkm2 = x;
  pkm1 = x + 1.0;
  qkm1 = z * x;
  ans = pkm1 / qkm1;
  do {
    c += 1.0;
    y += 1.0;
    z += 2.0;
    yc = y * c;
    pk = pkm1 * z - pkm2 * yc;
    qk = qkm1 * z - qkm2 * yc;
    if (qk != 0) {
      r = pk / qk;
      t = fabs((ans - r) / r);
      ans = r;
    } else {
      t = 1.0;
    }
    pkm2 = pkm1;
    pkm1 = pk;
    qkm2 = qkm1;
    qkm1 = qk;
    if (fabs(pk) > BIG) {
      pkm2 *= MACHEPF;
      pkm1 *= MACHEPF;
      qkm2 *= MACHEPF;
      qkm1 *= MACHEPF;
    }
  } while (t > MACHEPF);
  return ans * ax;
}
```