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

```asm
<thetag>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  movsd  %xmm0,-0x48(%rbp)
  mov    %rdi,-0x50(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  D32(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  D97(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x34(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x34(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  D45(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  cmpl   $0x9,-0x34(%rbp)
  jg     L24
  addl   $0x50,-0x34(%rbp)
L24:
  mov    -0x34(%rbp),%eax
  sub    $0x45,%eax
  lea    0x3(%rax),%edx
  test   %eax,%eax
  cmovs  %edx,%eax
  sar    $0x2,%eax
  mov    %eax,-0x30(%rbp)
  cmpl   $0x45,-0x34(%rbp)
  jg     L21
  mov    -0x34(%rbp),%eax
  sub    $0x48,%eax
  lea    0x3(%rax),%edx
  test   %eax,%eax
  cmovs  %edx,%eax
  sar    $0x2,%eax
  mov    %eax,-0x30(%rbp)
L21:
  mov    -0x34(%rbp),%eax
  sub    $0x46,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D93(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D98(%rip),%xmm0
  addsd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x30(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  addsd  -0x18(%rbp),%xmm0
  mov    -0x50(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D0(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D63(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  movsd  D23(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x2c(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x2c(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm1
  movsd  D23(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x10(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x28(%rbp),%xmm0
  jbe    L4
  movsd  -0x28(%rbp),%xmm1
  movsd  D23(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
L4:
  movsd  -0x28(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  pop    %rbp
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D32
type: byte[8]
name: disassemble_data
arguments: tag: D97
type: byte[8]
name: disassemble_data
arguments: tag: D45
type: f64
name: disassemble_data
arguments: tag: D93
type: f64
name: disassemble_data
arguments: tag: D98
type: f64
name: disassemble_data
arguments: tag: D0
type: f64
name: disassemble_data
arguments: tag: D63
type: f64
name: disassemble_data
arguments: tag: D23
type: f64

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

tool: {"D97": [0xfc,0xa9,0xf1,0xd2,0x4d,0x62,0x50,0x3f]}

tool: D45: 1000

tool: D93: 365

tool: D98: 7305

tool: D0: 6.3003880987

tool: D63: 1.72944494

tool: D23: 6.283185307179586

assistant: ```c
double thetag(double EP, double *DS50) {
  int JY, N, I;
  double YR, D, THETA, TEMP, THETAG;
  YR = (EP + 2.0E-7) * 1.0E-3;
  JY = (int)YR;
  YR = JY;
  D = EP - YR * 1.0E3;
  if (JY < 10)
    JY += 80;
  N = (JY - 69) / 4;
  if (JY < 70)
    N = (JY - 72) / 4;
  *DS50 = 7305.0 + 365.0 * (JY - 70) + N + D;
  THETA = 1.72944494 + 6.3003880987 * *DS50;
  TEMP = THETA / (6.2831853071795864770);
  I = (int)TEMP;
  TEMP = I;
  THETAG = THETA - TEMP * (6.2831853071795864770);
  if (THETAG < 0.0)
    THETAG += (6.2831853071795864770);
  return THETAG;
}
```