On Sat, Mar 21, 2015 at 7:32 PM, Mikael Mortensen <[email protected]> wrote: > 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. >
The difference (2.43 vs 2.32) is too small to be significant, at least without averaging over a number of runs. The options "different_nonzero_pattern" and "same” were removed from PETS 3.5. Garth > 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
