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

```asm
<kepler>:
  endbr64
  sub    $0x58,%rsp
  movsd  %xmm1,0x30(%rsp)
  call   <rev180@plt>
  lea    0x48(%rsp),%rdi
  lea    0x40(%rsp),%rsi
  movsd  D92(%rip),%xmm5
  movsd  %xmm0,0x20(%rsp)
  mulsd  %xmm0,%xmm5
  movapd %xmm5,%xmm0
  call   <sincos@plt>
  movsd  0x40(%rsp),%xmm1
  subsd  0x30(%rsp),%xmm1
  movsd  0x48(%rsp),%xmm0
  call   <atan2@plt>
  movsd  0x30(%rsp),%xmm7
  comisd D99(%rip),%xmm7
  movapd %xmm0,%xmm2
  movsd  D31(%rip),%xmm0
  mulsd  %xmm0,%xmm2
  jbe    L33
  mulsd  %xmm7,%xmm0
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm0,0x28(%rsp)
  movsd  D75(%rip),%xmm0
  divsd  %xmm7,%xmm0
  movsd  %xmm0,0x38(%rsp)
  movsd  D92(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  call   <cos@plt>
  mulsd  0x30(%rsp),%xmm0
  movsd  D39(%rip),%xmm4
  movsd  D70(%rip),%xmm3
  movsd  0x8(%rsp),%xmm2
  movsd  %xmm3,0x18(%rsp)
  subsd  %xmm0,%xmm4
  movsd  %xmm4,0x10(%rsp)
  nopw   0x0(%rax,%rax,1)
L4:
  movsd  D92(%rip),%xmm0
  movsd  %xmm2,0x8(%rsp)
  mulsd  %xmm2,%xmm0
  call   <sin@plt>
  movsd  0x28(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm2
  movsd  D9(%rip),%xmm6
  mulsd  %xmm0,%xmm1
  addsd  0x20(%rsp),%xmm1
  subsd  %xmm2,%xmm1
  divsd  0x10(%rsp),%xmm1
  addsd  %xmm1,%xmm2
  andpd  D51(%rip),%xmm1
  comisd %xmm1,%xmm6
  ja     L33
  comisd 0x18(%rsp),%xmm1
  jae    L33
  comisd 0x38(%rsp),%xmm1
  jbe    L27
  movsd  D92(%rip),%xmm0
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm2,0x8(%rsp)
  mulsd  %xmm2,%xmm0
  call   <cos@plt>
  mulsd  0x30(%rsp),%xmm0
  movsd  D39(%rip),%xmm4
  movsd  0x8(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm1
  subsd  %xmm0,%xmm4
  movsd  %xmm4,0x10(%rsp)
L27:
  movsd  %xmm1,0x18(%rsp)
  jmp    L4
  nopl   (%rax)
L33:
  movapd %xmm2,%xmm0
  add    $0x58,%rsp
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D92
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D99
type_guessed: f64[2]
name: unpack_data
arguments: data_label: D31
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D75
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D39
type_guessed: f64
name: unpack_data
arguments: data_label: D70
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D9
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D51
type_guessed: qword[2]

tool: {"D92": [0x39,0x9d,0x52,0xa2,0x46,0xdf,0x91,0x3f]}

tool: {"D99": [0.008,0.001]}

tool: {"D31": [0xf8,0xc1,0x63,0x1a,0xdc,0xa5,0x4c,0x40]}

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

tool: D39: 1

tool: {"D70": [0x00,0x00,0x00,0x20,0x5f,0xa0,0x02,0x42]}

tool: {"D9": [0x2b,0x69,0xa4,0x29,0x2b,0x1b,0x60,0x3e]}

tool: {"D51": [0x7fffffffffffffff,0x0000000000000000]}

assistant: ```c
double kepler(double m, double ex) {
  double m1, sinm, cosm, exd, exan, dexan, lim1, adko, adk, denom;
  int converged;
  m1 = rev180(m);
  sinm = sin(m1 * (3.14159265358979323846 / 180.0));
  cosm = cos(m1 * (3.14159265358979323846 / 180.0));
  exan = ((180.0 / 3.14159265358979323846) * atan2(sinm, cosm - ex));
  if (ex > 0.008) {
    exd = ex * (180.0 / 3.14159265358979323846);
    lim1 = 1E-3 / ex;
    adko = 1E+10;
    denom = 1.0 - ex * cos(exan * (3.14159265358979323846 / 180.0));
    do {
      dexan = (m1 + exd * sin(exan * (3.14159265358979323846 / 180.0)) - exan) / denom;
      exan = exan + dexan;
      adk = fabs(dexan);
      converged = adk < 3E-8 || adk >= adko;
      adko = adk;
      if (!converged && adk > lim1)
        denom = 1.0 - ex * cos(exan * (3.14159265358979323846 / 180.0));
    } while (!converged);
  }
  return exan;
}
```