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

```asm
<lomb_scargle>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r12
  push   %rbx
  sub    $0x1c0,%rsp
  mov    %edi,-0xb4(%rbp)
  mov    %esi,-0xb8(%rbp)
  mov    %edx,-0xbc(%rbp)
  mov    %ecx,-0xc0(%rbp)
  mov    %r8,-0xc8(%rbp)
  mov    %r9,-0xd0(%rbp)
  movsd  %xmm0,-0x140(%rbp)
  movsd  %xmm1,-0x148(%rbp)
  movsd  %xmm2,-0x150(%rbp)
  movsd  %xmm3,-0x158(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xd8(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xe0(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0xe8(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0xf0(%rbp)
  mov    0x30(%rbp),%rax
  mov    %rax,-0xf8(%rbp)
  mov    0x38(%rbp),%rax
  mov    %rax,-0x100(%rbp)
  mov    0x40(%rbp),%rax
  mov    %rax,-0x108(%rbp)
  mov    0x48(%rbp),%rax
  mov    %rax,-0x110(%rbp)
  mov    0x50(%rbp),%rax
  mov    %rax,-0x118(%rbp)
  mov    0x58(%rbp),%rax
  mov    %rax,-0x120(%rbp)
  mov    0x60(%rbp),%rax
  mov    %rax,-0x128(%rbp)
  mov    0x68(%rbp),%rax
  mov    %rax,-0x130(%rbp)
  mov    0x70(%rbp),%rax
  mov    %rax,-0x138(%rbp)
  mov    0x78(%rbp),%rax
  mov    %rax,-0x160(%rbp)
  mov    0x80(%rbp),%rax
  mov    %rax,-0x168(%rbp)
  mov    0x88(%rbp),%rax
  mov    %rax,-0x170(%rbp)
  mov    0x90(%rbp),%rax
  mov    %rax,-0x178(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x28(%rbp)
  xor    %eax,%eax
  mov    %rsp,%rax
  mov    %rax,%r12
  mov    -0xbc(%rbp),%eax
  add    %eax,%eax
  mov    %eax,-0xa4(%rbp)
  movsd  -0x148(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0xa8(%rbp)
  movq   $0x0,-0x88(%rbp)
  mov    -0xa8(%rbp),%edx
  mov    -0x178(%rbp),%rax
  mov    %edx,(%rax)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  mov    -0xb4(%rbp),%edx
  movslq %edx,%rax
  sub    $0x1,%rax
  mov    %rax,-0x70(%rbp)
  movslq %edx,%rax
  mov    %rax,%r10
  mov    $0x0,%r11d
  movslq %edx,%rax
  mov    %rax,%r14
  mov    $0x0,%r15d
  movslq %edx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%ebx
  mov    $0x0,%edx
  div    %rbx
  imul   $0x10,%rax,%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L11:
  cmp    %rcx,%rsp
  je     L54
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L11
L54:
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  mov    %rdx,%rax
  and    $0xfff,%eax
  test   %rax,%rax
  je     L71
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L71:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x68(%rbp)
  mov    -0xb4(%rbp),%edx
  movslq %edx,%rax
  sub    $0x1,%rax
  mov    %rax,-0x60(%rbp)
  movslq %edx,%rax
  mov    %rax,-0x190(%rbp)
  movq   $0x0,-0x188(%rbp)
  movslq %edx,%rax
  mov    %rax,-0x1a0(%rbp)
  movq   $0x0,-0x198(%rbp)
  movslq %edx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%ebx
  mov    $0x0,%edx
  div    %rbx
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L102:
  cmp    %rdx,%rsp
  je     L49
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L102
L49:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L104
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L104:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x58(%rbp)
  mov    -0xb4(%rbp),%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x50(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x1b0(%rbp)
  movq   $0x0,-0x1a8(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x1c0(%rbp)
  movq   $0x0,-0x1b8(%rbp)
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%ebx
  mov    $0x0,%edx
  div    %rbx
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L169:
  cmp    %rdx,%rsp
  je     L184
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L169
L184:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L53
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L53:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x48(%rbp)
  mov    -0xb4(%rbp),%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x40(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x1d0(%rbp)
  movq   $0x0,-0x1c8(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x1e0(%rbp)
  movq   $0x0,-0x1d8(%rbp)
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%ebx
  mov    $0x0,%edx
  div    %rbx
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L7:
  cmp    %rdx,%rsp
  je     L166
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L7
L166:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L129
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L129:
  mov    %rsp,%rax
  add    $0x7,%rax
  shr    $0x3,%rax
  shl    $0x3,%rax
  mov    %rax,-0x38(%rbp)
  movq   $0x0,-0x90(%rbp)
  jmp    L137
L76:
  mov    -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  lea    (%rdx,%rax,1),%rbx
  mov    -0xd8(%rbp),%r8
  mov    -0xe8(%rbp),%rdi
  mov    -0xe0(%rbp),%rcx
  mov    -0xd0(%rbp),%rdx
  mov    -0xc0(%rbp),%esi
  mov    -0xb4(%rbp),%eax
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <do_lomb@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x78(%rbp),%xmm0
  jbe    L105
  pxor   %xmm0,%xmm0
  ucomisd -0x80(%rbp),%xmm0
  jp     L105
  pxor   %xmm0,%xmm0
  ucomisd -0x80(%rbp),%xmm0
  jne    L105
  mov    -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  mov    -0x38(%rbp),%rdi
  mov    -0x48(%rbp),%rcx
  mov    -0xe8(%rbp),%rdx
  mov    -0xe0(%rbp),%rsi
  mov    -0xb4(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <copy_sincos@plt>
  mov    -0x90(%rbp),%rax
  mov    %rax,-0x88(%rbp)
L105:
  mov    -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x150(%rbp),%xmm0
  jbe    L56
  mov    -0x148(%rbp),%rdi
  mov    -0x68(%rbp),%r9
  mov    -0xe8(%rbp),%r10
  mov    -0xe0(%rbp),%rcx
  mov    -0xd0(%rbp),%rdx
  mov    -0xc0(%rbp),%esi
  mov    -0xb4(%rbp),%eax
  sub    $0x8,%rsp
  lea    -0xa8(%rbp),%r8
  push   %r8
  push   -0xd8(%rbp)
  push   -0x118(%rbp)
  push   -0x110(%rbp)
  push   -0x108(%rbp)
  push   -0x100(%rbp)
  push   -0x58(%rbp)
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %eax,%edi
  call   <do_lomb_zoom@plt>
  add    $0x40,%rsp
  movq   %xmm0,%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x160(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  mov    -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  lea    (%rdx,%rax,1),%rbx
  movsd  -0x158(%rbp),%xmm0
  mov    -0x140(%rbp),%rdi
  mov    -0x128(%rbp),%r9
  mov    -0x130(%rbp),%r10
  mov    -0x120(%rbp),%rcx
  mov    -0xc0(%rbp),%edx
  mov    -0xbc(%rbp),%esi
  mov    -0xb4(%rbp),%eax
  sub    $0x8,%rsp
  push   $0x0
  lea    -0xa0(%rbp),%r8
  push   %r8
  push   -0x168(%rbp)
  lea    -0x98(%rbp),%r8
  push   %r8
  push   -0x138(%rbp)
  push   -0xd0(%rbp)
  push   -0xd8(%rbp)
  push   -0x58(%rbp)
  push   -0x68(%rbp)
  movapd %xmm0,%xmm1
  movq   %rdi,%xmm0
  mov    %r10,%r8
  mov    %eax,%edi
  call   <refine_psd@plt>
  add    $0x50,%rsp
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x80(%rbp),%xmm0
  jbe    L56
  mov    -0x38(%rbp),%rdi
  mov    -0x48(%rbp),%rcx
  mov    -0x58(%rbp),%rdx
  mov    -0x68(%rbp),%rsi
  mov    -0xb4(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <copy_sincos@plt>
  mov    -0x90(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  mov    -0xa8(%rbp),%edx
  mov    -0x178(%rbp),%rax
  mov    %edx,(%rax)
  mov    -0x90(%rbp),%rax
  mov    %rax,-0x88(%rbp)
L56:
  mov    -0xe8(%rbp),%rdi
  mov    -0xe0(%rbp),%rcx
  mov    -0xf8(%rbp),%rdx
  mov    -0xf0(%rbp),%rsi
  mov    -0xb4(%rbp),%eax
  mov    $0x0,%r9d
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <update_sincos@plt>
  addq   $0x1,-0x90(%rbp)
L137:
  mov    -0xb8(%rbp),%eax
  cltq
  cmp    %rax,-0x90(%rbp)
  jb     L76
  mov    -0x88(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xc8(%rbp),%rax
  lea    (%rdx,%rax,1),%rbx
  movsd  -0x158(%rbp),%xmm0
  mov    -0x140(%rbp),%rdi
  mov    -0x128(%rbp),%r9
  mov    -0x130(%rbp),%r8
  mov    -0x120(%rbp),%rcx
  mov    -0xc0(%rbp),%edx
  mov    -0xbc(%rbp),%esi
  mov    -0xb4(%rbp),%eax
  sub    $0x8,%rsp
  push   $0x1
  push   -0x170(%rbp)
  push   -0x168(%rbp)
  push   -0x160(%rbp)
  push   -0x138(%rbp)
  push   -0xd0(%rbp)
  push   -0xd8(%rbp)
  push   -0x38(%rbp)
  push   -0x48(%rbp)
  movapd %xmm0,%xmm1
  movq   %rdi,%xmm0
  mov    %eax,%edi
  call   <refine_psd@plt>
  add    $0x50,%rsp
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    %r12,%rsp
  nop
  mov    -0x28(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L27
  call   <__stack_chk_fail@plt>
L27:
  lea    -0x20(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

assistant: ```c
void lomb_scargle(int numt, int numf, int nharm, int detrend_order, double psd[], double cn[], double wth[], double sinx[], double cosx[], double sinx_step[], double cosx_step[], double sinx_back[], double cosx_back[], double sinx_smallstep[], double cosx_smallstep[], double hat_matr[], double hat_hat[], double hat0[], double soln[], double chi0, double freq_zoom, double psdmin, double tone_control, double lambda0[], double lambda0_range[], double Tr[], int ifreq[]) {
  int i, k, npar = 2 * (int)nharm, ifr = (int)(freq_zoom) / 2;
  unsigned long j;
  unsigned long jmax = 0;
  *ifreq = ifr;
  double psdmax = 0., psd0max = 0., Trace, sinx1[numt], cosx1[numt], sinx2[numt], cosx2[numt], lambda, px, pxmax;
  for (j = 0; j < numf; j++) {
    psd[j] = do_lomb(numt, detrend_order, cn, sinx, cosx, wth);
    if (psd[j] > psd0max && psdmax == 0) {
      psd0max = psd[j];
      copy_sincos(numt, sinx, cosx, sinx2, cosx2);
      jmax = j;
    }
    if (psd[j] > (double)psdmin) {
      px = do_lomb_zoom(numt, detrend_order, cn, sinx, cosx, sinx1, cosx1, sinx_back, cosx_back, sinx_smallstep, cosx_smallstep, wth, freq_zoom, &ifr);
      lambda = *lambda0;
      psd[j] = refine_psd(numt, nharm, detrend_order, hat_matr, hat0, hat_hat, sinx1, cosx1, wth, cn, soln, &lambda, lambda0_range, chi0, tone_control, &Trace, 0);
      if (psd[j] > psdmax) {
        copy_sincos(numt, sinx1, cosx1, sinx2, cosx2);
        psdmax = psd[j];
        *ifreq = ifr;
        jmax = j;
      }
    }
    update_sincos(numt, sinx_step, cosx_step, sinx, cosx, 0);
  }
  psd[jmax] = refine_psd(numt, nharm, detrend_order, hat_matr, hat0, hat_hat, sinx2, cosx2, wth, cn, soln, lambda0, lambda0_range, chi0, tone_control, Tr, 1);
}
```