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

```asm
<interpolating>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  mov    %rdx,%rsi
  mov    %rcx,%rdx
  push   %r14
  push   %r13
  mov    %r9,%r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %r8,%rbx
  sub    $0x98,%rsp
  movsd  (%r15),%xmm1
  movsd  (%rsi),%xmm0
  mov    %rdi,0x18(%rsp)
  lea    0x60(%rsp),%r8
  lea    0x40(%rsp),%rdi
  mov    0xd0(%rsp),%r14
  subsd  %xmm1,%xmm0
  mov    0xd8(%rsp),%r12
  mov    0xe0(%rsp),%rbp
  mov    %fs:0x28,%rcx
  mov    %rcx,0x88(%rsp)
  xor    %ecx,%ecx
  mov    %r8,0x10(%rsp)
  movsd  %xmm0,0x40(%rsp)
  movsd  (%rdx),%xmm0
  movq   $0x0,0x60(%rsp)
  movq   $0x0,0x68(%rsp)
  subsd  %xmm1,%xmm0
  movsd  0x8(%r15),%xmm1
  movq   $0x0,0x70(%rsp)
  movq   $0x0,0x78(%rsp)
  movsd  %xmm0,0x48(%rsp)
  movsd  0x8(%rsi),%xmm0
  mov    %r8,%rsi
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x50(%rsp)
  movsd  0x8(%rdx),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x58(%rsp)
  call   <mat22Invert@plt>
  mov    0x18(%rsp),%rax
  mov    0x10(%rsp),%rdi
  lea    0x30(%rsp),%rdx
  lea    0x20(%rsp),%rsi
  movsd  %xmm0,0x8(%rsp)
  movq   $0x0,0x30(%rsp)
  movsd  (%rax),%xmm1
  movq   $0x0,0x38(%rsp)
  subsd  (%r15),%xmm1
  movsd  %xmm1,0x20(%rsp)
  movsd  0x8(%rax),%xmm1
  subsd  0x8(%r15),%xmm1
  movsd  %xmm1,0x28(%rsp)
  call   <mat221Multiply@plt>
  movsd  0x8(%rsp),%xmm0
  ucomisd D8(%rip),%xmm0
  jp     L3
  je     L21
L3:
  movsd  0x0(%r13),%xmm4
  movsd  (%r14),%xmm0
  movsd  0x30(%rsp),%xmm2
  movsd  (%r12),%xmm3
  subsd  %xmm4,%xmm0
  movsd  0x38(%rsp),%xmm1
  subsd  %xmm4,%xmm3
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm3
  addsd  %xmm4,%xmm0
  addsd  %xmm3,%xmm0
  mulsd  (%rbx),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movsd  0x8(%r13),%xmm4
  movsd  0x8(%r14),%xmm0
  movsd  0x8(%r12),%xmm3
  subsd  %xmm4,%xmm0
  subsd  %xmm4,%xmm3
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm3
  addsd  %xmm4,%xmm0
  addsd  %xmm3,%xmm0
  mulsd  0x8(%rbx),%xmm0
  movsd  %xmm0,0x8(%rbp)
  movsd  0x10(%r13),%xmm3
  movsd  0x10(%r14),%xmm0
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  movsd  0x10(%r12),%xmm2
  subsd  %xmm3,%xmm2
  mulsd  %xmm2,%xmm1
  addsd  %xmm3,%xmm0
  addsd  %xmm1,%xmm0
  mulsd  0x10(%rbx),%xmm0
  movsd  %xmm0,0x10(%rbp)
L22:
  mov    0x88(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L4
  add    $0x98,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L21:
  lea    D18(%rip),%rdi
  call   <printf@plt>
  jmp    L22
L4:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_tag: D8
data_type: f64
name: parse_data
arguments: data_tag: D18
data_type: string

tool: D8: 0

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


assistant: ```c
void interpolating(double x[2], double a[2], double b[2], double c[2], double rgb[3], double alpha[3], double beta[3], double gamma[3], double new_rgb[3]) {
  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) {
    new_rgb[0] = (alpha[0] + pq[0] * (beta[0] - alpha[0]) + pq[1] * (gamma[0] - alpha[0])) * rgb[0];
    new_rgb[1] = (alpha[1] + pq[0] * (beta[1] - alpha[1]) + pq[1] * (gamma[1] - alpha[1])) * rgb[1];
    new_rgb[2] = (alpha[2] + pq[0] * (beta[2] - alpha[2]) + pq[1] * (gamma[2] - alpha[2])) * rgb[2];
  } else {
    printf("The matrix doesn't have an Inverse, something is wrong here\n");
  }
}
```