Problem with ScaLapack PCLAPRNT (RH Linux, ifort 11.1, openmpi 1.2.4, mkl 10.2.x.x)
I have run into a problem with a fortran code when calling scalapack routine PCLAPRNT. As soon as this subroutine is called the code will die on a seg fault. If I use a hand built scalapack and blacs (using source from netlib.org) library with the same code then it will run correctly. It only seg faults when I link the MKL scalapack.
My machine is Linux running RH5 and I trie several flavors of MKL 10.2 including 10.2.6.038. I also use MKL for the blacs and blas. The compiler is ifort 11.1 and the mpi library I am using is a hand-compiler ifort based openmpi v1.2.4.
I could not reproduce your error, but that is understandable because you are using a number of custom-built components.
The example uses a couple of big arrays which are allocated on the stack. That could cause problems. You may try reducing LWORK and LRWORK to, say, 10000; you could make sure that you have raised the stack limit adequately; you could use the -heap-arrays compiler option; you could declare some/all of RWORK, WORK, A, Z with the SAVE attribute.
If none of these suggestions work, you may need to replace your custom-built libraries one at a time by a standard version and try running.