user: 请显示与下面列出的汇编代码相匹配的C语言实现:

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