user: 下面的汇编指令对应的C源代码是什么:

```asm
<solar_azimuth>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  mov    %edi,%r12d
  push   %rbp
  mov    %esi,%ebp
  push   %rbx
  sub    $0x40,%rsp
  lea    0x38(%rsp),%r13
  lea    0x30(%rsp),%r14
  movsd  %xmm2,(%rsp)
  movsd  %xmm0,0x28(%rsp)
  mov    %r14,%rsi
  mov    %r13,%rdi
  mulsd  D23(%rip),%xmm0
  movsd  %xmm1,0x10(%rsp)
  divsd  D19(%rip),%xmm0
  call   <sincos@plt>
  mov    %ebp,%esi
  movsd  (%rsp),%xmm2
  movsd  0x10(%rsp),%xmm5
  movsd  0x38(%rsp),%xmm3
  mov    %r12d,%edi
  movsd  0x30(%rsp),%xmm7
  movapd %xmm2,%xmm1
  movsd  %xmm2,0x18(%rsp)
  movapd %xmm5,%xmm0
  movsd  %xmm3,0x20(%rsp)
  movsd  %xmm7,0x8(%rsp)
  call   <LAT@plt>
  mov    %ebp,%edi
  mov    %eax,%ebx
  call   <declination@plt>
  mov    %ebx,%edi
  movsd  %xmm0,(%rsp)
  call   <hour_angle@plt>
  movsd  0x10(%rsp),%xmm5
  mov    %ebp,%esi
  mov    %r12d,%edi
  movsd  0x28(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm2
  movapd %xmm5,%xmm1
  movapd %xmm4,%xmm0
  call   <solar_zenith@plt>
  mulsd  D59(%rip),%xmm0
  mov    %r14,%rsi
  mov    %r13,%rdi
  call   <sincos@plt>
  movsd  D59(%rip),%xmm6
  mulsd  (%rsp),%xmm6
  movapd %xmm6,%xmm0
  call   <sin@plt>
  movsd  0x20(%rsp),%xmm3
  mulsd  0x30(%rsp),%xmm3
  movapd %xmm0,%xmm1
  movapd %xmm3,%xmm0
  subsd  %xmm1,%xmm0
  divsd  0x38(%rsp),%xmm0
  divsd  0x8(%rsp),%xmm0
  comisd D84(%rip),%xmm0
  jbe    L19
  movsd  D47(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L19
  mov    D84(%rip),%rax
  movq   %rax,%xmm0
L8:
  call   <acos@plt>
  cmp    $0xa8c0,%ebx
  ja     L21
  xorpd  D69(%rip),%xmm0
L21:
  mulsd  D88(%rip),%xmm0
  add    $0x40,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopw   0x0(%rax,%rax,1)
L19:
  movsd  D57(%rip),%xmm2
  comisd %xmm0,%xmm2
  jbe    L8
  movsd  D74(%rip),%xmm1
  cmpnltsd %xmm0,%xmm1
  andpd  %xmm1,%xmm0
  andnpd %xmm2,%xmm1
  orpd   %xmm1,%xmm0
  jmp    L8
```

tool_calls: name: read_data
arguments: data_label: D23
type: f64
name: read_data
arguments: data_label: D19
type: f64
name: read_data
arguments: data_label: D59
type: byte[8]
name: read_data
arguments: data_label: D84
type: f64[2]
name: read_data
arguments: data_label: D47
type: byte[8]
name: read_data
arguments: data_label: D69
type: f32[4]
name: read_data
arguments: data_label: D88
type: byte[8]
name: read_data
arguments: data_label: D57
type: f32[2]
name: read_data
arguments: data_label: D74
type: byte[8]

tool: D23: 3.1415926

tool: D19: 180

tool: {"D59": [0x18,0x39,0x35,0x9d,0x46,0xdf,0x91,0x3f]}

tool: {"D84": [1.0,3.1415926]}

tool: {"D47": [0x0b,0x7a,0x6f,0x0c,0x01,0x00,0xf0,0x3f]}

tool: {"D69": [0.0,-0.0,0.0,0.0]}

tool: {"D88": [0xaa,0x9f,0x96,0x22,0xdc,0xa5,0x4c,0x40]}

tool: {"D57": [0.0,-1.875]}

tool: {"D74": [0x0b,0x7a,0x6f,0x0c,0x01,0x00,0xf0,0xbf]}

assistant: ```c
double solar_azimuth(int time, int day, double latitude, double longitude, double long_std) {
  double theta = 0, omega = 0, delta = 0, phi = latitude * 3.1415926 / 180.0, psi = 0, cospsi = 0;
  int lat = LAT(time, day, longitude, long_std);
  delta = 3.1415926 / 180.0 * declination(day);
  omega = 3.1415926 / 180.0 * hour_angle(lat);
  theta = 3.1415926 / 180.0 * solar_zenith(time, day, latitude, longitude, long_std);
  cospsi = (cos(theta) * sin(phi) - sin(delta)) / sin(theta) / cos(phi);
  if (cospsi > 1.0 && cospsi < 1.0 + 1e-6)
    cospsi = 1.0;
  if (cospsi < -1.0 && cospsi > -1.0 - 1e-6)
    cospsi = -1.0;
  psi = -acos(cospsi);
  if (lat > 43200 || lat < 0)
    psi = -psi;
  psi *= (180.0 / 3.1415926);
  return psi;
}
```