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

```asm
<sp_coletree>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,0x18(%rsp)
  mov    %rsi,%rbp
  mov    %rdx,%r13
  mov    %ecx,%ebx
  mov    %r8d,%r14d
  mov    %r8d,0x14(%rsp)
  mov    %r9,0x8(%rsp)
  mov    %r8d,%edi
  call   <mxCallocInt@plt>
  mov    %rax,%r15
  mov    %r14d,%edi
  call   <initialize_disjoint_sets@plt>
  mov    %ebx,%edi
  call   <mxCallocInt@plt>
  mov    %rax,%r12
  test   %ebx,%ebx
  jle    L51
  mov    %ebx,%ebx
  lea    (%rax,%rbx,4),%rdx
L18:
  mov    0x14(%rsp),%ebx
  mov    %ebx,(%rax)
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L18
L51:
  mov    0x14(%rsp),%eax
  test   %eax,%eax
  jle    L32
  mov    %rbp,0x20(%rsp)
  mov    %eax,%eax
  mov    %rax,0x28(%rsp)
  mov    %rbp,%r9
  mov    $0x0,%r8d
L46:
  mov    %r8d,%esi
  mov    0x18(%rsp),%rax
  mov    (%rax,%r8,4),%edx
  mov    %r9,%rdi
  cmp    (%r9),%edx
  jge    L38
  movslq %edx,%rdx
L105:
  movslq 0x0(%r13,%rdx,4),%rax
  lea    (%r12,%rax,4),%rcx
  mov    (%rcx),%eax
  cmp    %esi,%eax
  cmovg  %esi,%eax
  mov    %eax,(%rcx)
  add    $0x1,%rdx
  cmp    %edx,(%rdi)
  jg     L105
L38:
  add    $0x1,%r8
  add    $0x4,%r9
  cmp    0x28(%rsp),%r8
  jne    L46
  mov    $0x0,%ebp
  jmp    L120
L8:
  add    $0x1,%rbx
  cmp    %ebx,(%r14)
  jle    L70
L57:
  movslq 0x0(%r13,%rbx,4),%rax
  mov    (%r12,%rax,4),%edi
  cmp    %ebp,%edi
  jge    L8
  call   <find@plt>
  mov    %eax,%esi
  cltq
  mov    (%r15,%rax,4),%eax
  cmp    %ebp,%eax
  je     L8
  cltq
  mov    0x8(%rsp),%rdx
  mov    0x10(%rsp),%ecx
  mov    %ecx,(%rdx,%rax,4)
  mov    0x4(%rsp),%edi
  call   <link@plt>
  mov    %eax,0x4(%rsp)
  cltq
  mov    0x10(%rsp),%ecx
  mov    %ecx,(%r15,%rax,4)
  jmp    L8
L70:
  add    $0x1,%rbp
  addq   $0x4,0x20(%rsp)
  cmp    0x28(%rsp),%rbp
  je     L32
L120:
  mov    %ebp,%edi
  mov    %ebp,0x10(%rsp)
  call   <make_set@plt>
  mov    %eax,0x4(%rsp)
  cltq
  mov    %ebp,(%r15,%rax,4)
  mov    0x8(%rsp),%rax
  mov    0x14(%rsp),%esi
  mov    %esi,(%rax,%rbp,4)
  mov    0x18(%rsp),%rax
  mov    (%rax,%rbp,4),%ebx
  mov    0x20(%rsp),%rax
  mov    %rax,%r14
  cmp    (%rax),%ebx
  jge    L70
  movslq %ebx,%rbx
  jmp    L57
L32:
  mov    %r15,%rdi
  call   <superlu_free@plt>
  mov    %r12,%rdi
  call   <superlu_free@plt>
  mov    $0x0,%eax
  call   <finalize_disjoint_sets@plt>
  mov    $0x0,%eax
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
int sp_coletree(int *acolst, int *acolend, int *arow, int nr, int nc, int *parent) {
  int *root;
  int *firstcol;
  int rset, cset;
  int row, col;
  int rroot;
  int p;
  root = mxCallocInt(nc);
  initialize_disjoint_sets(nc);
  firstcol = mxCallocInt(nr);
  for (row = 0; row < nr; firstcol[row++] = nc)
    ;
  for (col = 0; col < nc; col++)
    for (p = acolst[col]; p < acolend[col]; p++) {
      row = arow[p];
      firstcol[row] = ((firstcol[row]) < (col) ? (firstcol[row]) : (col));
    }
  for (col = 0; col < nc; col++) {
    cset = make_set(col);
    root[cset] = col;
    parent[col] = nc;
    for (p = acolst[col]; p < acolend[col]; p++) {
      row = firstcol[arow[p]];
      if (row >= col)
        continue;
      rset = find(row);
      rroot = root[rset];
      if (rroot != col) {
        parent[rroot] = col;
        cset = link(cset, rset);
        root[cset] = col;
      }
    }
  }
  superlu_free(root);
  superlu_free(firstcol);
  finalize_disjoint_sets();
  return 0;
}
```