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

```asm
<caldat>:
  endbr64
  mov    %rdx,%r8
  cmp    $0x231518,%rdi
  jle    L1
  lea    -0x1c7dd0(%rdi),%rax
  pxor   %xmm0,%xmm0
  movsd  D60(%rip),%xmm2
  cvtsi2sd %rax,%xmm0
  subsd  %xmm2,%xmm0
  divsd  D18(%rip),%xmm0
  cvttsd2si %xmm0,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  lea    0x1(%rdx,%rdi,1),%rax
  mulsd  %xmm2,%xmm0
  cvttsd2si %xmm0,%rdx
  sub    %rdx,%rax
L40:
  lea    0x5f4(%rax),%rdx
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  sub    $0x2534ca,%rax
  cvtsi2sd %rax,%xmm0
  subsd  D9(%rip),%xmm0
  divsd  D65(%rip),%xmm0
  addsd  D25(%rip),%xmm0
  cvttsd2si %xmm0,%rax
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm1
  lea    (%rax,%rax,8),%r9
  lea    (%rax,%r9,8),%r9
  lea    (%r9,%r9,4),%r9
  cvtsi2sd %r9,%xmm0
  mov    %rdx,%r9
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  movsd  D48(%rip),%xmm1
  cvttsd2si %xmm0,%r10
  pxor   %xmm0,%xmm0
  sub    %r10,%r9
  sub    %r10d,%edx
  cvtsi2sd %r9,%xmm0
  divsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%r9
  pxor   %xmm0,%xmm0
  cvtsi2sd %r9,%xmm0
  mulsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%r10
  sub    %r10d,%edx
  mov    %edx,(%r8)
  lea    -0x1(%r9),%edx
  sub    $0xd,%r9d
  cmp    $0xd,%edx
  cmovge %r9d,%edx
  mov    %edx,(%rsi)
  lea    -0x126b(%rax),%edx
  mov    %edx,(%rcx)
  cmpl   $0x2,(%rsi)
  jle    L42
  lea    -0x126c(%rax),%edx
  mov    %edx,(%rcx)
L42:
  test   %edx,%edx
  jg     L12
  sub    $0x1,%edx
  mov    %edx,(%rcx)
L12:
  test   %rdi,%rdi
  js     L2
  ret
  nopl   (%rax)
L1:
  test   %rdi,%rdi
  js     L55
  movsd  D60(%rip),%xmm2
  mov    %rdi,%rax
  jmp    L40
  nopl   (%rax)
L2:
  movabs $0xe5aae5e0b9ee9a9b,%rdx
  mov    %rdi,%rax
  imul   %rdx
  add    %rdi,%rdx
  sar    $0x3f,%rdi
  sar    $0xf,%rdx
  sub    %rdx,%rdi
  add    $0x1,%rdi
  imul   $0xffffff9c,%edi,%edi
  add    %edi,(%rcx)
  ret
  nopl   0x0(%rax,%rax,1)
L55:
  mov    %rdi,%rax
  movsd  D60(%rip),%xmm2
  movabs $0xe5aae5e0b9ee9a9b,%rdx
  imul   %rdx
  mov    %rdi,%rax
  sar    $0x3f,%rax
  add    %rdi,%rdx
  sar    $0xf,%rdx
  sub    %rdx,%rax
  add    $0x1,%rax
  imul   $0x8ead,%rax,%rax
  add    %rdi,%rax
  jmp    L40
```

tool_calls: name: extract_symbolic_data
arguments: tag: D60
type_infered: f64
name: extract_symbolic_data
arguments: tag: D18
type_infered: f64
name: extract_symbolic_data
arguments: tag: D9
type_infered: f64
name: extract_symbolic_data
arguments: tag: D65
type_infered: f64
name: extract_symbolic_data
arguments: tag: D25
type_infered: f64
name: extract_symbolic_data
arguments: tag: D48
type_infered: f64

tool: D60: 0.25

tool: D18: 36524.25

tool: D9: 122.1

tool: D65: 365.25

tool: D25: 6680

tool: D48: 30.6001

assistant: ```c
void caldat(long julian, int *mm, int *id, int *iyyy) {
  long ja, jalpha, jb, jc, jd, je;
  if (julian >= 2299161) {
    jalpha = (long)(((double)(julian - 1867216) - 0.25) / 36524.25);
    ja = julian + 1 + jalpha - (long)(0.25 * jalpha);
  } else if (julian < 0) {
    ja = julian + 36525 * (1 - julian / 36525);
  } else
    ja = julian;
  jb = ja + 1524;
  jc = (long)(6680.0 + ((double)(jb - 2439870) - 122.1) / 365.25);
  jd = (long)(365 * jc + (0.25 * jc));
  je = (long)((jb - jd) / 30.6001);
  *id = jb - jd - (long)(30.6001 * je);
  *mm = je - 1;
  if (*mm > 12)
    *mm -= 12;
  *iyyy = jc - 4715;
  if (*mm > 2)
    --(*iyyy);
  if (*iyyy <= 0)
    --(*iyyy);
  if (julian < 0)
    *iyyy -= 100 * (1 - julian / 36525);
}
```