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