Re: [ublas] A x = b; binding for large sparse matrix A

From: Gunter Winkler (guwi17_at_[hidden])
Date: 2007-11-07 18:40:04

Am Dienstag, 6. November 2007 14:05 schrieb Per Abrahamsen:
> Hi ublas,
>
> We currently use lu_factorize and ublas::matrix A to solve an
> equation system A x = b, and it works fine for the small matrixes we
> have tested it with.
>
> What binding / library do you suggest?

The AFAIK simplest way is to try 

http://guwi17.de/ublas/ublas_csparse.hpp

together with 

http://www.cise.ufl.edu/research/sparse/CSparse/

Alternatively you can use the boost-bindings library from the 
boost-sandbox repository which provided UMFPACK binding. (@Karl: 
correct me if I'm wrong)

in both cases you should prepare you matrix as 

compressed_matrix<double, row_major> for C-libraries

or

compressed_matrix<double, column_major, 1> for FORTRAN libraries

mfg
Gunter