user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<svd_topk>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,-0x68(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    %rdx,-0x78(%rbp)
  mov    %rcx,-0x80(%rbp)
  mov    %r8,-0x88(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  mov    -0x80(%rbp),%rax
  add    %rax,%rax
  mov    %rax,-0x28(%rbp)
  movq   $0x2,-0x20(%rbp)
  movq   $0x0,-0x40(%rbp)
  movq   $0x0,-0x38(%rbp)
  movq   $0x0,-0x60(%rbp)
  movq   $0x0,-0x58(%rbp)
  movq   $0x0,-0x50(%rbp)
  movq   $0x0,-0x30(%rbp)
  mov    -0x78(%rbp),%rdx
  mov    -0x70(%rbp),%rax
  cmp    %rax,%rdx
  cmovbe %rdx,%rax
  cmp    %rax,-0x28(%rbp)
  jae    L12
  mov    -0x28(%rbp),%rax
  imul   -0x78(%rbp),%rax
  mov    $0x4,%esi
  mov    %rax,%rdi
  call   <mem_alloc@plt>
  mov    %rax,-0x40(%rbp)
  mov    -0x78(%rbp),%rdx
  mov    -0x70(%rbp),%rax
  cmp    %rax,%rdx
  cmovae %rdx,%rax
  imul   -0x28(%rbp),%rax
  mov    $0x4,%esi
  mov    %rax,%rdi
  call   <mem_alloc@plt>
  mov    %rax,-0x38(%rbp)
  mov    -0x78(%rbp),%rdx
  mov    -0x70(%rbp),%rax
  cmp    %rax,%rdx
  cmovae %rdx,%rax
  imul   -0x28(%rbp),%rax
  mov    $0x4,%esi
  mov    %rax,%rdi
  call   <mem_alloc@plt>
  mov    %rax,-0x30(%rbp)
  cmpq   $0x0,-0x40(%rbp)
  je     L42
  cmpq   $0x0,-0x38(%rbp)
  je     L42
  cmpq   $0x0,-0x30(%rbp)
  je     L42
  movq   $0x0,-0x48(%rbp)
  jmp    L19
L50:
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x40(%rbp),%rax
  lea    (%rdx,%rax,1),%rbx
  mov    $0x0,%eax
  call   <rnorm@plt>
  movd   %xmm0,%eax
  mov    %eax,(%rbx)
  addq   $0x1,-0x48(%rbp)
L19:
  mov    -0x78(%rbp),%rax
  imul   -0x28(%rbp),%rax
  cmp    %rax,-0x48(%rbp)
  jb     L50
  mov    -0x30(%rbp),%r8
  mov    -0x28(%rbp),%rdi
  mov    -0x78(%rbp),%rcx
  mov    -0x70(%rbp),%rdx
  mov    -0x40(%rbp),%rsi
  mov    -0x68(%rbp),%rax
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <dotnn@plt>
  mov    -0x38(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    -0x70(%rbp),%rsi
  mov    -0x30(%rbp),%rax
  mov    %rax,%rdi
  call   <qr@plt>
  movq   $0x0,-0x48(%rbp)
  jmp    L41
L108:
  mov    -0x30(%rbp),%r8
  mov    -0x28(%rbp),%rdi
  mov    -0x78(%rbp),%rcx
  mov    -0x70(%rbp),%rdx
  mov    -0x38(%rbp),%rsi
  mov    -0x68(%rbp),%rax
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <dottn@plt>
  mov    -0x38(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    -0x78(%rbp),%rsi
  mov    -0x30(%rbp),%rax
  mov    %rax,%rdi
  call   <qr@plt>
  mov    -0x30(%rbp),%r8
  mov    -0x28(%rbp),%rdi
  mov    -0x78(%rbp),%rcx
  mov    -0x70(%rbp),%rdx
  mov    -0x38(%rbp),%rsi
  mov    -0x68(%rbp),%rax
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <dotnn@plt>
  mov    -0x38(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    -0x70(%rbp),%rsi
  mov    -0x30(%rbp),%rax
  mov    %rax,%rdi
  call   <qr@plt>
  addq   $0x1,-0x48(%rbp)
L41:
  mov    -0x48(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  jb     L108
  mov    -0x40(%rbp),%r8
  mov    -0x78(%rbp),%rdi
  mov    -0x28(%rbp),%rcx
  mov    -0x70(%rbp),%rdx
  mov    -0x68(%rbp),%rsi
  mov    -0x38(%rbp),%rax
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <dottn@plt>
  mov    -0x78(%rbp),%rdx
  mov    -0x28(%rbp),%rcx
  mov    -0x40(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <transpose@plt>
  lea    -0x50(%rbp),%r8
  lea    -0x60(%rbp),%rdi
  lea    -0x58(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    -0x78(%rbp),%rsi
  mov    -0x40(%rbp),%rax
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <svd_full@plt>
  mov    -0x60(%rbp),%rsi
  mov    -0x50(%rbp),%rax
  mov    -0x40(%rbp),%r8
  mov    -0x80(%rbp),%rdi
  mov    -0x28(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <trim@plt>
  mov    -0x30(%rbp),%r8
  mov    -0x80(%rbp),%rdi
  mov    -0x28(%rbp),%rcx
  mov    -0x70(%rbp),%rdx
  mov    -0x40(%rbp),%rsi
  mov    -0x38(%rbp),%rax
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <dotnn@plt>
  jmp    L90
L12:
  nop
  jmp    L77
L42:
  nop
L77:
  mov    -0x30(%rbp),%rax
  mov    %rax,%rdi
  call   <mem_free@plt>
  movq   $0x0,-0x30(%rbp)
L90:
  mov    -0x88(%rbp),%rax
  mov    -0x30(%rbp),%rdx
  mov    %rdx,(%rax)
  mov    -0x58(%rbp),%rax
  mov    %rax,%rdi
  call   <mem_free@plt>
  mov    -0x60(%rbp),%rax
  mov    %rax,%rdi
  call   <mem_free@plt>
  mov    -0x50(%rbp),%rax
  mov    %rax,%rdi
  call   <mem_free@plt>
  mov    -0x40(%rbp),%rax
  mov    %rax,%rdi
  call   <mem_free@plt>
  mov    -0x38(%rbp),%rax
  mov    %rax,%rdi
  call   <mem_free@plt>
  nop
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L99
  call   <__stack_chk_fail@plt>
L99:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
void svd_topk(const float *a, size_t m, size_t n, size_t k, float **ou) {
  const size_t l = 2 * k;
  const size_t p = 2;
  size_t i;
  float *b = NULL;
  float *q = NULL;
  float *s = NULL;
  float *u = NULL;
  float *v = NULL;
  float *y = NULL;
  if (l >= (((m) < (n)) ? (m) : (n)))
    goto error;
  b = mem_alloc(l * n, sizeof(float));
  q = mem_alloc(l * (((m) > (n)) ? (m) : (n)), sizeof(float));
  y = mem_alloc(l * (((m) > (n)) ? (m) : (n)), sizeof(float));
  if (b == NULL || q == NULL || y == NULL)
    goto error;
  for (i = 0; i < n * l; i++)
    b[i] = rnorm();
  dotnn(a, b, m, n, l, y);
  qr(y, m, l, q);
  for (i = 0; i < p; i++) {
    dottn(a, q, m, n, l, y);
    qr(y, n, l, q);
    dotnn(a, q, m, n, l, y);
    qr(y, m, l, q);
  }
  dottn(q, a, m, l, n, b);
  transpose(b, l, n);
  svd_full(b, n, l, &u, &s, &v);
  trim(v, s, l, l, k, b);
  dotnn(q, b, m, l, k, y);
  goto done;
error:
  do {
    mem_free(y);
    y = NULL;
  } while (0);
done:
  *ou = y;
  mem_free(u);
  mem_free(s);
  mem_free(v);
  mem_free(b);
  mem_free(q);
}
```