I notice that in PetscLinearSolver::solve we have:
--
#if PETSC_RELEASE_LESS_THAN(3,5,0)
ierr = KSPSetOperators(_ksp, submat, subprecond,
this->same_preconditioner ?
SAME_PRECONDITIONER : DIFFERENT_NONZERO_PATTERN);
#else
ierr = KSPSet
Total guess:
http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPSetReusePreconditioner.html
Sounds promising though.
On Wed, Jan 7, 2015 at 11:16 PM, David Knezevic
wrote:
> I notice that in PetscLinearSolver::solve we have:
>
> --
>
> #if PETSC_RELEASE_LESS
If you want to skip the new factorization you can simply call
KSPSetReusePreconditioner() appropriate in the #else case below
Barry
> On Jan 7, 2015, at 10:16 PM, David Knezevic
> wrote:
>
> I notice that in PetscLinearSolver::solve we have:
>
> --
>
> #if PETSC
