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

```asm
<jerasure_invertible_matrix>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %rdi,0x20(%rsp)
  mov    %esi,0x18(%rsp)
  mov    %edx,0x40(%rsp)
  test   %esi,%esi
  jle    L112
  mov    %esi,%eax
  movslq %esi,%rdx
  mov    %rdx,0x58(%rsp)
  lea    0x4(,%rdx,4),%rsi
  mov    %rsi,0x60(%rsp)
  lea    -0x1(%rax),%esi
  mov    %rsi,0x28(%rsp)
  add    $0x2,%rsi
  mov    %rsi,0x68(%rsp)
  mov    %eax,%eax
  mov    %rax,0x50(%rsp)
  lea    0x0(,%rdx,4),%rax
  mov    %rax,0x8(%rsp)
  mov    %rdi,0x48(%rsp)
  movq   $0x0,0x10(%rsp)
  movl   $0x0,0x44(%rsp)
  movq   $0x1,0x38(%rsp)
  lea    0x4(%rdi),%rax
  mov    %rax,0x30(%rsp)
  jmp    L18
L206:
  mov    0x38(%rsp),%rcx
  mov    %ecx,%eax
  mov    0x18(%rsp),%r9d
  cmp    %ecx,%r9d
  jle    L46
  mov    0x44(%rsp),%edx
  add    %r9d,%edx
  movslq %edx,%rdx
  lea    -0x1(%rcx,%rdx,1),%rdx
  mov    0x20(%rsp),%rcx
  lea    (%rcx,%rdx,4),%rdx
  mov    0x8(%rsp),%r8
L24:
  mov    (%rdx),%ecx
  test   %ecx,%ecx
  jne    L46
  add    $0x1,%eax
  add    %r8,%rdx
  cmp    %eax,%r9d
  jne    L24
L194:
  mov    %ecx,%eax
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L46:
  mov    0x18(%rsp),%edx
  cmp    %eax,%edx
  je     L204
  imul   %edx,%eax
  mov    %eax,%edx
  mov    0x20(%rsp),%rax
  mov    0x10(%rsp),%rdi
  lea    (%rax,%rdi,4),%rax
  mov    %rdi,%rcx
  add    0x28(%rsp),%rcx
  mov    0x30(%rsp),%rbx
  lea    (%rbx,%rcx,4),%r8
  movslq %edx,%rdx
  sub    %rdi,%rdx
L111:
  mov    (%rax),%ecx
  mov    (%rax,%rdx,4),%edi
  mov    %edi,(%rax)
  mov    %ecx,(%rax,%rdx,4)
  add    $0x4,%rax
  cmp    %r8,%rax
  jne    L111
  jmp    L167
L53:
  mov    0x40(%rsp),%r13d
  mov    %r13d,%edx
  mov    $0x1,%edi
  call   <galois_single_divide@plt>
  mov    %eax,%ebp
  mov    0x20(%rsp),%rdi
  mov    0x10(%rsp),%rax
  lea    (%rdi,%rax,4),%rbx
  add    0x28(%rsp),%rax
  mov    0x30(%rsp),%rdi
  lea    (%rdi,%rax,4),%r12
L186:
  mov    %r13d,%edx
  mov    %ebp,%esi
  mov    (%rbx),%edi
  call   <galois_single_multiply@plt>
  mov    %eax,(%rbx)
  add    $0x4,%rbx
  cmp    %r12,%rbx
  jne    L186
  jmp    L20
L50:
  movslq 0x4(%rsp),%rdx
  mov    0x20(%rsp),%rax
  lea    (%rax,%rdx,4),%rax
  mov    0x28(%rsp),%rdi
  lea    (%rdx,%rdi,1),%rcx
  mov    0x30(%rsp),%rdi
  lea    (%rdi,%rcx,4),%rsi
  mov    0x10(%rsp),%rcx
  sub    %rdx,%rcx
L56:
  mov    (%rax,%rcx,4),%edx
  xor    %edx,(%rax)
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L56
L140:
  add    $0x1,%r12d
  add    0x8(%rsp),%r13
  mov    0x18(%rsp),%ecx
  add    %ecx,0x4(%rsp)
  cmp    %r12d,%ecx
  je     L51
L201:
  mov    0x0(%r13),%ebp
  test   %ebp,%ebp
  je     L140
  cmp    $0x1,%ebp
  je     L50
  movslq 0x4(%rsp),%rax
  mov    0x20(%rsp),%rdi
  lea    (%rdi,%rax,4),%rbx
  mov    0x28(%rsp),%rdi
  lea    (%rax,%rdi,1),%rdx
  mov    0x30(%rsp),%rdi
  lea    (%rdi,%rdx,4),%r15
  mov    0x10(%rsp),%r14
  sub    %rax,%r14
  mov    %r12d,0x1c(%rsp)
  mov    0x40(%rsp),%r12d
L70:
  mov    (%rbx,%r14,4),%esi
  mov    %r12d,%edx
  mov    %ebp,%edi
  call   <galois_single_multiply@plt>
  xor    %eax,(%rbx)
  add    $0x4,%rbx
  cmp    %r15,%rbx
  jne    L70
  mov    0x1c(%rsp),%r12d
  jmp    L140
L51:
  addq   $0x1,0x38(%rsp)
  mov    0x38(%rsp),%rax
  mov    0x60(%rsp),%rdx
  add    %rdx,0x48(%rsp)
  mov    0x58(%rsp),%rdx
  add    %rdx,0x10(%rsp)
  cmp    0x68(%rsp),%rax
  je     L153
L18:
  mov    0x48(%rsp),%rax
  mov    %rax,%rsi
  mov    (%rax),%edi
  test   %edi,%edi
  je     L206
L167:
  mov    (%rsi),%esi
  cmp    $0x1,%esi
  jne    L53
L20:
  mov    0x38(%rsp),%rax
  mov    %eax,%r12d
  mov    0x50(%rsp),%rdi
  cmp    %rdi,%rax
  je     L232
  mov    0x8(%rsp),%r13
  add    0x48(%rsp),%r13
  mov    0x18(%rsp),%edi
  add    %edi,0x44(%rsp)
  mov    0x44(%rsp),%eax
  mov    %eax,0x4(%rsp)
  jmp    L201
L153:
  mov    $0x1,%ecx
  jmp    L194
L112:
  mov    $0x1,%ecx
  jmp    L194
L204:
  mov    %edi,%ecx
  jmp    L194
L232:
  mov    $0x1,%ecx
  jmp    L194
```

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