qr_decomp QR decomposition using Householder reflections on C++, Fortran, including using BLAS and LAPACK