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

```asm
<combust>:
  endbr64
  movsd  (%rdi),%xmm1
  movsd  0x8(%rdi),%xmm0
  movsd  D16(%rip),%xmm2
  mulsd  %xmm1,%xmm0
  addsd  %xmm1,%xmm0
  movsd  D23(%rip),%xmm1
  mulsd  0x20(%rdi),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rsi)
  movsd  (%rdi),%xmm1
  movsd  0x8(%rdi),%xmm4
  movsd  0x10(%rdi),%xmm3
  movapd %xmm1,%xmm0
  addsd  %xmm1,%xmm0
  mulsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm0
  movsd  D47(%rip),%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm4,%xmm1
  addsd  %xmm1,%xmm0
  movapd %xmm4,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm1,%xmm0
  movapd %xmm4,%xmm1
  mulsd  %xmm2,%xmm1
  mulsd  %xmm3,%xmm1
  movapd %xmm4,%xmm3
  addsd  %xmm1,%xmm0
  movsd  D60(%rip),%xmm1
  mulsd  %xmm1,%xmm3
  mulsd  0x18(%rdi),%xmm3
  addsd  %xmm3,%xmm0
  movsd  D72(%rip),%xmm3
  mulsd  %xmm3,%xmm4
  addsd  %xmm4,%xmm0
  movsd  D103(%rip),%xmm4
  mulsd  0x20(%rdi),%xmm4
  subsd  %xmm4,%xmm0
  movsd  %xmm0,0x8(%rsi)
  movsd  0x8(%rdi),%xmm5
  movsd  0x10(%rdi),%xmm4
  movapd %xmm5,%xmm0
  addsd  %xmm5,%xmm0
  mulsd  %xmm2,%xmm5
  mulsd  %xmm4,%xmm0
  mulsd  %xmm4,%xmm5
  mulsd  %xmm4,%xmm0
  addsd  %xmm5,%xmm0
  movsd  D59(%rip),%xmm5
  mulsd  %xmm4,%xmm5
  mulsd  %xmm4,%xmm5
  addsd  %xmm5,%xmm0
  movsd  D112(%rip),%xmm5
  mulsd  %xmm5,%xmm4
  addsd  %xmm4,%xmm0
  movsd  D12(%rip),%xmm4
  mulsd  0x20(%rdi),%xmm4
  subsd  %xmm4,%xmm0
  movsd  %xmm0,0x10(%rsi)
  movsd  0x8(%rdi),%xmm0
  movsd  0x18(%rdi),%xmm4
  mulsd  %xmm1,%xmm0
  movapd %xmm4,%xmm6
  addsd  %xmm4,%xmm6
  mulsd  %xmm4,%xmm0
  mulsd  %xmm6,%xmm4
  addsd  %xmm4,%xmm0
  movsd  D78(%rip),%xmm4
  mulsd  0x20(%rdi),%xmm4
  subsd  %xmm4,%xmm0
  movsd  %xmm0,0x18(%rsi)
  movsd  (%rdi),%xmm8
  movsd  0x8(%rdi),%xmm4
  movsd  0x10(%rdi),%xmm7
  movsd  0x18(%rdi),%xmm6
  mulsd  %xmm4,%xmm2
  movapd %xmm8,%xmm0
  mulsd  %xmm4,%xmm0
  mulsd  %xmm4,%xmm1
  mulsd  %xmm7,%xmm2
  addsd  %xmm8,%xmm0
  movsd  D38(%rip),%xmm8
  mulsd  %xmm6,%xmm1
  mulsd  %xmm4,%xmm8
  mulsd  %xmm6,%xmm6
  mulsd  %xmm4,%xmm8
  addsd  %xmm8,%xmm0
  movapd %xmm4,%xmm8
  mulsd  %xmm7,%xmm8
  mulsd  %xmm3,%xmm4
  mulsd  %xmm7,%xmm8
  addsd  %xmm8,%xmm0
  addsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  D89(%rip),%xmm1
  mulsd  %xmm7,%xmm1
  addsd  %xmm4,%xmm0
  mulsd  %xmm7,%xmm1
  mulsd  %xmm5,%xmm7
  addsd  %xmm1,%xmm0
  addsd  %xmm7,%xmm0
  addsd  %xmm6,%xmm0
  subsd  D18(%rip),%xmm0
  movsd  %xmm0,0x20(%rsi)
  ret
```

tool_calls: name: disassemble_data
arguments: data_tag: D16
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D23
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D47
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D60
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D72
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D103
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D59
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D112
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D12
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D78
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D38
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D89
type_guessed: f64
name: disassemble_data
arguments: data_tag: D18
type_guessed: f64

tool: {"D16": [0xb8,0x17,0x59,0xb1,0x46,0xdd,0x41,0x3f]}

tool: {"D23": [0x00,0x00,0x00,0x00,0x00,0x00,0x08,0x40]}

tool: {"D47": [0x9a,0x33,0x29,0xa5,0x6a,0x32,0xc8,0x3e]}

tool: {"D60": [0xb4,0x9a,0x97,0x8b,0x44,0xdd,0x01,0x3f]}

tool: {"D72": [0x7b,0x3b,0x1a,0x6b,0xa5,0x2e,0x9e,0x3e]}

tool: {"D103": [0x00,0x00,0x00,0x00,0x00,0x00,0x24,0x40]}

tool: {"D59": [0x4e,0x62,0x10,0x58,0x39,0xb4,0xd8,0x3f]}

tool: {"D112": [0x0d,0xe2,0x9e,0x0f,0x18,0xe9,0x3a,0x3f]}

tool: {"D12": [0x00,0x00,0x00,0x00,0x00,0x00,0x20,0x40]}

tool: {"D78": [0x00,0x00,0x00,0x00,0x00,0x00,0x44,0x40]}

tool: {"D38": [0xbc,0x77,0x1b,0x6e,0x9c,0x21,0xb0,0x3e]}

tool: D89: 0.193

tool: D18: 1

assistant: ```c
void combust(double *p, double *x, int m, int n, void *data) {
  double R, R5, R6, R7, R8, R9, R10;
  R = 10;
  R5 = 0.193;
  R6 = 4.10622 * 1e-4;
  R7 = 5.45177 * 1e-4;
  R8 = 4.4975 * 1e-7;
  R9 = 3.40735 * 1e-5;
  R10 = 9.615 * 1e-7;
  x[0] = p[0] * p[1] + p[0] - 3 * p[4];
  x[1] = 2 * p[0] * p[1] + p[0] + 3 * R10 * p[1] * p[1] + p[1] * p[2] * p[2] + R7 * p[1] * p[2] + R9 * p[1] * p[3] + R8 * p[1] - R * p[4];
  x[2] = 2 * p[1] * p[2] * p[2] + R7 * p[1] * p[2] + 2 * R5 * p[2] * p[2] + R6 * p[2] - 8 * p[4];
  x[3] = R9 * p[1] * p[3] + 2 * p[3] * p[3] - 4 * R * p[4];
  x[4] = p[0] * p[1] + p[0] + R10 * p[1] * p[1] + p[1] * p[2] * p[2] + R7 * p[1] * p[2] + R9 * p[1] * p[3] + R8 * p[1] + R5 * p[2] * p[2] + R6 * p[2] + p[3] * p[3] - 1.0;
}
```