OONumerics User : |
From: bmclean (bmclean_at_[hidden])
Date: 2006-02-07 01:55:12
Particularly I am looking right now at freePOOMA and PETSc... both seem viable, if difficult to set up. Does anyone have subjective experience with these two, in terms of install-ability, ease of use, completeness, active maintenance, and execution speed?
thanks, Ben.
> -------Original Message-------
> From: Ben McLean <bmclean_at_[hidden]>
> Subject: [oon-list] fast linalg using arrays of matrices
> Sent: 05 Feb '06 16:23
>
> Hi all,
>
>
> I am (re)developing a successful seismic modeling program that I
> developed with my own flexible but slow (uses temporaries)
> arrays-of-matrices classes, and am looking for an appropriate library to
> handle the linear algebra. I have tried several, some of which are
> exceptional (blitz++, tvmet, thanks guys), but none of which I have had
> simple success in achieving all of what I would like to.
>
>
> My requirements, in sort-of order of importance or necesssity, are
>
> - 3D arrays with stencil operations (like blitz)
>
> - extensive tensor/matrix/vector types (like TVMET)
>
> - the array cells can be the forementioned matrix types, and still
> support chained expressions (havent seen this work well yet)
>
> - the result must be fast
>
> - automatically or easily run in parallel mode on SMP machines (like
> POOMA).
>
>
> I have had limited success integrating TVMET as an acceptible cell type
> for blitz arrays, though I have to manually code up several arg and
> return type combinations for every operator. I wouldnt know where to
> start parallelizing the resulting expressions. Does anyone know how
> fast/stable/complete freepooma is, and whether that would be ideal for my
> purpose? Other ideas?
>
>
> Thanks for all comments and suggestions,
>
> Ben McLean.
>