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

```asm
<solar_coords_and_times>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  add    $0xffffffffffffff80,%rsp
  movsd  %xmm0,-0x58(%rbp)
  mov    %rdi,-0x60(%rbp)
  mov    %rsi,-0x68(%rbp)
  mov    %rdx,-0x70(%rbp)
  movsd  -0x58(%rbp),%xmm0
  movsd  D64(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D16(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D50(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  D213(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x38(%rbp),%rax
  mov    $0x6,%edi
  movq   %rax,%xmm0
  call   <fmod@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x38(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D32(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D167(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  D213(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x30(%rbp),%rax
  mov    $0x6,%edi
  movq   %rax,%xmm0
  call   <fmod@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x30(%rbp)
  movsd  -0x30(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sin@plt>
  movsd  D210(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x30(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sin@plt>
  movsd  D51(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  addsd  -0x78(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm1
  movsd  D213(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x38(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D164(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  D194(%rip),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D213(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x18(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movsd  %xmm0,-0x78(%rbp)
  mov    -0x20(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movapd %xmm0,%xmm4
  mulsd  -0x78(%rbp),%xmm4
  movsd  %xmm4,-0x78(%rbp)
  mov    -0x18(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x20(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sin@plt>
  movapd %xmm0,%xmm5
  mulsd  -0x80(%rbp),%xmm5
  movq   %xmm5,%rax
  movsd  -0x78(%rbp),%xmm1
  movq   %rax,%xmm0
  call   <atan2@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x48(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x48(%rbp),%xmm0
  jb     L1
  movsd  -0x48(%rbp),%xmm1
  movsd  D87(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
L1:
  movsd  -0x18(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sin@plt>
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x20(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sin@plt>
  movapd %xmm0,%xmm6
  mulsd  -0x78(%rbp),%xmm6
  movq   %xmm6,%rax
  movq   %rax,%xmm0
  call   <asin@plt>
  movq   %xmm0,%rax
  mov    -0x60(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    -0x58(%rbp),%rax
  movq   %rax,%xmm0
  call   <calc_gmst@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movsd  D132(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x48(%rbp),%xmm1
  movsd  D87(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x68(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x68(%rbp),%rax
  mov    (%rax),%rax
  mov    $0x1,%edi
  movq   %rax,%xmm0
  call   <fmod@plt>
  movq   %xmm0,%rax
  mov    -0x68(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    -0x68(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  jbe    L28
  mov    -0x68(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D9(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x68(%rbp),%rax
  movsd  %xmm0,(%rax)
L28:
  movsd  -0x20(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sin@plt>
  movsd  D110(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  subsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  D84(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sin@plt>
  movsd  D93(%rip),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x78(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  D182(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x68(%rbp),%rax
  movsd  (%rax),%xmm0
  addsd  -0x8(%rbp),%xmm0
  mov    -0x70(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x70(%rbp),%rax
  mov    (%rax),%rax
  mov    $0x1,%edi
  movq   %rax,%xmm0
  call   <fmod@plt>
  movq   %xmm0,%rax
  mov    -0x70(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    -0x70(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  ja     L36
  jmp    L29
L36:
  mov    -0x70(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D9(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  movsd  %xmm0,(%rax)
L29:
  nop
  leave
  ret
```

tool_calls: name: read_data
arguments: tag: D64
datatype: f64
name: read_data
arguments: tag: D16
datatype: byte[8]
name: read_data
arguments: tag: D50
datatype: f64
name: read_data
arguments: tag: D213
datatype: byte[8]
name: read_data
arguments: tag: D32
datatype: byte[8]
name: read_data
arguments: tag: D167
datatype: f64
name: read_data
arguments: tag: D210
datatype: f64
name: read_data
arguments: tag: D51
datatype: byte[8]
name: read_data
arguments: tag: D164
datatype: byte[8]
name: read_data
arguments: tag: D194
datatype: f64
name: read_data
arguments: tag: D87
datatype: byte[8]
name: read_data
arguments: tag: D132
datatype: f32[2]
name: read_data
arguments: tag: D9
datatype: f64
name: read_data
arguments: tag: D110
datatype: f64
name: read_data
arguments: tag: D84
datatype: f32[2]
name: read_data
arguments: tag: D93
datatype: byte[8]
name: read_data
arguments: tag: D182
datatype: f64

tool: D64: 2451545

tool: {"D16": [0xd4,0x65,0x8c,0x6a,0x6c,0x8a,0xef,0x3f]}

tool: D50: 280.466

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

tool: {"D32": [0x72,0xb0,0xed,0xa3,0x09,0x8a,0xef,0x3f]}

tool: D167: 357.528

tool: D210: 1.915

tool: {"D51": [0x7b,0x14,0xae,0x47,0xe1,0x7a,0x94,0x3f]}

tool: {"D164": [0x48,0xaf,0xbc,0x9a,0xf2,0xd7,0x9a,0x3e]}

tool: D194: 23.44

tool: {"D87": [0x18,0x2d,0x44,0x54,0xfb,0x21,0x19,0x40]}

tool: {"D132": [0.0,1.75]}

tool: D9: 1

tool: D110: 2.466

tool: {"D84": [0.0,2.25]}

tool: {"D93": [0x89,0x41,0x60,0xe5,0xd0,0x22,0xab,0x3f]}

tool: D182: 360

assistant: ```c
void solar_coords_and_times(double jtime, double *delta, double *gw_mean_sol_time, double *gw_appar_sol_time) {
  double a;
  double jdelta;
  double L;
  double g;
  double lambda;
  double epsilon;
  double alpha;
  double gmst;
  double E;
  jdelta = jtime - 2451545.;
  L = (280.466 + .9856474 * jdelta) * (3.14159265358979323846264338327 / 180.0);
  L = fmod(L, (3.14159265358979323846264338327 * 2.));
  g = (357.528 + .9856003 * jdelta) * (3.14159265358979323846264338327 / 180.0);
  g = fmod(g, (3.14159265358979323846264338327 * 2.));
  a = 1.915 * sin(g) + .020 * sin(2.0 * g);
  lambda = L + a * (3.14159265358979323846264338327 / 180.0);
  epsilon = (23.44 - .0000004 * jdelta) * (3.14159265358979323846264338327 / 180.0);
  alpha = atan2(cos(epsilon) * sin(lambda), cos(epsilon) * cos(lambda));
  if (alpha <= 0.)
    alpha += (3.14159265358979323846264338327 * 2.);
  *delta = asin(sin(epsilon) * sin(lambda));
  gmst = calc_gmst(jtime);
  *gw_mean_sol_time = (gmst - .5 - alpha / (3.14159265358979323846264338327 * 2.));
  if ((*gw_mean_sol_time = fmod(*gw_mean_sol_time, 1.)) < 0.)
    *gw_mean_sol_time += 1.;
  E = (-a + 2.466 * sin(2 * lambda) - .053 * sin(4 * lambda)) / 360.;
  *gw_appar_sol_time = *gw_mean_sol_time + E;
  if ((*gw_appar_sol_time = fmod(*gw_appar_sol_time, 1.)) < 0.)
    *gw_appar_sol_time += 1.;
}
```