user: 请从给定的汇编指令推导出C源代码:

```asm
<solar_coords_and_times>:
  endbr64
  push   %r14
  movapd %xmm0,%xmm2
  movq   %xmm0,%r14
  push   %r13
  push   %r12
  mov    %rdi,%r12
  mov    $0x6,%edi
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %rdx,%rbx
  sub    $0x30,%rsp
  subsd  D84(%rip),%xmm2
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x10(%rsp)
  mulsd  D142(%rip),%xmm0
  addsd  D125(%rip),%xmm0
  mulsd  D34(%rip),%xmm0
  call   <fmod@plt>
  movsd  0x10(%rsp),%xmm2
  mov    $0x6,%edi
  movsd  %xmm0,0x8(%rsp)
  movsd  D151(%rip),%xmm0
  movapd %xmm2,%xmm6
  movsd  %xmm2,0x20(%rsp)
  mulsd  D102(%rip),%xmm6
  addsd  %xmm6,%xmm0
  mulsd  D34(%rip),%xmm0
  call   <fmod@plt>
  cvttsd2si %xmm0,%edi
  movsd  %xmm0,0x18(%rsp)
  call   <sin@plt>
  movsd  0x18(%rsp),%xmm1
  movsd  %xmm0,0x10(%rsp)
  addsd  %xmm1,%xmm1
  cvttsd2si %xmm1,%edi
  call   <sin@plt>
  movsd  0x10(%rsp),%xmm3
  mulsd  D93(%rip),%xmm3
  mulsd  D49(%rip),%xmm0
  movsd  0x20(%rsp),%xmm2
  mulsd  D47(%rip),%xmm2
  movsd  D34(%rip),%xmm7
  movsd  0x8(%rsp),%xmm5
  movapd %xmm3,%xmm1
  addsd  %xmm0,%xmm1
  movapd %xmm2,%xmm0
  movsd  D153(%rip),%xmm2
  mulsd  %xmm1,%xmm7
  subsd  %xmm0,%xmm2
  movsd  %xmm1,0x18(%rsp)
  mulsd  D34(%rip),%xmm2
  addsd  %xmm7,%xmm5
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm5,0x8(%rsp)
  call   <cos@plt>
  movsd  %xmm0,0x10(%rsp)
  movsd  0x8(%rsp),%xmm0
  call   <cos@plt>
  movsd  0x10(%rsp),%xmm1
  cvttsd2si 0x8(%rsp),%r13d
  mulsd  %xmm0,%xmm1
  mov    %r13d,%edi
  movsd  %xmm1,0x20(%rsp)
  call   <sin@plt>
  mulsd  0x10(%rsp),%xmm0
  movsd  0x20(%rsp),%xmm1
  call   <atan2@plt>
  pxor   %xmm5,%xmm5
  movsd  D61(%rip),%xmm3
  movsd  0x28(%rsp),%xmm2
  comisd %xmm0,%xmm5
  movapd %xmm0,%xmm1
  jb     L14
  addsd  %xmm3,%xmm1
L14:
  cvttsd2si %xmm2,%edi
  movsd  %xmm3,0x28(%rsp)
  movsd  %xmm1,0x20(%rsp)
  call   <sin@plt>
  mov    %r13d,%edi
  movsd  %xmm0,0x10(%rsp)
  call   <sin@plt>
  mulsd  0x10(%rsp),%xmm0
  call   <asin@plt>
  movsd  %xmm0,(%r12)
  movq   %r14,%xmm0
  call   <calc_gmst@plt>
  movsd  0x20(%rsp),%xmm1
  movsd  0x28(%rsp),%xmm3
  mov    $0x1,%edi
  subsd  D167(%rip),%xmm0
  divsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%rbp)
  call   <fmod@plt>
  pxor   %xmm5,%xmm5
  comisd %xmm0,%xmm5
  jbe    L15
  addsd  D184(%rip),%xmm0
L15:
  movsd  0x8(%rsp),%xmm4
  movsd  %xmm0,0x0(%rbp)
  movapd %xmm4,%xmm0
  addsd  %xmm4,%xmm0
  cvttsd2si %xmm0,%edi
  call   <sin@plt>
  movsd  D118(%rip),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  0x8(%rsp),%xmm0
  mulsd  D53(%rip),%xmm0
  subsd  0x18(%rsp),%xmm1
  movsd  %xmm1,0x10(%rsp)
  cvttsd2si %xmm0,%edi
  call   <sin@plt>
  movsd  0x10(%rsp),%xmm1
  mov    $0x1,%edi
  movsd  D20(%rip),%xmm2
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  divsd  D140(%rip),%xmm0
  addsd  0x0(%rbp),%xmm0
  movsd  %xmm0,(%rbx)
  call   <fmod@plt>
  pxor   %xmm4,%xmm4
  comisd %xmm0,%xmm4
  jbe    L25
  addsd  D184(%rip),%xmm0
L25:
  movsd  %xmm0,(%rbx)
  add    $0x30,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D84
data_type: f64
name: extract_symbolic_data
arguments: label: D142
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D125
data_type: f64
name: extract_symbolic_data
arguments: label: D34
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D151
data_type: f64
name: extract_symbolic_data
arguments: label: D102
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D93
data_type: f64
name: extract_symbolic_data
arguments: label: D49
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D47
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D153
data_type: f64
name: extract_symbolic_data
arguments: label: D61
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D167
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D184
data_type: f64
name: extract_symbolic_data
arguments: label: D118
data_type: f64
name: extract_symbolic_data
arguments: label: D53
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D20
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D140
data_type: f64

tool: D84: 2451545

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

tool: D125: 280.466

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

tool: D151: 357.528

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

tool: D93: 1.915

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

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

tool: D153: 23.44

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

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

tool: D184: 1

tool: D118: 2.466

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

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

tool: D140: 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.;
}
```