user: 请将这段汇编代码片段转换成等效的C代码：

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

tool_calls: name: disassemble_data
arguments: data_tag: D20
data_type: byte[4]
name: disassemble_data
arguments: data_tag: D24
data_type: byte[4]

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

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

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);
}
```