user: Please show me the C implementation that matches the assembly code listed below:

```asm
<inv_mix_columns>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  movl   $0x0,-0x38(%rbp)
  jmp    L17
L31:
  movl   $0x0,-0x34(%rbp)
  jmp    L11
L16:
  mov    -0x38(%rbp),%eax
  lea    0x0(,%rax,4),%edx
  mov    -0x34(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x34(%rbp),%eax
  cltq
  mov    %edx,-0x30(%rbp,%rax,4)
  addl   $0x1,-0x34(%rbp)
L11:
  cmpl   $0x3,-0x34(%rbp)
  jle    L16
  mov    -0x30(%rbp),%eax
  mov    $0xe,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%ebx
  mov    -0x24(%rbp),%eax
  mov    $0x9,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x28(%rbp),%eax
  mov    $0xd,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x2c(%rbp),%eax
  mov    $0xb,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%edx
  mov    -0x38(%rbp),%eax
  shl    $0x2,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  xor    %ebx,%edx
  mov    %edx,(%rax)
  mov    -0x2c(%rbp),%eax
  mov    $0xe,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%ebx
  mov    -0x30(%rbp),%eax
  mov    $0x9,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x24(%rbp),%eax
  mov    $0xd,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x28(%rbp),%eax
  mov    $0xb,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%edx
  mov    -0x38(%rbp),%eax
  shl    $0x2,%eax
  add    $0x1,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  xor    %ebx,%edx
  mov    %edx,(%rax)
  mov    -0x28(%rbp),%eax
  mov    $0xe,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%ebx
  mov    -0x2c(%rbp),%eax
  mov    $0x9,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x30(%rbp),%eax
  mov    $0xd,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x24(%rbp),%eax
  mov    $0xb,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%edx
  mov    -0x38(%rbp),%eax
  shl    $0x2,%eax
  add    $0x2,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  xor    %ebx,%edx
  mov    %edx,(%rax)
  mov    -0x24(%rbp),%eax
  mov    $0xe,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%ebx
  mov    -0x28(%rbp),%eax
  mov    $0x9,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x2c(%rbp),%eax
  mov    $0xd,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  xor    %eax,%ebx
  mov    -0x30(%rbp),%eax
  mov    $0xb,%esi
  mov    %eax,%edi
  call   <gmul@plt>
  mov    %eax,%edx
  mov    -0x38(%rbp),%eax
  shl    $0x2,%eax
  add    $0x3,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  xor    %ebx,%edx
  mov    %edx,(%rax)
  addl   $0x1,-0x38(%rbp)
L17:
  cmpl   $0x3,-0x38(%rbp)
  jle    L31
  nop
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L4
  call   <__stack_chk_fail@plt>
L4:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
void inv_mix_columns(uint8_t *state) {
  uint8_t a[4];
  for (uint8_t k = 0; k < 4; k++) {
    for (uint8_t i = 0; i < 4; i++) {
      a[i] = state[i + 4 * k];
    }
    state[4 * k] = gmul(a[0], 14) ^ gmul(a[3], 9) ^ gmul(a[2], 13) ^ gmul(a[1], 11);
    state[1 + 4 * k] = gmul(a[1], 14) ^ gmul(a[0], 9) ^ gmul(a[3], 13) ^ gmul(a[2], 11);
    state[2 + 4 * k] = gmul(a[2], 14) ^ gmul(a[1], 9) ^ gmul(a[0], 13) ^ gmul(a[3], 11);
    state[3 + 4 * k] = gmul(a[3], 14) ^ gmul(a[2], 9) ^ gmul(a[1], 13) ^ gmul(a[0], 11);
  }
}
```