← Back to team overview

dolfin team mailing list archive

Re: Linear algebra

 

I forgot to include Trilinos in the list. Rob has mentioned it a few
times.

How does PETSc compare to Trilinos? (I know your answer may be biased,
but nevermind.)

Do you think C++ wrappers for PETSc will happen any time soon? One
option would be to do this as a FEniCS project.

/Anders

On Sun, Oct 24, 2004 at 04:09:37PM -0500, Matthew Knepley wrote:
> Anders Logg <logg@xxxxxxxxx> writes:
> > There are a couple of different options we could consider:
>     I can comment on the ones I know about.
> > 
> > 1. PETSc: http://www-unix.mcs.anl.gov/petsc/petsc-2/
> > 
> > Pros: Satisfies B1-4
> > Cons: Does not seem to satisfy B5. Maybe a future version of PETSc will?
>      PETSc proper will never move away from the C interface. However, I
> have already written Python wrappers with a fully OO interface (not hard).
> I used SIDL, so the same method could generate C++ wrappers as well.
> > 
> > 2. MTL: http://www.osl.iu.edu/research/mtl/
> > 
> > Pros: Satisfies B1 (?) and B5
> > Cons: Does not satisfy B2-4?
>      I believe they are working on parallelizing it, and Andy does maintain
> this, although I think it does not really have a large customer base yet. It
> is definitely open source.
> > 
> > 3. uBLAS: http://www.boost.org/libs/numeric/ublas/doc/index.htm
> > 
> > Pros: Satisfies B1-3, B5
> > Cons: Does not satisfy B4
>      I don't think parallelism will even happen here. It is much more likely
> from Andy Lumsdaine and MTL.
> > 
> > 4. POOMA: http://www.codesourcery.com/pooma/
> > 
> > Pros: ?
> > Cons: ?
>      I would not even consider using this.
> > 
> > 5. Sparselib++: http://math.nist.gov/sparselib++/
> > 
> > Pros: ?
> > Cons: ?
>      Or this.
> > 
> > 6. TNT: http://math.nist.gov/tnt/
> > 
> > Pros: ?
> > Cons: ?
>      Or this.
> > 
> > 7. Blitz++: http://www.oonumerics.org/blitz/
> > 
> > Pros: ?
> > Cons: ?
>     I like this project since it is so innovative (and Todd is cool). However,
> I did not think it was ever intended to be parallel, and I am not sure
> what the support for sparse matrices is like, but I think it might even be
> missing.
> 
>     I would note that you have not considered Trilinos, which is a library
> very similar to PETSc from SNL.
> 
>         Matt
> > 
> > /Anders



References