LAPACK and openmp

Discussion created by yakktr on May 31, 2008
Latest reply on Jun 11, 2008 by yakktr
Lapack routines in acml cannot be called from a pure subroutine


I am using dgetrf and dgetri for matrix inversion in a paralellized program by openmp. When I tag the subroutine with "PURE" the compiler gives error saying dgetrf and dgetri isnot pure. It is a simple parallel program, each core running the same code with different parameters so there is no message passing between cores.

 Does dgetri and dgetrf save variables which may  cause trouble in openmp?

While searching the net, I read that Intel MKL Lapack implementation is pure, should I use it or can I go with ACML