Hi,
In the attached test case (a.cpp), I'm trying to use the scalapack PZGESVD routine, but it crashes for a reason I cannot figure out. I'm calling from C++. If I call PDGESVD on a real (double precision) matrix, then everything works fine. I've attached the screen output from both a real compile (real.png) and a complex compile (complex.png). The only difference between the two compiles is that the real version defines dtype=1 at line 94 and the complex version defines dtype=2 at line 94. Also, the PZLAPRNT routine produces weird results unless you comment out line 118 (the PZGESVD call still crashes though).
I'm using MKL version 11.1.0 and intel MPI version 4.1 on a linux system.
Am I calling the routine incorrectly, or is this a library issue? Thanks in advance.
John
