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

```asm
<calcEquationOfTime>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  movsd  %xmm0,-0x68(%rbp)
  mov    -0x68(%rbp),%rax
  movq   %rax,%xmm0
  call   <calcObliquityCorrection@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x58(%rbp)
  mov    -0x68(%rbp),%rax
  movq   %rax,%xmm0
  call   <calcGeomMeanLongSun@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x50(%rbp)
  mov    -0x68(%rbp),%rax
  movq   %rax,%xmm0
  call   <calcEccentricityEarthOrbit@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x48(%rbp)
  mov    -0x68(%rbp),%rax
  movq   %rax,%xmm0
  call   <calcGeomMeanAnomalySun@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x58(%rbp),%rax
  movq   %rax,%xmm0
  call   <degToRad@plt>
  movq   %xmm0,%rax
  movsd  D56(%rip),%xmm0
  movq   %rax,%xmm3
  divsd  %xmm0,%xmm3
  movq   %xmm3,%rax
  movq   %rax,%xmm0
  call   <tan@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <degToRad@plt>
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x40(%rbp),%rax
  movq   %rax,%xmm0
  call   <degToRad@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x28(%rbp)
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <degToRad@plt>
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x20(%rbp)
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <degToRad@plt>
  movsd  D16(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x40(%rbp),%rax
  movq   %rax,%xmm0
  call   <degToRad@plt>
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x30(%rbp),%xmm1
  movsd  -0x48(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0x28(%rbp),%xmm0
  subsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm2
  movsd  D16(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  -0x38(%rbp),%xmm0
  mulsd  -0x28(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x38(%rbp),%xmm2
  movsd  D31(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  mulsd  -0x38(%rbp),%xmm1
  mulsd  -0x18(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x48(%rbp),%xmm2
  movsd  D40(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  mulsd  -0x48(%rbp),%xmm1
  mulsd  -0x10(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x8(%rbp),%rax
  movq   %rax,%xmm0
  call   <radToDeg@plt>
  movsd  D16(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_label: D56
type_guessed: f64
name: parse_data
arguments: data_label: D16
type_guessed: f64
name: parse_data
arguments: data_label: D31
type_guessed: f64
name: parse_data
arguments: data_label: D40
type_guessed: f64

tool: D56: 2

tool: D16: 4

tool: D31: 0.5

tool: D40: 1.25

assistant: ```c
double calcEquationOfTime(double t) {
  double epsilon = calcObliquityCorrection(t);
  double l0 = calcGeomMeanLongSun(t);
  double e = calcEccentricityEarthOrbit(t);
  double m = calcGeomMeanAnomalySun(t);
  double y = tan(degToRad(epsilon) / 2.0);
  y *= y;
  double sin2l0 = sin(2.0 * degToRad(l0));
  double sinm = sin(degToRad(m));
  double cos2l0 = cos(2.0 * degToRad(l0));
  double sin4l0 = sin(4.0 * degToRad(l0));
  double sin2m = sin(2.0 * degToRad(m));
  double Etime = y * sin2l0 - 2.0 * e * sinm + 4.0 * e * y * sinm * cos2l0 - 0.5 * y * y * sin4l0 - 1.25 * e * e * sin2m;
  return radToDeg(Etime) * 4.0;
}
```