user: Please show me the C implementation that matches the assembly code listed below:

```asm
<frlufact>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %ecx,%r13d
  push   %r12
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x38,%rsp
  cmp    %rdx,%rdi
  je     L46
  mov    %rdx,%rsi
  mov    %ecx,%edx
  imul   %ecx,%edx
  shl    $0x2,%edx
  call   <memcpy@plt>
L46:
  mov    %r13d,%eax
  sub    $0x1,%eax
  mov    %eax,0x28(%rsp)
  je     L117
  lea    0x1(%r13),%eax
  mov    %rbp,0x20(%rsp)
  mov    %r13d,%r14d
  mov    %eax,0x2c(%rsp)
  movl   $0x0,0x14(%rsp)
  movl   $0x0,0x1c(%rsp)
  movl   $0x0,0x18(%rsp)
  movl   $0x0,0x10(%rsp)
L62:
  cmp    0x10(%rsp),%r13d
  jbe    L37
  mov    0x10(%rsp),%r12d
  mov    0x14(%rsp),%ebp
  pxor   %xmm1,%xmm1
  mov    %r12d,%r15d
  nopl   0x0(%rax,%rax,1)
L18:
  mov    %ebp,%eax
  movss  %xmm1,0xc(%rsp)
  movss  (%rbx,%rax,4),%xmm0
  call   <fabs@plt>
  movss  0xc(%rsp),%xmm1
  comiss %xmm1,%xmm0
  maxss  %xmm1,%xmm0
  cmova  %r12d,%r15d
  add    $0x1,%r12d
  movaps %xmm0,%xmm1
  add    %r13d,%ebp
  cmp    %r12d,%r13d
  jne    L18
  mov    0x20(%rsp),%rax
  mov    %r15d,%edx
  imul   %r13d,%edx
  sub    0x18(%rsp),%edx
  mov    %r15d,(%rax)
  mov    0x14(%rsp),%eax
  nopl   0x0(%rax)
L21:
  mov    %eax,%ecx
  lea    (%rbx,%rcx,4),%rsi
  lea    (%rdx,%rax,1),%ecx
  add    $0x1,%eax
  lea    (%rbx,%rcx,4),%rcx
  movss  (%rsi),%xmm0
  movss  (%rcx),%xmm1
  movss  %xmm1,(%rsi)
  movss  %xmm0,(%rcx)
  cmp    %eax,%r14d
  jne    L21
L39:
  mov    0x1c(%rsp),%eax
  pxor   %xmm2,%xmm2
  addl   $0x1,0x10(%rsp)
  movss  (%rbx,%rax,4),%xmm0
  ucomiss %xmm2,%xmm0
  jp     L135
  lea    0x0(%r13,%r14,1),%r15d
  jne    L135
L96:
  mov    0x2c(%rsp),%edi
  add    %r13d,0x18(%rsp)
  mov    %r15d,%r14d
  add    %edi,0x1c(%rsp)
  add    %edi,0x14(%rsp)
  mov    0x10(%rsp),%edi
  addq   $0x4,0x20(%rsp)
  cmp    %edi,0x28(%rsp)
  jne    L62
L117:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L135:
  movss  D4(%rip),%xmm1
  lea    0x0(%r13,%r14,1),%r15d
  divss  %xmm0,%xmm1
  cmp    0x10(%rsp),%r13d
  jbe    L96
  mov    0x14(%rsp),%eax
  mov    0x10(%rsp),%r9d
  lea    0x0(%r13,%rax,1),%r10d
  mov    %r10d,%eax
  nopl   0x0(%rax)
L57:
  mov    %eax,%edx
  mov    %r9d,%r11d
  add    $0x1,%r9d
  add    %r13d,%eax
  lea    (%rbx,%rdx,4),%rdx
  movss  (%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,(%rdx)
  cmp    %r9d,%r13d
  jne    L57
  mov    0x18(%rsp),%edi
  lea    0x0(%r13,%r14,1),%r15d
  mov    %r15d,%esi
  sub    %r14d,%edi
  mov    0x10(%rsp),%r14d
  nopl   (%rax)
L59:
  mov    %r10d,%eax
  lea    (%rbx,%rax,4),%r8
  lea    0x1(%r10),%eax
  nopl   0x0(%rax,%rax,1)
L29:
  lea    (%rax,%rdi,1),%ecx
  mov    %eax,%edx
  add    $0x1,%eax
  movss  (%rbx,%rcx,4),%xmm1
  mulss  (%r8),%xmm1
  lea    (%rbx,%rdx,4),%rdx
  movss  (%rdx),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rdx)
  cmp    %eax,%esi
  jne    L29
  lea    0x1(%r14),%eax
  add    %r9d,%r10d
  sub    %r9d,%edi
  add    %r9d,%esi
  cmp    %r14d,%r11d
  je     L96
  mov    %eax,%r14d
  jmp    L59
L37:
  mov    0x20(%rsp),%rax
  mov    0x10(%rsp),%edi
  mov    %edi,(%rax)
  jmp    L39
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D4
type_guessed: byte[4]

tool: {"D4": [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];
    }
  }
}
```