<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Jul 7, 2014 at 5:36 PM, Tobin Isaac <span dir="ltr"><<a href="mailto:tisaac@ices.utexas.edu" target="_blank">tisaac@ices.utexas.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Hi,<br>
<br>
If I have a pointwise Jacobian function f, I know that I can call<br>
DMGetDS() and pass f to PetscDSSetJacobian(), and that f will be used<br>
by PetscFEIntegrate() and thus by DMPlexSNESComputeJacobianFEM().  It<br>
looks like PetscFEIntegrate() is only used by plex and not da.  Is<br>
there anyway that I can use f to construct a Jacobian with a da<br>
object?<br></blockquote><div><br></div><div>There is a completely untested way. My idea is that the DA is able to give back</div><div>cell closures the same way that Plex does, making all the FEM code work just</div><div>
as it does in the unstructured case. There are a few toy tests of this functionality</div><div>in src/dm/impls/plex/examples/tests/ex3.c. I have not made a thoroughgoing</div><div>example which switches between Plex and DA, but it should be possible. Let</div>
<div>me know if it fails.</div><div><br></div><div>All this DA code which does not assume a vertex discretization is activated by</div><div>giving it a default section.</div><div><br></div><div>  Thanks,</div><div><br></div>
<div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  Toby<br></blockquote></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener
</div></div>