user: Demonstrate the C code that generates the assembly output below:

```asm
<mark_cluster_numbers>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,%r13
  mov    %rsi,0x10(%rsp)
  mov    %rdx,0x18(%rsp)
  mov    %ecx,%ebx
  mov    %ecx,0x20(%rsp)
  mov    %r8d,%edx
  mov    %r8d,0xc(%rsp)
  mov    %r9,0x38(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x78(%rsp)
  xor    %eax,%eax
  mov    $0xffffffff,%ecx
  mov    %ebx,%esi
  call   <fill_matrix@plt>
  test   %ebx,%ebx
  jle    L56
  mov    %ebx,%eax
  mov    %rax,0x30(%rsp)
  mov    $0x0,%r15d
  movl   $0x0,0x8(%rsp)
  lea    0x50(%rsp),%rax
  mov    %rax,0x28(%rsp)
  jmp    L13
L95:
  add    $0x1,%r12d
  add    $0x8,%rbp
  cmp    %r12d,0x44(%rsp)
  jle    L93
L37:
  movslq 0x0(%rbp),%rax
  movslq 0x4(%rbp),%rdx
  mov    0x0(%r13,%rax,8),%rax
  cmpl   $0xffffffff,(%rax,%rdx,4)
  jne    L95
  sub    $0x8,%rsp
  mov    0x10(%rsp),%eax
  push   %rax
  mov    %rbp,%r9
  mov    0x1c(%rsp),%r8d
  mov    0x30(%rsp),%ecx
  mov    0x28(%rsp),%rdx
  mov    0x20(%rsp),%rsi
  mov    %r13,%rdi
  call   <mark_one_cluster_aux@plt>
  add    $0x10,%rsp
  jmp    L95
L93:
  addl   $0x1,0x8(%rsp)
L51:
  add    $0x1,%rbx
  cmp    %rbx,%r14
  je     L67
L12:
  mov    0x0(%r13,%r15,8),%rax
  cmpl   $0x0,(%rax,%rbx,4)
  jns    L51
  mov    0x24(%rsp),%eax
  mov    %eax,0x48(%rsp)
  mov    %ebx,0x4c(%rsp)
  sub    $0x8,%rsp
  lea    0x4c(%rsp),%rax
  push   %rax
  mov    0x38(%rsp),%r9
  lea    0x58(%rsp),%r8
  mov    0x1c(%rsp),%ecx
  mov    0x30(%rsp),%edx
  mov    0x28(%rsp),%rsi
  mov    0x20(%rsp),%rdi
  call   <get_bonded_neighbors@plt>
  mov    0x0(%r13,%r15,8),%rax
  mov    0x18(%rsp),%ecx
  mov    %ecx,(%rax,%rbx,4)
  add    $0x10,%rsp
  cmpl   $0x0,0x44(%rsp)
  jle    L93
  lea    0x50(%rsp),%rbp
  mov    $0x0,%r12d
  jmp    L37
L67:
  add    $0x1,%r15
  cmp    0x30(%rsp),%r15
  je     L72
L13:
  mov    %r15d,0x24(%rsp)
  mov    0xc(%rsp),%eax
  mov    %eax,%r14d
  mov    $0x0,%ebx
  test   %eax,%eax
  jg     L12
  jmp    L67
L56:
  movl   $0x0,0x8(%rsp)
L72:
  mov    0x38(%rsp),%rax
  test   %rax,%rax
  je     L39
  mov    0x8(%rsp),%edi
  mov    %edi,(%rax)
L39:
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L47
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L47:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void mark_cluster_numbers(int **site_marks, int **vbonds, int **hbonds, int M, int N, int *pnum_clusters) {
  int neighbors[4][2];
  int numnei;
  int i, j, m;
  int A[2];
  int cluster_number = 0;
  fill_matrix(site_marks, M, N, -1);
  for (i = 0; i < M; i++) {
    for (j = 0; j < N; j++) {
      if (site_marks[i][j] >= 0)
        continue;
      A[0] = i;
      A[1] = j;
      get_bonded_neighbors(vbonds, hbonds, M, N, A, neighbors, &numnei);
      site_marks[i][j] = cluster_number;
      for (m = 0; m < numnei; m++) {
        int Bi = neighbors[m][0];
        int Bj = neighbors[m][1];
        if (site_marks[Bi][Bj] == -1) {
          mark_one_cluster_aux(site_marks, vbonds, hbonds, M, N, neighbors[m], cluster_number);
        }
      }
      cluster_number++;
    }
  }
  if (pnum_clusters)
    *pnum_clusters = cluster_number;
}
```