user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<Aberration>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xb8,%rsp
  mov    %rdi,%rbp
  mov    %rsi,%rbx
  mov    %fs:0x28,%rax
  mov    %rax,0xa8(%rsp)
  xor    %eax,%eax
  mov    D143(%rip),%r15
  test   %edx,%edx
  js     L77
L41:
  mov    $0x0,%eax
  call   <LongitudeSun@plt>
  movsd  %xmm0,0x10(%rsp)
  mov    $0x0,%eax
  call   <Eccentricity@plt>
  movsd  %xmm0,0x8(%rsp)
  mov    $0x0,%eax
  call   <LongitudePerihelion@plt>
  movsd  %xmm0,0x48(%rsp)
  mov    $0x0,%eax
  call   <MeanObliquity@plt>
  movsd  %xmm0,0x30(%rsp)
  movsd  0x0(%rbp),%xmm6
  movsd  %xmm6,0x28(%rsp)
  mov    (%rbx),%r12
  movq   %r12,%xmm0
  andpd  D132(%rip),%xmm0
  comisd D105(%rip),%xmm0
  ja     L50
  lea    0x78(%rsp),%r13
  lea    0x70(%rsp),%r14
  movsd  D80(%rip),%xmm0
  mulsd  0x28(%rsp),%xmm0
  divsd  D91(%rip),%xmm0
  mov    %r14,%rsi
  mov    %r13,%rdi
  call   <sincos@plt>
  movsd  0x70(%rsp),%xmm4
  movsd  %xmm4,0x18(%rsp)
  movsd  0x78(%rsp),%xmm6
  movsd  %xmm6,0x20(%rsp)
  movsd  D80(%rip),%xmm0
  mulsd  0x10(%rsp),%xmm0
  divsd  D241(%rip),%xmm0
  mov    %r14,%rsi
  mov    %r13,%rdi
  call   <sincos@plt>
  movsd  0x70(%rsp),%xmm7
  movsd  %xmm7,0x38(%rsp)
  movsd  0x78(%rsp),%xmm2
  movsd  %xmm2,0x40(%rsp)
  movsd  D80(%rip),%xmm0
  mulsd  0x30(%rsp),%xmm0
  divsd  D241(%rip),%xmm0
  movsd  %xmm0,0x58(%rsp)
  call   <cos@plt>
  movsd  %xmm0,0x10(%rsp)
  movsd  0x18(%rsp),%xmm4
  mulsd  0x8(%rsp),%xmm4
  movsd  %xmm4,0x30(%rsp)
  movsd  D80(%rip),%xmm11
  mulsd  0x48(%rsp),%xmm11
  movapd %xmm11,%xmm0
  divsd  D241(%rip),%xmm0
  mov    %r14,%rsi
  mov    %r13,%rdi
  call   <sincos@plt>
  movsd  0x70(%rsp),%xmm8
  movsd  %xmm8,0x48(%rsp)
  movsd  0x78(%rsp),%xmm12
  movsd  %xmm12,0x50(%rsp)
  movq   %r12,%xmm13
  mulsd  D80(%rip),%xmm13
  movapd %xmm13,%xmm0
  divsd  D241(%rip),%xmm0
  mov    %r14,%rsi
  mov    %r13,%rdi
  call   <sincos@plt>
  mov    0x70(%rsp),%r14
  mov    0x78(%rsp),%r13
  movsd  0x18(%rsp),%xmm5
  xorpd  D166(%rip),%xmm5
  movapd %xmm5,%xmm0
  mulsd  0x38(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm3
  mulsd  %xmm3,%xmm0
  movsd  0x20(%rsp),%xmm6
  movapd %xmm6,%xmm1
  mulsd  0x40(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  0x30(%rsp),%xmm1
  mulsd  0x48(%rsp),%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm1,%xmm0
  movapd %xmm6,%xmm7
  mulsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm1
  mulsd  0x50(%rsp),%xmm1
  addsd  %xmm1,%xmm0
  movq   %r14,%xmm11
  divsd  %xmm11,%xmm0
  movq   %r15,%xmm9
  mulsd  D99(%rip),%xmm9
  movsd  %xmm9,0x60(%rsp)
  mulsd  %xmm9,%xmm0
  movsd  %xmm0,0x68(%rsp)
  movsd  0x58(%rsp),%xmm0
  call   <tan@plt>
  movq   %r14,%xmm14
  mulsd  %xmm0,%xmm14
  movapd %xmm14,%xmm2
  movq   %r13,%xmm15
  mulsd  0x20(%rsp),%xmm15
  subsd  %xmm15,%xmm2
  movsd  0x38(%rsp),%xmm0
  xorpd  D166(%rip),%xmm0
  movsd  0x10(%rsp),%xmm3
  mulsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  movq   %r13,%xmm7
  mulsd  0x18(%rsp),%xmm7
  movapd %xmm7,%xmm1
  mulsd  0x40(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  0x48(%rsp),%xmm6
  mulsd  0x8(%rsp),%xmm6
  movapd %xmm6,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  movq   %r13,%xmm5
  mulsd  0x30(%rsp),%xmm5
  movapd %xmm5,%xmm1
  mulsd  0x50(%rsp),%xmm1
  addsd  %xmm1,%xmm0
  mulsd  0x60(%rsp),%xmm0
  movsd  0x68(%rsp),%xmm1
  divsd  D68(%rip),%xmm1
  addsd  0x28(%rsp),%xmm1
  movsd  %xmm1,0x88(%rsp)
  divsd  D317(%rip),%xmm0
  movq   %r12,%xmm5
  addsd  %xmm5,%xmm0
  movsd  %xmm0,0x90(%rsp)
L31:
  movsd  0x88(%rsp),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movsd  0x90(%rsp),%xmm0
  movsd  %xmm0,(%rbx)
  mov    0xa8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L4
  add    $0xb8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L77:
  mov    D6(%rip),%r15
  jmp    L41
L50:
  lea    0xa0(%rsp),%rsi
  lea    0x98(%rsp),%rdi
  movq   %r12,%xmm1
  movapd %xmm6,%xmm0
  call   <CelestialToEcliptical@plt>
  movsd  0x98(%rsp),%xmm6
  movsd  0xa0(%rsp),%xmm4
  movsd  %xmm4,0x20(%rsp)
  lea    0x78(%rsp),%r12
  lea    0x70(%rsp),%r13
  movsd  0x48(%rsp),%xmm0
  movsd  %xmm6,0x18(%rsp)
  subsd  %xmm6,%xmm0
  mulsd  D80(%rip),%xmm0
  divsd  D241(%rip),%xmm0
  mov    %r13,%rsi
  mov    %r12,%rdi
  call   <sincos@plt>
  mov    0x70(%rsp),%r14
  movsd  0x78(%rsp),%xmm7
  movsd  %xmm7,0x28(%rsp)
  movsd  0x10(%rsp),%xmm0
  subsd  0x18(%rsp),%xmm0
  mulsd  D80(%rip),%xmm0
  divsd  D241(%rip),%xmm0
  mov    %r13,%rsi
  mov    %r12,%rdi
  call   <sincos@plt>
  movsd  0x78(%rsp),%xmm7
  movsd  %xmm7,0x30(%rsp)
  movq   %r14,%xmm0
  mulsd  0x8(%rsp),%xmm0
  subsd  0x70(%rsp),%xmm0
  movq   %xmm0,%r14
  movsd  D80(%rip),%xmm0
  mulsd  0x20(%rsp),%xmm0
  divsd  D241(%rip),%xmm0
  mov    %r13,%rsi
  mov    %r12,%rdi
  call   <sincos@plt>
  movsd  0x70(%rsp),%xmm1
  movsd  0x78(%rsp),%xmm0
  movapd %xmm1,%xmm3
  movsd  D126(%rip),%xmm2
  comisd %xmm1,%xmm2
  jbe    L32
  movapd %xmm2,%xmm3
L32:
  movq   %r15,%xmm2
  mulsd  D99(%rip),%xmm2
  movq   %r14,%xmm1
  divsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D317(%rip),%xmm3
  divsd  %xmm3,%xmm1
  addsd  0x18(%rsp),%xmm1
  movsd  %xmm1,0x10(%rsp)
  xorpd  D166(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  movsd  0x28(%rsp),%xmm2
  mulsd  0x8(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm5
  subsd  %xmm2,%xmm5
  mulsd  %xmm5,%xmm0
  divsd  %xmm3,%xmm0
  addsd  0x20(%rsp),%xmm0
  call   <Map180@plt>
  movsd  %xmm0,0x8(%rsp)
  comisd D117(%rip),%xmm0
  jbe    L35
  movsd  D241(%rip),%xmm0
  movapd %xmm0,%xmm6
  subsd  0x90(%rsp),%xmm6
  movsd  %xmm6,0x8(%rsp)
  addsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x10(%rsp)
L35:
  movsd  D215(%rip),%xmm0
  movsd  0x8(%rsp),%xmm5
  comisd %xmm5,%xmm0
  jbe    L70
  movsd  D185(%rip),%xmm0
  subsd  %xmm5,%xmm0
  movsd  %xmm0,0x8(%rsp)
  movsd  0x10(%rsp),%xmm4
  addsd  D241(%rip),%xmm4
  movsd  %xmm4,0x10(%rsp)
L70:
  movsd  0x10(%rsp),%xmm0
  call   <Map360@plt>
  lea    0x90(%rsp),%rsi
  lea    0x88(%rsp),%rdi
  movsd  0x8(%rsp),%xmm1
  call   <EclipticalToCelestial@plt>
  jmp    L31
L4:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: label: D143
data_type: f64
name: retrieve_data
arguments: label: D132
data_type: qword[2]
name: retrieve_data
arguments: label: D105
data_type: f64[2]
name: retrieve_data
arguments: label: D80
data_type: f64
name: retrieve_data
arguments: label: D91
data_type: f64
name: retrieve_data
arguments: label: D241
data_type: f64
name: retrieve_data
arguments: label: D166
data_type: byte[16]
name: retrieve_data
arguments: label: D99
data_type: f64
name: retrieve_data
arguments: label: D68
data_type: byte[8]
name: retrieve_data
arguments: label: D317
data_type: f64
name: retrieve_data
arguments: label: D6
data_type: byte[8]
name: retrieve_data
arguments: label: D126
data_type: f64
name: retrieve_data
arguments: label: D117
data_type: f64[2]
name: retrieve_data
arguments: label: D215
data_type: byte[8]
name: retrieve_data
arguments: label: D185
data_type: byte[8]

tool: D143: 1

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

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

tool: D80: 3.14159265358

tool: D91: 12

tool: D241: 180

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

tool: D99: 20.49552

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

tool: D317: 3600

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

tool: D126: 0.0001

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

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

tool: {"D185": [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;
}
```