user: 下面的汇编指令对应的C源代码是什么:

```asm
<gauss_jordan>:
  endbr64
  push   %r15
  mov    %rdx,%rax
  push   %r14
  sub    $0x1,%rax
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdx,-0x70(%rsp)
  mov    %rax,-0x18(%rsp)
  je     L323
  mov    %edx,%esi
  mov    %edx,-0x1c(%rsp)
  mov    %rdi,%r12
  sub    $0x1,%esi
  mov    %esi,-0x34(%rsp)
  js     L131
  sub    $0x1,%eax
  movl   $0xffffffff,-0x74(%rsp)
  lea    0x0(,%rdx,8),%r13
  pxor   %xmm3,%xmm3
  cmp    %esi,%eax
  cmovg  %esi,%eax
  xor    %r15d,%r15d
  xor    %r14d,%r14d
  mov    %r15d,-0x78(%rsp)
  mov    %edx,%r15d
  mov    %eax,-0x38(%rsp)
  movslq %esi,%rax
  lea    -0x2(%rdx),%esi
  mov    %esi,-0x68(%rsp)
  lea    -0x8(%rdi),%rsi
  mov    %rax,-0x30(%rsp)
  shl    $0x3,%rax
  mov    %rsi,-0x40(%rsp)
  lea    -0x3(%rdx),%esi
  mov    %rax,-0x10(%rsp)
  mov    %esi,-0x20(%rsp)
  jmp    L147
  nopl   0x0(%rax,%rax,1)
L246:
  lea    0x1(%r14),%eax
  addl   $0x1,-0x74(%rsp)
  mov    -0x58(%rsp),%r15d
  add    $0x1,%r14
  cmp    %r14d,-0x38(%rsp)
  jl     L278
L147:
  movslq -0x78(%rsp),%rax
  imul   -0x70(%rsp),%rax
  lea    (%rax,%r14,1),%rdx
  shl    $0x3,%rdx
  lea    (%r12,%rdx,1),%rcx
  movsd  (%rcx),%xmm0
  ucomisd %xmm3,%xmm0
  jp     L389
  lea    -0x1(%r15),%ebx
  mov    %ebx,-0x58(%rsp)
  je     L246
L389:
  mov    %r14d,-0x54(%rsp)
  cmp    -0x34(%rsp),%r14d
  jg     L355
  mov    -0x30(%rsp),%rsi
  add    %rax,%rsi
  shl    $0x3,%rsi
  lea    (%r12,%rsi,1),%rdi
  movsd  (%rdi),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rdi)
  cmp    -0x68(%rsp),%r14d
  jg     L355
  lea    -0x1(%r15),%edi
  mov    -0x20(%rsp),%ebx
  mov    -0x74(%rsp),%r9d
  xor    %r10d,%r10d
  mov    %edi,-0x58(%rsp)
  neg    %rdi
  lea    0x8(%rdx),%r8
  cmp    %r9d,%ebx
  lea    0x8(,%rdi,8),%rdi
  cmovl  %r10,%rdi
  lea    -0x8(%rsi,%rdi,1),%rdi
  cmp    %rdi,%r8
  setle  %dil
  cmp    %rdx,%rsi
  setle  %dl
  or     %dil,%dl
  je     L384
  cmp    %r9d,%ebx
  setge  %dil
  cmp    $0x2,%r15d
  setne  %dl
  test   %dl,%dil
  je     L384
  cmp    %r9d,-0x20(%rsp)
  mov    $0x1,%edi
  cmovge -0x58(%rsp),%edi
  lea    -0x10(%r12,%rsi,1),%rdx
  movsd  (%rcx),%xmm1
  mov    %edi,%esi
  mov    %rdx,%r8
  shr    %esi
  unpcklpd %xmm1,%xmm1
  sub    $0x1,%esi
  shl    $0x4,%rsi
  sub    %rsi,%r8
L236:
  movupd (%rdx),%xmm0
  mov    %rdx,%rsi
  sub    $0x10,%rdx
  shufpd $0x1,%xmm0,%xmm0
  divpd  %xmm1,%xmm0
  shufpd $0x1,%xmm0,%xmm0
  movups %xmm0,0x10(%rdx)
  cmp    %rsi,%r8
  jne    L236
  mov    -0x68(%rsp),%edx
  mov    %edi,%esi
  and    $0xfffffffe,%esi
  sub    %esi,%edx
  cmp    %esi,%edi
  je     L303
  movslq %edx,%rdx
  add    %rax,%rdx
  lea    (%r12,%rdx,8),%rdx
  movsd  (%rdx),%xmm0
  divsd  (%rcx),%xmm0
  movsd  %xmm0,(%rdx)
L303:
  mov    -0x78(%rsp),%edx
  test   %edx,%edx
  je     L282
  mov    %rax,%rcx
  mov    -0x74(%rsp),%ebx
  mov    $0x0,%esi
  mov    -0x34(%rsp),%r9d
  neg    %rcx
  mov    %r14,-0x28(%rsp)
  lea    0x0(,%r14,8),%rdx
  mov    -0x54(%rsp),%r14d
  shl    $0x3,%rcx
  mov    %r15d,-0x54(%rsp)
  sub    %r12,%rcx
  mov    %rcx,-0x60(%rsp)
  mov    %r15d,%ecx
  neg    %rcx
  cmp    %ebx,-0x68(%rsp)
  mov    $0x1,%ebx
  lea    0x8(,%rcx,8),%rcx
  cmovge %r15d,%ebx
  mov    %rcx,%rdi
  cmovl  %rsi,%rdi
  mov    -0x30(%rsp),%rsi
  mov    %ebx,-0x48(%rsp)
  xor    %ebp,%ebp
  lea    (%rax,%rsi,1),%rcx
  mov    %rdi,-0x50(%rsp)
  xor    %esi,%esi
  mov    %r12,%rdi
  lea    -0x8(%r12,%rcx,8),%r11
  mov    %ebx,%ecx
  and    $0xfffffffe,%ebx
  sub    %ebx,%r9d
  shr    %ecx
  mov    %ebx,-0x44(%rsp)
  mov    -0x10(%rsp),%rbx
  movslq %r9d,%r9
  lea    -0x1(%rcx),%r8d
  lea    (%rax,%r9,1),%rcx
  not    %r8
  lea    (%r12,%rax,8),%rax
  lea    (%r12,%rcx,8),%r10
  mov    -0x60(%rsp),%rcx
  shl    $0x4,%r8
  mov    %rax,%r15
  add    $0x8,%rcx
  mov    %rcx,-0x60(%rsp)
  nopl   0x0(%rax,%rax,1)
L3:
  mov    -0x74(%rsp),%ecx
  cmp    %ecx,-0x68(%rsp)
  setge  %al
  cmpl   $0x1,-0x54(%rsp)
  setne  %cl
  and    %ecx,%eax
  mov    -0x60(%rsp),%rcx
  add    %rdi,%rcx
  setne  %cl
  test   %cl,%al
  je     L112
  mov    -0x50(%rsp),%rax
  lea    0x8(%rdx),%rcx
  add    %rbx,%rax
  cmp    %rax,%rcx
  lea    0x8(%rbx),%rax
  setle  %cl
  cmp    %rdx,%rax
  setle  %al
  or     %cl,%al
  je     L112
  mov    -0x40(%rsp),%rax
  movsd  (%r12,%rdx,1),%xmm2
  lea    (%rax,%rbx,1),%rcx
  unpcklpd %xmm2,%xmm2
  xor    %eax,%eax
  nopl   0x0(%rax)
L321:
  movupd (%r11,%rax,1),%xmm1
  movupd (%rcx,%rax,1),%xmm0
  shufpd $0x1,%xmm1,%xmm1
  mulpd  %xmm2,%xmm1
  shufpd $0x1,%xmm0,%xmm0
  subpd  %xmm1,%xmm0
  shufpd $0x1,%xmm0,%xmm0
  movups %xmm0,(%rcx,%rax,1)
  sub    $0x10,%rax
  cmp    %r8,%rax
  jne    L321
  mov    -0x44(%rsp),%ecx
  cmp    %ecx,-0x48(%rsp)
  je     L17
  movsd  (%r12,%rdx,1),%xmm1
  mulsd  (%r10),%xmm1
  lea    (%r9,%rsi,1),%rax
  lea    (%r12,%rax,8),%rax
  movsd  (%rax),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L17:
  add    $0x1,%ebp
  add    %r13,%rbx
  add    %r13,%rdx
  add    -0x70(%rsp),%rsi
  add    %r13,%rdi
  cmp    -0x78(%rsp),%ebp
  jne    L3
  mov    -0x28(%rsp),%r14
L282:
  addl   $0x1,-0x78(%rsp)
  jmp    L246
  nopw   0x0(%rax,%rax,1)
L112:
  mov    -0x30(%rsp),%rax
  nopl   (%rax)
L150:
  movsd  (%r12,%rdx,1),%xmm1
  mulsd  (%r15,%rax,8),%xmm1
  movsd  (%rdi,%rax,8),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdi,%rax,8)
  sub    $0x1,%rax
  cmp    %eax,%r14d
  jle    L150
  jmp    L17
L278:
  mov    -0x78(%rsp),%r15d
  cmp    -0x18(%rsp),%r14
  jb     L123
L323:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L131:
  xor    %eax,%eax
  xor    %r15d,%r15d
  pxor   %xmm3,%xmm3
  xor    %r14d,%r14d
L123:
  movslq -0x34(%rsp),%rbx
  movslq %eax,%r8
  lea    -0x1(%rax),%r9d
  add    $0x1,%r8
  mov    %rbx,-0x68(%rsp)
  mov    -0x1c(%rsp),%ebx
  lea    -0x3(%rbx),%ebp
  lea    -0x2(%rbx),%esi
  mov    -0x70(%rsp),%rbx
  mov    %ebp,-0x78(%rsp)
  mov    -0x18(%rsp),%rbp
  mov    %esi,-0x74(%rsp)
  jmp    L107
  nopl   0x0(%rax,%rax,1)
L387:
  add    $0x1,%eax
  lea    0x1(%r8),%rdx
  add    $0x1,%r9d
  mov    %r8,%r14
  cmp    %r8,%rbp
  jbe    L323
  mov    %rdx,%r8
L107:
  movslq %r15d,%rdx
  imul   %rbx,%rdx
  lea    (%rdx,%r14,1),%rcx
  shl    $0x3,%rcx
  lea    (%r12,%rcx,1),%rsi
  movsd  (%rsi),%xmm0
  ucomisd %xmm3,%xmm0
  jp     L385
  je     L387
L385:
  mov    -0x34(%rsp),%r14d
  cmp    %r14d,%eax
  jg     L345
  mov    -0x68(%rsp),%rdi
  add    %rdx,%rdi
  shl    $0x3,%rdi
  lea    (%r12,%rdi,1),%r10
  movsd  (%r10),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%r10)
  cmp    %eax,-0x74(%rsp)
  jl     L345
  sub    %eax,%r14d
  xor    %r13d,%r13d
  lea    0x8(%rcx),%r11
  mov    %r14d,%r10d
  neg    %r10
  cmp    -0x78(%rsp),%r9d
  lea    0x8(,%r10,8),%r10
  cmovg  %r13,%r10
  lea    -0x8(%rdi,%r10,1),%r10
  cmp    %r10,%r11
  setle  %r10b
  cmp    %rcx,%rdi
  setle  %cl
  or     %r10b,%cl
  je     L199
  mov    -0x1c(%rsp),%ecx
  sub    %eax,%ecx
  cmp    $0x2,%ecx
  setne  %r10b
  cmp    -0x78(%rsp),%r9d
  setle  %cl
  test   %cl,%r10b
  je     L199
  cmp    -0x78(%rsp),%r9d
  mov    %r14d,%r10d
  mov    $0x1,%ecx
  movsd  (%rsi),%xmm1
  cmovg  %ecx,%r10d
  lea    -0x10(%r12,%rdi,1),%rdi
  lea    -0x10(%rdi),%rcx
  unpcklpd %xmm1,%xmm1
  mov    %r10d,%r11d
  shr    %r11d
  lea    -0x1(%r11),%r14d
  mov    %rcx,%r11
  shl    $0x4,%r14
  sub    %r14,%r11
  jmp    L353
  nopl   (%rax)
L58:
  sub    $0x10,%rcx
L353:
  movupd (%rdi),%xmm0
  shufpd $0x1,%xmm0,%xmm0
  divpd  %xmm1,%xmm0
  shufpd $0x1,%xmm0,%xmm0
  movups %xmm0,(%rdi)
  mov    %rcx,%rdi
  cmp    %r11,%rcx
  jne    L58
  mov    -0x74(%rsp),%ecx
  mov    %r10d,%edi
  and    $0xfffffffe,%edi
  sub    %edi,%ecx
  cmp    %edi,%r10d
  je     L345
  movslq %ecx,%rcx
  add    %rdx,%rcx
  lea    (%r12,%rcx,8),%rdx
  movsd  (%rdx),%xmm0
  divsd  (%rsi),%xmm0
  movsd  %xmm0,(%rdx)
L345:
  add    $0x1,%r15d
  jmp    L387
L355:
  lea    -0x1(%r15),%ebx
  mov    %ebx,-0x58(%rsp)
  jmp    L303
L199:
  movslq -0x74(%rsp),%rcx
  lea    (%r12,%rdx,8),%rdx
  nop
L96:
  movsd  (%rdx,%rcx,8),%xmm0
  divsd  (%rsi),%xmm0
  movsd  %xmm0,(%rdx,%rcx,8)
  sub    $0x1,%rcx
  cmp    %ecx,%eax
  jle    L96
  jmp    L345
L384:
  movslq -0x68(%rsp),%rdx
  lea    (%r12,%rax,8),%rsi
  mov    %r14d,%edi
  nopl   0x0(%rax)
L324:
  movsd  (%rsi,%rdx,8),%xmm0
  divsd  (%rcx),%xmm0
  movsd  %xmm0,(%rsi,%rdx,8)
  sub    $0x1,%rdx
  cmp    %edx,%edi
  jle    L324
  jmp    L303
```

assistant: ```c
void gauss_jordan(double *a, const size_t rows, const size_t cols) {
  int row = 0;
  for (int col = 0; col < cols - 1; ++col) {
    if (a[row * cols + col] != 0) {
      for (int i = cols - 1; i > col - 1; --i) {
        a[row * cols + i] /= a[row * cols + col];
      }
      for (int i = 0; i < row; ++i) {
        for (int j = cols - 1; j > col - 1; --j) {
          a[i * cols + j] -= a[i * cols + col] * a[row * cols + j];
        }
      }
      row++;
    }
  }
}
```