I suspect there are some randomness from external matsolverpackage. I use a parallel matrix as preconditioner, and the preconditioner is exactly LU factored, then I saw the KSP (gmres) residuals are different from different runs. Is it true of randomness from matsolverpackage? Thanks a lot!