<div class="gmail_quote">On Wed, Oct 12, 2011 at 01:10, Nathan Collier <span dir="ltr"><<a href="mailto:nathaniel.collier@gmail.com">nathaniel.collier@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Jed,<br>
<br>
I would talk to Lisandro about this. In our PETSC-nurbs code he<br>
implemented a few additional routines that let us approximate<br>
jacobians by coloring. I am not sure of the difference between this<br>
and what you want but maybe he has it done. The extra routines are in<br>
fdcolor.c and the usage is in NavierStokesKorteweg2D.c. Basically we<br>
have to get the SNES out of the TS and then set a few things.<br></blockquote><div><br></div><div>Thanks Nathan. I should clarify that my question was not about how the user can set this up (indeed, I made it work this way and there are a bunch of examples now), but rather about how to make it happen automatically and especially, how to compose correctly with geometric multigrid.</div>
<div><br></div></div>