user: 下面的汇编指令对应的C源代码是什么：

```asm
<pluto_ell>:
  endbr64
  push   %r12
  movapd %xmm0,%xmm1
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x60,%rsp
  movsd  D66(%rip),%xmm4
  movsd  D96(%rip),%xmm3
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  movsd  %xmm0,0x8(%rsp)
  lea    0x38(%rsp),%rdx
  movsd  D87(%rip),%xmm2
  lea    0x40(%rsp),%rsi
  lea    0x30(%rsp),%rdi
  mov    $0x5,%eax
  movsd  D62(%rip),%xmm0
  lea    0x28(%rsp),%rbp
  lea    0x20(%rsp),%r12
  call   <reduce_elements@plt>
  movsd  0x8(%rsp),%xmm1
  lea    0x48(%rsp),%rsi
  subsd  D58(%rip),%xmm1
  movsd  D91(%rip),%xmm0
  mov    D102(%rip),%rax
  lea    0x50(%rsp),%rdi
  mulsd  %xmm1,%xmm0
  movq   %rax,%xmm1
  mov    $0x2,%eax
  mulsd  D37(%rip),%xmm0
  divsd  D25(%rip),%xmm0
  call   <anomaly@plt>
  movsd  0x48(%rsp),%xmm0
  call   <cos@plt>
  mulsd  D102(%rip),%xmm0
  mov    %r12,%rsi
  mov    %rbp,%rdi
  movsd  D6(%rip),%xmm1
  subsd  %xmm0,%xmm1
  mulsd  D9(%rip),%xmm1
  movsd  0x40(%rsp),%xmm0
  addsd  0x50(%rsp),%xmm0
  movsd  %xmm1,0x10(%rbx)
  call   <sincos@plt>
  movsd  0x20(%rsp),%xmm1
  mov    %r12,%rsi
  mov    %rbp,%rdi
  movsd  0x28(%rsp),%xmm5
  movsd  0x30(%rsp),%xmm0
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm5,0x8(%rsp)
  call   <sincos@plt>
  movsd  0x8(%rsp),%xmm0
  mulsd  0x28(%rsp),%xmm0
  movsd  0x20(%rsp),%xmm6
  movsd  %xmm6,0x10(%rsp)
  call   <asin@plt>
  movsd  0x10(%rsp),%xmm6
  mulsd  0x8(%rsp),%xmm6
  movsd  0x18(%rsp),%xmm1
  movsd  %xmm0,0x8(%rbx)
  movapd %xmm6,%xmm0
  call   <atan2@plt>
  addsd  0x38(%rsp),%xmm0
  movsd  %xmm0,(%rbx)
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L4
  add    $0x60,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
L4:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: label: D66
type_guessed: byte[8]
name: unpack_data
arguments: label: D96
type_guessed: byte[8]
name: unpack_data
arguments: label: D87
type_guessed: byte[8]
name: unpack_data
arguments: label: D62
type_guessed: f64
name: unpack_data
arguments: label: D58
type_guessed: byte[8]
name: unpack_data
arguments: label: D91
type_guessed: byte[8]
name: unpack_data
arguments: label: D102
type_guessed: f64
name: unpack_data
arguments: label: D37
type_guessed: f64
name: unpack_data
arguments: label: D25
type_guessed: f64
name: unpack_data
arguments: label: D6
type_guessed: f64
name: unpack_data
arguments: label: D9
type_guessed: f64

tool: {"D66": [0x73,0xaf,0x54,0xd7,0xb3,0xcd,0xfe,0x3f]}

tool: {"D96": [0x36,0x29,0xfd,0xf5,0x1f,0xc5,0xff,0x3f]}

tool: {"D87": [0x30,0x2e,0xb0,0xaa,0x43,0x25,0xd3,0x3f]}

tool: D62: 36525

tool: {"D58": [0x80,0xed,0x7c,0x3f,0x31,0x20,0xe0,0x40]}

tool: {"D91": [0x21,0x00,0x11,0x8d,0x93,0x41,0x70,0x3f]}

tool: D102: 0.249

tool: D37: 3.141592653589793

tool: D25: 180

tool: D6: 1

tool: D9: 39.543

assistant: ```c
void pluto_ell(double mj, double *ret) {
  double a = 39.543, e = 0.2490, inc0 = 17.140, Om0 = 110.307, omeg0 = 113.768, mjp = 2448045.539 - 2415020.0, mjeq = (2451545.0 - 2415020.0), n = 144.9600 / 36525.;
  double inc, Om, omeg;
  double ma, ea, nu;
  double lo, slo, clo;
  reduce_elements(mjeq, mj, ((inc0)*3.141592653589793 / 180.), ((omeg0)*3.141592653589793 / 180.), ((Om0)*3.141592653589793 / 180.), &inc, &omeg, &Om);
  ma = (((mj - mjp) * n) * 3.141592653589793 / 180.);
  anomaly(ma, e, &nu, &ea);
  ret[2] = a * (1.0 - e * cos(ea));
  lo = omeg + nu;
  slo = sin(lo);
  clo = cos(lo);
  ret[1] = asin(slo * sin(inc));
  ret[0] = atan2(slo * cos(inc), clo) + Om;
}
```