<div dir="ltr">Jed,<div><br></div><div>Today I fought for several hours against a strange behaviour occuring when reusing a TS for multiple solves. You can reproduce it by applying the attached patch.<div><br></div><div>The issue is with DMCoarsenHookAdd called in SNESSetUpMatrices, and git blames you ;-).</div><div><br></div><div>Every time you call TSGetRHSMats_Private, TSGetIJacobian gets called, which in turns calls SNESSetUpMatrices, that adds the entry in the linked list of hooks via DMCoarsenHookAdd. This causes a considerable slow down when using the RHS interface with an implicit solver, as you can see from running the patched code. <br><div><br></div><div>can you please take a look?</div><div>-- <br><div class="gmail_signature">Stefano</div>
</div></div></div></div>