Add Polar decomposition with QDWH algorithm

The main objective of this PR is to add a Polar Decomposition algorithm (A = U_p H ) with U_p orthonormal, and H, Hermitian. To enable this, this PR adds:

• Add subroutines to compute the QR on the full matrix, as well as the Q generation.
• Fix issue #36 (closed)