On 10/09/2017 06:03 AM, Ines Arous wrote: > Hi GSL team, > > I am trying to use the QR function that is provided by GSL (which is very > effective compared to other QR implementation I found online). My problem is > that it computes the full decomposition while I need the thin QR > decomposition. > Is there a way to compute the thin QR decomposition using GSL? If so can you > please point me on how I can do that? > > Best regards, > Thanks! > Inès Arous
What exactly do you need to do? Are you trying to form the explicit Q1 and Q2 matrices? You almost never need to do this. If you just want the product Q^T b, you can use gsl_linalg_QR_QTvec(). If you want to compute the residual norm || y - A*x ||, this can be done without explicitly forming Q1 and Q2. If you really do need Q1, Q2 and R(NxN), you can use gsl_linalg_QR_unpack() and just take submatrix views to get each factor. Patrick
