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

```asm
<expint>:
  endbr64
  push   %r14
  push   %rbx
  sub    $0x38,%rsp
  movapd %xmm0,%xmm7
  mov    %fs:0x28,%rax
  mov    %rax,0x28(%rsp)
  xor    %eax,%eax
  pxor   %xmm0,%xmm0
  comisd %xmm7,%xmm0
  ja     L156
  movsd  %xmm1,0x8(%rsp)
  comisd %xmm1,%xmm0
  ja     L156
  lea    0x20(%rsp),%rdi
  movsd  %xmm7,0x10(%rsp)
  movapd %xmm7,%xmm0
  call   <modf@plt>
  test   %rax,%rax
  jne    L156
  movsd  0x8(%rsp),%xmm1
  ucomisd D130(%rip),%xmm1
  jp     L15
  jne    L15
  movsd  D160(%rip),%xmm0
  comisd 0x10(%rsp),%xmm0
  ja     L156
  movsd  0x10(%rsp),%xmm2
  ucomisd D130(%rip),%xmm2
  jp     L69
  je     L51
L69:
  movsd  D87(%rip),%xmm0
  movsd  0x10(%rsp),%xmm1
  subsd  %xmm0,%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rbx
  jmp    L60
L156:
  mov    D30(%rip),%rbx
L60:
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L259
  movq   %rbx,%xmm0
  add    $0x38,%rsp
  pop    %rbx
  pop    %r14
  ret
L51:
  movsd  0x8(%rsp),%xmm0
  xorpd  D113(%rip),%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm2
  divsd  0x8(%rsp),%xmm2
  movq   %xmm2,%rbx
  jmp    L60
L83:
  cmp    $0x14d,%eax
  je     L26
L260:
  pxor   %xmm4,%xmm4
  cvtsi2sd %eax,%xmm4
  addsd  0x10(%rsp),%xmm4
  mulsd  %xmm4,%xmm2
  addsd  %xmm1,%xmm2
  mulsd  %xmm4,%xmm3
  addsd  %xmm0,%xmm3
  add    $0x1,%eax
  pxor   %xmm4,%xmm4
  cvtsi2sd %eax,%xmm4
  mulsd  %xmm4,%xmm1
  movsd  0x8(%rsp),%xmm8
  movapd %xmm8,%xmm5
  mulsd  %xmm2,%xmm5
  addsd  %xmm5,%xmm1
  mulsd  %xmm4,%xmm0
  movapd %xmm8,%xmm4
  mulsd  %xmm3,%xmm4
  addsd  %xmm4,%xmm0
  movq   %rbx,%xmm4
  movapd %xmm1,%xmm5
  divsd  %xmm0,%xmm5
  movq   %xmm5,%rbx
  ucomisd %xmm4,%xmm5
  jp     L7
  je     L26
L7:
  comisd %xmm7,%xmm0
  jb     L83
  divsd  %xmm6,%xmm2
  divsd  %xmm6,%xmm1
  divsd  %xmm6,%xmm3
  divsd  %xmm6,%xmm0
  jmp    L83
L26:
  movsd  0x8(%rsp),%xmm0
  xorpd  D113(%rip),%xmm0
  call   <exp@plt>
  movq   %rbx,%xmm6
  mulsd  %xmm0,%xmm6
  movq   %xmm6,%rbx
  jmp    L60
L189:
  movsd  0x8(%rsp),%xmm0
  call   <log@plt>
  movsd  D6(%rip),%xmm1
  subsd  %xmm0,%xmm1
  movq   %xmm1,%rbx
  movsd  D87(%rip),%xmm1
  movapd %xmm1,%xmm0
  pxor   %xmm2,%xmm2
  movq   D113(%rip),%xmm4
  movapd %xmm1,%xmm6
  movsd  D71(%rip),%xmm5
  jmp    L44
L131:
  movq   %rbx,%xmm2
L44:
  xorpd  %xmm4,%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  %xmm1,%xmm0
  movapd %xmm0,%xmm3
  divsd  %xmm1,%xmm3
  movq   %rbx,%xmm7
  subsd  %xmm3,%xmm7
  movq   %xmm7,%rbx
  ucomisd %xmm2,%xmm7
  jp     L199
  je     L270
L199:
  addsd  %xmm6,%xmm1
  comisd %xmm1,%xmm5
  ja     L131
L270:
  movsd  D87(%rip),%xmm6
  movsd  %xmm6,0x18(%rsp)
  movsd  0x8(%rsp),%xmm6
  xorpd  D113(%rip),%xmm6
  movq   %xmm6,%r14
  movsd  0x10(%rsp),%xmm6
  comisd D87(%rip),%xmm6
  jbe    L60
L70:
  movq   %r14,%xmm0
  call   <exp@plt>
  movq   %rbx,%xmm1
  mulsd  0x8(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  0x18(%rsp),%xmm6
  divsd  %xmm6,%xmm0
  movq   %xmm0,%rbx
  addsd  D87(%rip),%xmm6
  movsd  %xmm6,0x18(%rsp)
  movsd  0x10(%rsp),%xmm2
  comisd %xmm6,%xmm2
  ja     L70
  jmp    L60
L15:
  movsd  0x10(%rsp),%xmm6
  ucomisd D130(%rip),%xmm6
  jp     L181
  je     L51
L181:
  movsd  0x8(%rsp),%xmm6
  comisd D120(%rip),%xmm6
  jbe    L189
  movapd %xmm6,%xmm0
  mov    $0x0,%eax
  mov    D87(%rip),%rbx
  movq   %rbx,%xmm3
  movq   %rbx,%xmm1
  pxor   %xmm2,%xmm2
  movsd  D136(%rip),%xmm7
  movapd %xmm7,%xmm6
  jmp    L260
L259:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_label: D130
type_infered: f64
name: unpack_data
arguments: data_label: D160
type_infered: byte[8]
name: unpack_data
arguments: data_label: D87
type_infered: f64[2]
name: unpack_data
arguments: data_label: D30
type_infered: byte[8]
name: unpack_data
arguments: data_label: D113
type_infered: byte[16]
name: unpack_data
arguments: data_label: D6
type_infered: byte[8]
name: unpack_data
arguments: data_label: D71
type_infered: byte[8]
name: unpack_data
arguments: data_label: D120
type_infered: byte[16]
name: unpack_data
arguments: data_label: D136
type_infered: byte[8]

tool: D130: 0

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

tool: {"D87": [1.0,-1.0]}

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

tool: {"D113": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: {"D6": [0x19,0xb6,0x6f,0xfc,0x8c,0x78,0xe2,0xbf]}

tool: {"D71": [0x00,0x00,0x00,0x00,0x00,0xd0,0x74,0x40]}

tool: {"D120": [0x00,0x00,0x00,0x00,0x00,0x00,0x08,0x40,0x1b,0x69,0x57,0x43,0xb8,0x17,0x9e,0x47]}

tool: {"D136": [0x1b,0x69,0x57,0x43,0xb8,0x17,0x9e,0x47]}

assistant: ```c
double expint(double n, double z) {
  double y;
  {
    double junk;
    if (n < 0 || z < 0 || modf(n, &junk))
      return -1.0;
    if (z == 0 && n < 2)
      return -1.0;
  }
  if (n == 0)
    return exp(-z) / z;
  if (z == 0)
    return 1 / (n - 1);
  if (z > 3) {
    double n0, n1, n2, n3, d0, d1, d2, d3, y_prev = 1;
    int i;
    n0 = 0;
    n1 = 1;
    d0 = 1;
    d1 = z;
    for (i = 0; i < 333; i++) {
      n2 = n0 * (n + i) + n1;
      d2 = d0 * (n + i) + d1;
      n3 = n1 * (1 + i) + n2 * z;
      d3 = d1 * (1 + i) + d2 * z;
      y = n3 / d3;
      if (y == y_prev)
        break;
      y_prev = y;
      n0 = n2;
      n1 = n3;
      d0 = d2;
      d1 = d3;
      if (d3 >= 1.0E+37) {
        n0 /= 1.0E+37;
        n1 /= 1.0E+37;
        d0 /= 1.0E+37;
        d1 /= 1.0E+37;
      }
    }
    y = exp(-z) * y;
  } else {
    double euler_constant = .577215664901532860606512;
    double y_prev = 0;
    double t, m;
    y = -euler_constant - log(z);
    t = 1;
    for (m = 1; m < 333; m++) {
      t = -t * z / m;
      y = y - t / m;
      if (y == y_prev)
        break;
      y_prev = y;
    }
    for (m = 1; m < n; m++)
      y = (exp(-z) - z * y) / m;
  }
  return y;
}
```