Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.

Pardiso Output = -1.#IND00

rditz
Beginner
770 Views
0 Kudos
5 Replies
Alexander_K_Intel2
770 Views

Hi,

It would be great if you provide small testcase to reproduce this issue on our side. Maybe there are some problems in matrix structure, maybe some input parameters like iparm are incorrect, and maybe condition number of input matrix is really huge. To choose correct answer we need to run testcase.

With best regards,

Alexander Kalinkin


0 Kudos
rditz
Beginner
770 Views
Hi, sorry but I am new to this.What do you exactly mean by testcase ? Should I post my code or should I send the .cpp file.Also sparse matrix solver is a part of a bigger FEM code.Will it be possible for you to pinpoint the error ?
Also I think there cannot be a mistake with iparm as it is solving correctly for lesser(magnitudewise) coefficients of A matrix.
Thank you for the reply.
0 Kudos
Alexander_K_Intel2
770 Views
Hi,
You could send .cpp example or file with ia, ja, a, rhs, iparm and correct solution arrays. Any variant would be convenientfor us if we reproduce your issue.
With best regards,
Alexander Kalinkin
0 Kudos
Gennady_F_Intel
Moderator
770 Views
please use this post to see how attach the code to the forum.
0 Kudos
rditz
Beginner
770 Views
Hey thanks for the reply.As I said earlier Pardiso solver is part of a larger code.I have attached the code and 3 txt files.Open each files through freopen.Note the output of each file.

1) 3bricks......smaller magnitude coefficient + small size of matrix
2) 3bricks_diff._coefficients...larger magnitude coefficient + small size of matrix
3) Increased_size_of_matrix...smaller magnitude coefficient + larger size of matrix

Once again thanks for the reply and the help
0 Kudos
Reply