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

```asm
<cdflib_fpser>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %rcx,-0x50(%rbp)
  movsd  D0(%rip),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm2
  movsd  D20(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jae    L39
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm3
  movsd  %xmm3,-0x58(%rbp)
  mov    -0x48(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <log@plt>
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  lea    D11(%rip),%rax
  mov    %rax,%rdi
  call   <cdflib_exparg@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm4
  comisd -0x18(%rbp),%xmm4
  jbe    L7
  movsd  -0x30(%rbp),%xmm0
  jmp    L42
L7:
  mov    -0x18(%rbp),%rax
  movq   %rax,%xmm0
  call   <exp@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x30(%rbp)
  jmp    L3
L39:
  nop
L3:
  mov    -0x40(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x30(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D0(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0x48(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x18(%rbp),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
L18:
  movsd  -0x28(%rbp),%xmm1
  movsd  D0(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0x48(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  -0x18(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x18(%rbp),%xmm0
  divsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x20(%rbp),%xmm0
  addsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movq   D69(%rip),%xmm1
  andpd  %xmm1,%xmm0
  comisd -0x10(%rbp),%xmm0
  jbe    L43
  jmp    L18
L43:
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  movsd  D0(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x30(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x30(%rbp),%xmm0
L42:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_tag: D0
data_type: f64
name: parse_data
arguments: data_tag: D20
data_type: byte[8]
name: parse_data
arguments: data_tag: D11
data_type: string
name: parse_data
arguments: data_tag: D69
data_type: qword

tool: D0: 1

tool: {"D20": [0xfc,0xa9,0xf1,0xd2,0x4d,0x62,0x50,0x3f]}

tool: D11: 

tool: D69: 9223372036854775807

assistant: ```c
double cdflib_fpser(double *a, double *b, double *x, double *eps) {
  static const int K1 = 1;
  double fpser, an, c, s, t, tol;
  fpser = 1.0e0;
  if (*a <= 1.e-3 * *eps) {
    goto S10;
  }
  fpser = 0.0e0;
  t = *a * log(*x);
  if (t < cdflib_exparg(&K1)) {
    return fpser;
  }
  fpser = exp(t);
S10:
  fpser = *b / *a * fpser;
  tol = *eps / *a;
  an = *a + 1.0e0;
  t = *x;
  s = t / an;
S20:
  an += 1.0e0;
  t = *x * t;
  c = t / an;
  s += c;
  if (fabs(c) > tol) {
    goto S20;
  }
  fpser *= (1.0e0 + *a * s);
  return fpser;
}
```