[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Most efficient and accurate method to solve linear system
From: |
Brian Gough |
Subject: |
Re: [Help-gsl] Most efficient and accurate method to solve linear system using GSL? |
Date: |
Fri, 05 Dec 2008 17:16:20 +0000 |
User-agent: |
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???
Hello,
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
manual.
--
Brian Gough
(GSL Maintainer)
Support freedom by joining the FSF!
http://www.fsf.org/news/fall-2008-fundraiser