[Libmesh-devel] solve_with_same_matrix?

2015-04-06 Thread David Knezevic
PETSc will automatically avoid recomputing preconditioners if you call solve multiple times with the same matrix and different right-hand-sides, but we don't currently get this benefit because in PetscLinearSolver we re-initialize the solver's matrix on each call to solve. I tested out adding solv

Re: [Libmesh-devel] solve_with_same_matrix?

2015-04-06 Thread Roy Stogner
On Mon, 6 Apr 2015, David Knezevic wrote: > PETSc will automatically avoid recomputing preconditioners if you call solve > multiple times with the same matrix and different right-hand-sides, > but we don't currently get this benefit because in PetscLinearSolver we > re-initialize the solver's m

Re: [Libmesh-devel] solve_with_same_matrix?

2015-04-06 Thread David Knezevic
On Mon, Apr 6, 2015 at 11:10 AM, Roy Stogner wrote: > > On Mon, 6 Apr 2015, David Knezevic wrote: > > PETSc will automatically avoid recomputing preconditioners if you call >> solve multiple times with the same matrix and different right-hand-sides, >> but we don't currently get this benefit bec