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

```asm
<ludecomp>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%r8
  mov    %rsi,%r10
  mov    %rdx,-0x8(%rsp)
  mov    %ecx,%edi
  cmp    %rsi,%r8
  je     L26
  mov    %ecx,%edx
  imul   %ecx,%edx
  test   %edx,%edx
  jle    L26
  mov    %edx,%edx
  mov    $0x0,%eax
L129:
  movss  (%r8,%rax,4),%xmm0
  movss  %xmm0,(%r10,%rax,4)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L129
L26:
  test   %edi,%edi
  jle    L162
  lea    -0x1(%rdi),%eax
  mov    %rax,-0x18(%rsp)
  mov    %edi,%edx
  mov    $0x0,%eax
  mov    -0x8(%rsp),%rcx
L156:
  mov    %eax,(%rcx,%rax,4)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L156
  lea    0x1(%rdi),%eax
  mov    %eax,-0xc(%rsp)
  movslq %edi,%rcx
  shl    $0x2,%rcx
  mov    %r10,%rbx
  movl   $0xffffffff,-0x1c(%rsp)
  movl   $0x0,-0x10(%rsp)
  movl   $0x0,-0x20(%rsp)
  mov    $0x0,%r11d
  lea    0x4(%r10),%rbp
  pxor   %xmm3,%xmm3
  movss  D0(%rip),%xmm2
  jmp    L231
L68:
  mov    %r12d,%r9d
  jmp    L191
L221:
  movslq %r13d,%rdx
  lea    (%r10,%rdx,4),%rax
  add    %r14,%rdx
  lea    0x0(%rbp,%rdx,4),%rsi
  mov    %rbx,%rdx
  pxor   %xmm1,%xmm1
L79:
  movss  (%rdx),%xmm0
  mulss  (%rax),%xmm0
  addss  %xmm0,%xmm1
  add    %rcx,%rdx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L79
L6:
  movss  (%r8),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%r8)
  add    $0x1,%r9d
  add    %rcx,%r8
  add    %edi,%r13d
  cmp    %r9d,%edi
  je     L84
L116:
  pxor   %xmm1,%xmm1
  test   %r12d,%r12d
  jg     L221
  jmp    L6
L84:
  lea    0x1(%r15),%r9d
  cmp    -0x18(%rsp),%r11
  je     L162
  movslq -0x10(%rsp),%rax
  lea    (%r10,%rax,4),%r14
  movslq -0x20(%rsp),%rax
  add    $0x1,%r11
  lea    (%rax,%r11,1),%rdx
  lea    (%r10,%rdx,4),%r8
  mov    %rax,%r13
  neg    %r13
  shl    $0x2,%r13
  lea    (%r10,%rax,4),%r15
  mov    -0x1c(%rsp),%edx
  add    %rax,%rdx
  lea    0x0(%rbp,%rdx,4),%rsi
  jmp    L14
L139:
  lea    (%r8,%r13,1),%rdx
  mov    %r15,%rax
  pxor   %xmm1,%xmm1
L187:
  movss  (%rdx),%xmm0
  mulss  (%rax),%xmm0
  addss  %xmm0,%xmm1
  add    %rcx,%rdx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L187
L70:
  movss  (%r8),%xmm0
  subss  %xmm1,%xmm0
  divss  (%r14),%xmm0
  movss  %xmm0,(%r8)
  add    $0x1,%r9d
  add    $0x4,%r8
  cmp    %r9d,%edi
  je     L157
L14:
  pxor   %xmm1,%xmm1
  test   %r12d,%r12d
  jg     L139
  jmp    L70
L157:
  add    %edi,-0x20(%rsp)
  mov    -0xc(%rsp),%esi
  add    %esi,-0x10(%rsp)
  add    $0x4,%rbx
  addl   $0x1,-0x1c(%rsp)
L231:
  mov    %r11d,%r15d
  mov    %r11d,%r9d
  mov    %r11d,%r12d
  cmp    %r11d,%edi
  jle    L1
  movslq -0x20(%rsp),%rax
  lea    (%r10,%rax,4),%rsi
  mov    %r11,%rax
  mov    $0xffffffff,%edx
  movaps %xmm3,%xmm1
L238:
  movss  (%rsi,%rax,4),%xmm0
  andps  %xmm2,%xmm0
  comiss %xmm1,%xmm0
  cmova  %eax,%edx
  maxss  %xmm1,%xmm0
  movaps %xmm0,%xmm1
  add    $0x1,%rax
  cmp    %eax,%edi
  jg     L238
  cmp    $0xffffffff,%edx
  je     L1
  cmp    %r12d,%edx
  je     L68
  mov    %rbx,%rax
  mov    $0x0,%esi
  movslq %edx,%rdx
  mov    %rdx,%r8
  sub    %r11,%r8
L55:
  movss  (%rax),%xmm0
  movss  (%rax,%r8,4),%xmm1
  movss  %xmm1,(%rax)
  movss  %xmm0,(%rax,%r8,4)
  add    $0x1,%esi
  add    %rcx,%rax
  cmp    %esi,%edi
  jne    L55
  mov    -0x8(%rsp),%r14
  mov    (%r14,%r11,4),%eax
  lea    (%r14,%rdx,4),%rdx
  mov    (%rdx),%esi
  mov    %esi,(%r14,%r11,4)
  mov    %eax,(%rdx)
L1:
  cmp    %r9d,%edi
  jle    L84
L191:
  mov    %edi,%eax
  imul   %r9d,%eax
  cltq
  add    %r11,%rax
  lea    (%r10,%rax,4),%r8
  mov    %r9d,%r13d
  imul   %edi,%r13d
  mov    -0x1c(%rsp),%r14d
  jmp    L116
L162:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: label: D0
type_infered: dword

tool: D0: 2147483647

assistant: ```c
int ludecomp(float A[], float LU[], int P[], int N) {
  int i, j, k;
  float max, tmp;
  int who, tmpi;
  float sum;
  if (A != LU) {
    for (i = 0; i < N * N; i++) {
      LU[i] = A[i];
    }
  }
  for (i = 0; i < N; i++) {
    P[i] = i;
  }
  for (i = 0; i < N; i++) {
    max = 0;
    who = -1;
    for (j = i; j < N; j++) {
      tmp = fabsf(LU[(i) * (N) + (j)]);
      if (tmp > max) {
        max = tmp;
        who = j;
      }
    }
    if (who == -1) {
      who = i;
    }
    if (who != i) {
      for (j = 0; j < N; j++) {
        tmp = LU[(j) * (N) + (i)];
        LU[(j) * (N) + (i)] = LU[(j) * (N) + (who)];
        LU[(j) * (N) + (who)] = tmp;
      }
      tmpi = P[i];
      P[i] = P[who];
      P[who] = tmpi;
    }
    for (j = i; j < N; j++) {
      sum = 0;
      for (k = 0; k < i; k++) {
        sum += LU[(k) * (N) + (i)] * LU[(j) * (N) + (k)];
      }
      LU[(j) * (N) + (i)] -= sum;
    }
    for (j = i + 1; j < N; j++) {
      sum = 0;
      for (k = 0; k < i; k++) {
        sum += LU[(k) * (N) + (j)] * LU[(i) * (N) + (k)];
      }
      LU[(i) * (N) + (j)] = (LU[(i) * (N) + (j)] - sum) / LU[(i) * (N) + (i)];
    }
  }
}
```