user: Translate the below assembly code into C programming language:

```asm
<babl_matrix_invert>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x68(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0x68(%rbp),%rcx
  lea    -0x50(%rbp),%rax
  mov    $0x48,%edx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x30(%rbp),%xmm2
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  movsd  -0x18(%rbp),%xmm3
  movsd  -0x28(%rbp),%xmm2
  mulsd  %xmm3,%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x48(%rbp),%xmm2
  movsd  -0x38(%rbp),%xmm3
  movsd  -0x10(%rbp),%xmm1
  mulsd  %xmm3,%xmm1
  movsd  -0x28(%rbp),%xmm4
  movsd  -0x20(%rbp),%xmm3
  mulsd  %xmm4,%xmm3
  subsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  -0x40(%rbp),%xmm2
  movsd  -0x38(%rbp),%xmm3
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm3,%xmm0
  movsd  -0x30(%rbp),%xmm4
  movsd  -0x20(%rbp),%xmm3
  mulsd  %xmm4,%xmm3
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  D7(%rip),%xmm0
  divsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x30(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm2
  movsd  -0x28(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  mov    -0x70(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x40(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x48(%rbp),%xmm2
  movsd  -0x10(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x8,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x48(%rbp),%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x40(%rbp),%xmm2
  movsd  -0x30(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x10,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x28(%rbp),%xmm1
  movsd  -0x20(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x38(%rbp),%xmm2
  movsd  -0x10(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x18,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x40(%rbp),%xmm2
  movsd  -0x20(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x20,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x38(%rbp),%xmm1
  movsd  -0x40(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x50(%rbp),%xmm2
  movsd  -0x28(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x28,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x38(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x20(%rbp),%xmm2
  movsd  -0x30(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x30,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x20(%rbp),%xmm1
  movsd  -0x48(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x50(%rbp),%xmm2
  movsd  -0x18(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x38,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x30(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x38(%rbp),%xmm2
  movsd  -0x48(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  add    $0x40,%rax
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L4
  call   <__stack_chk_fail@plt>
L4:
  leave
  ret
```

tool_calls: name: unpack_data
arguments: label: D7
data_type: f64

tool: D7: 1

assistant: ```c
void babl_matrix_invert(const double *in, double *out) {
  double mat[9];
  double det, invdet;
  memcpy(mat, in, sizeof(mat));
  det = mat[0 * 3 + 0] * (mat[1 * 3 + 1] * mat[2 * 3 + 2] - mat[2 * 3 + 1] * mat[1 * 3 + 2]) - mat[0 * 3 + 1] * (mat[1 * 3 + 0] * mat[2 * 3 + 2] - mat[1 * 3 + 2] * mat[2 * 3 + 0]) + mat[0 * 3 + 2] * (mat[1 * 3 + 0] * mat[2 * 3 + 1] - mat[1 * 3 + 1] * mat[2 * 3 + 0]);
  invdet = 1.0 / det;
  out[0 * 3 + 0] = (mat[1 * 3 + 1] * mat[2 * 3 + 2] - mat[2 * 3 + 1] * mat[1 * 3 + 2]) * invdet;
  out[0 * 3 + 1] = (mat[0 * 3 + 2] * mat[2 * 3 + 1] - mat[0 * 3 + 1] * mat[2 * 3 + 2]) * invdet;
  out[0 * 3 + 2] = (mat[0 * 3 + 1] * mat[1 * 3 + 2] - mat[0 * 3 + 2] * mat[1 * 3 + 1]) * invdet;
  out[1 * 3 + 0] = (mat[1 * 3 + 2] * mat[2 * 3 + 0] - mat[1 * 3 + 0] * mat[2 * 3 + 2]) * invdet;
  out[1 * 3 + 1] = (mat[0 * 3 + 0] * mat[2 * 3 + 2] - mat[0 * 3 + 2] * mat[2 * 3 + 0]) * invdet;
  out[1 * 3 + 2] = (mat[1 * 3 + 0] * mat[0 * 3 + 2] - mat[0 * 3 + 0] * mat[1 * 3 + 2]) * invdet;
  out[2 * 3 + 0] = (mat[1 * 3 + 0] * mat[2 * 3 + 1] - mat[2 * 3 + 0] * mat[1 * 3 + 1]) * invdet;
  out[2 * 3 + 1] = (mat[2 * 3 + 0] * mat[0 * 3 + 1] - mat[0 * 3 + 0] * mat[2 * 3 + 1]) * invdet;
  out[2 * 3 + 2] = (mat[0 * 3 + 0] * mat[1 * 3 + 1] - mat[1 * 3 + 0] * mat[0 * 3 + 1]) * invdet;
}
```