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>
