I see, but the call to solve also calls set_operator, so the code on the Q&A should work as far as I can tell? Why is "different_nonzero_pattern" faster than "same”? It used to (and should) be the other way around.
Mikael > On 21 Mar 2015, at 20:09, Garth N. Wells <[email protected]> wrote: > > Call ‘set_operator' when you want the preconditioner to be updated. > > Garth > > >> On 21 Mar 2015, at 18:54, Mikael Mortensen <[email protected]> >> wrote: >> >> Hi, >> >> Seems like there is a problem with preconditioner reuse in dolfin 1.5: >> http://fenicsproject.org/qa/6817/reuse-of-krylov-preconditioner?show=6838#c6838 >> >> >> Is this a bug? I can only see this issue that is marked as resolved: >> https://bitbucket.org/fenics-project/dolfin/issue/269/update-logic-on-petsc-preconditioner-re >> . The function set_petsc_operators that used to set the preconditioner >> structure is commented out in PETScKrylovSolver.cpp... >> >> Regards >> >> Mikael >> >> >> _______________________________________________ >> fenics mailing list >> [email protected] >> http://fenicsproject.org/mailman/listinfo/fenics > _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
