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

```asm
<esmQuaternionToMatrix>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x38(%rbp)
  movss  %xmm0,-0x3c(%rbp)
  movss  %xmm1,-0x40(%rbp)
  movss  %xmm2,-0x44(%rbp)
  movss  %xmm3,-0x48(%rbp)
  movss  -0x3c(%rbp),%xmm0
  mulss  %xmm0,%xmm0
  movss  %xmm0,-0x24(%rbp)
  movss  -0x40(%rbp),%xmm0
  mulss  %xmm0,%xmm0
  movss  %xmm0,-0x20(%rbp)
  movss  -0x44(%rbp),%xmm0
  mulss  %xmm0,%xmm0
  movss  %xmm0,-0x1c(%rbp)
  movss  -0x3c(%rbp),%xmm0
  mulss  -0x40(%rbp),%xmm0
  movss  %xmm0,-0x18(%rbp)
  movss  -0x3c(%rbp),%xmm0
  mulss  -0x44(%rbp),%xmm0
  movss  %xmm0,-0x14(%rbp)
  movss  -0x40(%rbp),%xmm0
  mulss  -0x44(%rbp),%xmm0
  movss  %xmm0,-0x10(%rbp)
  movss  -0x48(%rbp),%xmm0
  mulss  -0x3c(%rbp),%xmm0
  movss  %xmm0,-0xc(%rbp)
  movss  -0x48(%rbp),%xmm0
  mulss  -0x40(%rbp),%xmm0
  movss  %xmm0,-0x8(%rbp)
  movss  -0x48(%rbp),%xmm0
  mulss  -0x44(%rbp),%xmm0
  movss  %xmm0,-0x4(%rbp)
  movss  -0x20(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  movss  D19(%rip),%xmm0
  movaps %xmm0,%xmm1
  subss  %xmm2,%xmm1
  movss  -0x1c(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  subss  %xmm2,%xmm1
  movaps %xmm1,%xmm0
  mov    -0x38(%rbp),%rax
  movss  %xmm0,(%rax)
  movss  -0x18(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  -0x4(%rbp),%xmm0
  addss  %xmm0,%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x4,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x14(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  -0x8(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  subss  %xmm2,%xmm1
  movaps %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0xc,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x18(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  -0x4(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x10,%rax
  subss  %xmm2,%xmm1
  movaps %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x24(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  movss  D19(%rip),%xmm0
  movaps %xmm0,%xmm1
  subss  %xmm2,%xmm1
  movss  -0x1c(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x14,%rax
  subss  %xmm2,%xmm1
  movaps %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x10(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  -0xc(%rbp),%xmm0
  addss  %xmm0,%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x18,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x1c,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x14(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  -0x8(%rbp),%xmm0
  addss  %xmm0,%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x20,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x10(%rbp),%xmm0
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  movss  -0xc(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x24,%rax
  subss  %xmm2,%xmm1
  movaps %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x24(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  movss  D19(%rip),%xmm0
  movaps %xmm0,%xmm1
  subss  %xmm2,%xmm1
  movss  -0x20(%rbp),%xmm0
  movaps %xmm0,%xmm2
  addss  %xmm0,%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x28,%rax
  subss  %xmm2,%xmm1
  movaps %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  D19(%rip),%xmm0
  movss  %xmm0,(%rax)
  nop
  pop    %rbp
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D19
data_type: byte[4]

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

assistant: ```c
void esmQuaternionToMatrix(float *M, float x, float y, float z, float w) {
  float xsq;
  float ysq;
  float zsq;
  float xy;
  float xz;
  float yz;
  float wx;
  float wy;
  float wz;
  xsq = x * x;
  ysq = y * y;
  zsq = z * z;
  xy = x * y;
  xz = x * z;
  yz = y * z;
  wx = w * x;
  wy = w * y;
  wz = w * z;
  M[0] = 1 - (2 * ysq) - (2 * zsq);
  M[1] = (2 * xy) + (2 * wz);
  M[2] = (2 * xz) - (2 * wy);
  M[3] = 0;
  M[4] = (2 * xy) - (2 * wz);
  M[5] = 1 - (2 * xsq) - (2 * zsq);
  M[6] = (2 * yz) + (2 * wx);
  M[7] = 0;
  M[8] = (2 * xz) + (2 * wy);
  M[9] = (2 * yz) - (2 * wx);
  M[10] = 1 - (2 * xsq) - (2 * ysq);
  M[11] = 0;
  M[12] = 0;
  M[13] = 0;
  M[14] = 0;
  M[15] = 1;
}
```