tags:

views:

89

answers:

1

I wrote a program in c using MPI (Message Passing Inteface) that compute recursively the inverse of a lower triangular matrix. Every cpu sends 2 submatrices to other two cpus, they compute them and they give them back to the cpu caller. When the cpu caller has its submatrices it has to perform a matrix multiplication. In the recurrence equation the bottle neck is matrix multiplication. I implemented parallel multiplication with mpi in c but i'm not able to embed it into a function. Is it possible?

thanks, Simone

A: 

You can use PDGEMM/PSGEMM functions for parallel mpi matrix multiplication from SCALAPACK/PBLAS.

http://www.netlib.org/scalapack/pblas_qref.html#PvGEMM

download here http://www.netlib.org/scalapack/

osgx