Re: [ublas] How to solve systems of linear equations?

Subject: Re: [ublas] How to solve systems of linear equations?
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2010-03-09 15:39:01

Tim Odenthal wrote:

> 
>> You could give the MUMPS part of Boost.Numeric_bindings a shot. For
>> Mumps, see here http://mumps.enseeiht.fr/. For the "stable"
>> numeric_bindings, see
>> here http://svn.boost.org/svn/boost/sandbox/numeric_bindings-v1/  Support
>> for MUMPS in the development version of the bindings is expected within a
>> fortnight.
>> [...]
>> HtH,
>> Cheers,

Erm... a bit of a long fortnight, but we've taken a couple of steps already. 
Should not be long now :-)

> In the meantime, I implemented a highly specialized version of a
> ConjugateGradient solver myself, since this seemed easier than to adapt
> our code to use ublas::vectors/matrices everywhere - we use special
> vectors...

I would be interested in using your custom vectors as an example of how to 
bind foreign containers in the current numeric_bindings -- if you're 
interested, please let me know.

Regards,

Rutger