[petsc-users] MKL Pardiso returning NAN values

[Chris Hewson]

Hi There,

We have a matrix that is singular and trying to solve it. We first use an
iterative solve with KSPBCGS, the solution vector is nan values and the
converged reason from PETSc of KSP_DIVERGED_NANORINF, that's great and what
I would expect.

Sometimes in our program we redo a failed solve using the MKL Pardiso
solver, when the same matrix and vectors get put into that solver which is
a KSPPREONLY, I get KSP_CONVERGED_ITS as a converged reason and solution
vector with nan values in it.

Stepping through the PETSc calls, I see that the external call to Pardiso
doesn't return an error for this, so not really the fault of PETSc, but
curious if y'all have seen this before or a solution/workaround to this?

*Chris Hewson*
Senior Reservoir Simulation Engineer
ResFrac
+1.587.575.9792
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240710/da79904f/attachment.html>