user: Converte this assembly snippet into its equivalent C code:

```asm
<mat_inverse_matrix_d3>:
  endbr64
  push   %rbp
  mov    %rdi,%rbp
  mov    %rsi,%rdi
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x78,%rsp
  movsd  %xmm0,0x18(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x68(%rsp)
  xor    %eax,%eax
  call   <mat_get_determinant_d3@plt>
  movapd %xmm0,%xmm4
  movsd  %xmm0,0x10(%rsp)
  unpcklpd %xmm4,%xmm4
  movaps %xmm4,(%rsp)
  call   <mat_Dabs@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  movapd (%rsp),%xmm4
  comisd %xmm0,%xmm2
  ja     L21
  movsd  0x40(%rbx),%xmm9
  movsd  0x20(%rbx),%xmm5
  lea    0x20(%rsp),%rsi
  mov    %rbp,%rdi
  movsd  0x28(%rbx),%xmm11
  movupd 0x30(%rbx),%xmm6
  movupd (%rbx),%xmm0
  movsd  0x38(%rbx),%xmm8
  movapd %xmm5,%xmm7
  movapd %xmm9,%xmm2
  movsd  0x10(%rbx),%xmm13
  movapd %xmm6,%xmm12
  movapd %xmm11,%xmm14
  movsd  0x30(%rbx),%xmm10
  unpcklpd %xmm8,%xmm7
  unpcklpd %xmm9,%xmm14
  unpckhpd %xmm0,%xmm12
  movupd 0x10(%rbx),%xmm3
  unpcklpd %xmm13,%xmm2
  unpcklpd %xmm10,%xmm8
  unpcklpd %xmm0,%xmm6
  mulpd  %xmm14,%xmm12
  mulpd  %xmm7,%xmm2
  mulpd  %xmm7,%xmm6
  subpd  %xmm12,%xmm2
  movapd %xmm10,%xmm12
  unpcklpd %xmm12,%xmm12
  divpd  %xmm4,%xmm2
  movaps %xmm2,0x20(%rsp)
  movapd %xmm0,%xmm2
  unpckhpd %xmm12,%xmm2
  movapd %xmm11,%xmm12
  unpcklpd %xmm12,%xmm12
  mulpd  %xmm12,%xmm2
  movapd %xmm5,%xmm12
  unpcklpd %xmm9,%xmm12
  mulpd  %xmm3,%xmm12
  subpd  %xmm12,%xmm2
  movapd %xmm3,%xmm12
  movsd  %xmm9,%xmm12
  movapd %xmm3,%xmm9
  unpcklpd %xmm0,%xmm9
  divpd  %xmm4,%xmm2
  movaps %xmm2,0x30(%rsp)
  movapd %xmm0,%xmm2
  unpcklpd %xmm13,%xmm2
  mulpd  %xmm12,%xmm2
  movapd %xmm10,%xmm12
  unpcklpd %xmm11,%xmm12
  mulpd  %xmm12,%xmm9
  subpd  %xmm9,%xmm2
  divpd  %xmm4,%xmm2
  movaps %xmm2,0x40(%rsp)
  movapd %xmm3,%xmm2
  unpckhpd %xmm3,%xmm3
  unpckhpd %xmm0,%xmm2
  mulpd  %xmm8,%xmm2
  subpd  %xmm6,%xmm2
  divpd  %xmm4,%xmm2
  movaps %xmm2,0x50(%rsp)
  movapd %xmm0,%xmm2
  unpckhpd %xmm0,%xmm0
  mulsd  %xmm5,%xmm2
  mulsd  %xmm3,%xmm0
  subsd  %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  movsd  %xmm2,0x60(%rsp)
  call   <mat_copy_matrix_d3@plt>
  mov    $0x1,%eax
L16:
  mov    0x68(%rsp),%rdx
  sub    %fs:0x28,%rdx
  jne    L5
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  ret
  nop
L21:
  movapd %xmm1,%xmm0
  lea    D4(%rip),%rdi
  call   <warning_print@plt>
  xor    %eax,%eax
  jmp    L16
L5:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_tag: D4
type_guessed: string

tool: D4: spglib: No inverse matrix (det=%f)


assistant: ```c
int mat_inverse_matrix_d3(double m[3][3], double a[3][3], const double precision) {
  double det;
  double c[3][3];
  det = mat_get_determinant_d3(a);
  if (mat_Dabs(det) < precision) {
    warning_print("spglib: No inverse matrix (det=%f)\n", det);
    return 0;
  }
  c[0][0] = (a[1][1] * a[2][2] - a[1][2] * a[2][1]) / det;
  c[1][0] = (a[1][2] * a[2][0] - a[1][0] * a[2][2]) / det;
  c[2][0] = (a[1][0] * a[2][1] - a[1][1] * a[2][0]) / det;
  c[0][1] = (a[2][1] * a[0][2] - a[2][2] * a[0][1]) / det;
  c[1][1] = (a[2][2] * a[0][0] - a[2][0] * a[0][2]) / det;
  c[2][1] = (a[2][0] * a[0][1] - a[2][1] * a[0][0]) / det;
  c[0][2] = (a[0][1] * a[1][2] - a[0][2] * a[1][1]) / det;
  c[1][2] = (a[0][2] * a[1][0] - a[0][0] * a[1][2]) / det;
  c[2][2] = (a[0][0] * a[1][1] - a[0][1] * a[1][0]) / det;
  mat_copy_matrix_d3(m, c);
  return 1;
}
```