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

```asm
<td4RTransform>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x40,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %rdx,-0x38(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x10,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x20,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x30,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x14(%rbp)
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x14,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x24,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x34,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x10(%rbp)
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x18,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x28,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x38,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0xc(%rbp)
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x30(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x1c,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x2c,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm2
  mov    -0x30(%rbp),%rax
  add    $0x3c,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  mov    -0x38(%rbp),%rax
  add    $0xc,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  movss  -0x14(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x4,%rax
  movss  -0x10(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movss  -0xc(%rbp),%xmm0
  movss  %xmm0,(%rax)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L2
  call   <__stack_chk_fail@plt>
L2:
  leave
  ret
```

assistant: ```c
void td4RTransform(float *pin, float *f, float *pout) {
  float ptmp[3];
  ptmp[0] = pin[0] * f[0] + pin[1] * f[4] + pin[2] * f[8] + pin[3] * f[12];
  ptmp[1] = pin[0] * f[1] + pin[1] * f[5] + pin[2] * f[9] + pin[3] * f[13];
  ptmp[2] = pin[0] * f[2] + pin[1] * f[6] + pin[2] * f[10] + pin[3] * f[14];
  pout[3] = pin[0] * f[3] + pin[1] * f[7] + pin[2] * f[11] + pin[3] * f[15];
  pout[0] = ptmp[0];
  pout[1] = ptmp[1];
  pout[2] = ptmp[2];
}
```