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

```asm
<Matrix_inversion>:
  endbr64
  push   %r15
  movslq %edx,%r15
  push   %r14
  push   %r13
  mov    %r15,%r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  mov    %r15d,%edi
  imul   %r15d,%edi
  shl    $0x3,%r15
  sub    $0x48,%rsp
  movslq %edi,%rdi
  mov    %rsi,0x38(%rsp)
  shl    $0x3,%rdi
  call   <malloc@plt>
  mov    %r15,%rdi
  mov    %rax,%rbp
  call   <malloc@plt>
  mov    %rax,%r14
  test   %r13d,%r13d
  jle    L71
  mov    %r13d,%eax
  mov    %r13d,%esi
  mov    %r13d,%ecx
  mov    %r14,%r10
  shr    %eax
  mov    %rsi,0x30(%rsp)
  and    $0xfffffffe,%ecx
  xor    %r9d,%r9d
  lea    -0x1(%rax),%r12d
  lea    0x0(,%rsi,8),%rax
  movl   $0x1,0xc(%rsp)
  mov    %rax,0x28(%rsp)
  add    $0x1,%r12
  lea    (%r14,%rax,1),%r8
  mov    %ecx,%r14d
  shl    $0x4,%r12
  nop
L25:
  mov    %rbx,%rdi
  mov    %rbp,%rsi
  xor    %edx,%edx
  xor    %r11d,%r11d
  nopl   0x0(%rax,%rax,1)
L56:
  xor    %eax,%eax
  cmp    $0x1,%r13d
  je     L58
  nopl   0x0(%rax)
L9:
  movupd (%rdi,%rax,1),%xmm1
  movups %xmm1,(%rsi,%rax,1)
  add    $0x10,%rax
  cmp    %r12,%rax
  jne    L9
  mov    %r14d,%eax
  cmp    %r13d,%r14d
  je     L37
L79:
  add    %edx,%eax
  cltq
  movsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,0x0(%rbp,%rax,8)
L37:
  add    $0x1,%r11d
  add    %r13d,%edx
  add    %r15,%rdi
  add    %r15,%rsi
  cmp    %r11d,%r13d
  jne    L56
  mov    0x28(%rsp),%rdx
  mov    %r10,%rdi
  xor    %esi,%esi
  mov    %r8,0x20(%rsp)
  mov    %r9,0x10(%rsp)
  call   <memset@plt>
  mov    0x10(%rsp),%r9
  mov    %r13d,%edx
  mov    %rbp,%rdi
  mov    %rax,%r10
  mov    D1(%rip),%rax
  mov    %r10,%rsi
  mov    %r9,0x18(%rsp)
  mov    %rax,(%r10,%r9,8)
  xor    %eax,%eax
  mov    %r10,0x10(%rsp)
  call   <Cholesky_decomposition@plt>
  imul   0xc(%rsp),%eax
  mov    0x18(%rsp),%r9
  mov    0x10(%rsp),%r10
  mov    0x20(%rsp),%r8
  mov    %eax,0xc(%rsp)
  mov    0x38(%rsp),%rax
  lea    (%rax,%r9,8),%rdx
  mov    %r10,%rax
  nopl   (%rax)
L75:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,(%rdx)
  add    %r15,%rdx
  cmp    %rax,%r8
  jne    L75
  add    $0x1,%r9
  cmp    0x30(%rsp),%r9
  jne    L25
  mov    %r10,%r14
L83:
  mov    %r14,%rdi
  call   <free@plt>
  mov    %rbp,%rdi
  call   <free@plt>
  mov    0xc(%rsp),%eax
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L58:
  xor    %eax,%eax
  jmp    L79
L71:
  movl   $0x1,0xc(%rsp)
  jmp    L83
```

tool_calls: name: parse_data
arguments: data_tag: D1
datatype: f64

tool: D1: 1

assistant: ```c
int Matrix_inversion(double *ata, double *ata_inv, int n) {
  double *tmpa, *b;
  int *l;
  int r_value;
  register int i, j, k;
  tmpa = malloc(n * n * sizeof(double));
  b = malloc(n * sizeof(double));
  l = malloc(n * sizeof(int));
  r_value = 1;
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      for (k = 0; k < n; k++)
        tmpa[j * n + k] = ata[j * n + k];
      b[j] = 0.0;
    }
    b[i] = 1.0;
    r_value *= Cholesky_decomposition(tmpa, b, n);
    for (j = 0; j < n; j++)
      ata_inv[j * n + i] = b[j];
  }
  free(l);
  free(b);
  free(tmpa);
  return r_value;
}
```