user: Converte this assembly snippet into its equivalent C code:

```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  D45(%rip),%xmm0
  comiss %xmm2,%xmm0
  jbe    L13
  pxor   %xmm0,%xmm0
  add    $0x38,%rsp
  ret
  nopw   0x0(%rax,%rax,1)
L13:
  movaps %xmm2,%xmm0
  call   <expf@plt>
  movss  0xc(%rsp),%xmm1
  pxor   %xmm8,%xmm8
  movss  D16(%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)
L71:
  addss  %xmm9,%xmm8
  addss  %xmm9,%xmm7
  movaps %xmm5,%xmm10
  addss  D34(%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     L64
  jne    L64
  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 D92(%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  D16(%rip),%xmm9
  ja     L8
  movaps %xmm5,%xmm1
  movaps %xmm6,%xmm3
  movaps %xmm10,%xmm5
  movaps %xmm2,%xmm6
  jmp    L71
  nopl   0x0(%rax,%rax,1)
L8:
  movaps %xmm9,%xmm11
L12:
  movss  D94(%rip),%xmm3
  movss  D94(%rip),%xmm1
  mulss  %xmm6,%xmm3
  movss  D94(%rip),%xmm6
  mulss  %xmm5,%xmm1
  movss  D94(%rip),%xmm5
  mulss  %xmm2,%xmm6
  mulss  %xmm10,%xmm5
L26:
  comiss D94(%rip),%xmm11
  ja     L71
  movss  0x2c(%rsp),%xmm0
  mulss  0x20(%rsp),%xmm0
  add    $0x38,%rsp
  ret
  nopl   0x0(%rax)
L64:
  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 D92(%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  D16(%rip),%xmm9
  ja     L12
  movaps %xmm5,%xmm1
  movaps %xmm6,%xmm3
  movaps %xmm10,%xmm5
  movaps %xmm2,%xmm6
  jmp    L26
```

tool_calls: name: retrieve_data
arguments: label: D45
type: byte[4]
name: retrieve_data
arguments: label: D16
type: f32
name: retrieve_data
arguments: label: D34
type: f32
name: retrieve_data
arguments: label: D92
type: f32
name: retrieve_data
arguments: label: D94
type: byte[4]

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

tool: D16: 1

tool: D34: 2

tool: D92: 16777216

tool: {"D94": [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;
}
```