PARDISO, _DOUBLE_COMPLEX_t, exception thrown

[Perpf2000]

Hi!
I'm trying to replace my current sparse solver in a c++ code by MKL/PARDISO. When I run the code, factorization goes well. However, on phase=33 (solve) calling PARDISO(...) throws an exception. As no source is available for debugging I',m stuck.
I suspect that it might have to do with my using the complex-format, which is poorly documented in the manual...

Is there anyone who can point me the right direction?

thanks & best regards!

[Gennady_F_Intel]

Hello Nils,

1. We have the example of program to show the use of the "PARDISO" routine for complex unsymmetric linear systems:

please see MKLROOTexamplessolversourcepardiso_unsym_complex_c.c

may be it will help you to understand how to use complex datatype for PARDISO.

2. Can you get us the test case, If it will not help you?

Feel free to reply in Private if needed.

--Gennady

[Perpf2000]

Quoting - Gennady Fedorov (Intel)

Hello Nils,

1. We have the example of program to show the use of the "PARDISO" routine for complex unsymmetric linear systems:

please see MKLROOTexamplessolversourcepardiso_unsym_complex_c.c

may be it will help you to understand how to use complex datatype for PARDISO.

2. Can you get us the test case, If it will not help you?

Feel free to reply in Private if needed.

--Gennady

Wow! You guys are fast!!

I just came back to post: "Once again it was my own stupidity that gave me a hard time..."
I found the problem, it had nothing to do with the complex-issue.

Thnaks again for your quick response!