user: What is the c source code of the assembly code below?

```asm
<mat4_rotation_quaternion>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    -0x40(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x24(%rbp)
  mov    -0x40(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x20(%rbp)
  mov    -0x40(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x1c(%rbp)
  mov    -0x40(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x18(%rbp)
  mov    -0x40(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x14(%rbp)
  mov    -0x40(%rbp),%rax
  add    $0x20,%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x10(%rbp)
  mov    -0x40(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0xc(%rbp)
  mov    -0x40(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x8(%rbp)
  mov    -0x40(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x40(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%rbp)
  movss  -0x20(%rbp),%xmm0
  subss  -0x1c(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  D29(%rip),%xmm0
  subss  %xmm1,%xmm0
  mov    -0x38(%rbp),%rax
  movss  %xmm0,(%rax)
  movss  -0x18(%rbp),%xmm0
  addss  -0x14(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x4,%rax
  addss  %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x10(%rbp),%xmm0
  subss  -0xc(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  addss  %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0xc,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x18(%rbp),%xmm0
  subss  -0x14(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x10,%rax
  addss  %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x24(%rbp),%xmm0
  subss  -0x1c(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x14,%rax
  movss  D29(%rip),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x8(%rbp),%xmm0
  addss  -0x4(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x18,%rax
  addss  %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x1c,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x10(%rbp),%xmm0
  addss  -0xc(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x20,%rax
  addss  %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x8(%rbp),%xmm0
  subss  -0x4(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x24,%rax
  addss  %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x24(%rbp),%xmm0
  subss  -0x20(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  mov    -0x38(%rbp),%rax
  add    $0x28,%rax
  movss  D29(%rip),%xmm0
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x2c,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x30,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x34,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x38,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x3c,%rax
  movss  D29(%rip),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  pop    %rbp
  ret
```

tool_calls: name: unpack_data
arguments: label: D29
type_infered: f32

tool: D29: 1

assistant: ```c
float *mat4_rotation_quaternion(float *result, float *q) {
  float xx = q[0] * q[0];
  float yy = q[1] * q[1];
  float zz = q[2] * q[2];
  float xy = q[0] * q[1];
  float zw = q[2] * q[3];
  float xz = q[8] * q[0];
  float yw = q[1] * q[3];
  float yz = q[1] * q[2];
  float xw = q[0] * q[3];
  result[0] = 1.0f - 2.0f * (yy - zz);
  result[1] = 2.0f * (xy + zw);
  result[2] = 2.0f * (xz - yw);
  result[3] = 0.0f;
  result[4] = 2.0f * (xy - zw);
  result[5] = 1.0f - 2.0f * (xx - zz);
  result[6] = 2.0f * (yz + xw);
  result[7] = 0.0f;
  result[8] = 2.0f * (xz + yw);
  result[9] = 2.0f * (yz - xw);
  result[10] = 1.0f - 2.0f * (xx - yy);
  result[11] = 0.0f;
  result[12] = 0.0f;
  result[13] = 0.0f;
  result[14] = 0.0f;
  result[15] = 1.0f;
  return result;
}
```