Internal consistency check failure when using DSS

[Tianyu_D_]

Hi there,

I am working on a project to use DSS to solve large sparse linear equations. However I met with internal consistency check failure. I could see no problem with my codes. Could anyone help me to figure the problem out? Thank you so much. The attached is my project file.

316170

[Steven_L_Intel1]

For some reason, your attachment did not make it through. Is this an MKL solver you are calling?

[Tianyu_D_]

Yes. I am calling mkl_dss solver.

Here is my project file.

316430

[Tianyu_D_]

Here is my project file.

[Steven_L_Intel1]

Ok - I have moved this to the MKL forum.

[Gennady_F_Intel]

ok, we will look into the problem.

[Gennady_F_Intel · Accepted solution]

this is the problem with the test.

From MKL documentation:  DSS_REORDER:...

If opt contains the option MKL_DSS_MY_ORDER, then you must supply a permutation vector in the array perm. In this case, the array perm is of length nRows, where nRows is the number of rows in the matrix as defined by the previous call to dss_define_structure.

If opt contains the option MKL_DSS_GET_ORDER, then the permutation vector computed during the dss_reorder call is copied to the array perm. In this case you must allocate the array perm beforehand. The permutation vector is computed in the same way as if the option MKL_DSS_AUTO_ORDER is set.....

finally, making the following changes, should help:

!          error = dss_reorder( handle, MKL_DSS_GET_ORDER, Idum)
           error = dss_reorder( handle, MKL_DSS_AUTO_ORDER, Idum)

regards, Gennady

[Tianyu_D_]

Thank you so much!

Problem solved.