Great! gsl_linalg_QR_lssolve has worked without errors, now to test the results.
Thank you very much! On Tue, Apr 24, 2012 at 6:15 PM, Rhys Ulerich <[email protected]>wrote: > >> > Suppose we have A of size 9x2 and b of size 9x1 > >> > How would I solve Ax = b for x using GSL? Which method to use? > > >> Probably gsl_linalg_QR_lssolve. You can read up on the QR routines at > >> http://www.gnu.org/software/gsl/manual/html_node/QR-Decomposition.html. > > > That's what I'm looking into and although it seems I've supposedly > followed > > the correct steps, it still complains: > > > > gsl: qr.c:112: ERROR: QR matrix must be square > > qr.c:112 is within the routine gsl_linalg_QR_solve and not the routine > gsl_linalg_QR_lssolve. I expect the former to complain about > nonsquare matrices. Try calling the latter. > > - Rhys > -- alkar
