topic Input error message - help in Intel® oneAPI Math Kernel Library

Input error message - help

Brad_Doudican — Mon, 14 Nov 2011 18:46:04 GMT

Hi - I'm having trouble deciphering an error code. When I run Phase=11 in Pardiso I get the following error messages. However, I've checked the input arrays and they seem to be reasonable. Any thoughts on what I'm missing?

When iparm(27)=0:

*** Error in PARDISO ( sequence_ido,parameters) error_num= 8
*** Input check: ia(neqns+1)_new -1 _old 0 are incompatible
*** Input parameters: inconsistent error= 8 max_fac_store_in: 1
matrix_number_in : 1 matrix_type_in : 2
ido_in : 11 neqns_in : 14751
ia(neqns_in+1)-1 : -1 nb_in : 1

And when iparm(27)=1:

*** Error in PARDISO ( sequence_ido,parameters) error_num= 18
*** Input check: unexpected error with working arrays ia and/or ja

To derive the input arrays I've used the following code to extract the sparse storage form from a full upper triangular matrix GLKF(NEQ,NEQ) where nomenclature JA = COLUMNS, IA = ROWINDEX, A = VALUES:

KOUNT=COUNT(GLKF/=0d0)

ALLOCATE ( VALUES(KOUNT), COLUMNS(KOUNT), ROWINDEX(NEQ+1) )

KOUNTER = KOUNT
KOUNT = 0
ROWINDEX(1)=1

DO J=1,NEQ
ROWINDEX(J+1)=ROWINDEX(J)
DO I=1,NEQ
IF (GLKF(J,I)/=0.OR.I==J) THEN
ROWINDEX(J+1)=ROWINDEX(J+1)+1
KOUNT=KOUNT+1
VALUES(KOUNT)=GLKF(J,I)
COLUMNS(KOUNT)=I
END IF
END DO
END DO

Many thanks!

Input error message - help

mecej4 — Mon, 14 Nov 2011 19:12:51 GMT

The values of KOUNT from the first line of code shown above will not agree with the that after the nested loops since you handle zero-valued diagonal elements differently in the two places.

Depending on how you call PARDISO, a consistency check would catch this discrepancy.

Input error message - help

Brad_Doudican — Mon, 14 Nov 2011 19:28:06 GMT

Thanks mecej4, but earlier in the code I check for zeros on the diagonal - there are none if we get to this point. Any thoughts on the specific error messages?

Input error message - help

mecej4 — Mon, 14 Nov 2011 21:25:42 GMT

If there are no zeros on the diagonal, why do you have this statement:

IF (GLKF(J,I)/=0.OR.I==J) THEN

If there are some zero values on the principal diagonal when this statement is executed, the statements in the IF block will be executed even for those zero elements.

Input error message - help

Brad_Doudican — Mon, 14 Nov 2011 23:15:03 GMT

I was worried that would be a hang up in getting my question answered. The .OR.I==J code functions but doesn't have an effect - it's residual from a previous attempt at remedy. As I stated previously, no zeros on the diagonal. That's taken care of previously.

Do you know where I can find out more about the various error codes cited in my original post?

Input error message - help

Alexander_K_Intel2 — Tue, 15 Nov 2011 01:03:48 GMT

Hi,

One can see the next line: "ia(neqns_in+1)-1 : -1" so its seem that you have some problem with ia array. Could you print the value of neqns and ia(neqns_in+1) before first pardiso call?

With best regards,

Alexander Kalinkin