OONumerics User :
Date view Thread view Subject view Author view

From: Ben McLean (bmclean_at_[hidden])
Date: 2006-02-05 11:23:08

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.

Date view Thread view Subject view Author view