[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Most efficient and accurate method to solve linear system
Re: [Help-gsl] Most efficient and accurate method to solve linear system using GSL?
Fri, 05 Dec 2008 17:16:20 +0000
Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI)
At Wed, 03 Dec 2008 21:26:09 -0700,
Evan Ruzanski wrote:
> I'm using the GSL C library to solve a linear system, Y = HX, with a
> design matrix, H, of real numbers sized with many more rows than
> columns, i.e., 107910 X 27. I'm looking for advice on the best function
> and method to use in terms of speed, stability, and accuracy. Can you
> please help???
For a general matrix try gsl_linalg_SV_decomp_mod which is designed
for the case of M>>N. The reference for it can be found in the GSL
Support freedom by joining the FSF!