Page 1 of 1

mkl+scalapack issue

Posted: Thu Jan 30, 2014 11:45 pm
by Jordan
Dear all,

I am experiencing some issue using mkl+scalapack.
I compilte the code with -O0 -g -traceback options and with the ifort 14.0.1 compiler. The MKL is the version 11.0.1 and openmpi 1.7.2

Everything goes fine except the segfault that occures at the very beginning of the code.
The functions stack is
abi_linalg_init
init_scalapack
build_grid_scalapack
BLACS_GRIDINIT which is in th MKL library and can not be debugged.

I removed the scalapack linal flavor and it works.

I am not sure I really need scalapack but it is always good to report errors we find.

Best,

Jordan

EDIT : Scalapack of the MKL library, architecture Marenostrum III at the BSC (Intel Xeon E5-2670). Abinit is my 7.7.1-private, I can try the released 7.6

Re: mkl+scalapack issue

Posted: Tue Feb 04, 2014 5:34 pm
by pouillon
As usual: which version of Abinit? Which version of ScaLAPACK? Which architecture? ...

A config.log would also be welcome.

Re: mkl+scalapack issue

Posted: Tue Feb 04, 2014 6:25 pm
by gmatteo
I've encountered similar problems with MKL+Scalapack in older versions of Abinit (I guess Jordan is using v7.6 or 7.7.1).
I didn't have time to investigate the problem and I decided to disable Scalapack in the GW code.
Of course I'm not happy with this solution, and I would prefer to have a suite of unit tests for the Scalapack wrappers.
pFunit seems an interesting framework (http://pfunit.sourceforge.net/) as it provides supports for MPI and Omp

Re: mkl+scalapack issue

Posted: Tue Feb 04, 2014 6:41 pm
by Jordan
Yes Matteo.
I did almost the same as you -> disable ScaLAPACK in the whole abinit since the segfault arrives at the very beginning of the code....

Re: mkl+scalapack issue  [SOLVED]

Posted: Fri Feb 07, 2014 4:11 pm
by Jordan
I finally found out what the problem was.
I use openmpi and to use mkl+scalapack, one should link with mkl_blacs_openmpi !
Everything works smoothly now.

Jordan