complex arrays with PARDISO

[ringlery]

Greetings,
I'm trying to determine how to use complex arrays with the PARDISO solver. Does one simply change the declaration in mkl_pardiso.h so that you can pass an array of structs that represent a complex number, say _MKL_Complex16? I can't find any information in documentation other than where it states that it can be done. Any information would be appreciated.

[Sergey_K_Intel1]

If you are calling PARDISO from C/C++ program,please take a look at file mkl_dss.h to learn which complex types you can use for complex arrays. For example you can describe the array of non-zeros as follows:

_DOUBLE_COMPLEX_t cValues[your size];

For Fortran program, array of values in the case ofcomplex number should be described as follows

DOUBLE COMPLEX VALUES(NONZEROS)

All the best

Sergey