Set operators also for unrestricted solve & without preconditioner.

libMesh · Jun 30, 2023 · f2ddeeb · f2ddeeb
1 parent dc964ef
commit f2ddeeb
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/src/solvers/petsc_linear_solver.C b/src/solvers/petsc_linear_solver.C
@@ -632,6 +632,11 @@ PetscLinearSolver<T>::solve_base (SparseMatrix<T> * matrix,
           ierr = KSPSetOperators(_ksp, mat, const_cast<PetscMatrix<T> *>(precond)->mat());
           LIBMESH_CHKERR(ierr);
         }
+      else
+        {
+          ierr = KSPSetOperators(_ksp, mat, mat);
+          LIBMESH_CHKERR(ierr);
+        }
 
       if (this->_preconditioner)
         {