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

```asm
<jerasure_invertible_matrix>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rdi,0x28(%rsp)
  mov    %esi,0x4(%rsp)
  test   %esi,%esi
  jle    L133
  movslq 0x4(%rsp),%rbx
  mov    %edx,%r13d
  movq   $0x0,0x18(%rsp)
  movq   $0x0,0x20(%rsp)
  mov    %rbx,%rax
  lea    0x4(,%rbx,4),%rdi
  mov    %rbx,0x48(%rsp)
  mov    %eax,(%rsp)
  mov    0x28(%rsp),%rax
  mov    %rdi,0x50(%rsp)
  lea    -0x1(%rbx),%edi
  mov    %rdi,0x10(%rsp)
  mov    %ebx,%edi
  shl    $0x2,%rbx
  mov    %rax,0x38(%rsp)
  add    $0x4,%rax
  mov    %rdi,0x58(%rsp)
  mov    %rbx,0x8(%rsp)
  mov    %rax,0x30(%rsp)
L15:
  mov    (%rsp),%eax
  sub    0x4(%rsp),%eax
  mov    %eax,0x44(%rsp)
  mov    0x38(%rsp),%rax
  mov    (%rax),%esi
  mov    0x20(%rsp),%eax
  lea    0x1(%rax),%ebx
  test   %esi,%esi
  je     L208
  cmp    $0x1,%esi
  jne    L51
L43:
  mov    0x10(%rsp),%rdi
  cmp    %rdi,0x20(%rsp)
  je     L133
  mov    0x8(%rsp),%rbp
  add    0x38(%rsp),%rbp
  nopl   0x0(%rax)
L45:
  mov    0x0(%rbp),%r12d
  test   %r12d,%r12d
  je     L143
  cmp    $0x1,%r12d
  je     L197
  movslq (%rsp),%rax
  mov    0x28(%rsp),%rdi
  mov    %ebx,0x40(%rsp)
  lea    (%rdi,%rax,4),%r14
  mov    0x10(%rsp),%rdi
  lea    (%rax,%rdi,1),%rdx
  mov    0x30(%rsp),%rdi
  lea    (%rdi,%rdx,4),%rdx
  mov    0x18(%rsp),%rdi
  mov    %rdx,%r15
  sub    %rax,%rdi
  mov    %rdi,%rbx
  nopl   0x0(%rax)
L47:
  mov    (%r14,%rbx,4),%esi
  mov    %r13d,%edx
  mov    %r12d,%edi
  call   <galois_single_multiply@plt>
  xor    %eax,(%r14)
  add    $0x4,%r14
  cmp    %r14,%r15
  jne    L47
  mov    0x40(%rsp),%ebx
L143:
  mov    0x4(%rsp),%ecx
  add    $0x1,%ebx
  add    %ecx,(%rsp)
  add    0x8(%rsp),%rbp
  cmp    %ebx,%ecx
  jne    L45
  mov    0x50(%rsp),%rbx
  mov    0x4(%rsp),%edi
  add    %rbx,0x38(%rsp)
  mov    0x44(%rsp),%ebx
  addq   $0x1,0x20(%rsp)
  mov    0x20(%rsp),%rax
  lea    (%rbx,%rdi,2),%ebx
  mov    0x48(%rsp),%rdi
  add    %rdi,0x18(%rsp)
  mov    %ebx,(%rsp)
  cmp    0x58(%rsp),%rax
  jne    L15
L133:
  mov    $0x1,%esi
L171:
  add    $0x68,%rsp
  mov    %esi,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L197:
  movslq (%rsp),%rdx
  mov    0x10(%rsp),%rdi
  mov    0x28(%rsp),%rax
  lea    (%rdx,%rdi,1),%rcx
  mov    0x30(%rsp),%rdi
  lea    (%rax,%rdx,4),%rax
  lea    (%rdi,%rcx,4),%rsi
  mov    0x18(%rsp),%rcx
  sub    %rdx,%rcx
  nopw   0x0(%rax,%rax,1)
L46:
  mov    (%rax,%rcx,4),%edx
  xor    %edx,(%rax)
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L46
  jmp    L143
  nopl   0x0(%rax,%rax,1)
L208:
  mov    0x4(%rsp),%r8d
  cmp    %ebx,%r8d
  jle    L23
  movslq (%rsp),%rax
  mov    0x28(%rsp),%rdx
  mov    0x8(%rsp),%r9
  mov    %rax,%rdi
  add    0x20(%rsp),%rax
  lea    (%rdx,%rax,4),%rcx
  mov    %edi,%eax
  mov    %ebx,%edi
  jmp    L170
  cs nopw 0x0(%rax,%rax,1)
L123:
  add    $0x1,%edi
  add    %r8d,%eax
  add    %r9,%rcx
  cmp    %edi,%r8d
  je     L171
L170:
  mov    (%rcx),%esi
  mov    %eax,%edx
  test   %esi,%esi
  je     L123
L86:
  mov    0x18(%rsp),%rsi
  mov    0x10(%rsp),%rcx
  movslq %edx,%rdx
  mov    0x28(%rsp),%rax
  mov    0x30(%rsp),%rdi
  add    %rsi,%rcx
  sub    %rsi,%rdx
  lea    (%rax,%rsi,4),%rax
  lea    (%rdi,%rcx,4),%rdi
  nop
L87:
  mov    (%rax),%ecx
  mov    (%rax,%rdx,4),%esi
  mov    %esi,(%rax)
  mov    %ecx,(%rax,%rdx,4)
  add    $0x4,%rax
  cmp    %rdi,%rax
  jne    L87
  mov    0x38(%rsp),%rax
  mov    (%rax),%esi
  cmp    $0x1,%esi
  je     L43
L51:
  mov    $0x1,%edi
  mov    %r13d,%edx
  call   <galois_single_divide@plt>
  mov    0x28(%rsp),%rdi
  mov    %eax,%ebp
  mov    0x18(%rsp),%rax
  lea    (%rdi,%rax,4),%r12
  mov    0x30(%rsp),%rdi
  add    0x10(%rsp),%rax
  lea    (%rdi,%rax,4),%r14
  xchg   %ax,%ax
L140:
  mov    (%r12),%edi
  mov    %r13d,%edx
  mov    %ebp,%esi
  add    $0x4,%r12
  call   <galois_single_multiply@plt>
  mov    %eax,-0x4(%r12)
  cmp    %r14,%r12
  jne    L140
  jmp    L43
L23:
  je     L171
  mov    0x4(%rsp),%edx
  imul   %ebx,%edx
  jmp    L86
```

assistant: ```c
int jerasure_invertible_matrix(int *mat, int rows, int w) {
  int cols, i, j, k, x, rs2;
  int row_start, tmp, inverse;
  cols = rows;
  for (i = 0; i < cols; i++) {
    row_start = cols * i;
    if (mat[row_start + i] == 0) {
      for (j = i + 1; j < rows && mat[cols * j + i] == 0; j++)
        ;
      if (j == rows)
        return 0;
      rs2 = j * cols;
      for (k = 0; k < cols; k++) {
        tmp = mat[row_start + k];
        mat[row_start + k] = mat[rs2 + k];
        mat[rs2 + k] = tmp;
      }
    }
    tmp = mat[row_start + i];
    if (tmp != 1) {
      inverse = galois_single_divide(1, tmp, w);
      for (j = 0; j < cols; j++) {
        mat[row_start + j] = galois_single_multiply(mat[row_start + j], inverse, w);
      }
    }
    k = row_start + i;
    for (j = i + 1; j != cols; j++) {
      k += cols;
      if (mat[k] != 0) {
        if (mat[k] == 1) {
          rs2 = cols * j;
          for (x = 0; x < cols; x++) {
            mat[rs2 + x] ^= mat[row_start + x];
          }
        } else {
          tmp = mat[k];
          rs2 = cols * j;
          for (x = 0; x < cols; x++) {
            mat[rs2 + x] ^= galois_single_multiply(tmp, mat[row_start + x], w);
          }
        }
      }
    }
  }
  return 1;
}
```