Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Highlighted
Beginner
5 Views

does pardiso work for a 1 element matrix?

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
Highlighted
Beginner
5 Views

Nevermind, the problem was a dumb mistake on my part.
0 Kudos