I would like to do a QR factorization using LAPACK. From the documentation available here https://software.intel.com/en-us/node/521003#E832D468-0891-40EC-9468-6868C8EC9AD0 , I've decided to use geqrf for the factorization.
As I need to solve a Least square problem, I need to solve R.x = (Q1)^T b as explained in the documentation. You can apply matrix Q with the function ormqr. But, which routine should I use to solve R.x = y ?
Velvia, can you explain what you are attempting to do? If all you want is the least-squares solution to overdetermined equations with a full-rank matrix, you can use the "driver" routine ?GELS().