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

```asm
<denseGETRF>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %rdi,0x10(%rsp)
  mov    %rdx,0x18(%rsp)
  mov    %rcx,0x38(%rsp)
  test   %rdx,%rdx
  je     L90
  lea    (%rdi,%rdx,8),%rax
  mov    %rsi,%r15
  xor    %ebp,%ebp
  movq   $0x0,0x30(%rsp)
  mov    %rax,0x28(%rsp)
  nopl   0x0(%rax)
L192:
  mov    0x10(%rsp),%rax
  mov    (%rax,%rbp,1),%r13
  mov    0x30(%rsp),%rax
  lea    0x1(%rax),%rbx
  mov    %rbx,0x20(%rsp)
  cmp    %rbx,%r15
  jbe    L182
  mov    %rax,0x8(%rsp)
  nopl   (%rax)
L215:
  movsd  0x0(%r13,%rbx,8),%xmm0
  call   <fabs@plt>
  mov    0x8(%rsp),%rdx
  mov    %rax,(%rsp)
  lea    0x0(,%rdx,8),%r14
  lea    0x0(%r13,%r14,1),%r12
  movsd  (%r12),%xmm0
  call   <fabs@plt>
  cmp    %rax,(%rsp)
  jle    L214
  mov    %rbx,0x8(%rsp)
  lea    0x0(,%rbx,8),%r14
  lea    0x0(%r13,%r14,1),%r12
L214:
  add    $0x1,%rbx
  cmp    %rbx,%r15
  ja     L215
  pxor   %xmm4,%xmm4
  ucomisd (%r12),%xmm4
  mov    0x38(%rsp),%rax
  mov    0x8(%rsp),%rsi
  mov    %rsi,(%rax,%rbp,1)
  jp     L82
  je     L53
L82:
  mov    0x30(%rsp),%rsi
  cmp    %rsi,0x8(%rsp)
  jne    L55
L116:
  mov    0x20(%rsp),%rax
  movsd  D4(%rip),%xmm1
  divsd  0x0(%r13,%rbp,1),%xmm1
  nopw   0x0(%rax,%rax,1)
L191:
  movsd  0x0(%r13,%rax,8),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%r15
  ja     L191
L51:
  mov    0x18(%rsp),%rdx
  cmp    %rdx,0x20(%rsp)
  je     L90
  mov    0x20(%rsp),%rdi
  jmp    L72
  nopw   0x0(%rax,%rax,1)
L212:
  add    $0x1,%rdi
  cmp    %rdi,0x18(%rsp)
  jbe    L26
L72:
  mov    0x10(%rsp),%rax
  pxor   %xmm3,%xmm3
  mov    (%rax,%rdi,8),%rcx
  movsd  (%rcx,%rbp,1),%xmm2
  ucomisd %xmm3,%xmm2
  jp     L76
  je     L212
L76:
  mov    0x20(%rsp),%rax
  cmp    %rax,%r15
  jbe    L212
  nop
L30:
  movsd  0x0(%r13,%rax,8),%xmm1
  movsd  (%rcx,%rax,8),%xmm0
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%r15
  ja     L30
  add    $0x1,%rdi
  cmp    %rdi,0x18(%rsp)
  ja     L72
L26:
  mov    0x20(%rsp),%rax
  add    $0x8,%rbp
  mov    %rax,0x30(%rsp)
  jmp    L192
L55:
  mov    0x10(%rsp),%rcx
  nopl   (%rax)
L199:
  mov    (%rcx),%rax
  add    $0x8,%rcx
  lea    (%rax,%r14,1),%r8
  add    %rbp,%rax
  movsd  (%r8),%xmm0
  movsd  (%rax),%xmm1
  movsd  %xmm1,(%r8)
  movsd  %xmm0,(%rax)
  cmp    %rcx,0x28(%rsp)
  jne    L199
  jmp    L116
L90:
  movq   $0x0,0x20(%rsp)
L53:
  mov    0x20(%rsp),%rax
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L182:
  pxor   %xmm5,%xmm5
  ucomisd 0x0(%r13,%rbp,1),%xmm5
  mov    0x38(%rsp),%rax
  mov    0x30(%rsp),%rdx
  mov    %rdx,(%rax,%rbp,1)
  jp     L51
  jne    L51
  jmp    L53
```

tool_calls: name: unpack_data
arguments: data_tag: D4
type_infered: f64

tool: D4: 1

assistant: ```c
long int denseGETRF(double **a, const size_t m, const size_t n, long int *p) {
  long int i, j;
  long int l;
  double *col_j, *col_k;
  double temp, mult, a_kj;
  for (long int k = 0; k < n; k++) {
    col_k = a[k];
    l = k;
    for (i = k + 1; i < m; i++)
      if (fabs(col_k[i]) > fabs(col_k[l]))
        l = i;
    p[k] = l;
    if (col_k[l] == 0.0)
      return (k + 1);
    if (l != k) {
      for (i = 0; i < n; i++) {
        temp = a[i][l];
        a[i][l] = a[i][k];
        a[i][k] = temp;
      }
    }
    mult = 1.0 / col_k[k];
    for (i = k + 1; i < m; i++)
      col_k[i] *= mult;
    for (j = k + 1; j < n; j++) {
      col_j = a[j];
      a_kj = col_j[k];
      if (a_kj != 0.0) {
        for (i = k + 1; i < m; i++)
          col_j[i] -= a_kj * col_k[i];
      }
    }
  }
  return (0);
}
```