Change number of threads for pardiso

[Dong_Zheng]

I'm using MKL10.0.1.014 and testing MKL 10.1.0.006 on Linux. It seems that when evern I use "omp_set_num_threads" and/or "mkl_set_num_threads" to change the number of threads, MKL pardiso crashes.

The only way I can make it work is by doing the following

Use "mkl_set_num_threads" to change the number of threads
Unload pardiso with "phase=-1"
Reinitiate pardiso with "phase=11"
Re-do factoriation.

The draw back of the above method is that I have to re-do the factoriation. Is there a better way to do it?

DZ

[basel]

Dong,

>The draw back of the above method is that I have to re-do the factoriation. Is there a better way to do it?

No - this functionality is very useful but it is not yet implemented in the current version of PARDISO.

Regards,
Olaf

[Sergey_K_Intel1]

Dong

Would you please provide more details about the failure? In particular I'd like to take a look at iparm settings.

It looks like you are trying to use out-of-core version with iparm(60)=1. For this particular case mkl_num_threads has to be 1 since our out-of-core version has not been parallelized yet.

Thanks in advance

All the best

Sergey