Q1=randn(4,1);
Q1=Q1/norm(Q1);
R1=quatern2rotMat_nolo(Q1);
Q2=randn(4,1);
Q2=Q2/norm(Q2);
R2=quatern2rotMat_nolo(Q2);
R12=R1*R2;
Q12=rotMat2qRichard_nolo(R12)
quaternProd(Q1,Q2)
R1t2=R1'*R2;
Q1t2=rotMat2qRichard_nolo(R1t2)
quaternProd([Q1(1);-Q1(2:4)],Q2)
R12t=R1*R2';
Q12t=rotMat2qRichard_nolo(R12t)
quaternProd(Q1,[Q2(1);-Q2(2:4)])