QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Fix quaternion conjugation order on vector rotation. 

Now it does not crash on C4v and the projectors seem to be indeed
projectors, but the translation operator decomposition is still wrong.


Former-commit-id: 6273c5a76b4e424f77b12514df209b995999b24b
This commit is contained in:
Marek Nečada 2019-03-12 22:22:02 +02:00
parent 208e05d811
commit 02ec8804d0
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
qpms/wigner.h
View File

@ -152,8 +152,8 @@ static inline cart3_t qpms_quat_rot_cart3(qpms_quat_t q, const cart3_t v) {
//const qpms_quat_t qc = qpms_quat_normalise(qpms_quat_pow(q, -1)); // implementation of _pow wrong! //const qpms_quat_t qc = qpms_quat_normalise(qpms_quat_pow(q, -1)); // implementation of _pow wrong!
const qpms_quat_t qc = qpms_quat_conj(q); const qpms_quat_t qc = qpms_quat_conj(q);
const qpms_quat_t vv = qpms_quat_from_cart3(v); const qpms_quat_t vv = qpms_quat_from_cart3(v);
return qpms_quat_to_cart3(qpms_quat_mult(qc, return qpms_quat_to_cart3(qpms_quat_mult(q,
qpms_quat_mult(vv, q))); qpms_quat_mult(vv, qc)));
} }
/// Wigner D matrix element from a rotator quaternion for integer \a l. /// Wigner D matrix element from a rotator quaternion for integer \a l.