OONumerics User : |
From: Dominik Heide (dheide_at_[hidden])
Date: 2005-03-24 05:31:35
On Wednesday 23 March 2005 21:32, Nick Cummings wrote:
> I'd like the post a follow-up question to one I posted a few months ago.
> First, let me thank those who responded for their helpful comments. I
> gave MTL a serious look but had some issues with solving complex linear
> systems. In the end, I've been using GMM++ and have the main linear
> algebra portion of what I wanted to do complete.
>
> My program requires a matrix that is the result of a rather lengthy matrix
> expression, which will be changed later for other problems. I would like
> to design things so that this matrix expression may be input in a fairly
> transparent way, for which I'd really like to use operator overloading.
> ...
I wanted to do the same half a year agoo. The main problem when dooing this is
to avoid the creation of temporary variables. TBCI
(http://plasimo.phys.tue.nl/TBCI/) should do what you want to do with while
avoiding temoraries but when I tried it the classes I used were buggy and the
documentation is very poor.
I ended up using PetSc (http://www-unix.mcs.anl.gov/petsc/petsc-2/) which is
written in C ( C++ compatible) which is implemented completly paralell and
has very nice features, like different iterative solvers a nonlinear solver,
debugging facilities for paralell programms ,...
If you find a good, fast, bugfree and well documented objectoriented C++
library with overloaded operators for matrices and vectors please let us
know.
Dominik