<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-"><br>
>> have you thought about how to fix this?<br>
>><br>
>> <a href="https://lists.mcs.anl.gov/pipermail/petsc-dev/2017-October/021387.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>pipermail/petsc-dev/2017-<wbr>October/021387.html</a><br>
<br>
</span>I feel like I'm missing context.  I know I failed to follow up in a<br>
thread at some time in the past, but I don't really understand the<br>
problem or what is being attempted.  What is the failing test case?<br>
<span class="gmail-"><br></span></blockquote><div><br></div><div>There's no failing test; however, you can see considerable slow downs when calling TSSolve on the same ts multiple times. Just copying and paste the original mail message below: it seems clear to me the code path. However, what it's not clear is why we need to call DMCoarsenHookAdd so many times. Cannot be set once? </div><div><br></div><div>Every time you call TSGetRHSMats_Private, TSGetIJacobian gets called,<br></div><div><div>which in turns calls SNESSetUpMatrices, that adds the entry in the linked</div><div>list of hooks via DMCoarsenHookAdd. This causes a considerable slow down</div><div>when using the RHS interface with an implicit solver, as you can see from</div><div>running the patched code.</div></div><div><br></div></div>
</div></div>