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

```asm
<atof>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  xor    %r12d,%r12d
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x28,%rsp
  jmp    L12
  nopl   (%rax)
L67:
  add    $0x1,%r12d
L12:
  movsbl (%rbx),%edi
  mov    %rbx,%r13
  add    $0x1,%rbx
  call   <isspace@plt>
  test   %rax,%rax
  jne    L67
  movzbl 0x0(%r13),%eax
  cmp    $0x2d,%al
  je     L130
  mov    $0x1,%r14d
  cmp    $0x2b,%al
  je     L181
L25:
  movslq %r12d,%r13
  movq   $0x0,(%rsp)
  add    %rbp,%r13
  mov    %r13,%r15
  jmp    L23
  cs nopw 0x0(%rax,%rax,1)
L163:
  movsbl (%r15),%eax
  movsd  D19(%rip),%xmm0
  pxor   %xmm1,%xmm1
  add    $0x1,%r15
  mulsd  (%rsp),%xmm0
  sub    $0x30,%eax
  cvtsi2sd %eax,%xmm1
  addsd  %xmm0,%xmm1
  movsd  %xmm1,(%rsp)
L23:
  mov    %r12d,%ebx
  movsbl (%r15),%edi
  sub    %r13d,%ebx
  add    %r15d,%ebx
  call   <isdigit@plt>
  test   %rax,%rax
  jne    L163
  xor    %eax,%eax
  cmpb   $0x2e,(%r15)
  sete   %al
  add    %eax,%ebx
  mov    D75(%rip),%rax
  movslq %ebx,%r13
  add    %rbp,%r13
  mov    %rax,0x8(%rsp)
  mov    %r13,%r15
  jmp    L151
  nopw   0x0(%rax,%rax,1)
L68:
  movsbl (%r15),%eax
  pxor   %xmm1,%xmm1
  add    $0x1,%r15
  movsd  D19(%rip),%xmm0
  mulsd  (%rsp),%xmm0
  movsd  D19(%rip),%xmm3
  mulsd  0x8(%rsp),%xmm3
  sub    $0x30,%eax
  cvtsi2sd %eax,%xmm1
  addsd  %xmm0,%xmm1
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm1,(%rsp)
L151:
  mov    %ebx,%r12d
  movsbl (%r15),%edi
  sub    %r13d,%r12d
  add    %r15d,%r12d
  call   <isdigit@plt>
  test   %rax,%rax
  jne    L68
  movzbl (%r15),%eax
  and    $0xffffffdf,%eax
  cmp    $0x45,%al
  jne    L50
  lea    0x1(%r12),%eax
  cltq
  movzbl 0x0(%rbp,%rax,1),%edx
  lea    -0x2b(%rdx),%ecx
  and    $0xfd,%ecx
  je     L88
L85:
  add    %rax,%rbp
  pxor   %xmm1,%xmm1
  jmp    L185
  nop
L97:
  movsbl 0x0(%rbp),%eax
  mulsd  D19(%rip),%xmm1
  pxor   %xmm0,%xmm0
  add    $0x1,%rbp
  sub    $0x30,%eax
  cvtsi2sd %eax,%xmm0
  addsd  %xmm0,%xmm1
L185:
  movsbl 0x0(%rbp),%edi
  movsd  %xmm1,0x10(%rsp)
  call   <isdigit@plt>
  movsd  0x10(%rsp),%xmm1
  test   %rax,%rax
  jne    L97
  pxor   %xmm2,%xmm2
  mov    D75(%rip),%rax
  ucomisd %xmm2,%xmm1
  movq   %rax,%xmm2
  jp     L159
  je     L89
L159:
  movsd  D19(%rip),%xmm0
  nopl   0x0(%rax)
L83:
  subsd  D75(%rip),%xmm1
  mulsd  %xmm0,%xmm2
  pxor   %xmm7,%xmm7
  ucomisd %xmm7,%xmm1
  jp     L83
  jne    L83
L89:
  pxor   %xmm0,%xmm0
  cmpl   $0x2d,0x1c(%rsp)
  cvtsi2sd %r14d,%xmm0
  mulsd  (%rsp),%xmm0
  je     L129
  mulsd  %xmm2,%xmm0
  divsd  0x8(%rsp),%xmm0
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  cs nopw 0x0(%rax,%rax,1)
L130:
  mov    $0xffffffff,%r14d
L181:
  add    $0x1,%r12d
  jmp    L25
  nop
L88:
  movsbl %dl,%eax
  mov    %eax,0x1c(%rsp)
  lea    0x2(%r12),%eax
  cltq
  jmp    L85
  nopl   0x0(%rax,%rax,1)
L129:
  divsd  %xmm2,%xmm0
  divsd  0x8(%rsp),%xmm0
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L50:
  pxor   %xmm0,%xmm0
  cvtsi2sd %r14d,%xmm0
  mulsd  (%rsp),%xmm0
  divsd  0x8(%rsp),%xmm0
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_tag: D19
type: f64
name: read_data
arguments: data_tag: D75
type: f64

tool: D19: 10

tool: D75: 1

assistant: ```c
double atof(char s[]) {
  double val, power, pot, e = 1.0;
  int i, sign, sign_pot;
  for (i = 0; isspace(s[i]); i++)
    ;
  sign = (s[i] == '-') ? -1 : 1;
  if (s[i] == '+' || s[i] == '-')
    i++;
  for (val = 0.0; isdigit(s[i]); i++)
    val = 10.0 * val + (s[i] - '0');
  if (s[i] == '.')
    i++;
  for (power = 1.0; isdigit(s[i]); i++) {
    val = 10.0 * val + (s[i] - '0');
    power *= 10;
  }
  if (s[i] == 'e' || s[i] == 'E')
    i++;
  else
    return sign * val / power;
  if (s[i] == '-' || s[i] == '+')
    sign_pot = s[i++];
  for (pot = 0.0; isdigit(s[i]); i++)
    pot = 10 * pot + (s[i] - '0');
  while (pot != 0) {
    e *= 10;
    pot--;
  }
  if (sign_pot == '-')
    return sign * val / e / power;
  else
    return sign * val * e / power;
}
```