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

[behzad baghapour]

Thank you. How do you handle the sparse matrix?

On Mon, Aug 20, 2012 at 12:24 PM, Gong Ding <gdiso at ustc.edu> wrote:

> 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/863513e8/attachment.html>