user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<getSTcoordinates>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xd0,%rsp
  mov    %rdi,-0x88(%rbp)
  mov    %rsi,-0x90(%rbp)
  mov    %rdx,-0x98(%rbp)
  mov    %rcx,-0xa0(%rbp)
  mov    %r8,-0xa8(%rbp)
  mov    %r9,-0xb0(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xb8(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xc0(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0xc8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0x98(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0xa0(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x98(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0xa0(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  lea    -0x30(%rbp),%rdx
  lea    -0x50(%rbp),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <mat22Invert@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x78(%rbp)
  mov    -0x88(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0x88(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  lea    -0x60(%rbp),%rdx
  lea    -0x70(%rbp),%rcx
  lea    -0x30(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <mat221Multiply@plt>
  pxor   %xmm0,%xmm0
  ucomisd -0x78(%rbp),%xmm0
  jp     L31
  pxor   %xmm0,%xmm0
  ucomisd -0x78(%rbp),%xmm0
  je     L22
L31:
  mov    -0xb0(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x60(%rbp),%xmm2
  mov    -0xb8(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0xb0(%rbp),%rax
  movsd  (%rax),%xmm3
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x58(%rbp),%xmm2
  mov    -0xc0(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0xb0(%rbp),%rax
  movsd  (%rax),%xmm3
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mov    -0xc8(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0xb0(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  movsd  -0x60(%rbp),%xmm2
  mov    -0xb8(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0xb0(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm3
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x58(%rbp),%xmm2
  mov    -0xc0(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0xb0(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm3
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  mov    -0xc8(%rbp),%rax
  add    $0x8,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L17
L22:
  lea    D1(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
L17:
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L3
  call   <__stack_chk_fail@plt>
L3:
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D1
data_type: string

tool: D1: The matrix doesn't have an Inverse, something is wrong here


assistant: ```c
void getSTcoordinates(double x[2], double a[2], double b[2], double c[2], double rgb[3], double alpha[2], double beta[2], double gamma[2], double coordinates[2]) {
  double m[2][2] = {{b[0] - a[0], c[0] - a[0]}, {b[1] - a[1], c[1] - a[1]}};
  double Inv[2][2] = {{0.0, 0.0}, {0.0, 0.0}};
  double det = mat22Invert(m, Inv);
  double vec[2] = {x[0] - a[0], x[1] - a[1]};
  double pq[2] = {0.0, 0.0};
  mat221Multiply(Inv, vec, pq);
  if (det != 0) {
    coordinates[0] = alpha[0] + pq[0] * (beta[0] - alpha[0]) + pq[1] * (gamma[0] - alpha[0]);
    coordinates[1] = alpha[1] + pq[0] * (beta[1] - alpha[1]) + pq[1] * (gamma[1] - alpha[1]);
  } else {
    printf("The matrix doesn't have an Inverse, something is wrong here\n");
  }
}
```