user: Translate the below assembly code into C programming language:

```asm
<lib3ds_matrix_rotate_quat>:
  endbr64
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x68,%rsp
  movss  (%rsi),%xmm0
  movss  0x4(%rsi),%xmm3
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  movss  0x8(%rsi),%xmm2
  movss  0xc(%rsi),%xmm1
  mulss  %xmm0,%xmm0
  mulss  %xmm3,%xmm3
  mulss  %xmm2,%xmm2
  mulss  %xmm1,%xmm1
  addss  %xmm3,%xmm0
  addss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,0xc(%rsp)
  call   <fabs@plt>
  test   %eax,%eax
  jle    L22
  movss  D7(%rip),%xmm1
  movss  0xc(%rsp),%xmm0
  movss  D9(%rip),%xmm2
  divss  %xmm0,%xmm1
L2:
  movss  (%rbx),%xmm3
  movss  0x4(%rbx),%xmm5
  mov    %rbp,%rsi
  mov    %rbp,%rdi
  movss  0x8(%rbx),%xmm8
  movss  0xc(%rbx),%xmm4
  movabs $0x3f80000000000000,%rax
  lea    0x10(%rsp),%rdx
  movaps %xmm3,%xmm7
  movaps %xmm5,%xmm6
  movaps %xmm3,%xmm10
  mov    %rax,0x48(%rsp)
  mulss  %xmm1,%xmm6
  movaps %xmm4,%xmm9
  movl   $0x0,0x2c(%rsp)
  mulss  %xmm1,%xmm7
  movl   $0x0,0x1c(%rsp)
  mulss  %xmm8,%xmm1
  movl   $0x0,0x3c(%rsp)
  movq   $0x0,0x40(%rsp)
  mulss  %xmm6,%xmm9
  mulss  %xmm6,%xmm10
  mulss  %xmm5,%xmm6
  movaps %xmm1,%xmm0
  movaps %xmm4,%xmm1
  mulss  %xmm7,%xmm1
  mulss  %xmm0,%xmm4
  mulss  %xmm3,%xmm7
  mulss  %xmm0,%xmm5
  mulss  %xmm0,%xmm3
  movaps %xmm6,%xmm11
  mulss  %xmm8,%xmm0
  movaps %xmm2,%xmm8
  addss  %xmm0,%xmm11
  addss  %xmm7,%xmm0
  addss  %xmm6,%xmm7
  subss  %xmm11,%xmm8
  movss  %xmm8,0x10(%rsp)
  movaps %xmm10,%xmm8
  subss  %xmm4,%xmm8
  addss  %xmm10,%xmm4
  movss  %xmm4,0x14(%rsp)
  movaps %xmm2,%xmm4
  subss  %xmm7,%xmm2
  subss  %xmm0,%xmm4
  movss  %xmm8,0x20(%rsp)
  movaps %xmm5,%xmm0
  movaps %xmm9,%xmm8
  addss  %xmm3,%xmm8
  subss  %xmm1,%xmm0
  subss  %xmm9,%xmm3
  addss  %xmm5,%xmm1
  movss  %xmm2,0x38(%rsp)
  movss  %xmm4,0x24(%rsp)
  movss  %xmm8,0x30(%rsp)
  movss  %xmm0,0x34(%rsp)
  movss  %xmm3,0x18(%rsp)
  movss  %xmm1,0x28(%rsp)
  call   <lib3ds_matrix_mult@plt>
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L13
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  ret
  nopl   0x0(%rax)
L22:
  movss  D9(%rip),%xmm2
  movaps %xmm2,%xmm1
  jmp    L2
L13:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_tag: D7
type_guessed: byte[4]
name: read_data
arguments: data_tag: D9
type_guessed: byte[4]

tool: {"D7": [0x00,0x00,0x00,0x40]}

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

assistant: ```c
void lib3ds_matrix_rotate_quat(float m[4][4], float q[4]) {
  float s, xs, ys, zs, wx, wy, wz, xx, xy, xz, yy, yz, zz, l;
  float R[4][4];
  l = q[0] * q[0] + q[1] * q[1] + q[2] * q[2] + q[3] * q[3];
  if (fabs(l) < (1e-5)) {
    s = 1.0f;
  } else {
    s = 2.0f / l;
  }
  xs = q[0] * s;
  ys = q[1] * s;
  zs = q[2] * s;
  wx = q[3] * xs;
  wy = q[3] * ys;
  wz = q[3] * zs;
  xx = q[0] * xs;
  xy = q[0] * ys;
  xz = q[0] * zs;
  yy = q[1] * ys;
  yz = q[1] * zs;
  zz = q[2] * zs;
  R[0][0] = 1.0f - (yy + zz);
  R[1][0] = xy - wz;
  R[2][0] = xz + wy;
  R[0][1] = xy + wz;
  R[1][1] = 1.0f - (xx + zz);
  R[2][1] = yz - wx;
  R[0][2] = xz - wy;
  R[1][2] = yz + wx;
  R[2][2] = 1.0f - (xx + yy);
  R[3][0] = R[3][1] = R[3][2] = R[0][3] = R[1][3] = R[2][3] = 0.0f;
  R[3][3] = 1.0f;
  lib3ds_matrix_mult(m, m, R);
}
```