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)
Edited by Mathieu Faverge