Re: [Libmesh-users] Recycling preconditioned in system_of_equations7.C

2016-11-15 Thread Roy Stogner
On Tue, 15 Nov 2016, Salazar De Troya, Miguel wrote: What would be a good place to insert PetscDiffSolver::init()? Maybe in DifferentiableSystem::init() after the user initialization?  I believe so. --- Roy--

Re: [Libmesh-users] Recycling preconditioned in system_of_equations7.C

2016-11-15 Thread Salazar De Troya, Miguel
What would be a good place to insert PetscDiffSolver::init()? Maybe in DifferentiableSystem::init() after the user initialization? Miguel From: "Kong, Fande" mailto:fande.k...@inl.gov>> Date: Tuesday, November 15, 2016 at 11:44 AM To: Roy Stogner mailto:royst...@ices.utexas.edu>> Cc: Miguel Salaz

Re: [Libmesh-users] Recycling preconditioned in system_of_equations7.C

2016-11-15 Thread Kong, Fande
In petsc_nonlinear_solver, there is the same issue. After each solve, SNES will get destroyed.Unfortunately, there are already a lot of applications built on this wrong interface. However, we should fix this in some point. Fande, On Tue, Nov 15, 2016 at 12:36 PM, Roy Stogner wrote: > > On

Re: [Libmesh-users] Recycling preconditioned in system_of_equations7.C

2016-11-15 Thread Roy Stogner
On Tue, 15 Nov 2016, Salazar De Troya, Miguel wrote: I am modifying PetscDiffSolve to be able to reuse the preconditioner. However, every time PetscDiffSolver::solve() is called, a new SNES is created by calling PetscDiffSolver::init(), why is this? I'm not sure, but if I had to bet money the

Re: [Libmesh-users] Recycling preconditioned in system_of_equations7.C

2016-11-15 Thread Salazar De Troya, Miguel
I am modifying PetscDiffSolve to be able to reuse the preconditioner. However, every time PetscDiffSolver::solve() is called, a new SNES is created by calling PetscDiffSolver::init(), why is this? I don’t understand this coupling. I believe the solver should be initialized elsewhere or maybe inc