user: Converte this assembly snippet into its equivalent C code:

```asm
<_libtcas_filter_MitchellNetravali>:
  endbr64
  movapd %xmm0,%xmm3
  pxor   %xmm0,%xmm0
  movsd  D18(%rip),%xmm5
  movsd  D41(%rip),%xmm4
  comisd %xmm3,%xmm0
  movapd %xmm1,%xmm6
  movapd %xmm2,%xmm8
  mulsd  %xmm5,%xmm6
  mulsd  %xmm4,%xmm8
  jbe    L26
  xorpd  D51(%rip),%xmm3
L26:
  movsd  D78(%rip),%xmm7
  comisd %xmm3,%xmm7
  ja     L16
  movsd  D95(%rip),%xmm5
  comisd %xmm3,%xmm5
  jbe    L1
  movq   D51(%rip),%xmm7
  movapd %xmm1,%xmm0
  movapd %xmm1,%xmm5
  mulsd  %xmm4,%xmm5
  mulsd  D112(%rip),%xmm1
  xorpd  %xmm7,%xmm0
  xorpd  %xmm7,%xmm6
  subsd  %xmm8,%xmm0
  movsd  D72(%rip),%xmm8
  mulsd  %xmm2,%xmm8
  divsd  %xmm4,%xmm0
  addsd  %xmm8,%xmm5
  divsd  %xmm4,%xmm5
  mulsd  %xmm3,%xmm0
  addsd  %xmm5,%xmm0
  movsd  D87(%rip),%xmm5
  mulsd  %xmm2,%xmm5
  mulsd  D5(%rip),%xmm2
  mulsd  %xmm3,%xmm0
  subsd  %xmm5,%xmm6
  addsd  %xmm2,%xmm1
  divsd  %xmm4,%xmm6
  divsd  %xmm4,%xmm1
  addsd  %xmm6,%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm1,%xmm0
L1:
  ret
  nopl   0x0(%rax)
L16:
  movsd  D63(%rip),%xmm0
  subsd  D13(%rip),%xmm6
  mulsd  %xmm1,%xmm0
  addsd  %xmm1,%xmm1
  addsd  %xmm8,%xmm6
  divsd  %xmm4,%xmm6
  subsd  %xmm0,%xmm5
  movapd %xmm4,%xmm0
  subsd  %xmm1,%xmm0
  movapd %xmm5,%xmm2
  subsd  %xmm8,%xmm2
  divsd  %xmm4,%xmm2
  divsd  %xmm4,%xmm0
  mulsd  %xmm3,%xmm2
  addsd  %xmm6,%xmm2
  mulsd  %xmm3,%xmm2
  mulsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm0
  ret
```

tool_calls: name: read_data
arguments: tag: D18
data_type: byte[8]
name: read_data
arguments: tag: D41
data_type: f64
name: read_data
arguments: tag: D51
data_type: byte[16]
name: read_data
arguments: tag: D78
data_type: byte[8]
name: read_data
arguments: tag: D95
data_type: byte[8]
name: read_data
arguments: tag: D112
data_type: byte[8]
name: read_data
arguments: tag: D72
data_type: byte[8]
name: read_data
arguments: tag: D87
data_type: byte[8]
name: read_data
arguments: tag: D5
data_type: byte[8]
name: read_data
arguments: tag: D63
data_type: byte[8]
name: read_data
arguments: tag: D13
data_type: byte[8]

tool: {"D18": [0x00,0x00,0x00,0x00,0x00,0x00,0x28,0x40]}

tool: D41: 6

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

tool: {"D78": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

tool: {"D95": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x40]}

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

tool: {"D72": [0x00,0x00,0x00,0x00,0x00,0x00,0x3e,0x40]}

tool: {"D87": [0x00,0x00,0x00,0x00,0x00,0x00,0x48,0x40]}

tool: {"D5": [0x00,0x00,0x00,0x00,0x00,0x00,0x38,0x40]}

tool: {"D63": [0x00,0x00,0x00,0x00,0x00,0x00,0x22,0x40]}

tool: {"D13": [0x00,0x00,0x00,0x00,0x00,0x00,0x32,0x40]}

assistant: ```c
double _libtcas_filter_MitchellNetravali(double x, double b, double c) {
  double p0, p2, p3, q0, q1, q2, q3;
  p0 = (6 - 2 * b) / 6.0;
  p2 = (-18 + 12 * b + 6 * c) / 6.0;
  p3 = (12 - 9 * b - 6 * c) / 6.0;
  q0 = (8 * b + 24 * c) / 6.0;
  q1 = (-12 * b - 48 * c) / 6.0;
  q2 = (6 * b + 30 * c) / 6.0;
  q3 = (-b - 6 * c) / 6.0;
  if (x < 0)
    x = -x;
  if (x < 1)
    return p0 + (p2 + p3 * x) * x * x;
  else if (x < 2)
    return q0 + (q1 + (q2 + q3 * x) * x) * x;
  else
    return 0;
}
```