Re: [Help-gsl] Linear System Solve

[Rhys Ulerich]

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