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

```asm
<ut>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %esi,-0x4c(%rbp)
  mov    %edx,-0x50(%rbp)
  mov    -0x50(%rbp),%eax
  cltq
  shl    $0x2,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0x18(%rbp)
  movq   $0x0,-0x30(%rbp)
  jmp    L31
L20:
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x30(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x30(%rbp),%rax
  add    $0x1,%rax
  mov    %rax,-0x28(%rbp)
  jmp    L79
L106:
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x28(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x50(%rbp),%edx
  mov    -0x8(%rbp),%rcx
  mov    -0x10(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <rowcmp@plt>
  test   %rax,%rax
  je     L148
  mov    -0x50(%rbp),%edx
  mov    -0x10(%rbp),%rcx
  mov    -0x18(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
  mov    -0x50(%rbp),%edx
  mov    -0x8(%rbp),%rcx
  mov    -0x10(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
  mov    -0x50(%rbp),%edx
  mov    -0x18(%rbp),%rcx
  mov    -0x8(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
L148:
  addq   $0x1,-0x28(%rbp)
L79:
  mov    -0x4c(%rbp),%eax
  cltq
  cmp    %rax,-0x28(%rbp)
  jl     L106
  mov    -0x30(%rbp),%rax
  add    $0x1,%rax
  mov    %rax,-0x28(%rbp)
  jmp    L67
L27:
  movq   $0x0,-0x20(%rbp)
  jmp    L115
L144:
  addq   $0x1,-0x20(%rbp)
L115:
  mov    -0x50(%rbp),%eax
  cltq
  cmp    %rax,-0x20(%rbp)
  jge    L143
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x30(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  ucomiss %xmm1,%xmm0
  jp     L143
  pxor   %xmm1,%xmm1
  ucomiss %xmm1,%xmm0
  je     L144
L143:
  mov    -0x50(%rbp),%eax
  cltq
  cmp    %rax,-0x20(%rbp)
  je     L87
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x28(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x30(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x34(%rbp)
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x28(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  jmp    L5
L132:
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x28(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x30(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mulss  -0x34(%rbp),%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  imul   -0x28(%rbp),%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
L5:
  addq   $0x1,-0x20(%rbp)
  mov    -0x50(%rbp),%eax
  cltq
  cmp    %rax,-0x20(%rbp)
  jl     L132
  addq   $0x1,-0x28(%rbp)
L67:
  mov    -0x4c(%rbp),%eax
  cltq
  cmp    %rax,-0x28(%rbp)
  jl     L27
  jmp    L71
L87:
  nop
L71:
  addq   $0x1,-0x30(%rbp)
L31:
  mov    -0x30(%rbp),%rax
  lea    0x1(%rax),%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  cmp    %rax,%rdx
  jl     L20
  mov    -0x18(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  nop
  leave
  ret
```

assistant: ```c
void ut(float *mat, int m, int n) {
  long i, j, k;
  float *current, *other, *tmp = malloc(n * sizeof(float));
  for (i = 0; i + 1 < m; ++i) {
    current = mat + i * n;
    for (j = i + 1; j < m; ++j) {
      other = mat + j * n;
      if (rowcmp(current, other, n)) {
        memcpy(tmp, current, n);
        memcpy(current, other, n);
        memcpy(other, tmp, n);
      }
    }
    for (j = i + 1; j < m; ++j) {
      for (k = 0; k < n && !mat[i * n + k]; ++k)
        ;
      if (k == n)
        break;
      float c = mat[j * n + k] / mat[i * n + k];
      mat[j * n + k] = 0.0;
      while (++k < n)
        mat[j * n + k] -= c * mat[i * n + k];
    }
  }
  free(tmp);
}
```