[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## Re: [Help-gsl] Linear System Solve

**From**: |
Rhys Ulerich |

**Subject**: |
Re: [Help-gsl] Linear System Solve |

**Date**: |
Tue, 24 Apr 2012 10:15:52 -0500 |

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