Boost logo

Ublas :

From: Jens Seidel (jensseidel_at_[hidden])
Date: 2008-03-06 09:11:58

On Thu, Mar 06, 2008 at 02:32:50PM +0100, Jonas wrote:
> Which library/solver (that has uBLAS bindings, preferrably) would you
> recommend for solving a very sparse, real system Ax = b with size on
> the order of 10^5? My initial guess from monitoring the list would be
> umfpack. Would that be a good candidate?

Mmh, this is sufficient large to require special algorithms. Please
explain in more detail the properties of your matrix. Is it symmetric,
positive definite, an M matrix, ...? Does it result from
discretisation of a (partial) differential equation, which one?
Does there exist a hierarchy of matrices?

It is nearly impossible for me to answer your question without these

You could also look at Dolfin (
which contains some simple standard solvers and has also interfaces to
more difficult ones (such as Algebraic Multigrid).

Do you have a typical matrix somewhere available for download (in text
format)? Maybe I could test my Algebraic Multigrid solver (no publically
available, but I could give you some references) with it?