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

```asm
<u2f>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  push   %r14
  mov    %rdi,%r14
  push   %r13
  mov    %r8,%r13
  push   %r12
  mov    %r8,%r12
  push   %rbp
  mov    %rcx,%rbp
  push   %rbx
  mov    %edx,%ebx
  sub    $0x538,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x528(%rsp)
  xor    %eax,%eax
  test   %r8,%r8
  je     L217
L111:
  cmp    $0x8,%ebx
  jle    L132
  lea    0x30(%rsp),%r9
  lea    0x10(%rsp),%rcx
  mov    %r14,%rdi
  mov    %ebx,%edx
  mov    %r9,%r8
  mov    %r15,%rsi
  mov    %r9,0x8(%rsp)
  mov    %rcx,(%rsp)
  call   <normu@plt>
  mov    0x8(%rsp),%r9
  mov    (%rsp),%r8
  mov    %ebx,%ecx
  mov    %r14,%rdi
  mov    %r15,%rdx
  mov    %r13,%rsi
  call   <lin_fmN@plt>
  mov    $0x9,%ecx
  mov    %ebx,%edx
  mov    %r13,%rsi
  lea    0x2a0(%rsp),%r14
  mov    %r14,%rdi
  call   <cov_mat@plt>
  lea    0x50(%rsp),%rsi
  mov    $0x9,%edx
  mov    %r14,%rdi
  call   <lap_eig@plt>
  mov    $0x9,%esi
  mov    %r14,%rdi
  call   <trnm@plt>
  movsd  0x50(%rsp),%xmm0
  movsd  0x58(%rsp),%xmm1
  mov    $0x1,%eax
  comisd %xmm1,%xmm0
  jbe    L86
  movsd  0x60(%rsp),%xmm0
  comisd %xmm0,%xmm1
  ja     L94
L143:
  movapd %xmm1,%xmm0
  movsd  0x68(%rsp),%xmm1
  comisd %xmm1,%xmm0
  ja     L210
L98:
  movapd %xmm0,%xmm1
  movsd  0x70(%rsp),%xmm0
  comisd %xmm0,%xmm1
  ja     L65
L158:
  movapd %xmm1,%xmm0
  movsd  0x78(%rsp),%xmm1
  comisd %xmm1,%xmm0
  ja     L103
L37:
  movapd %xmm0,%xmm1
  movsd  0x80(%rsp),%xmm0
  comisd %xmm0,%xmm1
  ja     L151
L12:
  movapd %xmm1,%xmm0
  movsd  0x88(%rsp),%xmm1
  comisd %xmm1,%xmm0
  ja     L87
L205:
  movapd %xmm0,%xmm1
L23:
  ucomisd 0x90(%rsp),%xmm1
  mov    $0x8,%edx
  cmova  %rdx,%rax
  lea    (%r14,%rax,8),%rax
  jmp    L150
  nopl   0x0(%rax)
L132:
  mov    %ebx,%ecx
  mov    %r15,%rdx
  mov    %r13,%rsi
  mov    %r14,%rdi
  call   <lin_fm@plt>
  lea    0x2a0(%rsp),%rdx
  mov    %r13,%rsi
  lea    0x50(%rsp),%rdi
  mov    $0x8,%r9d
  lea    0xa0(%rsp),%r8
  mov    $0x9,%ecx
  call   <svduv@plt>
  lea    0x2e0(%rsp),%rax
L150:
  movupd 0x40(%rax),%xmm2
  movupd 0xd0(%rax),%xmm3
  mov    %rbp,%rdi
  movupd 0x160(%rax),%xmm4
  movupd 0x1f0(%rax),%xmm5
  movsd  0x240(%rax),%xmm0
  movlpd (%rax),%xmm2
  movlpd 0x90(%rax),%xmm3
  movlpd 0x120(%rax),%xmm4
  movlpd 0x1b0(%rax),%xmm5
  movups %xmm2,0x0(%rbp)
  movups %xmm3,0x10(%rbp)
  movups %xmm4,0x20(%rbp)
  movups %xmm5,0x30(%rbp)
  movsd  %xmm0,0x40(%rbp)
  call   <singulF@plt>
  cmp    $0x8,%ebx
  jg     L203
  test   %r12,%r12
  je     L176
L221:
  mov    0x528(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L42
  add    $0x538,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L94:
  movsd  0x68(%rsp),%xmm1
  mov    $0x2,%eax
  comisd %xmm1,%xmm0
  jbe    L98
L210:
  movsd  0x70(%rsp),%xmm0
  mov    $0x3,%eax
  comisd %xmm0,%xmm1
  jbe    L158
L65:
  movsd  0x78(%rsp),%xmm1
  mov    $0x4,%eax
  comisd %xmm1,%xmm0
  jbe    L37
L103:
  movsd  0x80(%rsp),%xmm0
  mov    $0x5,%eax
  comisd %xmm0,%xmm1
  jbe    L12
L151:
  movsd  0x88(%rsp),%xmm1
  mov    $0x6,%eax
  comisd %xmm1,%xmm0
  jbe    L205
L87:
  mov    $0x7,%eax
  jmp    L23
  nopl   0x0(%rax)
L86:
  movapd %xmm0,%xmm1
  movsd  0x60(%rsp),%xmm0
  xor    %eax,%eax
  comisd %xmm0,%xmm1
  ja     L94
  jmp    L143
  nopl   0x0(%rax,%rax,1)
L203:
  lea    0x30(%rsp),%rdx
  lea    0x10(%rsp),%rsi
  mov    %rbp,%rdi
  call   <denormF@plt>
  test   %r12,%r12
  jne    L221
L176:
  mov    %r13,%rdi
  call   <free@plt>
  jmp    L221
L217:
  lea    (%rdx,%rdx,8),%edi
  shl    $0x3,%edi
  call   <malloc@plt>
  mov    %rax,%r13
  jmp    L111
L42:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void u2f(const double *u, const int *inl, int len, double *F, double *buffer) {
  double A1[3], A2[3];
  double *Z, V[9 * 9], U[8 * 8], D[9], *p;
  int i, j;
  if (buffer == NULL)
    Z = (double *)malloc(sizeof(double) * 9 * len);
  else
    Z = buffer;
  if (len > 8) {
    normu(u, inl, len, A1, A2);
    lin_fmN(u, Z, inl, len, A1, A2);
    cov_mat(V, Z, len, 9);
    lap_eig(V, D, 9);
    trnm(V, 9);
  } else {
    lin_fm(u, Z, inl, len);
    svduv(D, Z, V, 9, U, 8);
  }
  if (len > 8) {
    j = 0;
    for (i = 1; i < 9; i++)
      if (D[i] < D[j])
        j = i;
    p = V + j;
  } else
    p = V + 8;
  for (i = 0; i < 9; i++) {
    F[i] = *p;
    p += 9;
  }
  singulF(F);
  if (len > 8)
    denormF(F, A1, A2);
  if (buffer == NULL)
    free(Z);
}
```