user: Please show me the C implementation that matches the assembly code listed below:

```asm
<quat_multiply>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x8(%rbp)
  mov    %rsi,-0x10(%rbp)
  mov    %rdx,-0x18(%rbp)
  mov    -0x10(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm1
  mov    -0x18(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x10(%rbp),%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x10(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  mulss  %xmm2,%xmm1
  subss  %xmm1,%xmm0
  mov    -0x8(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x10(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x10(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm1
  mulss  %xmm2,%xmm1
  mov    -0x8(%rbp),%rax
  add    $0x4,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x10(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm1
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x10(%rbp),%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movss  (%rax),%xmm1
  mulss  %xmm2,%xmm1
  mov    -0x8(%rbp),%rax
  add    $0x8,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x10(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm1
  mov    -0x18(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  movss  (%rax),%xmm1
  mulss  %xmm2,%xmm1
  subss  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  mulss  %xmm2,%xmm1
  subss  %xmm1,%xmm0
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm1
  mulss  %xmm2,%xmm1
  mov    -0x8(%rbp),%rax
  add    $0xc,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x8(%rbp),%rax
  pop    %rbp
  ret
```

assistant: ```c
float *quat_multiply(float *result, float *a, float *b) {
  result[0] = a[3] * b[0] + a[0] * b[3] + a[1] * b[2] - a[2] * b[1];
  result[1] = a[3] * b[1] + a[1] * b[3] + a[2] * b[0] - a[0] * b[2];
  result[2] = a[3] * b[2] + a[2] * b[3] + a[0] * b[1] - a[1] * b[0];
  result[3] = a[3] * b[3] - a[0] * b[0] - a[1] * b[1] - a[2] * b[2];
  return result;
}
```