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

```asm
<asteroid>:
  endbr64
  push   %r15
  xor    %eax,%eax
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x18,%rsp
  call   <rand@plt>
  pxor   %xmm8,%xmm8
  movslq %eax,%rdx
  mov    %eax,%ecx
  imul   $0x68db8bad,%rdx,%rdx
  sar    $0x1f,%ecx
  sar    $0x2c,%rdx
  sub    %ecx,%edx
  imul   $0x2710,%edx,%edx
  sub    %edx,%eax
  movslq %eax,%rdx
  sar    $0x1f,%eax
  imul   $0x66666667,%rdx,%rdx
  sar    $0x22,%rdx
  sub    %eax,%edx
  xor    %eax,%eax
  cvtsi2sd %edx,%xmm8
  movsd  %xmm8,0x8(%rsp)
  call   <rand@plt>
  mov    %eax,%r13d
  xor    %eax,%eax
  call   <rand@plt>
  mov    %eax,%r12d
  xor    %eax,%eax
  call   <rand@plt>
  mov    %eax,%r14d
  xor    %eax,%eax
  call   <rand@plt>
  mov    %eax,%r15d
  xor    %eax,%eax
  call   <rand@plt>
  mov    %eax,0x4(%rsp)
  xor    %eax,%eax
  call   <rand@plt>
  mov    %eax,(%rsp)
  xor    %eax,%eax
  call   <rand@plt>
  mov    %eax,%ebx
  xor    %eax,%eax
  call   <rand@plt>
  mov    %eax,%ebp
  xor    %eax,%eax
  call   <rand@plt>
  movslq (%rsp),%rdx
  movslq %r14d,%r8
  mov    %r12d,%r9d
  pxor   %xmm0,%xmm0
  movsd  D20(%rip),%xmm1
  pxor   %xmm7,%xmm7
  pxor   %xmm6,%xmm6
  pxor   %xmm5,%xmm5
  pxor   %xmm4,%xmm4
  pxor   %xmm3,%xmm3
  mov    %rdx,%rdi
  imul   $0x51eb851f,%rdx,%rdx
  mov    %edi,%ecx
  pxor   %xmm2,%xmm2
  sar    $0x26,%rdx
  sar    $0x1f,%ecx
  imul   $0x68db8bad,%r8,%r8
  sub    %ecx,%edx
  sar    $0x1f,%r9d
  imul   $0xc8,%edx,%edx
  sar    $0x2c,%r8
  sub    %edx,%edi
  movslq 0x4(%rsp),%rdx
  sub    $0x64,%edi
  mov    %rdx,%rsi
  imul   $0x51eb851f,%rdx,%rdx
  mov    %esi,%ecx
  cvtsi2sd %edi,%xmm7
  sar    $0x1f,%ecx
  sar    $0x26,%rdx
  sub    %ecx,%edx
  mov    %r15d,%ecx
  imul   $0xc8,%edx,%edx
  sar    $0x1f,%ecx
  sub    %edx,%esi
  movslq %r15d,%rdx
  imul   $0x51eb851f,%rdx,%rdx
  sub    $0x64,%esi
  cvtsi2sd %esi,%xmm6
  sar    $0x26,%rdx
  sub    %ecx,%edx
  imul   $0xc8,%edx,%edx
  sub    %edx,%r15d
  mov    %r14d,%edx
  sar    $0x1f,%edx
  lea    -0x64(%r15),%ecx
  sub    %edx,%r8d
  cvtsi2sd %ecx,%xmm5
  imul   $0x2710,%r8d,%r8d
  sub    %r8d,%r14d
  movslq %r12d,%r8
  imul   $0x68db8bad,%r8,%r8
  lea    -0x1388(%r14),%edx
  cvtsi2sd %edx,%xmm4
  sar    $0x2c,%r8
  sub    %r9d,%r8d
  mov    %r13d,%r9d
  imul   $0x2710,%r8d,%r8d
  sub    %r8d,%r12d
  movslq %r13d,%r8
  sub    $0x1388,%r12d
  imul   $0x68db8bad,%r8,%r8
  sar    $0x1f,%r9d
  sub    $0x20,%rsp
  cvtsi2sd %r12d,%xmm3
  movsd  0x28(%rsp),%xmm8
  sar    $0x2c,%r8
  sub    %r9d,%r8d
  imul   $0x2710,%r8d,%r8d
  sub    %r8d,%r13d
  mov    %eax,%r8d
  sar    $0x1f,%r8d
  cvtsi2sd %r13d,%xmm2
  subsd  D5(%rip),%xmm2
  shr    $0x18,%r8d
  add    %r8d,%eax
  movzbl %al,%eax
  sub    %r8d,%eax
  cvtsi2sd %eax,%xmm0
  mov    %ebp,%eax
  sar    $0x1f,%eax
  shr    $0x18,%eax
  add    %eax,%ebp
  mulsd  %xmm1,%xmm0
  movzbl %bpl,%ebp
  sub    %eax,%ebp
  mov    %ebx,%eax
  sar    $0x1f,%eax
  shr    $0x18,%eax
  add    %eax,%ebx
  movsd  %xmm0,0x10(%rsp)
  pxor   %xmm0,%xmm0
  movzbl %bl,%ebx
  cvtsi2sd %ebp,%xmm0
  sub    %eax,%ebx
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rsp)
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm8,%xmm1
  divsd  D0(%rip),%xmm1
  movsd  %xmm0,(%rsp)
  movapd %xmm8,%xmm0
  call   <makeBody@plt>
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D20
datatype: byte[8]
name: extract_symbolic_data
arguments: tag: D5
datatype: byte[8]
name: extract_symbolic_data
arguments: tag: D0
datatype: f64

tool: {"D20": [0x00,0x00,0x00,0x00,0x00,0x00,0x70,0x3f]}

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

tool: D0: 10

assistant: ```c
void asteroid() {
  int factor = 10000;
  double m = rand() % factor / 10;
  double r = m / 10;
  double xpos = (double)(rand() % factor) - factor / 2;
  double ypos = (rand() % factor) - factor / 2;
  double zpos = (rand() % factor) - factor / 2;
  double i = (rand() % 200) - 100;
  double j = (rand() % 200) - 100;
  double k = (rand() % 200) - 100;
  double rc = (double)(rand() % 256) / 256.0;
  double g = (double)(rand() % 256) / 256.0;
  double b = (double)(rand() % 256) / 256.0;
  makeBody(m, r, xpos, ypos, zpos, i, j, k, rc, g, b);
}
```