<div dir="ltr">HI Xiangdong<div><br></div><div>We have use the Tapenade AD package to compute sparse matrix jacobians for CFD problems. We use the result from the AD to populate a PETSc sparse matrix and then use PETSc to solve the resulting linear system. </div><div><br></div><div>Hope that helps.</div><div><br></div><div>Gaetan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 1, 2014 at 4:23 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@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="HOEnZb"><div class="h5"><br>
On Oct 1, 2014, at 3:16 PM, Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br>
<br>
> Xiangdong <<a href="mailto:epscodes@gmail.com">epscodes@gmail.com</a>> writes:<br>
><br>
>> Is there any reason that AD is not supported in PETSc? Is it because the<br>
>> ADIFOR/ADIC library itself?<br>
><br>
> Those libraries are not maintained and have license problems.  You're<br>
> welcome to use an AD tool yourself.  The support in PETSc was mostly<br>
> just making simple examples cleaner; it wasn't widely used in real<br>
> applications.<br>
<br>
</div></div>   Wasn’t used at all in anything but trivial codes.<br>
<br>
<br>
<br>
</blockquote></div><br></div>