0-based CSR Format (3-Array Variation) with PARDISO

0-based CSR Format (3-Array Variation) with PARDISO

Daniel Vlasic's picture

Hi,

I'd like to be able to use 0-based sparse matrix format when solving sparse systems with PARDISO. MKL docummentation mentions that sparse routines work with either 0- or 1-based formats, but PARDISO interface says otherwise. Is this possible?

Thanks,

Daniel

4 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.
Gennady Fedorov (Intel)'s picture

Daniel, what do you mean by that?"PARDISO interface says otherwise". set iparm[34] == 1 and use C-Ctyle indexing for ia and ja arrays. --Gennady

Daniel Vlasic's picture

This is not mentioned in the documentation for my MKL (10.2.6.037).

From mklman.pdf page 2540: "Note that the row and columns numbers start from 1."

Gennady Fedorov (Intel)'s picture

Ok, I see. that's right. MKL's PARDISO started to supportzero-based (C-style) array indexing since 10.3. You can participate 10.3 beta version program here. You can find there the example "pardiso_sym_0_based.c" which shows how to call PARDISO with these arrays. --Gennady

Login to leave a comment.