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

```asm
<sp_coletree>:
  endbr64
  push   %r15
  mov    %r8d,%r15d
  push   %r14
  mov    %rdx,%r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %ecx,%ebx
  sub    $0x38,%rsp
  mov    %rdi,0x20(%rsp)
  mov    %r8d,%edi
  mov    %rsi,0x8(%rsp)
  mov    %r8d,0x2c(%rsp)
  mov    %r9,0x10(%rsp)
  call   <mxCallocInt@plt>
  mov    %r15d,%edi
  mov    %rax,%rbp
  call   <initialize_disjoint_sets@plt>
  mov    %ebx,%edi
  call   <mxCallocInt@plt>
  mov    %rax,%r15
  test   %ebx,%ebx
  jle    L13
  lea    -0x1(%rbx),%eax
  cmp    $0x2,%eax
  jbe    L8
  mov    %ebx,%ecx
  movd   0x2c(%rsp),%xmm1
  mov    %r15,%rax
  shr    $0x2,%ecx
  shl    $0x4,%rcx
  pshufd $0x0,%xmm1,%xmm0
  add    %r15,%rcx
  nopl   (%rax)
L37:
  movups %xmm0,(%rax)
  add    $0x10,%rax
  cmp    %rcx,%rax
  jne    L37
  mov    %ebx,%eax
  and    $0xfffffffc,%eax
  test   $0x3,%bl
  je     L13
L60:
  mov    0x2c(%rsp),%edx
  movslq %eax,%rcx
  lea    0x0(,%rcx,4),%rsi
  mov    %edx,(%r15,%rcx,4)
  lea    0x1(%rax),%ecx
  cmp    %ebx,%ecx
  jge    L13
  add    $0x2,%eax
  mov    %edx,0x4(%r15,%rsi,1)
  cmp    %eax,%ebx
  jle    L13
  mov    %edx,0x8(%r15,%rsi,1)
L13:
  movslq 0x2c(%rsp),%rbx
  mov    0x8(%rsp),%r9
  xor    %r10d,%r10d
  mov    %rbx,0x18(%rsp)
  test   %ebx,%ebx
  jle    L58
  nopl   0x0(%rax,%rax,1)
L35:
  mov    0x20(%rsp),%rax
  mov    %r10d,%edi
  movslq (%rax,%r10,4),%rcx
  cmp    (%r9),%ecx
  jge    L54
  nopl   0x0(%rax)
L102:
  movslq (%r14,%rcx,4),%rax
  lea    (%r15,%rax,4),%rsi
  mov    (%rsi),%eax
  cmp    %edi,%eax
  cmovg  %edi,%eax
  add    $0x1,%rcx
  mov    %eax,(%rsi)
  cmp    %ecx,(%r9)
  jg     L102
L54:
  add    $0x1,%r10
  add    $0x4,%r9
  cmp    %r10,0x18(%rsp)
  jne    L35
  xor    %r13d,%r13d
  nopl   0x0(%rax)
L146:
  mov    %r13d,%edi
  mov    %r13d,0x28(%rsp)
  call   <make_set@plt>
  mov    0x2c(%rsp),%edx
  mov    %eax,%ebx
  cltq
  mov    %r13d,0x0(%rbp,%rax,4)
  mov    0x10(%rsp),%rax
  mov    %edx,(%rax,%r13,4)
  mov    0x20(%rsp),%rax
  mov    0x8(%rsp),%rdx
  movslq (%rax,%r13,4),%r12
  mov    (%rdx,%r13,4),%esi
  cmp    %esi,%r12d
  jge    L104
  nopw   0x0(%rax,%rax,1)
L87:
  movslq (%r14,%r12,4),%rax
  mov    (%r15,%rax,4),%edi
  cmp    %r13d,%edi
  jge    L149
  call   <find@plt>
  mov    %eax,%esi
  cltq
  movslq 0x0(%rbp,%rax,4),%rax
  cmp    %r13d,%eax
  je     L121
  mov    0x28(%rsp),%edx
  mov    0x10(%rsp),%rcx
  mov    %ebx,%edi
  mov    %edx,(%rcx,%rax,4)
  call   <link@plt>
  mov    0x28(%rsp),%edx
  mov    %eax,%ebx
  cltq
  mov    %edx,0x0(%rbp,%rax,4)
L121:
  mov    0x8(%rsp),%rax
  mov    (%rax,%r13,4),%esi
L149:
  add    $0x1,%r12
  cmp    %r12d,%esi
  jg     L87
L104:
  add    $0x1,%r13
  cmp    %r13,0x18(%rsp)
  jne    L146
L58:
  mov    %rbp,%rdi
  call   <superlu_free@plt>
  mov    %r15,%rdi
  call   <superlu_free@plt>
  xor    %eax,%eax
  call   <finalize_disjoint_sets@plt>
  add    $0x38,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L8:
  xor    %eax,%eax
  jmp    L60
```

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;
}
```