implements the Householder transformation for quaternion matrices.

Translation of matlab code in arXiv:1203.6151v4. Not used yet.