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

does pardiso work for a 1 element matrix?

Brian_Murphy
New Contributor II
283 Views

I am finding that pardiso does not work for a 1 element matrix.  Calling pardiso with phase=11 returns error=0, which means success, but the pt(:) pointer array is all zero's.  So when I subsequently call pardiso with phase=23, it fails with error=-7.

Can someone please confirm that pardiso does or does not work for a 1 element matrix?

I want to run my code for a one element matrix because I am trying to run a validation case to compare to sample calculations in a finite element textbook.  The case is a transient integration of a single dof system, thus the system matrix is a single nonzero element.

0 Kudos
1 Reply
Brian_Murphy
New Contributor II
283 Views
Nevermind, the problem was a dumb mistake on my part.
0 Kudos
Reply