user: What would the C source look like for the assembly code provided here:

```asm
<lu_invert>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %esi,%ebx
  sub    $0x58,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x48(%rsp)
  xor    %eax,%eax
  cmp    $0xa,%esi
  jg     L41
  lea    0x20(%rsp),%rax
  lea    0x18(%rsp),%rcx
  mov    %rax,%rdx
  mov    %rax,0x8(%rsp)
  call   <lu_decomp@plt>
  test   %rax,%rax
  jne    L16
  lea    -0x1(%rbx),%r12d
  xor    %edx,%edx
  xor    %edi,%edi
  mov    %r12d,%ecx
  mov    %r12d,%esi
  call   <dmatrix@plt>
  mov    0x8(%rsp),%r14
  mov    %rax,%r15
  test   %ebx,%ebx
  jle    L26
L70:
  xor    %esi,%esi
  nopl   (%rax)
L0:
  mov    0x0(%rbp,%rsi,8),%rcx
  mov    (%r15,%rsi,8),%rdx
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L91:
  movsd  (%rcx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%ebx
  jg     L91
  add    $0x1,%rsi
  cmp    %esi,%ebx
  jg     L0
  xor    %r13d,%r13d
  nopl   (%rax)
L88:
  mov    0x0(%rbp,%r13,8),%rcx
  lea    0x0(,%r13,8),%rdx
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L71:
  movq   $0x0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%ebx
  jg     L71
  mov    D8(%rip),%rax
  mov    %ebx,%esi
  mov    %r15,%rdi
  add    $0x1,%r13
  mov    %rax,(%rcx,%rdx,1)
  mov    %r14,%rdx
  call   <lu_backsub@plt>
  cmp    %r13d,%ebx
  jg     L88
  mov    %r15,%rdi
  mov    %r12d,%r8d
  xor    %ecx,%ecx
  mov    %r12d,%edx
  xor    %esi,%esi
  xor    %r15d,%r15d
  call   <free_dmatrix@plt>
  cmp    0x8(%rsp),%r14
  je     L33
  mov    %r12d,%edx
  xor    %esi,%esi
  mov    %r14,%rdi
  call   <free_ivector@plt>
  jmp    L33
L41:
  lea    -0x1(%rsi),%r12d
  xor    %edi,%edi
  mov    %r12d,%esi
  call   <ivector@plt>
  lea    0x18(%rsp),%rcx
  mov    %ebx,%esi
  mov    %rbp,%rdi
  mov    %rax,%rdx
  mov    %rax,%r14
  call   <lu_decomp@plt>
  test   %rax,%rax
  je     L8
  lea    0x20(%rsp),%rax
  cmp    %rax,%r14
  je     L16
  mov    %r12d,%edx
  xor    %esi,%esi
  mov    %r14,%rdi
  call   <free_ivector@plt>
L16:
  mov    $0x1,%r15d
L33:
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L54
  add    $0x58,%rsp
  mov    %r15d,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L26:
  mov    %r12d,%r8d
  xor    %ecx,%ecx
  mov    %r12d,%edx
  xor    %esi,%esi
  mov    %rax,%rdi
  xor    %r15d,%r15d
  call   <free_dmatrix@plt>
  jmp    L33
L8:
  mov    %r12d,%ecx
  xor    %edx,%edx
  mov    %r12d,%esi
  xor    %edi,%edi
  call   <dmatrix@plt>
  mov    %rax,%r15
  lea    0x20(%rsp),%rax
  mov    %rax,0x8(%rsp)
  jmp    L70
L54:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D8
type_guessed: f64

tool: D8: 1

assistant: ```c
int lu_invert(double **a, int n) {
  int i, j;
  double rip;
  int *pivx, PIVX[10];
  double **y;
  if (n <= 10)
    pivx = PIVX;
  else
    pivx = ivector(0, n - 1);
  if (lu_decomp(a, n, pivx, &rip)) {
    if (pivx != PIVX)
      free_ivector(pivx, 0, n - 1);
    return 1;
  }
  y = dmatrix(0, n - 1, 0, n - 1);
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      y[i][j] = a[i][j];
    }
  }
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++)
      a[i][j] = 0.0;
    a[i][i] = 1.0;
    lu_backsub(y, n, pivx, a[i]);
  }
  free_dmatrix(y, 0, n - 1, 0, n - 1);
  if (pivx != PIVX)
    free_ivector(pivx, 0, n - 1);
  return 0;
}
```