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

```asm
<randlc>:
  endbr64
  movsd  D16(%rip),%xmm4
  movapd %xmm0,%xmm2
  pxor   %xmm1,%xmm1
  movsd  D12(%rip),%xmm5
  movsd  (%rdi),%xmm3
  mulsd  %xmm4,%xmm0
  cvttsd2si %xmm0,%eax
  cvtsi2sd %eax,%xmm1
  movapd %xmm1,%xmm0
  mulsd  %xmm5,%xmm0
  subsd  %xmm0,%xmm2
  movapd %xmm3,%xmm0
  mulsd  %xmm4,%xmm0
  cvttsd2si %xmm0,%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movapd %xmm0,%xmm6
  mulsd  %xmm5,%xmm6
  mulsd  %xmm2,%xmm0
  subsd  %xmm6,%xmm3
  mulsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm2
  addsd  %xmm0,%xmm1
  mulsd  %xmm1,%xmm4
  cvttsd2si %xmm4,%eax
  pxor   %xmm4,%xmm4
  cvtsi2sd %eax,%xmm4
  mulsd  %xmm5,%xmm4
  subsd  %xmm4,%xmm1
  movapd %xmm1,%xmm0
  mulsd  %xmm5,%xmm0
  addsd  %xmm2,%xmm0
  movsd  D30(%rip),%xmm2
  movapd %xmm0,%xmm1
  mulsd  %xmm2,%xmm1
  cvttsd2si %xmm1,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  D37(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdi)
  mulsd  %xmm2,%xmm0
  ret
```

tool_calls: name: parse_data
arguments: label: D16
type_guessed: byte[8]
name: parse_data
arguments: label: D12
type_guessed: byte[8]
name: parse_data
arguments: label: D30
type_guessed: byte[8]
name: parse_data
arguments: label: D37
type_guessed: byte[8]

tool: {"D16": [0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x3e]}

tool: {"D12": [0x00,0x00,0x00,0x00,0x00,0x00,0x60,0x41]}

tool: {"D30": [0x00,0x00,0x00,0x00,0x00,0x00,0x10,0x3d]}

tool: {"D37": [0x00,0x00,0x00,0x00,0x00,0x00,0xd0,0x42]}

assistant: ```c
double randlc(double *x, double a) {
  double t1, t2, t3, t4, a1, a2, x1, x2, z;
  t1 = (0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5) * a;
  a1 = (int)t1;
  a2 = a - (2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0) * a1;
  t1 = (0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5) * (*x);
  x1 = (int)t1;
  x2 = (*x) - (2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0) * x1;
  t1 = a1 * x2 + a2 * x1;
  t2 = (int)((0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5) * t1);
  z = t1 - (2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0) * t2;
  t3 = (2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0) * z + a2 * x2;
  t4 = (int)(((0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5) * (0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5)) * t3);
  (*x) = t3 - ((2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0) * (2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0 * 2.0)) * t4;
  return (((0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5) * (0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5 * 0.5)) * (*x));
}
```