OK. Does any of them consider a sparse function differentiation or one should define a dense matrix to recieve the result? If so, then is there any possiblity to face with  memory allocation problem in Large scale simulations?<br>
<br><br><br><div class="gmail_quote">On Mon, Aug 20, 2012 at 8:55 AM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
This is a choice for you to make. Operator overloading is usually less intrusive to use, but also more limited in capability.<div class="HOEnZb"><div class="h5"><br><br><div class="gmail_quote">On Mon, Aug 20, 2012 at 12:24 AM, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com" target="_blank">behzad.baghapour@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Jed,<br><br>Which AD tool you suggest in the case of fast evaluations? Operator Overlapping types like ADOL-C or source code transformation types like ADIC?<br>

<br>In addition, I prefer the minimum possible change to the my code structure :-)<br>
<br>Thanks a lot,<br>BehZad<div><div><br><br><div class="gmail_quote">On Sun, Aug 19, 2012 at 10:06 PM, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com" target="_blank">behzad.baghapour@gmail.com</a>></span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">OK. Thanks. I should try it.<div><div><br><br><div class="gmail_quote">On Sun, Aug 19, 2012 at 9:56 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>



<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote"><div>On Sun, Aug 19, 2012 at 1:19 PM, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com" target="_blank">behzad.baghapour@gmail.com</a>></span> wrote:<br>



<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear developers,<br><br>I need to evaluate Jacobian matrix using an AD tool. I read from the Petsc manual that it can only work with structured grid via DMMG interface. Is there a possible way to use it for unstructured grid and set it for KSP structure?<br>




</blockquote><div><br></div></div><div>ADIC has not been supported for a long time so the "automatic" AD support in PETSc is no longer supported either. Unfortunately, the successors to ADIC, while more "capable", are much more involved to build and to use.</div>




<div><br></div><div>You have always been able to use any AD tool internally to define the Jacobian (in assembled or unassembled form).</div></div>
</blockquote></div><br>
</div></div></blockquote></div><br>
</div></div></blockquote></div><br>
</div></div></blockquote></div><br>