user: Please show me the C implementation that matches the assembly code listed below:

```asm
<cdrift_kepu_guess>:
  endbr64
  sub    $0x68,%rsp
  pxor   %xmm7,%xmm7
  movapd %xmm0,%xmm6
  movsd  %xmm5,0x38(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  comisd %xmm7,%xmm3
  jbe    L86
  movapd %xmm0,%xmm5
  movsd  D36(%rip),%xmm0
  divsd  %xmm1,%xmm5
  comisd %xmm5,%xmm0
  jb     L11
  mulsd  %xmm6,%xmm6
  movapd %xmm1,%xmm0
  addsd  %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm4,%xmm6
  mulsd  %xmm1,%xmm0
  divsd  %xmm0,%xmm6
  subsd  %xmm6,%xmm5
  movapd %xmm5,%xmm0
L137:
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L112
  add    $0x68,%rsp
  ret
  nopl   0x0(%rax,%rax,1)
L11:
  movapd %xmm2,%xmm8
  divsd  %xmm3,%xmm8
  movapd %xmm8,%xmm0
  mulsd  %xmm8,%xmm0
  mulsd  %xmm8,%xmm0
  divsd  %xmm0,%xmm2
  ucomisd %xmm2,%xmm7
  ja     L7
  sqrtsd %xmm2,%xmm2
L119:
  movapd %xmm8,%xmm0
  movsd  D8(%rip),%xmm5
  mulsd  %xmm2,%xmm0
  divsd  %xmm8,%xmm1
  mulsd  %xmm8,%xmm0
  divsd  %xmm0,%xmm4
  subsd  %xmm1,%xmm5
  movapd %xmm5,%xmm0
  mulsd  %xmm5,%xmm0
  movapd %xmm4,%xmm1
  mulsd  %xmm4,%xmm1
  addsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm7
  ja     L60
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x8(%rsp)
L14:
  mulsd  %xmm6,%xmm2
  lea    0x50(%rsp),%rdx
  lea    0x48(%rsp),%rsi
  movsd  %xmm3,0x20(%rsp)
  lea    0x40(%rsp),%rdi
  movsd  %xmm5,0x18(%rsp)
  movsd  %xmm4,0x10(%rsp)
  movq   $0x0,0x48(%rsp)
  movq   $0x0,0x50(%rsp)
  subsd  %xmm4,%xmm2
  movsd  %xmm2,0x40(%rsp)
  call   <cmco_sine@plt>
  movsd  0x10(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm5
  pxor   %xmm7,%xmm7
  mulsd  0x50(%rsp),%xmm4
  movsd  D31(%rip),%xmm0
  mulsd  0x48(%rsp),%xmm5
  movsd  0x20(%rsp),%xmm3
  addsd  %xmm5,%xmm4
  comisd %xmm7,%xmm4
  jae    L129
  movsd  D22(%rip),%xmm0
L129:
  mulsd  0x8(%rsp),%xmm0
  ucomisd %xmm3,%xmm7
  movsd  0x40(%rsp),%xmm1
  addsd  %xmm0,%xmm1
  ja     L59
  sqrtsd %xmm3,%xmm3
L87:
  divsd  %xmm3,%xmm1
  movapd %xmm1,%xmm0
  jmp    L137
  nopl   0x0(%rax,%rax,1)
L86:
  lea    0x38(%rsp),%rdi
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm0,0x8(%rsp)
  call   <cdrift_kepu_p3solve@plt>
  movsd  0x8(%rsp),%xmm6
  movsd  0x10(%rsp),%xmm1
  test   %eax,%eax
  je     L104
  divsd  %xmm1,%xmm6
  movapd %xmm6,%xmm0
  jmp    L137
  nopl   0x0(%rax,%rax,1)
L104:
  movsd  0x38(%rsp),%xmm0
  jmp    L137
L112:
  call   <__stack_chk_fail@plt>
L60:
  movsd  %xmm3,0x30(%rsp)
  movsd  %xmm6,0x28(%rsp)
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm5,0x18(%rsp)
  movsd  %xmm2,0x10(%rsp)
  call   <sqrt@plt>
  movsd  0x30(%rsp),%xmm3
  movsd  0x28(%rsp),%xmm6
  movsd  0x20(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm5
  movsd  %xmm0,0x8(%rsp)
  movsd  0x10(%rsp),%xmm2
  jmp    L14
L59:
  movapd %xmm3,%xmm0
  movsd  %xmm1,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm1
  movapd %xmm0,%xmm3
  jmp    L87
L7:
  movapd %xmm2,%xmm0
  movsd  %xmm4,0x28(%rsp)
  movsd  %xmm3,0x20(%rsp)
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm6,0x10(%rsp)
  movsd  %xmm8,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x28(%rsp),%xmm4
  pxor   %xmm7,%xmm7
  movsd  0x20(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm6
  movapd %xmm0,%xmm2
  movsd  0x8(%rsp),%xmm8
  jmp    L119
```

tool_calls: name: unpack_data
arguments: tag: D36
type_infered: f64
name: unpack_data
arguments: tag: D8
type_infered: f64
name: unpack_data
arguments: tag: D31
type_infered: f64
name: unpack_data
arguments: tag: D22
type_infered: byte[8]

tool: D36: 0.4

tool: D8: 1

tool: D31: 0.85

tool: {"D22": [0x33,0x33,0x33,0x33,0x33,0x33,0xeb,0xbf]}

assistant: ```c
double cdrift_kepu_guess(const double dt, const double r0, const double mu, const double alpha, const double u, double s) {
  if (alpha > 0.0) {
    if (dt / r0 <= 0.4) {
      s = dt / r0 - (dt * dt * u) / (2.0 * r0 * r0 * r0);
      return s;
    } else {
      const double a = mu / alpha;
      const double en = sqrt(mu / (a * a * a));
      const double ec = 1.0 - r0 / a;
      const double es = u / (en * a * a);
      const double e = sqrt(ec * ec + es * es);
      double y = en * dt - es;
      double sy = 0.0, cy = 0.0;
      cmco_sine(&y, &sy, &cy);
      const double sigma = (es * cy + ec * sy) >= 0 ? fabs(1.0) : -fabs(1.0);
      const double x = y + sigma * 0.85 * e;
      s = x / sqrt(alpha);
    }
  } else {
    const int iflg = cdrift_kepu_p3solve(dt, r0, mu, alpha, u, &s);
    if (iflg) {
      s = dt / r0;
    }
  }
  return s;
}
```