topic Segmentation fault with complex ScaLAPACK in Intel® oneAPI Math Kernel Library

Segmentation fault with complex ScaLAPACK

asd__asdqwe — Tue, 23 Jun 2015 12:42:02 GMT

Hello,

I'm trying to understand why the following simple program is segfaulting. When I'm using double instead of std::complex<double>, everything runs smoothly.

Thanks.

are there any specific input

Gennady_F_Intel — Tue, 23 Jun 2015 14:16:24 GMT

are there any specific input parameters? the # of mpi processes?

Could you please also share

Evarist_F_Intel — Tue, 23 Jun 2015 14:42:12 GMT

Could you please also share your link line? (interface and compiler is the most interesting part)

As far as I understand you

Evarist_F_Intel — Tue, 23 Jun 2015 17:05:14 GMT

As far as I understand you don't use GNU Fortran in your project (even though you built MPICH with GNU compiler and GNU Fortran bindings). Hence you do not need to link against libmkl_gf_lp64. Could you please try to use libmkl_intel_lp64 instead of gf?

mpicxx -cxx=icpc -lmkl_intel_thread -lmkl_core -lmkl_intel_lp64 -liomp5 -lpthread -lmpifort -lmkl_blacs_intelmpi_lp64 -lmkl_scalapack_lp64 scalapack_z.cpp -std=c++11

I am sorry for late respond.

Evarist_F_Intel — Thu, 25 Jun 2015 13:28:02 GMT

I am sorry for late respond. Yes I reproduced the crash.

The problem is in incorrect number of input parameters for the function pzgeqpf (actually it requires 13 parameters, while you pass only 11). Please take a look at the function declaration either in mkl/include/mkl_scalapack.h, or in MKL documentation or in Netlib documentation (http://www.netlib.org/scalapack/explore-html/de/dbc/pzgeqpf_8f_source.html).

yes, the problem is escalated

Gennady_F_Intel — Fri, 03 Jul 2015 18:52:49 GMT

yes, the problem is escalated. we will change these description at the on the next updates.

the doc has been fixed into

Gennady_F_Intel — Thu, 25 Feb 2016 08:01:43 GMT

the doc has been fixed into MKL v.11.3.1. we are sorry for delaying into notification.