user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<_glp_qmd_genqmd>:
  endbr64
  push   %r15
  mov    %r8,%r11
  mov    %r9,%r15
  push   %r14
  push   %r13
  push   %r12
  mov    %rcx,%r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x98,%rsp
  mov    (%rdi),%r10d
  mov    0xd8(%rsp),%rax
  mov    %rdx,0x30(%rsp)
  mov    0xe8(%rsp),%r8
  mov    0xf8(%rsp),%rdx
  mov    %rdi,0x58(%rsp)
  mov    %rax,0x20(%rsp)
  mov    0xf0(%rsp),%r9
  mov    0xe0(%rsp),%rax
  mov    %rsi,0x28(%rsp)
  mov    %r8,0x68(%rsp)
  mov    0xd0(%rsp),%r14
  mov    %rax,0x38(%rsp)
  mov    %r9,0x48(%rsp)
  mov    %rdx,0x8(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x88(%rsp)
  xor    %eax,%eax
  movl   $0x0,(%rdx)
  mov    (%rdi),%edi
  movl   $0x1,0x80(%rsp)
  test   %edi,%edi
  jle    L154
  mov    %rsi,%rbp
  mov    $0x1,%eax
  mov    %rbx,%rsi
  mov    %rbp,%rcx
  nopl   0x0(%rax,%rax,1)
L147:
  mov    %eax,(%r12,%rax,4)
  mov    %eax,(%r11,%rax,4)
  movl   $0x0,(%r14,%rax,4)
  movl   $0x1,(%r8,%rax,4)
  movl   $0x0,(%r9,%rax,4)
  mov    0x4(%rcx,%rax,4),%edx
  sub    (%rcx,%rax,4),%edx
  mov    %edx,(%r15,%rax,4)
  cmp    %edx,%r10d
  mov    (%rsi),%edi
  cmovg  %edx,%r10d
  add    $0x1,%rax
  cmp    %eax,%edi
  jge    L147
  mov    0x8(%rsp),%rax
  mov    (%rax),%ecx
L120:
  lea    0x84(%rsp),%rax
  mov    $0x1,%ebp
  mov    %rax,0x40(%rsp)
  lea    0x80(%rsp),%rax
  mov    %rax,0x50(%rsp)
  lea    0x7c(%rsp),%rax
  mov    %rax,0x60(%rsp)
  nopl   0x0(%rax)
L194:
  mov    %edi,%r13d
  mov    $0x1,%ebx
L51:
  cmp    %ebp,%ebx
  cmovl  %ebp,%ebx
  cmp    %edi,%ebx
  jg     L170
  movslq %ebx,%rax
  nopl   0x0(%rax)
L179:
  movslq (%r12,%rax,4),%rdx
  mov    %eax,%ebx
  lea    (%r14,%rdx,4),%rsi
  mov    %rdx,%r8
  mov    (%rsi),%r9d
  test   %r9d,%r9d
  js     L150
  mov    (%r15,%rdx,4),%edx
  cmp    %edx,%r10d
  jge    L175
  cmp    %edx,%r13d
  cmovg  %edx,%r13d
L150:
  add    $0x1,%rax
  cmp    %eax,%edi
  jge    L179
L170:
  mov    %r13d,%r10d
  jmp    L194
  nopw   0x0(%rax,%rax,1)
L175:
  mov    0x8(%rsp),%rax
  mov    %r10d,0x1c(%rsp)
  add    %ecx,%edx
  sub    $0x8,%rsp
  mov    %r8d,0x88(%rsp)
  mov    %r15,%rcx
  mov    %r14,%r8
  mov    %r11,0x18(%rsp)
  mov    %edx,(%rax)
  movl   $0x1,(%rsi)
  push   0x40(%rsp)
  push   0x70(%rsp)
  push   0x38(%rsp)
  mov    0x60(%rsp),%r9
  mov    0x50(%rsp),%rdx
  mov    0x48(%rsp),%rsi
  mov    0x70(%rsp),%rdi
  call   <_glp_qmd_qmdrch@plt>
  mov    0xa0(%rsp),%ecx
  movslq %ebp,%rax
  mov    0x3c(%rsp),%r10d
  mov    0x30(%rsp),%r11
  add    $0x20,%rsp
  mov    %r13d,0x18(%rsp)
  mov    0x48(%rsp),%r13
  nopw   0x0(%rax,%rax,1)
L73:
  movslq %ecx,%rsi
  movslq (%r12,%rax,4),%rdx
  mov    %eax,%ebp
  lea    (%r11,%rsi,4),%rdi
  movslq (%rdi),%r9
  mov    %edx,(%r12,%r9,4)
  mov    %r9d,(%r11,%rdx,4)
  mov    %ecx,(%r12,%rax,4)
  mov    %eax,(%rdi)
  add    $0x1,%rax
  movl   $0xffffffff,(%r15,%rsi,4)
  mov    0x0(%r13,%rsi,4),%ecx
  test   %ecx,%ecx
  jg     L73
  movslq 0x84(%rsp),%rax
  mov    0x18(%rsp),%r13d
  test   %eax,%eax
  jg     L5
L3:
  mov    0x58(%rsp),%rax
  mov    (%rax),%edi
  cmp    %ebp,%edi
  jle    L115
L171:
  mov    0x8(%rsp),%rax
  add    $0x1,%ebp
  mov    (%rax),%ecx
  jmp    L51
L5:
  mov    0x38(%rsp),%rcx
  movslq 0x7c(%rsp),%rdx
  mov    %r11,0x10(%rsp)
  mov    %r15,%r8
  mov    %r10d,0x18(%rsp)
  lea    0x4(%rcx,%rdx,4),%rdx
  push   %rdx
  mov    0x28(%rsp),%rcx
  lea    0x4(%rcx,%rax,4),%rax
  push   %rax
  push   %r14
  push   0x60(%rsp)
  mov    0x60(%rsp),%rdx
  mov    0x88(%rsp),%r9
  mov    0x50(%rsp),%rsi
  mov    0x48(%rsp),%rdi
  mov    %rcx,0x40(%rsp)
  call   <_glp_qmd_qmdupd@plt>
  mov    0xa4(%rsp),%edx
  movslq 0xa0(%rsp),%rax
  add    $0x20,%rsp
  mov    0x18(%rsp),%r10d
  mov    0x10(%rsp),%r11
  test   %edx,%edx
  movl   $0x0,(%r14,%rax,4)
  jle    L93
  mov    0x20(%rsp),%rcx
  sub    $0x1,%edx
  lea    0x4(%rcx),%rax
  lea    0x8(%rcx,%rdx,4),%rsi
  jmp    L33
  nopl   (%rax)
L163:
  cmp    %ecx,%r13d
  cmovg  %ecx,%r13d
L108:
  add    $0x4,%rax
  cmp    %rax,%rsi
  je     L93
L33:
  movslq (%rax),%rdx
  lea    (%r14,%rdx,4),%rcx
  mov    (%rcx),%edi
  test   %edi,%edi
  js     L108
  movl   $0x0,(%rcx)
  mov    (%r15,%rdx,4),%ecx
  cmp    %r10d,%ecx
  jg     L163
  add    $0x4,%rax
  mov    %r10d,%r13d
  mov    (%r11,%rdx,4),%ebx
  mov    %ecx,%r10d
  cmp    %rax,%rsi
  jne    L33
L93:
  mov    0x7c(%rsp),%ecx
  test   %ecx,%ecx
  jle    L3
  mov    %r11,0x10(%rsp)
  sub    $0x8,%rsp
  mov    %r14,%rcx
  mov    %r10d,0x20(%rsp)
  push   0x40(%rsp)
  mov    0x60(%rsp),%rdi
  mov    0x30(%rsp),%r9
  mov    0x50(%rsp),%r8
  mov    0x40(%rsp),%rdx
  mov    0x38(%rsp),%rsi
  call   <_glp_qmd_qmdqt@plt>
  pop    %rax
  pop    %rdx
  mov    0x58(%rsp),%rax
  mov    0x10(%rsp),%r11
  mov    0x18(%rsp),%r10d
  mov    (%rax),%edi
  cmp    %ebp,%edi
  jg     L171
L115:
  mov    0x88(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L101
  add    $0x98,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L154:
  xor    %ecx,%ecx
  jmp    L120
L101:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void _glp_qmd_genqmd(int *_neqns, int xadj[], int adjncy[], int perm[], int invp[], int deg[], int marker[], int rchset[], int nbrhd[], int qsize[], int qlink[], int *_nofsub) {
  int inode, ip, irch, j, mindeg, ndeg, nhdsze, node, np, num, nump1, nxnode, rchsze, search, thresh;
  mindeg = (*_neqns);
  (*_nofsub) = 0;
  for (node = 1; node <= (*_neqns); node++) {
    perm[node] = node;
    invp[node] = node;
    marker[node] = 0;
    qsize[node] = 1;
    qlink[node] = 0;
    ndeg = xadj[node + 1] - xadj[node];
    deg[node] = ndeg;
    if (ndeg < mindeg)
      mindeg = ndeg;
  }
  num = 0;
s200:
  search = 1;
  thresh = mindeg;
  mindeg = (*_neqns);
s300:
  nump1 = num + 1;
  if (nump1 > search)
    search = nump1;
  for (j = search; j <= (*_neqns); j++) {
    node = perm[j];
    if (marker[node] >= 0) {
      ndeg = deg[node];
      if (ndeg <= thresh)
        goto s500;
      if (ndeg < mindeg)
        mindeg = ndeg;
    }
  }
  goto s200;
s500:
  search = j;
  (*_nofsub) += deg[node];
  marker[node] = 1;
  _glp_qmd_qmdrch(&node, xadj, adjncy, deg, marker, &rchsze, rchset, &nhdsze, nbrhd);
  nxnode = node;
s600:
  num++;
  np = invp[nxnode];
  ip = perm[num];
  perm[np] = ip;
  invp[ip] = np;
  perm[num] = nxnode;
  invp[nxnode] = num;
  deg[nxnode] = -1;
  nxnode = qlink[nxnode];
  if (nxnode > 0)
    goto s600;
  if (rchsze > 0) {
    _glp_qmd_qmdupd(xadj, adjncy, &rchsze, rchset, deg, qsize, qlink, marker, &rchset[rchsze + 1], &nbrhd[nhdsze + 1]);
    marker[node] = 0;
    for (irch = 1; irch <= rchsze; irch++) {
      inode = rchset[irch];
      if (marker[inode] >= 0) {
        marker[inode] = 0;
        ndeg = deg[inode];
        if (ndeg < mindeg)
          mindeg = ndeg;
        if (ndeg <= thresh) {
          mindeg = thresh;
          thresh = ndeg;
          search = invp[inode];
        }
      }
    }
    if (nhdsze > 0)
      _glp_qmd_qmdqt(&node, xadj, adjncy, marker, &rchsze, rchset, nbrhd);
  }
  if (num < (*_neqns))
    goto s300;
  return;
}
```