user: 请将以下的汇编代码转换成C语言代码:

```asm
<frlufact>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,%rbx
  mov    %rsi,0x28(%rsp)
  mov    %ecx,%r14d
  cmp    %rdx,%rdi
  je     L28
  mov    %rdx,%rsi
  mov    %ecx,%edx
  imul   %ecx,%edx
  shl    $0x2,%edx
  call   <memcpy@plt>
L28:
  cmp    $0x1,%r14d
  je     L60
  lea    0x1(%r14),%eax
  mov    %eax,0x1c(%rsp)
  lea    -0x1(%r14),%eax
  mov    %rax,0x20(%rsp)
  mov    %r14d,%r15d
  movl   $0x0,0x14(%rsp)
  movl   $0x0,0x10(%rsp)
  movq   $0x0,0x8(%rsp)
  jmp    L1
L2:
  addq   $0x1,0x8(%rsp)
  mov    0x8(%rsp),%rax
  mov    0x1c(%rsp),%esi
  add    %esi,0x10(%rsp)
  add    %r14d,0x14(%rsp)
  add    %r14d,%r15d
  cmp    0x20(%rsp),%rax
  je     L60
L1:
  mov    0x8(%rsp),%rax
  mov    %eax,0x18(%rsp)
  cmp    %eax,%r14d
  jbe    L50
  mov    0x10(%rsp),%r13d
  mov    0x8(%rsp),%rax
  mov    %eax,%r12d
  mov    %eax,%ebp
  pxor   %xmm4,%xmm4
  movss  %xmm4,0x4(%rsp)
L75:
  mov    %r13d,%eax
  movss  (%rbx,%rax,4),%xmm0
  call   <fabs@plt>
  movss  0x4(%rsp),%xmm2
  comiss %xmm2,%xmm0
  maxss  %xmm2,%xmm0
  movss  %xmm0,0x4(%rsp)
  cmova  %ebp,%r12d
  add    $0x1,%ebp
  add    %r14d,%r13d
  cmp    %ebp,%r14d
  jne    L75
  mov    0x28(%rsp),%rax
  mov    0x8(%rsp),%rsi
  mov    %r12d,(%rax,%rsi,4)
  imul   %ebp,%r12d
  mov    0x10(%rsp),%eax
  sub    0x14(%rsp),%r12d
L37:
  mov    %eax,%edx
  lea    (%rbx,%rdx,4),%rcx
  movss  (%rcx),%xmm0
  lea    (%r12,%rax,1),%edx
  lea    (%rbx,%rdx,4),%rdx
  movss  (%rdx),%xmm1
  movss  %xmm1,(%rcx)
  movss  %xmm0,(%rdx)
  add    $0x1,%eax
  cmp    %r15d,%eax
  jne    L37
L99:
  mov    0x10(%rsp),%eax
  movss  (%rbx,%rax,4),%xmm0
  pxor   %xmm3,%xmm3
  ucomiss %xmm3,%xmm0
  jp     L10
  je     L2
L10:
  movss  D1(%rip),%xmm1
  divss  %xmm0,%xmm1
  mov    0x18(%rsp),%ebp
  add    $0x1,%ebp
  cmp    %ebp,%r14d
  jbe    L2
  mov    0x10(%rsp),%eax
  lea    (%r14,%rax,1),%r10d
  mov    %r10d,%eax
  mov    %ebp,%r9d
L94:
  mov    %eax,%edx
  lea    (%rbx,%rdx,4),%rdx
  movaps %xmm1,%xmm0
  mulss  (%rdx),%xmm0
  movss  %xmm0,(%rdx)
  mov    %r9d,%r11d
  add    $0x1,%r9d
  add    %r14d,%eax
  cmp    %r9d,%r14d
  jne    L94
  mov    0x14(%rsp),%edi
  sub    %r15d,%edi
  lea    (%r9,%r15,1),%esi
L89:
  mov    %r10d,%eax
  lea    (%rbx,%rax,4),%r8
  lea    0x1(%r10),%eax
L134:
  mov    %eax,%edx
  lea    (%rbx,%rdx,4),%rdx
  lea    (%rax,%rdi,1),%ecx
  movss  (%rbx,%rcx,4),%xmm1
  mulss  (%r8),%xmm1
  movss  (%rdx),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rdx)
  add    $0x1,%eax
  cmp    %eax,%esi
  jne    L134
  lea    0x1(%rbp),%eax
  add    %r9d,%r10d
  sub    %r9d,%edi
  add    %r9d,%esi
  cmp    %ebp,%r11d
  je     L2
  mov    %eax,%ebp
  jmp    L89
L50:
  mov    0x28(%rsp),%rsi
  mov    0x8(%rsp),%rax
  mov    %eax,(%rsi,%rax,4)
  jmp    L99
L60:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D1
type_guessed: byte[4]

tool: {"D1": [0x00,0x00,0x80,0x3f]}

assistant: ```c
void frlufact(float *u, unsigned int *p, const float *a, unsigned int d) {
  unsigned int i, j, k, mu;
  float f1, f2;
  if (u != a)
    memcpy(u, a, d * d * sizeof(u[0]));
  for (k = 0; k < d - 1; k++) {
    for (f1 = 0, i = mu = k; i < d; i++) {
      f2 = fabs(u[i * d + k]);
      if (f2 > f1) {
        f1 = f2;
        mu = i;
      }
    }
    p[k] = mu;
    for (i = k; i < d; i++)
      do {
        float z;
        z = (u[k * d + i]);
        (u[k * d + i]) = (u[mu * d + i]);
        (u[mu * d + i]) = z;
      } while (0);
    f1 = u[k * d + k];
    if (f1 != 0) {
      f1 = 1 / f1;
      for (i = k + 1; i < d; i++)
        u[i * d + k] *= f1;
      for (i = k + 1; i < d; i++)
        for (j = k + 1; j < d; j++)
          u[i * d + j] -= u[i * d + k] * u[k * d + j];
    }
  }
}
```