Am 14.11.2008 um 23:00 schrieb Greg Lindahl:

On Thu, Nov 13, 2008 at 11:07:59PM -0500, Cesar Andres Polindara Lopez wrote:

I have developed a finite element model for the simulation of
poroelastic materials. The model requires the implementation of a time
integration scheme (also known as step by step algorithms) hence
iterations are performed. For each iteration I have to solve a linear
system of equations so I have to calculate the inverse of a matrix and
then execute a multiplication.

Good news: look up BLAS, and you will find a standard interface for

http://netlib.org/liblist.html

doing matrix computations. There are parallel implementations
available: Atlas, AMD ACML, Intel's math lib, etc etc.

If you want to extend it to mutilple machines, you can also look into PBLAS as part of ScaLAPACK (above link) and/or it's Intel implementation in Intel's MKL.

-- Reuti

Bad news: Are you really inverting the matrix? You should probably be
doing something like LU decomposition.

-- greg


_______________________________________________
Beowulf mailing list, Beowulf@beowulf.org
To change your subscription (digest mode or unsubscribe) visit http://www.beowulf.org/mailman/listinfo/beowulf

_______________________________________________
Beowulf mailing list, Beowulf@beowulf.org
To change your subscription (digest mode or unsubscribe) visit 
http://www.beowulf.org/mailman/listinfo/beowulf

Reply via email to