user: What is the c source code of the assembly code below?

```asm
<floats_xytcov_mul>:
  endbr64
  push   %r15
  mov    %r8,%r15
  push   %r14
  mov    %rdx,%r14
  push   %r13
  mov    %rdi,%r13
  push   %r12
  mov    %rcx,%r12
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %r9,%rbx
  sub    $0x38,%rsp
  movss  (%rdi),%xmm2
  movss  0x4(%rdi),%xmm4
  movss  0x8(%rdi),%xmm0
  movss  0x4(%rdx),%xmm1
  movss  %xmm2,0xc(%rsp)
  movss  (%rdx),%xmm2
  movss  %xmm1,0x1c(%rsp)
  movss  %xmm2,0x2c(%rsp)
  movss  %xmm4,0x10(%rsp)
  movss  %xmm0,0x8(%rsp)
  call   <sin@plt>
  pxor   %xmm4,%xmm4
  movss  0x8(%rsp),%xmm0
  cvtsi2ss %rax,%xmm4
  movss  %xmm4,0x4(%rsp)
  call   <cos@plt>
  movss  0x4(%rsp),%xmm4
  pxor   %xmm0,%xmm0
  movss  0x10(%rbp),%xmm1
  cvtsi2ss %rax,%xmm0
  movss  0x20(%r12),%xmm12
  movss  0x2c(%rsp),%xmm2
  movaps %xmm4,%xmm6
  xorps  D9(%rip),%xmm6
  movss  %xmm1,0x18(%rsp)
  movss  0x1c(%rsp),%xmm1
  movss  %xmm12,0x28(%rsp)
  movss  0x20(%rbp),%xmm8
  movss  (%r12),%xmm5
  movaps %xmm6,%xmm14
  movaps %xmm1,%xmm12
  movss  0x14(%r12),%xmm9
  movss  0x10(%r12),%xmm13
  mulss  %xmm2,%xmm14
  movss  %xmm5,0x1c(%rsp)
  movss  0x4(%rbp),%xmm3
  movss  0x8(%rbp),%xmm11
  mulss  %xmm0,%xmm12
  movss  %xmm9,0x24(%rsp)
  movss  0x8(%r12),%xmm15
  movss  0x14(%rbp),%xmm10
  mulss  %xmm4,%xmm1
  movss  %xmm3,0x14(%rsp)
  movss  0x4(%r12),%xmm3
  mulss  %xmm0,%xmm2
  movss  %xmm15,0x20(%rsp)
  movaps %xmm14,%xmm7
  subss  %xmm12,%xmm7
  subss  %xmm1,%xmm2
  movaps %xmm8,%xmm1
  mulss  %xmm7,%xmm1
  movaps %xmm1,%xmm9
  movaps %xmm5,%xmm1
  movaps %xmm6,%xmm5
  mulss  %xmm13,%xmm5
  mulss  %xmm0,%xmm1
  movss  %xmm5,0x8(%rsp)
  movaps %xmm11,%xmm5
  addss  %xmm11,%xmm5
  movss  %xmm1,0x4(%rsp)
  movaps %xmm3,%xmm1
  addss  %xmm3,%xmm1
  movaps %xmm5,%xmm15
  movaps %xmm9,%xmm5
  addss  %xmm11,%xmm9
  mulss  %xmm7,%xmm15
  mulss  %xmm7,%xmm5
  mulss  %xmm10,%xmm7
  movaps %xmm9,%xmm11
  mulss  %xmm2,%xmm11
  addss  %xmm15,%xmm5
  movss  0x4(%rsp),%xmm15
  mulss  %xmm0,%xmm15
  addss  %xmm5,%xmm15
  movaps %xmm1,%xmm5
  mulss  %xmm0,%xmm5
  mulss  %xmm4,%xmm1
  mulss  %xmm6,%xmm5
  mulss  %xmm0,%xmm1
  addss  %xmm15,%xmm5
  movss  0x8(%rsp),%xmm15
  mulss  %xmm6,%xmm15
  addss  %xmm15,%xmm5
  addss  0x0(%rbp),%xmm5
  movss  0x20(%rsp),%xmm15
  movss  %xmm5,(%rbx)
  movss  0x14(%rsp),%xmm5
  addss  %xmm11,%xmm5
  movaps %xmm6,%xmm11
  mulss  %xmm3,%xmm11
  addss  0x4(%rsp),%xmm11
  mulss  %xmm0,%xmm3
  addss  0x8(%rsp),%xmm3
  addss  %xmm5,%xmm7
  mulss  %xmm4,%xmm11
  mulss  %xmm0,%xmm3
  addss  %xmm7,%xmm11
  movaps %xmm10,%xmm7
  addss  %xmm10,%xmm7
  addss  %xmm3,%xmm11
  movaps %xmm15,%xmm3
  mulss  %xmm0,%xmm3
  mulss  %xmm2,%xmm7
  movss  %xmm11,0x4(%rbx)
  addss  %xmm9,%xmm3
  movss  0x24(%rsp),%xmm9
  movss  %xmm11,0xc(%rbx)
  mulss  %xmm9,%xmm6
  addss  %xmm3,%xmm6
  movaps %xmm8,%xmm3
  addss  0x28(%rsp),%xmm8
  mulss  %xmm2,%xmm3
  movss  %xmm6,0x8(%rbx)
  movss  %xmm6,0x18(%rbx)
  movss  %xmm8,0x20(%rbx)
  movaps %xmm3,%xmm5
  addss  %xmm10,%xmm3
  mulss  %xmm2,%xmm5
  addss  0xc(%rsp),%xmm2
  addss  %xmm7,%xmm5
  movss  0x1c(%rsp),%xmm7
  mulss  %xmm4,%xmm7
  mulss  %xmm4,%xmm7
  mulss  %xmm15,%xmm4
  addss  %xmm7,%xmm5
  movaps %xmm13,%xmm7
  mulss  %xmm0,%xmm7
  addss  %xmm3,%xmm4
  addss  %xmm5,%xmm1
  mulss  %xmm0,%xmm7
  mulss  %xmm9,%xmm0
  addss  %xmm7,%xmm1
  addss  0x18(%rsp),%xmm1
  addss  %xmm0,%xmm4
  movaps %xmm12,%xmm0
  subss  %xmm14,%xmm0
  addss  0x10(%rsp),%xmm0
  movss  %xmm1,0x10(%rbx)
  movss  %xmm4,0x14(%rbx)
  movss  %xmm4,0x1c(%rbx)
  movss  %xmm2,(%r15)
  movss  %xmm0,0x4(%r15)
  movss  0x8(%r13),%xmm0
  addss  0x8(%r14),%xmm0
  movss  %xmm0,0x8(%r15)
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: data_label: D9
datatype: byte[16]

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

assistant: ```c
void floats_xytcov_mul(const float xyta[3], const float Ca[9], const float xytb[3], const float Cb[9], float xytr[3], float Cr[9]) {
  float xa = xyta[0], ya = xyta[1], ta = xyta[2];
  float xb = xytb[0], yb = xytb[1];
  float sa = (float)sin(ta), ca = (float)cos(ta);
  float P11 = Ca[0], P12 = Ca[1], P13 = Ca[2];
  float P22 = Ca[4], P23 = Ca[5];
  float P33 = Ca[8];
  float Q11 = Cb[0], Q12 = Cb[1], Q13 = Cb[2];
  float Q22 = Cb[4], Q23 = Cb[5];
  float Q33 = Cb[8];
  float JA13 = -sa * xb - ca * yb;
  float JA23 = ca * xb - sa * yb;
  float JB11 = ca;
  float JB12 = -sa;
  float JB21 = sa;
  float JB22 = ca;
  Cr[0] = P33 * JA13 * JA13 + 2 * P13 * JA13 + Q11 * JB11 * JB11 + 2 * Q12 * JB11 * JB12 + Q22 * JB12 * JB12 + P11;
  Cr[1] = P12 + JA23 * (P13 + JA13 * P33) + JA13 * P23 + JB21 * (JB11 * Q11 + JB12 * Q12) + JB22 * (JB11 * Q12 + JB12 * Q22);
  Cr[2] = P13 + JA13 * P33 + JB11 * Q13 + JB12 * Q23;
  Cr[3] = Cr[1];
  Cr[4] = P33 * JA23 * JA23 + 2 * P23 * JA23 + Q11 * JB21 * JB21 + 2 * Q12 * JB21 * JB22 + Q22 * JB22 * JB22 + P22;
  Cr[5] = P23 + JA23 * P33 + JB21 * Q13 + JB22 * Q23;
  Cr[6] = Cr[2];
  Cr[7] = Cr[5];
  Cr[8] = P33 + Q33;
  xytr[0] = ca * xb - sa * yb + xa;
  xytr[1] = sa * xb + ca * yb + ya;
  xytr[2] = xyta[2] + xytb[2];
}
```