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

```asm
<Aberration>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0xc8,%rsp
  movsd  D281(%rip),%xmm3
  mov    %fs:0x28,%rax
  mov    %rax,0xb8(%rsp)
  xor    %eax,%eax
  movq   %xmm3,%r14
  test   %edx,%edx
  jns    L57
  movsd  D348(%rip),%xmm4
  movq   %xmm4,%r14
L57:
  xor    %eax,%eax
  call   <LongitudeSun@plt>
  xor    %eax,%eax
  movq   %xmm0,%r15
  call   <Eccentricity@plt>
  xor    %eax,%eax
  movsd  %xmm0,(%rsp)
  call   <LongitudePerihelion@plt>
  xor    %eax,%eax
  movsd  %xmm0,0x30(%rsp)
  call   <MeanObliquity@plt>
  movsd  (%rbx),%xmm1
  movsd  0x0(%rbp),%xmm3
  movapd %xmm0,%xmm9
  movapd %xmm1,%xmm0
  andpd  D43(%rip),%xmm0
  comisd D74(%rip),%xmm0
  ja     L53
  movsd  D239(%rip),%xmm0
  lea    0x88(%rsp),%r12
  lea    0x80(%rsp),%r13
  movsd  %xmm1,0x48(%rsp)
  mov    %r13,%rsi
  mov    %r12,%rdi
  movsd  %xmm9,0x18(%rsp)
  mulsd  %xmm3,%xmm0
  movsd  %xmm3,0x78(%rsp)
  divsd  D274(%rip),%xmm0
  call   <sincos@plt>
  mov    %r13,%rsi
  movsd  D239(%rip),%xmm0
  movq   %r15,%xmm6
  mov    %r12,%rdi
  movsd  0x88(%rsp),%xmm4
  movsd  0x80(%rsp),%xmm3
  mulsd  %xmm6,%xmm0
  movsd  %xmm4,0x10(%rsp)
  divsd  D186(%rip),%xmm0
  movsd  %xmm3,0x8(%rsp)
  call   <sincos@plt>
  movsd  0x18(%rsp),%xmm9
  movsd  0x80(%rsp),%xmm5
  movsd  0x88(%rsp),%xmm1
  mulsd  D239(%rip),%xmm9
  divsd  D186(%rip),%xmm9
  movapd %xmm9,%xmm0
  movsd  %xmm5,0x20(%rsp)
  movsd  %xmm9,0x70(%rsp)
  movsd  %xmm1,0x28(%rsp)
  call   <cos@plt>
  movsd  0x8(%rsp),%xmm3
  mulsd  (%rsp),%xmm3
  mov    %r13,%rsi
  movsd  %xmm0,0x40(%rsp)
  mov    %r12,%rdi
  movsd  D239(%rip),%xmm0
  mulsd  0x30(%rsp),%xmm0
  divsd  D186(%rip),%xmm0
  movsd  %xmm3,0x18(%rsp)
  call   <sincos@plt>
  movsd  0x48(%rsp),%xmm1
  mov    %r13,%rsi
  mov    %r12,%rdi
  movsd  D239(%rip),%xmm0
  movsd  0x80(%rsp),%xmm7
  movsd  0x88(%rsp),%xmm11
  movsd  %xmm1,0x68(%rsp)
  mulsd  %xmm1,%xmm0
  movsd  %xmm7,0x30(%rsp)
  movsd  %xmm11,0x38(%rsp)
  divsd  D186(%rip),%xmm0
  call   <sincos@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  0x40(%rsp),%xmm5
  movq   %r14,%xmm12
  movq   D162(%rip),%xmm8
  movsd  0x28(%rsp),%xmm0
  movsd  0x88(%rsp),%xmm4
  movsd  %xmm5,0x58(%rsp)
  movsd  0x80(%rsp),%xmm7
  xorpd  %xmm8,%xmm2
  movsd  0x70(%rsp),%xmm9
  mulsd  0x20(%rsp),%xmm2
  movsd  %xmm4,0x60(%rsp)
  movsd  0x10(%rsp),%xmm4
  mulsd  D33(%rip),%xmm12
  movsd  %xmm7,0x50(%rsp)
  mulsd  %xmm4,%xmm0
  mulsd  %xmm5,%xmm2
  movsd  %xmm12,0x48(%rsp)
  subsd  %xmm0,%xmm2
  movsd  0x18(%rsp),%xmm0
  mulsd  0x30(%rsp),%xmm0
  mulsd  %xmm5,%xmm0
  movsd  (%rsp),%xmm5
  mulsd  %xmm4,%xmm5
  addsd  %xmm0,%xmm2
  movsd  0x38(%rsp),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm2
  movapd %xmm9,%xmm0
  divsd  %xmm7,%xmm2
  mulsd  %xmm12,%xmm2
  movsd  %xmm2,0x40(%rsp)
  call   <tan@plt>
  movsd  0x50(%rsp),%xmm7
  movsd  0x60(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm6
  movsd  0x20(%rsp),%xmm5
  mulsd  %xmm0,%xmm7
  movsd  0x40(%rsp),%xmm2
  movq   D162(%rip),%xmm8
  divsd  D245(%rip),%xmm2
  mulsd  %xmm4,%xmm6
  movsd  0x78(%rsp),%xmm3
  xorpd  %xmm8,%xmm5
  movapd %xmm5,%xmm0
  movsd  0x58(%rsp),%xmm5
  addsd  %xmm3,%xmm2
  mulsd  %xmm5,%xmm0
  movapd %xmm7,%xmm9
  movsd  0x8(%rsp),%xmm7
  subsd  %xmm6,%xmm9
  movsd  0x48(%rsp),%xmm6
  mulsd  %xmm4,%xmm7
  mulsd  0x28(%rsp),%xmm7
  mulsd  0x18(%rsp),%xmm4
  mulsd  %xmm9,%xmm0
  mulsd  0x38(%rsp),%xmm4
  subsd  %xmm7,%xmm0
  movsd  0x30(%rsp),%xmm7
  mulsd  (%rsp),%xmm7
  mulsd  %xmm7,%xmm5
  mulsd  %xmm9,%xmm5
  addsd  %xmm5,%xmm0
  addsd  %xmm4,%xmm0
  mulsd  %xmm0,%xmm6
  divsd  D152(%rip),%xmm6
  movsd  0x68(%rsp),%xmm1
  addsd  %xmm1,%xmm6
L7:
  movsd  %xmm2,0x0(%rbp)
  movsd  %xmm6,(%rbx)
  mov    0xb8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L34
  add    $0xc8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L53:
  movapd %xmm3,%xmm0
  lea    0xb0(%rsp),%rsi
  lea    0xa8(%rsp),%rdi
  call   <CelestialToEcliptical@plt>
  movsd  0x30(%rsp),%xmm0
  movsd  0xa8(%rsp),%xmm4
  lea    0x88(%rsp),%r12
  movsd  0xb0(%rsp),%xmm3
  lea    0x80(%rsp),%r13
  mov    %r12,%rdi
  subsd  %xmm4,%xmm0
  mov    %r13,%rsi
  movsd  %xmm4,0x8(%rsp)
  mulsd  D239(%rip),%xmm0
  movsd  %xmm3,0x10(%rsp)
  divsd  D186(%rip),%xmm0
  call   <sincos@plt>
  mov    %r13,%rsi
  movsd  0x88(%rsp),%xmm6
  movsd  0x80(%rsp),%xmm2
  mov    %r12,%rdi
  movq   %r15,%xmm0
  subsd  0x8(%rsp),%xmm0
  mulsd  D239(%rip),%xmm0
  movsd  %xmm6,0x18(%rsp)
  movsd  %xmm2,0x20(%rsp)
  divsd  D186(%rip),%xmm0
  call   <sincos@plt>
  movsd  0x20(%rsp),%xmm2
  mulsd  (%rsp),%xmm2
  mov    %r13,%rsi
  subsd  0x80(%rsp),%xmm2
  movsd  D239(%rip),%xmm0
  mov    %r12,%rdi
  movsd  0x88(%rsp),%xmm1
  mulsd  0x10(%rsp),%xmm0
  divsd  D186(%rip),%xmm0
  movsd  %xmm1,0x28(%rsp)
  movsd  %xmm2,0x20(%rsp)
  call   <sincos@plt>
  movsd  0x20(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm1
  movq   %r14,%xmm7
  movsd  D347(%rip),%xmm0
  mulsd  D33(%rip),%xmm7
  maxsd  0x80(%rsp),%xmm0
  mulsd  (%rsp),%xmm1
  movsd  D152(%rip),%xmm4
  divsd  %xmm0,%xmm2
  movsd  0x88(%rsp),%xmm0
  xorpd  D162(%rip),%xmm0
  mulsd  %xmm7,%xmm0
  movapd %xmm1,%xmm3
  movsd  0x28(%rsp),%xmm1
  subsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm0
  divsd  %xmm4,%xmm0
  mulsd  %xmm7,%xmm2
  divsd  %xmm4,%xmm2
  addsd  0x10(%rsp),%xmm0
  addsd  0x8(%rsp),%xmm2
  movsd  %xmm2,0x8(%rsp)
  call   <Map180@plt>
  comisd D163(%rip),%xmm0
  movsd  0x8(%rsp),%xmm2
  movapd %xmm0,%xmm1
  jbe    L4
  movsd  D186(%rip),%xmm1
  addsd  D186(%rip),%xmm2
  subsd  0xa0(%rsp),%xmm1
L4:
  movsd  D134(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L46
  movsd  D173(%rip),%xmm0
  addsd  D186(%rip),%xmm2
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
L46:
  movapd %xmm2,%xmm0
  movsd  %xmm1,(%rsp)
  call   <Map360@plt>
  movsd  (%rsp),%xmm1
  lea    0xa0(%rsp),%rsi
  lea    0x98(%rsp),%rdi
  call   <EclipticalToCelestial@plt>
  movsd  0x98(%rsp),%xmm2
  movsd  0xa0(%rsp),%xmm6
  jmp    L7
L34:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_tag: D281
type_guessed: f64
name: retrieve_data
arguments: data_tag: D348
type_guessed: byte[8]
name: retrieve_data
arguments: data_tag: D43
type_guessed: qword[2]
name: retrieve_data
arguments: data_tag: D74
type_guessed: f64[2]
name: retrieve_data
arguments: data_tag: D239
type_guessed: f64
name: retrieve_data
arguments: data_tag: D274
type_guessed: f64
name: retrieve_data
arguments: data_tag: D186
type_guessed: f64
name: retrieve_data
arguments: data_tag: D162
type_guessed: byte[16]
name: retrieve_data
arguments: data_tag: D33
type_guessed: f64
name: retrieve_data
arguments: data_tag: D245
type_guessed: byte[8]
name: retrieve_data
arguments: data_tag: D152
type_guessed: f64
name: retrieve_data
arguments: data_tag: D347
type_guessed: f64
name: retrieve_data
arguments: data_tag: D163
type_guessed: f64[2]
name: retrieve_data
arguments: data_tag: D134
type_guessed: byte[8]
name: retrieve_data
arguments: data_tag: D173
type_guessed: byte[8]

tool: D281: 1

tool: {"D348": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: {"D43": [0x7fffffffffffffff,0x0000000000000000]}

tool: {"D74": [85.0,3.14159265358]}

tool: D239: 3.14159265358

tool: D274: 12

tool: D186: 180

tool: {"D162": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: D33: 20.49552

tool: {"D245": [0x00,0x00,0x00,0x00,0x00,0x5e,0xea,0x40]}

tool: D152: 3600

tool: D347: 0.0001

tool: {"D163": [90.0,-90.0]}

tool: {"D134": [0x00,0x00,0x00,0x00,0x00,0x80,0x56,0xc0]}

tool: {"D173": [0x00,0x00,0x00,0x00,0x00,0x80,0x66,0xc0]}

assistant: ```c
void Aberration(double *ra, double *dec, int dirflag) {
  double ra_in, dec_in, ra_out, dec_out;
  double lambda, beta;
  double elong, elat;
  double dlong, dlat;
  double ec, lp, glsun, ka, eps0;
  double c1, c2;
  double dra, ddec;
  double dir;
  dir = 1.0;
  if (dirflag < 0) {
    dir = -1.0;
  }
  ka = 20.49552;
  glsun = LongitudeSun();
  ec = Eccentricity();
  lp = LongitudePerihelion();
  eps0 = MeanObliquity();
  ra_in = *ra;
  dec_in = *dec;
  if (fabs(dec_in) > (double)85.) {
    CelestialToEcliptical(ra_in, dec_in, &lambda, &beta);
    elong = lambda;
    elat = beta;
    c1 = -cos((glsun - elong) * 3.14159265358 / 180.) + ec * cos((lp - elong) * 3.14159265358 / 180.);
    c2 = cos(elat * 3.14159265358 / 180.);
    if (c2 < (double)0.0001)
      c2 = 0.0001;
    dlong = dir * ka * (c1 / c2) / 3600.;
    c1 = -sin(elat * 3.14159265358 / 180.);
    c2 = sin((glsun - elong) * 3.14159265358 / 180.) - ec * sin((lp - elong) * 3.14159265358 / 180.);
    dlat = dir * ka * c1 * c2 / 3600.;
    elong = elong + dlong;
    elat = elat + dlat;
    elat = Map180(elat);
    if (elat > 90.) {
      elat = 180. - dec_out;
      elong = elong + 180.;
    }
    if (elat < -90.) {
      elat = -180. - elat;
      elong = elong + 180.;
    }
    elong = Map360(elong);
    EclipticalToCelestial(elong, elat, &ra_out, &dec_out);
  } else {
    dra = (-cos(ra_in * 3.14159265358 / 12.) * cos(glsun * 3.14159265358 / 180.) * cos(eps0 * 3.14159265358 / 180.) - sin(ra_in * 3.14159265358 / 12.) * sin(glsun * 3.14159265358 / 180.) + ec * cos(ra_in * 3.14159265358 / 12.) * cos(lp * 3.14159265358 / 180.) * cos(eps0 * 3.14159265358 / 180) + ec * sin(ra_in * 3.14159265358 / 12.) * sin(lp * 3.14159265358 / 180.)) / cos(dec_in * 3.14159265358 / 180.);
    dra = dir * ka * dra;
    dra = dra / (15. * 3600.);
    ddec = -cos(glsun * 3.14159265358 / 180) * cos(eps0 * 3.14159265358 / 180.) * (tan(eps0 * 3.14159265358 / 180) * cos(dec_in * 3.14159265358 / 180.) - sin(ra_in * 3.14159265358 / 12.) * sin(dec_in * 3.14159265358 / 180.)) - cos(ra_in * 3.14159265358 / 12.) * sin(dec_in * 3.14159265358 / 180.) * sin(glsun * 3.14159265358 / 180) + ec * cos(lp * 3.14159265358 / 180.) * cos(eps0 * 3.14159265358 / 180.) * (tan(eps0 * 3.14159265358 / 180.) * cos(dec_in * 3.14159265358 / 180.) - sin(ra_in * 3.14159265358 / 12.) * sin(dec_in * 3.14159265358 / 180.)) + ec * cos(ra_in * 3.14159265358 / 12.) * sin(dec_in * 3.14159265358 / 180.) * sin(lp * 3.14159265358 / 180.);
    ddec = dir * ka * ddec;
    ddec = ddec / 3600.;
    ra_out = ra_in + dra;
    dec_out = dec_in + ddec;
  }
  *ra = ra_out;
  *dec = dec_out;
}
```