user: 请从给定的汇编指令推导出C源代码：

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

tool_calls: name: retrieve_data
arguments: tag: D0
datatype: 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;
}
```