topic Re: PARDISO, _DOUBLE_COMPLEX_t, exception thrown in Intel® oneAPI Math Kernel Library

PARDISO, _DOUBLE_COMPLEX_t, exception thrown

Perpf2000 — Mon, 31 Aug 2009 14:35:20 GMT

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!

Re: PARDISO, _DOUBLE_COMPLEX_t, exception thrown

Gennady_F_Intel — Mon, 31 Aug 2009 15:29:47 GMT

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

Re: PARDISO, _DOUBLE_COMPLEX_t, exception thrown

Perpf2000 — Mon, 31 Aug 2009 15:37:37 GMT

Quoting - Gennady Fedorov (Intel)

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!