user: 请从给定的汇编指令推导出C源代码：

```asm
<math_matrix_rotate>:
  endbr64
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x68,%rsp
  lea    0x10(%rsp),%rbp
  movss  %xmm2,0x8(%rsp)
  movss  %xmm1,0xc(%rsp)
  mov    %rbp,%rsi
  movss  %xmm0,(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  call   <math_matrix_copy@plt>
  movss  (%rsp),%xmm4
  movaps %xmm4,%xmm0
  movss  %xmm4,0x4(%rsp)
  call   <cos@plt>
  movss  0x4(%rsp),%xmm4
  movss  %xmm0,(%rsp)
  movaps %xmm4,%xmm0
  call   <sin@plt>
  movss  (%rsp),%xmm3
  mov    %rbp,%rsi
  mov    %rbx,%rdi
  movss  0x20(%rsp),%xmm11
  movss  0x30(%rsp),%xmm7
  movss  0x24(%rsp),%xmm10
  movss  0x34(%rsp),%xmm6
  movaps %xmm11,%xmm4
  movaps %xmm7,%xmm5
  movss  0x28(%rsp),%xmm9
  mulss  %xmm3,%xmm4
  mulss  %xmm0,%xmm5
  mulss  %xmm0,%xmm11
  mulss  %xmm3,%xmm7
  addss  %xmm5,%xmm4
  movaps %xmm6,%xmm5
  mulss  %xmm0,%xmm5
  mulss  %xmm3,%xmm6
  subss  %xmm11,%xmm7
  movss  %xmm4,0x10(%rbx)
  movaps %xmm10,%xmm4
  mulss  %xmm3,%xmm4
  mulss  %xmm0,%xmm10
  movss  %xmm7,0x20(%rbx)
  addss  %xmm5,%xmm4
  movss  0x38(%rsp),%xmm5
  subss  %xmm10,%xmm6
  movaps %xmm5,%xmm8
  movss  %xmm4,0x14(%rbx)
  mulss  %xmm0,%xmm8
  movaps %xmm9,%xmm4
  mulss  %xmm3,%xmm4
  movss  %xmm6,0x24(%rbx)
  mulss  %xmm0,%xmm9
  mulss  %xmm3,%xmm5
  addss  %xmm8,%xmm4
  movss  0x2c(%rsp),%xmm8
  movaps %xmm8,%xmm12
  subss  %xmm9,%xmm5
  movss  %xmm4,0x18(%rbx)
  movss  0x3c(%rsp),%xmm4
  mulss  %xmm0,%xmm8
  mulss  %xmm3,%xmm12
  movaps %xmm4,%xmm13
  movss  %xmm5,0x28(%rbx)
  mulss  %xmm0,%xmm13
  mulss  %xmm3,%xmm4
  addss  %xmm13,%xmm12
  subss  %xmm8,%xmm4
  movss  %xmm12,0x1c(%rbx)
  movss  %xmm4,0x2c(%rbx)
  call   <math_matrix_copy@plt>
  movss  0xc(%rsp),%xmm1
  movaps %xmm1,%xmm0
  movss  %xmm1,0x4(%rsp)
  call   <cos@plt>
  movss  0x4(%rsp),%xmm1
  movss  %xmm0,(%rsp)
  movaps %xmm1,%xmm0
  call   <sin@plt>
  movss  (%rsp),%xmm3
  mov    %rbp,%rsi
  mov    %rbx,%rdi
  movss  0x10(%rsp),%xmm10
  movss  0x30(%rsp),%xmm6
  movss  0x14(%rsp),%xmm9
  movss  0x34(%rsp),%xmm5
  movaps %xmm10,%xmm1
  movaps %xmm6,%xmm4
  movss  0x18(%rsp),%xmm8
  mulss  %xmm3,%xmm1
  mulss  %xmm0,%xmm4
  mulss  %xmm0,%xmm10
  mulss  %xmm3,%xmm6
  subss  %xmm4,%xmm1
  movaps %xmm5,%xmm4
  mulss  %xmm0,%xmm4
  mulss  %xmm3,%xmm5
  addss  %xmm10,%xmm6
  movss  %xmm1,(%rbx)
  movaps %xmm9,%xmm1
  mulss  %xmm3,%xmm1
  mulss  %xmm0,%xmm9
  movss  %xmm6,0x20(%rbx)
  subss  %xmm4,%xmm1
  movss  0x38(%rsp),%xmm4
  addss  %xmm9,%xmm5
  movaps %xmm4,%xmm7
  movss  %xmm1,0x4(%rbx)
  mulss  %xmm0,%xmm7
  movaps %xmm8,%xmm1
  mulss  %xmm3,%xmm1
  movss  %xmm5,0x24(%rbx)
  mulss  %xmm0,%xmm8
  mulss  %xmm3,%xmm4
  subss  %xmm7,%xmm1
  movss  0x1c(%rsp),%xmm7
  movaps %xmm7,%xmm11
  addss  %xmm8,%xmm4
  movss  %xmm1,0x8(%rbx)
  movss  0x3c(%rsp),%xmm1
  mulss  %xmm0,%xmm7
  mulss  %xmm3,%xmm11
  movaps %xmm1,%xmm12
  movss  %xmm4,0x28(%rbx)
  mulss  %xmm0,%xmm12
  mulss  %xmm3,%xmm1
  subss  %xmm12,%xmm11
  addss  %xmm7,%xmm1
  movss  %xmm11,0xc(%rbx)
  movss  %xmm1,0x2c(%rbx)
  call   <math_matrix_copy@plt>
  movss  0x8(%rsp),%xmm2
  movaps %xmm2,%xmm0
  movss  %xmm2,0x4(%rsp)
  call   <cos@plt>
  movss  0x4(%rsp),%xmm2
  movss  %xmm0,(%rsp)
  movaps %xmm2,%xmm0
  call   <sin@plt>
  movss  0x20(%rsp),%xmm5
  movss  0x10(%rsp),%xmm9
  movss  (%rsp),%xmm1
  movss  0x14(%rsp),%xmm8
  movaps %xmm9,%xmm2
  movaps %xmm5,%xmm3
  movss  0x24(%rsp),%xmm4
  movss  0x18(%rsp),%xmm7
  mulss  %xmm0,%xmm3
  mulss  %xmm1,%xmm2
  mulss  %xmm1,%xmm5
  mulss  %xmm0,%xmm9
  addss  %xmm3,%xmm2
  movaps %xmm4,%xmm3
  mulss  %xmm0,%xmm3
  mulss  %xmm1,%xmm4
  subss  %xmm9,%xmm5
  movss  %xmm2,(%rbx)
  movaps %xmm8,%xmm2
  mulss  %xmm1,%xmm2
  mulss  %xmm0,%xmm8
  movss  %xmm5,0x10(%rbx)
  addss  %xmm3,%xmm2
  movss  0x28(%rsp),%xmm3
  subss  %xmm8,%xmm4
  movaps %xmm3,%xmm6
  movss  %xmm2,0x4(%rbx)
  mulss  %xmm0,%xmm6
  movaps %xmm7,%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm4,0x14(%rbx)
  mulss  %xmm1,%xmm3
  mulss  %xmm0,%xmm7
  addss  %xmm6,%xmm2
  movss  0x1c(%rsp),%xmm6
  movaps %xmm6,%xmm10
  subss  %xmm7,%xmm3
  movss  %xmm2,0x8(%rbx)
  movss  0x2c(%rsp),%xmm2
  mulss  %xmm1,%xmm10
  mulss  %xmm0,%xmm6
  movaps %xmm2,%xmm11
  movss  %xmm3,0x18(%rbx)
  mulss  %xmm0,%xmm11
  mulss  %xmm1,%xmm2
  addss  %xmm11,%xmm10
  subss  %xmm6,%xmm2
  movss  %xmm10,0xc(%rbx)
  movss  %xmm2,0x1c(%rbx)
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L3
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  ret
L3:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void math_matrix_rotate(float *matrix, const float x, const float y, const float z) {
  float cache[16];
  float cosine;
  float sine;
  math_matrix_copy(matrix, cache);
  cosine = cos(x);
  sine = sin(x);
  matrix[4] = cache[4] * cosine + cache[8] * sine;
  matrix[5] = cache[5] * cosine + cache[9] * sine;
  matrix[6] = cache[6] * cosine + cache[10] * sine;
  matrix[7] = cache[7] * cosine + cache[11] * sine;
  matrix[8] = cache[8] * cosine - cache[4] * sine;
  matrix[9] = cache[9] * cosine - cache[5] * sine;
  matrix[10] = cache[10] * cosine - cache[6] * sine;
  matrix[11] = cache[11] * cosine - cache[7] * sine;
  math_matrix_copy(matrix, cache);
  cosine = cos(y);
  sine = sin(y);
  matrix[0] = cache[0] * cosine - cache[8] * sine;
  matrix[1] = cache[1] * cosine - cache[9] * sine;
  matrix[2] = cache[2] * cosine - cache[10] * sine;
  matrix[3] = cache[3] * cosine - cache[11] * sine;
  matrix[8] = cache[8] * cosine + cache[0] * sine;
  matrix[9] = cache[9] * cosine + cache[1] * sine;
  matrix[10] = cache[10] * cosine + cache[2] * sine;
  matrix[11] = cache[11] * cosine + cache[3] * sine;
  math_matrix_copy(matrix, cache);
  cosine = cos(z);
  sine = sin(z);
  matrix[0] = cache[0] * cosine + cache[4] * sine;
  matrix[1] = cache[1] * cosine + cache[5] * sine;
  matrix[2] = cache[2] * cosine + cache[6] * sine;
  matrix[3] = cache[3] * cosine + cache[7] * sine;
  matrix[4] = cache[4] * cosine - cache[0] * sine;
  matrix[5] = cache[5] * cosine - cache[1] * sine;
  matrix[6] = cache[6] * cosine - cache[2] * sine;
  matrix[7] = cache[7] * cosine - cache[3] * sine;
}
```