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

```asm
<gfks_rotate_vector_by_quaternion>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xa0,%rsp
  mov    %rdi,-0x88(%rbp)
  movss  %xmm0,-0x8c(%rbp)
  movss  %xmm1,-0x90(%rbp)
  movss  %xmm2,-0x94(%rbp)
  movss  %xmm3,-0x98(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movss  -0x8c(%rbp),%xmm0
  mulss  -0x90(%rbp),%xmm0
  addss  %xmm0,%xmm0
  movss  %xmm0,-0x80(%rbp)
  movss  -0x8c(%rbp),%xmm0
  mulss  -0x94(%rbp),%xmm0
  addss  %xmm0,%xmm0
  movss  %xmm0,-0x7c(%rbp)
  movss  -0x8c(%rbp),%xmm0
  mulss  -0x98(%rbp),%xmm0
  addss  %xmm0,%xmm0
  movss  %xmm0,-0x78(%rbp)
  movss  -0x90(%rbp),%xmm0
  mulss  -0x98(%rbp),%xmm0
  addss  %xmm0,%xmm0
  movss  %xmm0,-0x74(%rbp)
  movss  -0x90(%rbp),%xmm0
  mulss  -0x94(%rbp),%xmm0
  addss  %xmm0,%xmm0
  movss  %xmm0,-0x70(%rbp)
  movss  -0x94(%rbp),%xmm0
  mulss  -0x98(%rbp),%xmm0
  addss  %xmm0,%xmm0
  movss  %xmm0,-0x6c(%rbp)
  movss  -0x8c(%rbp),%xmm0
  addss  %xmm0,%xmm0
  cvttss2si %xmm0,%eax
  mov    $0x2,%esi
  mov    %eax,%edi
  call   <powf@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x68(%rbp)
  movss  -0x90(%rbp),%xmm0
  addss  %xmm0,%xmm0
  cvttss2si %xmm0,%eax
  mov    $0x2,%esi
  mov    %eax,%edi
  call   <powf@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x64(%rbp)
  movss  -0x94(%rbp),%xmm0
  addss  %xmm0,%xmm0
  cvttss2si %xmm0,%eax
  mov    $0x2,%esi
  mov    %eax,%edi
  call   <powf@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x60(%rbp)
  movss  D26(%rip),%xmm0
  subss  -0x64(%rbp),%xmm0
  subss  -0x60(%rbp),%xmm0
  movss  %xmm0,-0x50(%rbp)
  movss  -0x80(%rbp),%xmm0
  addss  -0x6c(%rbp),%xmm0
  movss  %xmm0,-0x4c(%rbp)
  movss  -0x7c(%rbp),%xmm0
  subss  -0x74(%rbp),%xmm0
  movss  %xmm0,-0x48(%rbp)
  movss  -0x80(%rbp),%xmm0
  subss  -0x6c(%rbp),%xmm0
  movss  %xmm0,-0x44(%rbp)
  movss  D26(%rip),%xmm0
  subss  -0x68(%rbp),%xmm0
  subss  -0x60(%rbp),%xmm0
  movss  %xmm0,-0x40(%rbp)
  movss  -0x70(%rbp),%xmm0
  addss  -0x78(%rbp),%xmm0
  movss  %xmm0,-0x3c(%rbp)
  movss  -0x7c(%rbp),%xmm0
  addss  -0x74(%rbp),%xmm0
  movss  %xmm0,-0x38(%rbp)
  movss  -0x70(%rbp),%xmm0
  subss  -0x78(%rbp),%xmm0
  movss  %xmm0,-0x34(%rbp)
  movss  D26(%rip),%xmm0
  subss  -0x68(%rbp),%xmm0
  subss  -0x64(%rbp),%xmm0
  movss  %xmm0,-0x30(%rbp)
  mov    -0x88(%rbp),%rax
  movss  (%rax),%xmm1
  movss  -0x50(%rbp),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  movss  -0x44(%rbp),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  movss  -0x38(%rbp),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x5c(%rbp)
  mov    -0x88(%rbp),%rax
  movss  (%rax),%xmm1
  movss  -0x4c(%rbp),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  movss  -0x40(%rbp),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  movss  -0x34(%rbp),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x58(%rbp)
  mov    -0x88(%rbp),%rax
  movss  (%rax),%xmm1
  movss  -0x48(%rbp),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  movss  -0x3c(%rbp),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x88(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  movss  -0x30(%rbp),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x54(%rbp)
  movss  -0x5c(%rbp),%xmm0
  mov    -0x88(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x88(%rbp),%rax
  add    $0x4,%rax
  movss  -0x58(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x88(%rbp),%rax
  add    $0x8,%rax
  movss  -0x54(%rbp),%xmm0
  movss  %xmm0,(%rax)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L9
  call   <__stack_chk_fail@plt>
L9:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D26
type: byte[4]

tool: {"D26": [0x00,0x00,0x80,0x3f]}

assistant: ```c
void gfks_rotate_vector_by_quaternion(float vec[3], float x, float y, float z, float w) {
  float xy2 = x * y * 2;
  float xz2 = x * z * 2;
  float xw2 = x * w * 2;
  float yw2 = y * w * 2;
  float yz2 = y * z * 2;
  float zw2 = z * w * 2;
  float _2x2 = powf(2 * x, 2);
  float _2y2 = powf(2 * y, 2);
  float _2z2 = powf(2 * z, 2);
  float Q[16];
  Q[0] = 1.0f - _2y2 - _2z2;
  Q[1] = xy2 + zw2;
  Q[2] = xz2 - yw2;
  Q[3] = xy2 - zw2;
  Q[4] = 1.0f - _2x2 - _2z2;
  Q[5] = yz2 + xw2;
  Q[6] = xz2 + yw2;
  Q[7] = yz2 - xw2;
  Q[8] = 1.0f - _2x2 - _2y2;
  float result[3];
  result[0] = (vec[0] * Q[0]) + (vec[1] * Q[3]) + (vec[2] * Q[6]);
  result[1] = (vec[0] * Q[1]) + (vec[1] * Q[4]) + (vec[2] * Q[7]);
  result[2] = (vec[0] * Q[2]) + (vec[1] * Q[5]) + (vec[2] * Q[8]);
  vec[0] = result[0];
  vec[1] = result[1];
  vec[2] = result[2];
}
```