>> > 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
