user: What is the c source code of the assembly code below?

```asm
<svd_topk>:
  endbr64
  push   %r15
  push   %r14
  lea    (%rcx,%rcx,1),%r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,0x10(%rsp)
  mov    %rcx,0x18(%rsp)
  mov    %r8,0x8(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x48(%rsp)
  xor    %eax,%eax
  cmp    %rsi,%rdx
  mov    %rsi,%rax
  movq   $0x0,0x30(%rsp)
  movq   $0x0,0x38(%rsp)
  cmovbe %rdx,%rax
  movq   $0x0,0x40(%rsp)
  cmp    %r14,%rax
  ja     L91
  xor    %r13d,%r13d
  xor    %r15d,%r15d
  xor    %ebx,%ebx
L21:
  mov    %r13,%rdi
  xor    %r13d,%r13d
  call   <mem_free@plt>
L49:
  mov    0x8(%rsp),%rax
  mov    0x38(%rsp),%rdi
  mov    %r13,(%rax)
  call   <mem_free@plt>
  mov    0x30(%rsp),%rdi
  call   <mem_free@plt>
  mov    0x40(%rsp),%rdi
  call   <mem_free@plt>
  mov    %rbx,%rdi
  call   <mem_free@plt>
  mov    %r15,%rdi
  call   <mem_free@plt>
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L9
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L91:
  mov    %r14,%rax
  mov    %rdx,%rbp
  mov    %rsi,%r12
  mov    $0x4,%esi
  imul   %rdx,%rax
  mov    %r12,%r13
  mov    %rax,%rdi
  mov    %rax,0x20(%rsp)
  call   <mem_alloc@plt>
  cmp    %r12,%rbp
  mov    $0x4,%esi
  cmovae %rbp,%r13
  mov    %rax,%rbx
  imul   %r14,%r13
  mov    %r13,%rdi
  call   <mem_alloc@plt>
  mov    %r13,%rdi
  mov    $0x4,%esi
  mov    %rax,%r15
  call   <mem_alloc@plt>
  mov    %rax,%r13
  test   %rbx,%rbx
  je     L21
  test   %r15,%r15
  je     L21
  test   %rax,%rax
  je     L21
  mov    0x20(%rsp),%rax
  test   %rax,%rax
  je     L22
  xor    %edx,%edx
  mov    %r15,0x20(%rsp)
  mov    %rax,%r15
  mov    %rbp,0x28(%rsp)
  mov    %rdx,%rbp
  nopw   0x0(%rax,%rax,1)
L54:
  xor    %eax,%eax
  call   <rnorm@plt>
  movss  %xmm0,(%rbx,%rbp,4)
  add    $0x1,%rbp
  cmp    %rbp,%r15
  jne    L54
  mov    0x20(%rsp),%r15
  mov    0x28(%rsp),%rbp
L22:
  mov    0x10(%rsp),%rdi
  mov    %rbp,%rcx
  mov    %r12,%rdx
  mov    %rbx,%rsi
  mov    %r13,%r9
  mov    %r14,%r8
  call   <dotnn@plt>
  mov    %r12,%rsi
  mov    %r15,%rcx
  mov    %r14,%rdx
  mov    %r13,%rdi
  call   <qr@plt>
  mov    %rbp,%rax
  mov    %rbx,0x28(%rsp)
  mov    %r12,%rbp
  movq   $0x2,0x20(%rsp)
  mov    0x10(%rsp),%rbx
  mov    %rax,%r12
L90:
  mov    %r13,%r9
  mov    %r14,%r8
  mov    %r12,%rcx
  mov    %rbp,%rdx
  mov    %r15,%rsi
  mov    %rbx,%rdi
  call   <dottn@plt>
  mov    %r15,%rcx
  mov    %r14,%rdx
  mov    %r12,%rsi
  mov    %r13,%rdi
  call   <qr@plt>
  mov    %r13,%r9
  mov    %r14,%r8
  mov    %r12,%rcx
  mov    %rbp,%rdx
  mov    %r15,%rsi
  mov    %rbx,%rdi
  call   <dotnn@plt>
  mov    %r15,%rcx
  mov    %r14,%rdx
  mov    %rbp,%rsi
  mov    %r13,%rdi
  call   <qr@plt>
  cmpq   $0x1,0x20(%rsp)
  jne    L61
  mov    0x28(%rsp),%rbx
  mov    %r12,%rax
  mov    0x10(%rsp),%rsi
  mov    %rbp,%r12
  mov    %r14,%rcx
  mov    %rax,%rbp
  mov    %rax,%r8
  mov    %r12,%rdx
  mov    %rbx,%r9
  mov    %r15,%rdi
  call   <dottn@plt>
  mov    %rbp,%rdx
  mov    %r14,%rsi
  mov    %rbx,%rdi
  call   <transpose@plt>
  lea    0x38(%rsp),%rcx
  mov    %r14,%rdx
  mov    %rbp,%rsi
  lea    0x40(%rsp),%r9
  lea    0x30(%rsp),%r8
  mov    %rbx,%rdi
  call   <svd_full@plt>
  mov    0x18(%rsp),%rbp
  mov    %rbx,%r9
  mov    %r14,%rcx
  mov    0x30(%rsp),%rsi
  mov    0x40(%rsp),%rdi
  mov    %r14,%rdx
  mov    %rbp,%r8
  call   <trim@plt>
  mov    %r13,%r9
  mov    %rbp,%r8
  mov    %r14,%rcx
  mov    %r12,%rdx
  mov    %rbx,%rsi
  mov    %r15,%rdi
  call   <dotnn@plt>
  jmp    L49
  nop
L61:
  movq   $0x1,0x20(%rsp)
  jmp    L90
L9:
  call   <__stack_chk_fail@plt>
```

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