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

```asm
<passf>:
  endbr64
  push   %r15
  mov    %r8d,%r15d
  push   %r14
  mov    %esi,%r14d
  mov    %ecx,%esi
  push   %r13
  mov    %r14d,%eax
  push   %r12
  shr    $0x1f,%eax
  push   %rbp
  add    %r14d,%eax
  push   %rbx
  sar    %eax
  mov    %edx,%ebx
  sub    $0x20,%rsp
  mov    %edx,-0x44(%rsp)
  add    $0x1,%edx
  mov    0x68(%rsp),%rbp
  mov    %eax,0x14(%rsp)
  mov    %edx,%eax
  shr    $0x1f,%eax
  mov    %rdi,0x8(%rsp)
  add    %edx,%eax
  mov    %r9,-0x70(%rsp)
  sar    %eax
  mov    %eax,-0x74(%rsp)
  mov    %ebx,%eax
  imul   %r14d,%ebx
  mov    %ebx,0x4(%rsp)
  cmp    %ecx,%r14d
  jge    L584
  cmp    $0x2,%eax
  jle    L510
  mov    -0x44(%rsp),%eax
  mov    %ecx,-0x60(%rsp)
  movl   $0x1,-0x68(%rsp)
  sub    $0x1,%eax
  mov    %rbp,0x68(%rsp)
  mov    %eax,%ebx
  mov    %eax,-0x40(%rsp)
  imul   %r14d,%eax
  imul   %ecx,%ebx
  mov    %r15d,-0x28(%rsp)
  mov    %r14d,%r15d
  mov    %eax,-0x54(%rsp)
  mov    %ecx,%eax
  imul   %r14d,%eax
  mov    %ebx,-0x58(%rsp)
  movslq 0x4(%rsp),%rbx
  mov    %rbx,-0x38(%rsp)
  lea    0x0(,%rbx,8),%r8
  movslq %r14d,%rbx
  mov    %eax,-0x48(%rsp)
  lea    0x0(,%rbx,8),%r10
  mov    %eax,%r13d
  mov    %rbx,-0x50(%rsp)
L74:
  test   %r14d,%r14d
  jle    L238
  movslq -0x54(%rsp),%rax
  mov    -0x70(%rsp),%rbx
  movslq %r15d,%r12
  movslq %r13d,%rdx
  mov    -0x58(%rsp),%r9d
  lea    (%rbx,%rax,8),%r11
  sub    %rax,%r12
  mov    -0x60(%rsp),%eax
  mov    0x68(%rsp),%rbx
  imul   %r14d,%r9d
  shl    $0x3,%r12
  imul   %r14d,%eax
  lea    (%rbx,%rdx,8),%rbp
  xor    %ebx,%ebx
  movslq %r9d,%r9
  cltq
  sub    %rax,%r9
  cs nopw 0x0(%rax,%rax,1)
L470:
  test   %esi,%esi
  jle    L262
  lea    (%r11,%r12,1),%rcx
  mov    %rbp,%rdx
  mov    %r11,%rax
  xor    %edi,%edi
L394:
  movsd  (%rcx),%xmm0
  addsd  (%rax),%xmm0
  add    $0x1,%edi
  movsd  %xmm0,(%rdx)
  movsd  (%rcx),%xmm0
  subsd  (%rax),%xmm0
  add    %r8,%rcx
  add    %r8,%rax
  movsd  %xmm0,(%rdx,%r9,8)
  add    %r10,%rdx
  cmp    %edi,%esi
  jne    L394
L262:
  add    $0x1,%ebx
  add    $0x8,%r11
  add    $0x8,%rbp
  cmp    %ebx,%r14d
  jne    L470
L238:
  addl   $0x1,-0x68(%rsp)
  add    -0x48(%rsp),%r13d
  add    %r14d,%r15d
  add    %esi,-0x60(%rsp)
  mov    -0x68(%rsp),%eax
  sub    %esi,-0x58(%rsp)
  sub    %r14d,-0x54(%rsp)
  cmp    %eax,-0x74(%rsp)
  jg     L74
  mov    -0x28(%rsp),%r15d
  mov    0x68(%rsp),%rbp
  test   %r14d,%r14d
  jle    L488
L217:
  mov    -0x38(%rsp),%r8
  mov    -0x50(%rsp),%rdi
  mov    %rbp,%r10
  xor    %r11d,%r11d
  mov    -0x70(%rsp),%r9
  shl    $0x3,%r8
  shl    $0x3,%rdi
L6:
  mov    %r10,%rcx
  mov    %r9,%rdx
  xor    %eax,%eax
  test   %esi,%esi
  jle    L163
  xchg   %ax,%ax
L60:
  movsd  (%rdx),%xmm0
  add    $0x1,%eax
  add    %r8,%rdx
  movsd  %xmm0,(%rcx)
  add    %rdi,%rcx
  cmp    %eax,%esi
  jne    L60
L163:
  add    $0x1,%r11d
  add    $0x8,%r9
  add    $0x8,%r10
  cmp    %r11d,%r14d
  jne    L6
L210:
  mov    $0x2,%eax
  sub    %r14d,%eax
  mov    %eax,(%rsp)
  mov    -0x44(%rsp),%eax
  cmp    $0x2,%eax
  jle    L386
  movslq %r14d,%rbx
  mov    %rbx,-0x50(%rsp)
  lea    -0x2(%rax),%ebx
  sub    $0x1,%eax
  mov    %ebx,-0x28(%rsp)
  mov    %eax,-0x40(%rsp)
L559:
  mov    -0x50(%rsp),%rdi
  mov    -0x40(%rsp),%ebx
  mov    %r15d,-0x54(%rsp)
  mov    %r15d,-0x50(%rsp)
  mov    0x78(%rsp),%r12
  imul   %r15d,%ebx
  lea    0x0(,%rdi,8),%rax
  mov    %esi,0x1c(%rsp)
  mov    %rax,-0x30(%rsp)
  movslq (%rsp),%rax
  movl   $0x0,-0x60(%rsp)
  mov    %rax,%rcx
  add    %rdi,%rax
  mov    0x78(%rsp),%rdi
  mov    %ebx,-0x48(%rsp)
  mov    %ebx,0x10(%rsp)
  lea    -0x8(%rdi,%rax,8),%r13
  mov    0x70(%rsp),%rdi
  movslq %r15d,%rax
  mov    %ebx,-0x58(%rsp)
  movl   $0x1,-0x68(%rsp)
  lea    (%rdi,%rax,8),%rax
  mov    %rbp,0x68(%rsp)
  mov    0x4(%rsp),%ebp
  mov    %rax,-0x18(%rsp)
  movslq %ebx,%rax
  lea    (%rdi,%rax,8),%rax
  mov    %rax,-0x20(%rsp)
  lea    (%r15,%r15,1),%eax
  mov    %eax,-0xc(%rsp)
  mov    -0x28(%rsp),%eax
  mov    %r14d,-0x28(%rsp)
  imul   %r15d,%eax
  mov    %eax,-0x10(%rsp)
  lea    (%r14,%rcx,1),%eax
  mov    %rdi,%r14
  mov    %eax,-0x24(%rsp)
  lea    -0x1(%r15),%eax
  mov    %eax,0x18(%rsp)
  mov    %rax,-0x40(%rsp)
  mov    0x60(%rsp),%rax
  add    $0x8,%rax
  mov    %rax,-0x8(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L621:
  mov    -0x24(%rsp),%r8d
  add    -0x60(%rsp),%r8d
  test   %r15d,%r15d
  jle    L377
  movslq -0x54(%rsp),%rax
  mov    0x60(%rsp),%rbx
  pxor   %xmm1,%xmm1
  cvtsi2sdl 0x80(%rsp),%xmm1
  mov    -0x20(%rsp),%rdx
  mov    -0x18(%rsp),%rdi
  lea    (%rbx,%rax,8),%rsi
  movslq -0x58(%rsp),%rax
  mov    -0x40(%rsp),%r9
  lea    (%rbx,%rax,8),%rcx
  xor    %eax,%eax
  nopl   0x0(%rax)
L143:
  movsd  -0x8(%r13),%xmm0
  mulsd  (%rdi,%rax,8),%xmm0
  mov    %rax,%r10
  addsd  (%r14,%rax,8),%xmm0
  movsd  %xmm0,(%rsi,%rax,8)
  movsd  0x0(%r13),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  (%rdx,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %r10,%r9
  jne    L143
L377:
  mov    -0x28(%rsp),%ebx
  add    %ebx,-0x60(%rsp)
  cmpl   $0x4,-0x44(%rsp)
  jle    L235
  movslq -0x50(%rsp),%rax
  mov    0x60(%rsp),%rbx
  mov    %r13,-0x38(%rsp)
  mov    $0x2,%r9d
  movslq -0x48(%rsp),%rcx
  mov    -0x10(%rsp),%r11d
  lea    (%rbx,%rax,8),%rbx
  mov    -0xc(%rsp),%r10d
  mov    -0x60(%rsp),%r13d
  mov    %rbx,-0x70(%rsp)
  mov    -0x40(%rsp),%rbx
  sub    %rax,%rcx
  lea    (%rbx,%rax,1),%rdx
  mov    -0x8(%rsp),%rbx
  lea    (%rbx,%rdx,8),%rdi
  nopl   0x0(%rax,%rax,1)
L357:
  add    %r13d,%r8d
  mov    %r8d,%eax
  sub    %ebp,%eax
  cmp    %r8d,%ebp
  cmovl  %eax,%r8d
  movslq %r8d,%rax
  shl    $0x3,%rax
  movsd  -0x10(%r12,%rax,1),%xmm3
  movsd  -0x8(%r12,%rax,1),%xmm0
  test   %r15d,%r15d
  jle    L258
  pxor   %xmm1,%xmm1
  movslq %r11d,%rbx
  movslq %r10d,%rsi
  mov    -0x70(%rsp),%rax
  cvtsi2sdl 0x80(%rsp),%xmm1
  lea    (%r14,%rbx,8),%rdx
  sub    %rbx,%rsi
  mulsd  %xmm1,%xmm0
  nopl   0x0(%rax)
L71:
  movsd  (%rdx,%rsi,8),%xmm1
  add    $0x8,%rdx
  mulsd  %xmm3,%xmm1
  addsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  movsd  -0x8(%rdx),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  (%rax,%rcx,8),%xmm1
  movsd  %xmm1,(%rax,%rcx,8)
  add    $0x8,%rax
  cmp    %rax,%rdi
  jne    L71
L258:
  add    $0x1,%r9d
  add    %r15d,%r10d
  sub    %r15d,%r11d
  cmp    %r9d,-0x74(%rsp)
  jg     L357
  mov    -0x38(%rsp),%r13
L235:
  addl   $0x1,-0x68(%rsp)
  add    -0x30(%rsp),%r13
  add    %r15d,-0x50(%rsp)
  mov    -0x68(%rsp),%eax
  sub    %r15d,-0x48(%rsp)
  add    %r15d,-0x54(%rsp)
  sub    %r15d,-0x58(%rsp)
  cmp    %eax,-0x74(%rsp)
  jg     L621
  mov    0x18(%rsp),%eax
  mov    0x1c(%rsp),%esi
  mov    %r15d,%r8d
  mov    $0x1,%edi
  mov    0x70(%rsp),%rbx
  mov    -0x28(%rsp),%r14d
  mov    0x68(%rsp),%rbp
  mov    -0x74(%rsp),%r10d
  lea    0x8(%rbx,%rax,8),%rcx
  mov    %rbx,%r9
  nopl   0x0(%rax,%rax,1)
L443:
  mov    %r9,%rax
  movslq %r8d,%rdx
  test   %r15d,%r15d
  jle    L83
  nopl   0x0(%rax,%rax,1)
L144:
  movsd  (%rax),%xmm0
  addsd  (%rax,%rdx,8),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rcx
  jne    L144
L83:
  add    $0x1,%edi
  add    %r15d,%r8d
  cmp    %edi,%r10d
  jg     L443
  lea    -0x2(%r15),%eax
  mov    %esi,-0x68(%rsp)
  mov    0x10(%rsp),%r9d
  mov    %r15d,%r11d
  shr    %eax
  mov    %rbp,0x68(%rsp)
  mov    -0x74(%rsp),%r13d
  mov    $0x1,%r10d
  lea    (%rax,%rax,1),%r12
  mov    0x70(%rsp),%rax
  mov    0x60(%rsp),%rbp
  mov    %r14d,-0x70(%rsp)
  lea    0x10(%rax),%rbx
  mov    %rax,%r14
  cs nopw 0x0(%rax,%rax,1)
L558:
  cmp    $0x1,%r15d
  jle    L548
  movslq %r11d,%r8
  movslq %r9d,%rsi
  lea    0x0(,%r8,8),%rax
  lea    (%r12,%r8,1),%rdi
  lea    0x0(%rbp,%rax,1),%rcx
  lea    0x0(%rbp,%rsi,8),%rdx
  add    %r14,%rax
  sub    %r8,%rsi
  lea    (%rbx,%rdi,8),%rdi
  nopl   0x0(%rax)
L409:
  movsd  (%rcx),%xmm0
  subsd  0x8(%rdx),%xmm0
  add    $0x10,%rcx
  add    $0x10,%rdx
  movsd  %xmm0,(%rax)
  movsd  -0x10(%rcx),%xmm0
  addsd  -0x8(%rdx),%xmm0
  movsd  %xmm0,(%rax,%rsi,8)
  movsd  -0x8(%rcx),%xmm0
  addsd  -0x10(%rdx),%xmm0
  movsd  %xmm0,0x8(%rax)
  movsd  -0x8(%rcx),%xmm0
  subsd  -0x10(%rdx),%xmm0
  movsd  %xmm0,0x8(%rax,%rsi,8)
  add    $0x10,%rax
  cmp    %rax,%rdi
  jne    L409
L548:
  add    $0x1,%r10d
  add    %r15d,%r11d
  sub    %r15d,%r9d
  cmp    %r10d,%r13d
  jg     L558
  mov    -0x70(%rsp),%r14d
  mov    -0x68(%rsp),%esi
  mov    0x68(%rsp),%rbp
L386:
  mov    0x8(%rsp),%rax
  cmp    $0x2,%r14d
  je     L107
  movl   $0x0,(%rax)
  movslq %r15d,%rdx
  xor    %eax,%eax
  test   %r15d,%r15d
  jle    L37
  mov    0x60(%rsp),%rcx
  mov    0x70(%rsp),%rdi
  cs nopw 0x0(%rax,%rax,1)
L102:
  movsd  (%rdi,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L102
L37:
  mov    -0x44(%rsp),%r9d
  cmp    $0x1,%r9d
  jle    L155
  mov    %esi,%ebx
  movslq %r14d,%rax
  mov    0x58(%rsp),%rdi
  mov    $0x1,%r11d
  imul   %r14d,%ebx
  mov    %rax,-0x50(%rsp)
  lea    0x0(,%rax,8),%r8
  mov    %ebx,%ecx
  nopl   (%rax)
L360:
  movslq %ecx,%rax
  xor    %edx,%edx
  shl    $0x3,%rax
  test   %esi,%esi
  jle    L439
  nopl   (%rax)
L485:
  movsd  0x0(%rbp,%rax,1),%xmm0
  add    $0x1,%edx
  movsd  %xmm0,(%rdi,%rax,1)
  movsd  0x8(%rbp,%rax,1),%xmm0
  movsd  %xmm0,0x8(%rdi,%rax,1)
  add    %r8,%rax
  cmp    %edx,%esi
  jne    L485
L439:
  lea    0x1(%r11),%eax
  add    %ebx,%ecx
  cmp    %eax,%r9d
  je     L5
  mov    %eax,%r11d
  jmp    L360
L584:
  mov    -0x44(%rsp),%edi
  cmp    $0x2,%edi
  jle    L213
  lea    (%r14,%r14,1),%eax
  mov    %r14d,%ebx
  mov    %r14d,-0x54(%rsp)
  neg    %eax
  imul   %ecx,%ebx
  mov    %r8d,-0x30(%rsp)
  mov    %eax,-0x40(%rsp)
  lea    -0x2(%rdi),%eax
  mov    %eax,-0x28(%rsp)
  imul   %r14d,%eax
  mov    %ebx,-0x50(%rsp)
  mov    %ebx,-0x48(%rsp)
  lea    -0x1(%r14),%ebx
  mov    %eax,%r13d
  lea    (%rcx,%rcx,1),%eax
  mov    %rbx,-0x68(%rsp)
  mov    -0x70(%rsp),%rbx
  neg    %eax
  movl   $0x1,-0x58(%rsp)
  imul   %r14d,%eax
  add    $0x8,%rbx
  mov    %rbp,0x68(%rsp)
  mov    0x4(%rsp),%ebp
  mov    %rbx,-0x60(%rsp)
  mov    %eax,-0x38(%rsp)
  mov    %ecx,%eax
  imul   %r13d,%eax
  mov    %eax,%r15d
L23:
  mov    -0x48(%rsp),%r11d
  mov    -0x54(%rsp),%r10d
  xor    %r12d,%r12d
  test   %esi,%esi
  jle    L617
  nopl   (%rax)
L507:
  test   %r14d,%r14d
  jle    L420
  mov    0x68(%rsp),%rdi
  lea    (%r11,%r15,1),%r9d
  movslq %r10d,%rbx
  lea    (%r10,%r13,1),%ecx
  movslq %r9d,%r9
  mov    -0x60(%rsp),%r8
  mov    -0x70(%rsp),%rax
  movslq %ecx,%rcx
  lea    (%rdi,%r9,8),%rdx
  mov    -0x68(%rsp),%rdi
  sub    %rbx,%rcx
  lea    (%rax,%rbx,8),%rax
  add    %rbx,%rdi
  lea    (%r8,%rdi,8),%r8
  movslq %r11d,%rdi
  sub    %r9,%rdi
  nopw   0x0(%rax,%rax,1)
L629:
  movsd  (%rax),%xmm0
  addsd  (%rax,%rcx,8),%xmm0
  movsd  %xmm0,(%rdx,%rdi,8)
  movsd  (%rax),%xmm0
  subsd  (%rax,%rcx,8),%xmm0
  add    $0x8,%rax
  add    $0x8,%rdx
  movsd  %xmm0,-0x8(%rdx)
  cmp    %r8,%rax
  jne    L629
L420:
  add    $0x1,%r12d
  add    %ebp,%r10d
  add    %r14d,%r11d
  cmp    %r12d,%esi
  jne    L507
L617:
  addl   $0x1,-0x58(%rsp)
  mov    -0x50(%rsp),%edi
  add    %r14d,-0x54(%rsp)
  mov    -0x58(%rsp),%eax
  add    %edi,-0x48(%rsp)
  add    -0x40(%rsp),%r13d
  add    -0x38(%rsp),%r15d
  cmp    %eax,-0x74(%rsp)
  jg     L23
  mov    -0x30(%rsp),%r15d
  mov    0x68(%rsp),%rbp
  test   %esi,%esi
  jle    L549
L225:
  mov    -0x44(%rsp),%ebx
  mov    -0x70(%rsp),%r12
  xor    %r9d,%r9d
  xor    %r10d,%r10d
  xor    %r11d,%r11d
  lea    -0x1(%r14),%r8d
  imul   %r14d,%ebx
L215:
  test   %r14d,%r14d
  jle    L109
  movslq %r10d,%rax
  lea    (%r12,%rax,8),%rdi
  movslq %r9d,%rax
  lea    0x0(%rbp,%rax,8),%rcx
  xor    %eax,%eax
  nopl   0x0(%rax)
L366:
  movsd  (%rdi,%rax,8),%xmm0
  mov    %rax,%rdx
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %r8,%rdx
  jne    L366
L109:
  add    $0x1,%r11d
  add    %ebx,%r10d
  add    %r14d,%r9d
  cmp    %r11d,%esi
  jg     L215
  jmp    L210
L448:
  mov    -0x50(%rsp),%rbx
  mov    %esi,%r15d
  mov    $0x1,%r13d
  lea    -0x4(%r14),%r11d
  imul   %r14d,%r15d
  mov    %r13d,-0x74(%rsp)
  shr    %r11d
  mov    0x58(%rsp),%r13
  lea    0x0(,%rbx,8),%rax
  add    %r11,%r11
  mov    %rax,-0x70(%rsp)
  movslq (%rsp),%rax
  mov    %r15d,-0x68(%rsp)
  mov    %r15d,%r12d
  mov    0x80(%rsp),%r15d
  add    %rbx,%rax
  mov    0x78(%rsp),%rbx
  lea    (%rbx,%rax,8),%r10
  lea    0x10(%rbp),%rbx
L502:
  mov    %r12d,%r9d
  xor    %r8d,%r8d
  test   %esi,%esi
  jle    L268
  nop
L137:
  cmp    $0x3,%r14d
  jle    L542
  pxor   %xmm2,%xmm2
  movslq %r9d,%rdx
  cvtsi2sd %r15d,%xmm2
  lea    0x0(,%rdx,8),%rcx
  add    %r11,%rdx
  lea    0x0(%rbp,%rcx,1),%rax
  lea    (%rbx,%rdx,8),%rdi
  add    %r13,%rcx
  mov    %r10,%rdx
  nopl   0x0(%rax)
L264:
  movsd  (%rdx),%xmm0
  add    $0x10,%rax
  add    $0x10,%rdx
  add    $0x10,%rcx
  movsd  -0x8(%rdx),%xmm1
  mulsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  0x8(%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rcx)
  movsd  -0x8(%rdx),%xmm0
  movsd  -0x10(%rdx),%xmm1
  mulsd  0x8(%rax),%xmm1
  mulsd  %xmm2,%xmm0
  mulsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rcx)
  cmp    %rax,%rdi
  jne    L264
L542:
  add    $0x1,%r8d
  add    %r14d,%r9d
  cmp    %r8d,%esi
  jne    L137
L268:
  addl   $0x1,-0x74(%rsp)
  add    -0x70(%rsp),%r10
  mov    -0x74(%rsp),%eax
  add    -0x68(%rsp),%r12d
  cmp    %eax,-0x44(%rsp)
  jne    L502
L155:
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L5:
  cmp    %esi,0x14(%rsp)
  jg     L448
  lea    -0x4(%r14),%eax
  xor    %r12d,%r12d
  mov    %r11d,-0x74(%rsp)
  mov    %ebx,%r15d
  and    $0xfffffffe,%eax
  mov    %ebx,-0x70(%rsp)
  mov    0x58(%rsp),%r11
  mov    $0x1,%r13d
  mov    %eax,-0x68(%rsp)
  mov    0x80(%rsp),%ebx
  lea    0x2(%r12),%eax
  cmp    $0x3,%r14d
  jle    L370
L151:
  mov    0x78(%rsp),%rdi
  cltq
  mov    $0x3,%r10d
  lea    (%rdi,%rax,8),%rdx
  movslq %r15d,%rax
  lea    0x18(,%rax,8),%r9
  nopl   0x0(%rax)
L536:
  test   %esi,%esi
  jle    L241
  pxor   %xmm2,%xmm2
  lea    0x0(%rbp,%r9,1),%rax
  lea    (%r11,%r9,1),%rcx
  xor    %edi,%edi
  cvtsi2sd %ebx,%xmm2
  nopw   0x0(%rax,%rax,1)
L48:
  movsd  0x8(%rdx),%xmm1
  movsd  (%rdx),%xmm0
  add    $0x1,%edi
  mulsd  -0x8(%rax),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rcx)
  movsd  0x8(%rdx),%xmm0
  movsd  (%rdx),%xmm1
  mulsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm0
  mulsd  -0x8(%rax),%xmm0
  add    %r8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rcx)
  add    %r8,%rcx
  cmp    %edi,%esi
  jne    L48
L241:
  add    $0x2,%r10d
  add    $0x10,%rdx
  add    $0x10,%r9
  cmp    %r10d,%r14d
  jg     L536
  mov    -0x68(%rsp),%eax
  lea    0x4(%r12,%rax,1),%r12d
L227:
  lea    0x1(%r13),%eax
  add    -0x70(%rsp),%r15d
  cmp    %r13d,-0x74(%rsp)
  je     L155
  mov    %eax,%r13d
  lea    0x2(%r12),%eax
  cmp    $0x3,%r14d
  jg     L151
L370:
  mov    %eax,%r12d
  jmp    L227
L107:
  movl   $0x1,(%rax)
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L510:
  movslq 0x4(%rsp),%rax
  mov    %rax,-0x38(%rsp)
  movslq %r14d,%rax
  mov    %rax,-0x50(%rsp)
  test   %r14d,%r14d
  jg     L217
L446:
  mov    $0x2,%eax
  sub    %r14d,%eax
  mov    %eax,(%rsp)
  jmp    L386
L213:
  test   %ecx,%ecx
  jg     L225
  jmp    L446
L488:
  mov    $0x2,%eax
  sub    %r14d,%eax
  mov    %eax,(%rsp)
  mov    -0x44(%rsp),%eax
  sub    $0x2,%eax
  mov    %eax,-0x28(%rsp)
  jmp    L559
L549:
  mov    $0x2,%eax
  sub    %r14d,%eax
  mov    %eax,(%rsp)
  movslq %r14d,%rax
  mov    %rax,-0x50(%rsp)
  mov    -0x44(%rsp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x40(%rsp)
  jmp    L559
```

assistant: ```c
void passf(int *nac, int ido, int ip, int l1, int idl1, const double cc[], double c1[], double c2[], double ch[], double ch2[], const double wa[], int isign) {
  int idij, idlj, idot, ipph, i, j, k, l, jc, lc, ik, nt, idj, idl, inc, idp;
  double wai, war;
  idot = ido / 2;
  nt = ip * idl1;
  ipph = (ip + 1) / 2;
  idp = ip * ido;
  if (ido >= l1) {
    for (j = 1; j < ipph; j++) {
      jc = ip - j;
      for (k = 0; k < l1; k++) {
        for (i = 0; i < ido; i++) {
          ch[i + (k + j * l1) * ido] = cc[i + (j + k * ip) * ido] + cc[i + (jc + k * ip) * ido];
          ch[i + (k + jc * l1) * ido] = cc[i + (j + k * ip) * ido] - cc[i + (jc + k * ip) * ido];
        }
      }
    }
    for (k = 0; k < l1; k++)
      for (i = 0; i < ido; i++)
        ch[i + k * ido] = cc[i + k * ip * ido];
  } else {
    for (j = 1; j < ipph; j++) {
      jc = ip - j;
      for (i = 0; i < ido; i++) {
        for (k = 0; k < l1; k++) {
          ch[i + (k + j * l1) * ido] = cc[i + (j + k * ip) * ido] + cc[i + (jc + k * ip) * ido];
          ch[i + (k + jc * l1) * ido] = cc[i + (j + k * ip) * ido] - cc[i + (jc + k * ip) * ido];
        }
      }
    }
    for (i = 0; i < ido; i++)
      for (k = 0; k < l1; k++)
        ch[i + k * ido] = cc[i + k * ip * ido];
  }
  idl = 2 - ido;
  inc = 0;
  for (l = 1; l < ipph; l++) {
    lc = ip - l;
    idl += ido;
    for (ik = 0; ik < idl1; ik++) {
      c2[ik + l * idl1] = ch2[ik] + wa[idl - 2] * ch2[ik + idl1];
      c2[ik + lc * idl1] = isign * wa[idl - 1] * ch2[ik + (ip - 1) * idl1];
    }
    idlj = idl;
    inc += ido;
    for (j = 2; j < ipph; j++) {
      jc = ip - j;
      idlj += inc;
      if (idlj > idp)
        idlj -= idp;
      war = wa[idlj - 2];
      wai = wa[idlj - 1];
      for (ik = 0; ik < idl1; ik++) {
        c2[ik + l * idl1] += war * ch2[ik + j * idl1];
        c2[ik + lc * idl1] += isign * wai * ch2[ik + jc * idl1];
      }
    }
  }
  for (j = 1; j < ipph; j++)
    for (ik = 0; ik < idl1; ik++)
      ch2[ik] += ch2[ik + j * idl1];
  for (j = 1; j < ipph; j++) {
    jc = ip - j;
    for (ik = 1; ik < idl1; ik += 2) {
      ch2[ik - 1 + j * idl1] = c2[ik - 1 + j * idl1] - c2[ik + jc * idl1];
      ch2[ik - 1 + jc * idl1] = c2[ik - 1 + j * idl1] + c2[ik + jc * idl1];
      ch2[ik + j * idl1] = c2[ik + j * idl1] + c2[ik - 1 + jc * idl1];
      ch2[ik + jc * idl1] = c2[ik + j * idl1] - c2[ik - 1 + jc * idl1];
    }
  }
  *nac = 1;
  if (ido == 2)
    return;
  *nac = 0;
  for (ik = 0; ik < idl1; ik++)
    c2[ik] = ch2[ik];
  for (j = 1; j < ip; j++) {
    for (k = 0; k < l1; k++) {
      c1[(k + j * l1) * ido + 0] = ch[(k + j * l1) * ido + 0];
      c1[(k + j * l1) * ido + 1] = ch[(k + j * l1) * ido + 1];
    }
  }
  if (idot <= l1) {
    idij = 0;
    for (j = 1; j < ip; j++) {
      idij += 2;
      for (i = 3; i < ido; i += 2) {
        idij += 2;
        for (k = 0; k < l1; k++) {
          c1[i - 1 + (k + j * l1) * ido] = wa[idij - 2] * ch[i - 1 + (k + j * l1) * ido] - isign * wa[idij - 1] * ch[i + (k + j * l1) * ido];
          c1[i + (k + j * l1) * ido] = wa[idij - 2] * ch[i + (k + j * l1) * ido] + isign * wa[idij - 1] * ch[i - 1 + (k + j * l1) * ido];
        }
      }
    }
  } else {
    idj = 2 - ido;
    for (j = 1; j < ip; j++) {
      idj += ido;
      for (k = 0; k < l1; k++) {
        idij = idj;
        for (i = 3; i < ido; i += 2) {
          idij += 2;
          c1[i - 1 + (k + j * l1) * ido] = wa[idij - 2] * ch[i - 1 + (k + j * l1) * ido] - isign * wa[idij - 1] * ch[i + (k + j * l1) * ido];
          c1[i + (k + j * l1) * ido] = wa[idij - 2] * ch[i + (k + j * l1) * ido] + isign * wa[idij - 1] * ch[i - 1 + (k + j * l1) * ido];
        }
      }
    }
  }
}
```