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

```asm
<lomb_scargle>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x118,%rsp
  mov    %edi,%r13d
  mov    %edx,-0x60(%rbp)
  mov    %ecx,-0x5c(%rbp)
  mov    %r8,%r14
  mov    %r9,-0x68(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x70(%rbp)
  mov    0x18(%rbp),%r12
  mov    0x20(%rbp),%r15
  mov    0x28(%rbp),%rax
  mov    %rax,-0x80(%rbp)
  mov    0x30(%rbp),%rax
  mov    %rax,-0x88(%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,-0xb8(%rbp)
  mov    0x60(%rbp),%rax
  mov    %rax,-0xc0(%rbp)
  mov    0x68(%rbp),%rax
  mov    %rax,-0xc8(%rbp)
  mov    0x70(%rbp),%rax
  mov    %rax,-0xd0(%rbp)
  movsd  %xmm0,-0xd8(%rbp)
  movsd  %xmm1,-0x120(%rbp)
  movsd  %xmm2,-0x90(%rbp)
  movsd  %xmm3,-0xe0(%rbp)
  mov    0x78(%rbp),%rax
  mov    %rax,-0xe8(%rbp)
  mov    0x80(%rbp),%rax
  mov    %rax,-0xf0(%rbp)
  mov    0x88(%rbp),%rax
  mov    %rax,-0x138(%rbp)
  mov    0x90(%rbp),%rcx
  mov    %rcx,-0x140(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  cvttsd2si %xmm1,%edx
  mov    %edx,%eax
  shr    $0x1f,%eax
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x4c(%rbp)
  mov    %eax,(%rcx)
  movslq %edi,%rdx
  shl    $0x3,%rdx
  lea    0xf(%rdx),%rax
  mov    %rax,%rdi
  and    $0xfffffffffffffff0,%rdi
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L90:
  cmp    %rcx,%rsp
  je     L80
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L90
L80:
  mov    %rdi,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L121
  orq    $0x0,-0x8(%rsp,%rax,1)
L121:
  mov    %rsp,-0xa8(%rbp)
  lea    0xf(%rdx),%rax
  mov    %rax,%rdi
  and    $0xfffffffffffffff0,%rdi
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L50:
  cmp    %rcx,%rsp
  je     L156
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L50
L156:
  mov    %rdi,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L63
  orq    $0x0,-0x8(%rsp,%rax,1)
L63:
  mov    %rsp,-0xb0(%rbp)
  lea    0xf(%rdx),%rax
  mov    %rax,%rdi
  and    $0xfffffffffffffff0,%rdi
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L17:
  cmp    %rcx,%rsp
  je     L97
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L17
L97:
  mov    %rdi,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L91
  orq    $0x0,-0x8(%rsp,%rax,1)
L91:
  mov    %rsp,-0x128(%rbp)
  add    $0xf,%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffffff0,%rax
  and    $0xfffffffffffff000,%rdx
  mov    %rsp,%rcx
  sub    %rdx,%rcx
L78:
  cmp    %rcx,%rsp
  je     L77
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L78
L77:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  je     L46
  orq    $0x0,-0x8(%rsp,%rdx,1)
L46:
  mov    %rsp,-0x130(%rbp)
  movslq %esi,%rax
  mov    %rax,-0x98(%rbp)
  test   %rax,%rax
  je     L143
  pxor   %xmm1,%xmm1
  movsd  %xmm1,-0xa0(%rbp)
  movsd  %xmm1,-0x78(%rbp)
  movq   $0x0,-0xf8(%rbp)
  mov    $0x0,%ebx
  jmp    L111
L102:
  movsd  (%r14,%rbx,8),%xmm0
  comisd -0x90(%rbp),%xmm0
  ja     L70
L194:
  mov    $0x0,%r9d
  mov    %r15,%r8
  mov    %r12,%rcx
  mov    -0x88(%rbp),%rdx
  mov    -0x80(%rbp),%rsi
  mov    %r13d,%edi
  call   <update_sincos@plt>
  add    $0x1,%rbx
  cmp    -0x98(%rbp),%rbx
  je     L142
L111:
  mov    -0x70(%rbp),%r9
  mov    %r15,%r8
  mov    %r12,%rcx
  mov    -0x68(%rbp),%rdx
  mov    -0x5c(%rbp),%esi
  mov    %r13d,%edi
  call   <do_lomb@plt>
  movsd  %xmm0,-0x58(%rbp)
  movsd  %xmm0,(%r14,%rbx,8)
  movsd  -0x78(%rbp),%xmm5
  pxor   %xmm6,%xmm6
  ucomisd %xmm6,%xmm5
  setnp  %al
  mov    $0x0,%edx
  cmovne %edx,%eax
  test   %al,%al
  je     L102
  comisd -0xa0(%rbp),%xmm0
  jbe    L102
  mov    -0x130(%rbp),%r8
  mov    -0x128(%rbp),%rcx
  mov    %r15,%rdx
  mov    %r12,%rsi
  mov    %r13d,%edi
  call   <copy_sincos@plt>
  movsd  -0x58(%rbp),%xmm3
  movsd  %xmm3,-0xa0(%rbp)
  mov    %rbx,-0xf8(%rbp)
  jmp    L102
L70:
  sub    $0x8,%rsp
  lea    -0x4c(%rbp),%rax
  push   %rax
  push   -0x70(%rbp)
  push   -0x118(%rbp)
  push   -0x110(%rbp)
  push   -0x108(%rbp)
  push   -0x100(%rbp)
  push   -0xb0(%rbp)
  movsd  -0x120(%rbp),%xmm0
  mov    -0xa8(%rbp),%r9
  mov    %r15,%r8
  mov    %r12,%rcx
  mov    -0x68(%rbp),%rdx
  mov    -0x5c(%rbp),%esi
  mov    %r13d,%edi
  call   <do_lomb_zoom@plt>
  mov    -0xe8(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  add    $0x38,%rsp
  push   $0x0
  lea    -0x48(%rbp),%rax
  push   %rax
  push   -0xf0(%rbp)
  lea    -0x40(%rbp),%rax
  push   %rax
  push   -0xd0(%rbp)
  push   -0x68(%rbp)
  push   -0x70(%rbp)
  push   -0xb0(%rbp)
  push   -0xa8(%rbp)
  movsd  -0xe0(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  mov    -0xc0(%rbp),%r9
  mov    -0xc8(%rbp),%r8
  mov    -0xb8(%rbp),%rcx
  mov    -0x5c(%rbp),%edx
  mov    -0x60(%rbp),%esi
  mov    %r13d,%edi
  call   <refine_psd@plt>
  movsd  %xmm0,(%r14,%rbx,8)
  add    $0x50,%rsp
  comisd -0x78(%rbp),%xmm0
  jbe    L194
  mov    -0x130(%rbp),%r8
  mov    -0x128(%rbp),%rcx
  mov    -0xb0(%rbp),%rdx
  mov    -0xa8(%rbp),%rsi
  mov    %r13d,%edi
  call   <copy_sincos@plt>
  movsd  (%r14,%rbx,8),%xmm3
  movsd  %xmm3,-0x78(%rbp)
  mov    -0x4c(%rbp),%eax
  mov    -0x140(%rbp),%rsi
  mov    %eax,(%rsi)
  mov    %rbx,-0xf8(%rbp)
  jmp    L194
L143:
  mov    -0x98(%rbp),%rax
  mov    %rax,-0xf8(%rbp)
L142:
  mov    -0xf8(%rbp),%rax
  lea    (%r14,%rax,8),%rbx
  sub    $0x8,%rsp
  push   $0x1
  push   -0x138(%rbp)
  push   -0xf0(%rbp)
  push   -0xe8(%rbp)
  push   -0xd0(%rbp)
  push   -0x68(%rbp)
  push   -0x70(%rbp)
  push   -0x130(%rbp)
  push   -0x128(%rbp)
  movsd  -0xe0(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  mov    -0xc0(%rbp),%r9
  mov    -0xc8(%rbp),%r8
  mov    -0xb8(%rbp),%rcx
  mov    -0x5c(%rbp),%edx
  mov    -0x60(%rbp),%esi
  mov    %r13d,%edi
  call   <refine_psd@plt>
  movsd  %xmm0,(%rbx)
  add    $0x50,%rsp
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L191
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L191:
  call   <__stack_chk_fail@plt>
```

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);
}
```