user: What is the c source code of the assembly code below?

```asm
<mat_inverse_matrix_d3>:
  endbr64
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x60,%rsp
  mov    %rdi,%rbp
  mov    %rsi,%rbx
  movq   %xmm0,%r12
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  mov    %rsi,%rdi
  call   <mat_get_determinant_d3@plt>
  movsd  %xmm0,0x8(%rsp)
  call   <mat_Dabs@plt>
  movq   %r12,%xmm4
  comisd %xmm0,%xmm4
  ja     L3
  movsd  0x20(%rbx),%xmm0
  movsd  0x40(%rbx),%xmm3
  movsd  0x28(%rbx),%xmm5
  movsd  0x38(%rbx),%xmm6
  movapd %xmm0,%xmm1
  mulsd  %xmm3,%xmm1
  movapd %xmm5,%xmm2
  mulsd  %xmm6,%xmm2
  subsd  %xmm2,%xmm1
  movsd  0x8(%rsp),%xmm10
  divsd  %xmm10,%xmm1
  movsd  %xmm1,0x10(%rsp)
  movsd  0x30(%rbx),%xmm2
  movsd  0x18(%rbx),%xmm4
  movapd %xmm5,%xmm1
  mulsd  %xmm2,%xmm1
  movapd %xmm3,%xmm7
  mulsd  %xmm4,%xmm7
  subsd  %xmm7,%xmm1
  divsd  %xmm10,%xmm1
  movsd  %xmm1,0x28(%rsp)
  movapd %xmm6,%xmm1
  mulsd  %xmm4,%xmm1
  movapd %xmm0,%xmm7
  mulsd  %xmm2,%xmm7
  subsd  %xmm7,%xmm1
  divsd  %xmm10,%xmm1
  movsd  %xmm1,0x40(%rsp)
  movsd  0x10(%rbx),%xmm1
  movsd  0x8(%rbx),%xmm8
  movapd %xmm6,%xmm7
  mulsd  %xmm1,%xmm7
  movapd %xmm3,%xmm9
  mulsd  %xmm8,%xmm9
  subsd  %xmm9,%xmm7
  divsd  %xmm10,%xmm7
  movsd  %xmm7,0x18(%rsp)
  movsd  (%rbx),%xmm7
  mulsd  %xmm7,%xmm3
  movapd %xmm2,%xmm9
  mulsd  %xmm1,%xmm9
  subsd  %xmm9,%xmm3
  divsd  %xmm10,%xmm3
  movsd  %xmm3,0x30(%rsp)
  mulsd  %xmm8,%xmm2
  mulsd  %xmm7,%xmm6
  subsd  %xmm6,%xmm2
  divsd  %xmm10,%xmm2
  movsd  %xmm2,0x48(%rsp)
  movapd %xmm5,%xmm2
  mulsd  %xmm8,%xmm2
  movapd %xmm0,%xmm3
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  divsd  %xmm10,%xmm2
  movsd  %xmm2,0x20(%rsp)
  mulsd  %xmm4,%xmm1
  mulsd  %xmm7,%xmm5
  subsd  %xmm5,%xmm1
  divsd  %xmm10,%xmm1
  movsd  %xmm1,0x38(%rsp)
  mulsd  %xmm7,%xmm0
  mulsd  %xmm8,%xmm4
  subsd  %xmm4,%xmm0
  divsd  %xmm10,%xmm0
  movsd  %xmm0,0x50(%rsp)
  lea    0x10(%rsp),%rsi
  mov    %rbp,%rdi
  call   <mat_copy_matrix_d3@plt>
  mov    $0x1,%eax
L17:
  mov    0x58(%rsp),%rdx
  sub    %fs:0x28,%rdx
  jne    L7
  add    $0x60,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
L3:
  movsd  0x8(%rsp),%xmm0
  lea    D0(%rip),%rdi
  call   <warning_print@plt>
  mov    $0x0,%eax
  jmp    L17
L7:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_tag: D0
type_guessed: string

tool: D0: 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;
}
```