Add a testing for the QR factorization of two square matrices on top of each other, and the associated generation of the Q that is required in QDWH