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.