* Sparse matrix linear solver added to GSL
@ 2014-01-30 23:40 Patrick Alken
0 siblings, 0 replies; only message in thread
From: Patrick Alken @ 2014-01-30 23:40 UTC (permalink / raw)
To: gsl-discuss
Hello all,
The sparse matrix routines are in pretty good shape now and they are
all located in a branch called 'sparse' on the gsl git. Furthermore I've
added a sparse linear solver (for A x = b) based on the Generalized
Residual Minimum Method (GMRES) iterative solver. GMRES is considered a
very good/robust solver for a wide class of matrices, though it does
have convergence issues on some matrices. In these cases it helps to
precondition the system but I have not currently implemented any
preconditioners. Maybe one day I'll get around to that.
Everything is documented, with an example of solving the 1D Poisson
equation sparse system, and I've added lots of automated tests so things
are working well. Anyone who uses sparse matrices in their work I'd be
interested in any feedback. I plan to merge the sparse branch into the
master branch in the next week or so.
Patrick
^ permalink raw reply [flat|nested] only message in thread
only message in thread, other threads:[~2014-01-30 23:40 UTC | newest]
Thread overview: (only message) (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2014-01-30 23:40 Sparse matrix linear solver added to GSL Patrick Alken
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).