Re: [Help-gsl] Linear System Solve

[Dimitris K.]

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 <address@hidden>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