user: What would the C source look like for the assembly code provided here?

```asm
<sunpos>:
  endbr64
  push   %r13
  mov    $0x1,%eax
  push   %r12
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x38,%rsp
  movsd  %xmm0,0x8(%rsp)
  lea    0x28(%rsp),%r12
  lea    0x20(%rsp),%r13
  mulsd  D66(%rip),%xmm0
  addsd  D27(%rip),%xmm0
  call   <revolution@plt>
  pxor   %xmm3,%xmm3
  mov    %r13,%rsi
  mov    %r12,%rdi
  movsd  0x8(%rsp),%xmm1
  movsd  D64(%rip),%xmm0
  cvtsi2sd %eax,%xmm3
  movsd  D14(%rip),%xmm2
  mulsd  %xmm1,%xmm0
  addsd  D17(%rip),%xmm0
  mulsd  D8(%rip),%xmm1
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm0,0x8(%rsp)
  movsd  D125(%rip),%xmm0
  mulsd  %xmm3,%xmm0
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x10(%rsp)
  call   <sincos@plt>
  movsd  0x10(%rsp),%xmm2
  mov    %r13,%rsi
  movsd  D24(%rip),%xmm0
  movsd  0x20(%rsp),%xmm1
  movsd  0x18(%rsp),%xmm3
  mov    %r12,%rdi
  mulsd  %xmm2,%xmm0
  mulsd  0x28(%rsp),%xmm0
  mulsd  %xmm2,%xmm1
  addsd  D53(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  addsd  %xmm3,%xmm0
  mulsd  D125(%rip),%xmm0
  call   <sincos@plt>
  movsd  0x10(%rsp),%xmm2
  movsd  0x20(%rsp),%xmm1
  movsd  D53(%rip),%xmm0
  movsd  0x28(%rsp),%xmm3
  subsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm2
  subsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  ucomisd %xmm0,%xmm2
  ja     L44
  sqrtsd %xmm0,%xmm0
L37:
  mulsd  %xmm0,%xmm3
  movapd %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm3,%xmm4
  mulsd  %xmm3,%xmm4
  addsd  %xmm4,%xmm0
  ucomisd %xmm0,%xmm2
  ja     L1
  sqrtsd %xmm0,%xmm0
L22:
  movsd  %xmm0,0x0(%rbp)
  movapd %xmm3,%xmm0
  call   <atan2@plt>
  mulsd  D24(%rip),%xmm0
  movsd  D13(%rip),%xmm1
  addsd  0x8(%rsp),%xmm0
  comisd %xmm1,%xmm0
  jb     L40
  subsd  %xmm1,%xmm0
L40:
  movsd  %xmm0,(%rbx)
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
L44:
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm1
  pxor   %xmm2,%xmm2
  jmp    L37
L1:
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm1
  jmp    L22
```

tool_calls: name: extract_symbolic_data
arguments: label: D66
type: f64
name: extract_symbolic_data
arguments: label: D27
type: f64
name: extract_symbolic_data
arguments: label: D64
type: byte[8]
name: extract_symbolic_data
arguments: label: D14
type: f64
name: extract_symbolic_data
arguments: label: D17
type: f64
name: extract_symbolic_data
arguments: label: D8
type: byte[8]
name: extract_symbolic_data
arguments: label: D125
type: byte[8]
name: extract_symbolic_data
arguments: label: D24
type: byte[8]
name: extract_symbolic_data
arguments: label: D53
type: f64
name: extract_symbolic_data
arguments: label: D13
type: f64

tool: D66: 0.9856002585

tool: D27: 356.047

tool: {"D64": [0x46,0xc9,0xbd,0x56,0xc8,0xb0,0x08,0x3f]}

tool: D14: 0.016709

tool: D17: 282.9404

tool: {"D8": [0xd9,0x5b,0xf3,0xca,0x26,0xc6,0x13,0x3e]}

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

tool: {"D24": [0xf8,0xc1,0x63,0x1a,0xdc,0xa5,0x4c,0x40]}

tool: D53: 1

tool: D13: 360

assistant: ```c
void sunpos(double d, double *lon, double *r) {
  double M, w, e, E, x, y, v;
  M = revolution(356.0470 + 0.9856002585 * d);
  w = 282.9404 + 4.70935E-5 * d;
  e = 0.016709 - 1.151E-9 * d;
  E = M + e * (180.0 / 3.1415926535897932384) * sin((M) * (3.1415926535897932384 / 180.0)) * (1.0 + e * cos((M) * (3.1415926535897932384 / 180.0)));
  x = cos((E) * (3.1415926535897932384 / 180.0)) - e;
  y = sqrt(1.0 - e * e) * sin((E) * (3.1415926535897932384 / 180.0));
  *r = sqrt(x * x + y * y);
  v = ((180.0 / 3.1415926535897932384) * atan2(y, x));
  *lon = v + w;
  if (*lon >= 360.0)
    *lon -= 360.0;
}
```