I have a problem with QR-factorization using dgeqrf with my fortran code. I need to create Q and R matrices, so that A = Q * R and then compute err=||A - Q*R|| / ||A||, wher
