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

```asm
<glmCrossProduct4d>:
  endbr64
  movsd  0x10(%rsi),%xmm4
  movsd  0x8(%rdx),%xmm3
  movsd  0x18(%rcx),%xmm11
  movsd  0x8(%rsi),%xmm2
  movsd  0x10(%rdx),%xmm7
  movsd  0x18(%rsi),%xmm8
  movsd  0x8(%rcx),%xmm9
  movsd  0x18(%rdx),%xmm6
  movsd  0x10(%rcx),%xmm5
  movq   D9(%rip),%xmm0
  movapd %xmm4,%xmm1
  xorpd  %xmm0,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  %xmm11,%xmm1
  movapd %xmm2,%xmm10
  mulsd  %xmm7,%xmm10
  mulsd  %xmm11,%xmm10
  addsd  %xmm10,%xmm1
  mulsd  %xmm8,%xmm7
  mulsd  %xmm9,%xmm7
  subsd  %xmm7,%xmm1
  mulsd  %xmm6,%xmm4
  mulsd  %xmm9,%xmm4
  addsd  %xmm4,%xmm1
  mulsd  %xmm8,%xmm3
  mulsd  %xmm5,%xmm3
  addsd  %xmm3,%xmm1
  mulsd  %xmm6,%xmm2
  mulsd  %xmm5,%xmm2
  subsd  %xmm2,%xmm1
  movsd  %xmm1,(%rdi)
  movsd  0x10(%rsi),%xmm5
  movsd  (%rdx),%xmm4
  movsd  0x18(%rcx),%xmm11
  movsd  (%rsi),%xmm3
  movsd  0x10(%rdx),%xmm1
  movsd  0x18(%rsi),%xmm8
  movsd  (%rcx),%xmm9
  movsd  0x18(%rdx),%xmm7
  movsd  0x10(%rcx),%xmm6
  movapd %xmm5,%xmm2
  mulsd  %xmm4,%xmm2
  mulsd  %xmm11,%xmm2
  movapd %xmm3,%xmm10
  mulsd  %xmm1,%xmm10
  mulsd  %xmm11,%xmm10
  subsd  %xmm10,%xmm2
  mulsd  %xmm8,%xmm1
  mulsd  %xmm9,%xmm1
  addsd  %xmm2,%xmm1
  mulsd  %xmm7,%xmm5
  mulsd  %xmm9,%xmm5
  subsd  %xmm5,%xmm1
  mulsd  %xmm8,%xmm4
  mulsd  %xmm6,%xmm4
  subsd  %xmm4,%xmm1
  mulsd  %xmm7,%xmm3
  mulsd  %xmm6,%xmm3
  addsd  %xmm3,%xmm1
  movsd  %xmm1,0x8(%rdi)
  movsd  0x8(%rsi),%xmm3
  movsd  (%rdx),%xmm2
  movsd  0x18(%rcx),%xmm10
  movsd  (%rsi),%xmm1
  movsd  0x8(%rdx),%xmm6
  movsd  0x18(%rsi),%xmm7
  movsd  (%rcx),%xmm8
  movsd  0x18(%rdx),%xmm5
  movsd  0x8(%rcx),%xmm4
  xorpd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  %xmm10,%xmm0
  movapd %xmm1,%xmm9
  mulsd  %xmm6,%xmm9
  mulsd  %xmm10,%xmm9
  addsd  %xmm9,%xmm0
  mulsd  %xmm7,%xmm6
  mulsd  %xmm8,%xmm6
  subsd  %xmm6,%xmm0
  mulsd  %xmm5,%xmm3
  mulsd  %xmm8,%xmm3
  addsd  %xmm3,%xmm0
  mulsd  %xmm7,%xmm2
  mulsd  %xmm4,%xmm2
  addsd  %xmm2,%xmm0
  mulsd  %xmm5,%xmm1
  mulsd  %xmm4,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%rdi)
  movsd  0x10(%rsi),%xmm4
  movsd  0x8(%rdx),%xmm1
  movsd  (%rcx),%xmm10
  movsd  0x8(%rsi),%xmm2
  movsd  0x10(%rdx),%xmm3
  movsd  (%rdx),%xmm8
  movsd  0x8(%rcx),%xmm9
  movsd  (%rsi),%xmm7
  movsd  0x10(%rcx),%xmm6
  movapd %xmm4,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm10,%xmm0
  movapd %xmm2,%xmm5
  mulsd  %xmm3,%xmm5
  mulsd  %xmm10,%xmm5
  subsd  %xmm5,%xmm0
  mulsd  %xmm8,%xmm4
  mulsd  %xmm9,%xmm4
  subsd  %xmm4,%xmm0
  mulsd  %xmm7,%xmm3
  mulsd  %xmm9,%xmm3
  addsd  %xmm3,%xmm0
  mulsd  %xmm8,%xmm2
  mulsd  %xmm6,%xmm2
  addsd  %xmm2,%xmm0
  mulsd  %xmm7,%xmm1
  mulsd  %xmm6,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%rdi)
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D9
type_infered: byte[8]

tool: {"D9": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

assistant: ```c
void glmCrossProduct4d(double result[4], const double v1[4], const double v2[4], const double v3[4]) {
  result[0] = -v1[2] * v2[1] * v3[3] + v1[1] * v2[2] * v3[3] - v1[3] * v2[2] * v3[1] + v1[2] * v2[3] * v3[1] + v1[3] * v2[1] * v3[2] - v1[1] * v2[3] * v3[2];
  result[1] = v1[2] * v2[0] * v3[3] - v1[0] * v2[2] * v3[3] + v1[3] * v2[2] * v3[0] - v1[2] * v2[3] * v3[0] - v1[3] * v2[0] * v3[2] + v1[0] * v2[3] * v3[2];
  result[2] = -v1[1] * v2[0] * v3[3] + v1[0] * v2[1] * v3[3] - v1[3] * v2[1] * v3[0] + v1[1] * v2[3] * v3[0] + v1[3] * v2[0] * v3[1] - v1[0] * v2[3] * v3[1];
  result[3] = v1[2] * v2[1] * v3[0] - v1[1] * v2[2] * v3[0] - v1[2] * v2[0] * v3[1] + v1[0] * v2[2] * v3[1] + v1[1] * v2[0] * v3[2] - v1[0] * v2[1] * v3[2];
}
```