[petsc-users] How to use AD for unstructured grid

[Gong Ding]

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.
 
You can find some information and source code here.
https://github.com/cogenda/Genius-TCAD-Open
The AD tool is just a head file include/math/adolc.h (stable) and its improved version include/math/adsmt.h (not fully tested).
 
BTW,implement AD by operator overloading is also possible in python.
 
Dear developers,
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?
Thanks,
BehZad
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120820/975d0157/attachment.html>