Thank you. How do you handle the sparse matrix?<br><br><div class="gmail_quote">On Mon, Aug 20, 2012 at 12:24 PM, Gong Ding <span dir="ltr"><<a href="mailto:gdiso@ustc.edu" target="_blank">gdiso@ustc.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>My code genius (2D/3D unstructured nonlinear solver for semiconductor simulation ) use adol-c for quite a long time. This AD tool is implemented by c++ operator overloading. Using an AD tool makes the development of Jacobian much much easier than handwriting. AD tool usually runs 2x slower than handwriting code, but build the Jacobian matrix is usually not a bottleneck for large scale nonlinear problems. <br>
 <br></div><div>You can find some information and source code here. <br><a href="https://github.com/cogenda/Genius-TCAD-Open" target="_blank">https://github.com/cogenda/Genius-TCAD-Open</a><br>The AD tool is just a head file include/math/adolc.h (stable) and its improved version include/math/adsmt.h (not fully tested).<br>
 <br><br>BTW,implement AD by operator overloading is also possible in python. <br></div><div class="HOEnZb"><div class="h5"><div> </div><blockquote style="padding-left:5px;margin-left:5px;border-left:#a0c6e5 2px solid;margin-right:0px">
<br>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>

<br>Thanks,<br>BehZad<br>
</blockquote><span></span></div></div></blockquote></div><br>