Don't assume commutativity for multiplication in QR calculations.

I tried to QR factorize some matrices of quaternions and found out than I had assumed commutativity of multiplication. That shouldn't be the case after this pr.