user: What is the c source code of the assembly code below:

```asm
<solar_azimuth>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %edi,-0x44(%rbp)
  mov    %esi,-0x48(%rbp)
  movsd  %xmm0,-0x50(%rbp)
  movsd  %xmm1,-0x58(%rbp)
  movsd  %xmm2,-0x60(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  D137(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D22(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x60(%rbp),%xmm0
  mov    -0x58(%rbp),%rcx
  mov    -0x48(%rbp),%edx
  mov    -0x44(%rbp),%eax
  movapd %xmm0,%xmm1
  movq   %rcx,%xmm0
  mov    %edx,%esi
  mov    %eax,%edi
  call   <LAT@plt>
  mov    %eax,-0x34(%rbp)
  mov    -0x48(%rbp),%eax
  mov    %eax,%edi
  call   <declination@plt>
  movsd  D10(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x34(%rbp),%eax
  mov    %eax,%edi
  call   <hour_angle@plt>
  movsd  D10(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x60(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mov    -0x50(%rbp),%rcx
  mov    -0x48(%rbp),%edx
  mov    -0x44(%rbp),%eax
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rcx,%xmm0
  mov    %edx,%esi
  mov    %eax,%edi
  call   <solar_zenith@plt>
  movsd  D10(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x20(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movsd  %xmm0,-0x68(%rbp)
  mov    -0x8(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movapd %xmm0,%xmm3
  mulsd  -0x68(%rbp),%xmm3
  movsd  %xmm3,-0x68(%rbp)
  mov    -0x10(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x68(%rbp),%xmm3
  subsd  %xmm0,%xmm3
  movsd  %xmm3,-0x68(%rbp)
  mov    -0x20(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x68(%rbp),%xmm3
  divsd  %xmm0,%xmm3
  movsd  %xmm3,-0x68(%rbp)
  mov    -0x8(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movapd %xmm0,%xmm1
  movsd  -0x68(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm0
  movsd  D20(%rip),%xmm1
  comisd %xmm1,%xmm0
  jbe    L46
  movsd  D108(%rip),%xmm0
  comisd -0x28(%rbp),%xmm0
  jbe    L46
  movsd  D20(%rip),%xmm0
  movsd  %xmm0,-0x28(%rbp)
L46:
  movsd  D1(%rip),%xmm0
  comisd -0x28(%rbp),%xmm0
  jbe    L11
  movsd  -0x28(%rbp),%xmm0
  comisd D104(%rip),%xmm0
  jbe    L11
  movsd  D1(%rip),%xmm0
  movsd  %xmm0,-0x28(%rbp)
L11:
  mov    -0x28(%rbp),%rax
  movq   %rax,%xmm0
  call   <acos@plt>
  movq   %xmm0,%rax
  movq   D17(%rip),%xmm0
  movq   %rax,%xmm4
  xorpd  %xmm0,%xmm4
  movapd %xmm4,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  cmpl   $0xa8c0,-0x34(%rbp)
  jg     L21
  cmpl   $0x0,-0x34(%rbp)
  jns    L32
L21:
  movsd  -0x30(%rbp),%xmm0
  movq   D17(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
L32:
  movsd  -0x30(%rbp),%xmm1
  movsd  D11(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: read_data
arguments: tag: D137
type_infered: f64
name: read_data
arguments: tag: D22
type_infered: f64
name: read_data
arguments: tag: D10
type_infered: byte[8]
name: read_data
arguments: tag: D20
type_infered: f64
name: read_data
arguments: tag: D108
type_infered: byte[8]
name: read_data
arguments: tag: D1
type_infered: f32[2]
name: read_data
arguments: tag: D104
type_infered: byte[16]
name: read_data
arguments: tag: D17
type_infered: f32[2]
name: read_data
arguments: tag: D11
type_infered: byte[8]

tool: D137: 3.1415926

tool: D22: 180

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

tool: D20: 1

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

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

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

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

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

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;
}
```