Doug,
let me add another thought I had: While the vertex locations of the mesh are 
hard-coded as doubles (and you won't be able to change that), there is a class 
MappingEulerian that allows you to provide a displacement field. This field is 
parameterized by a vector of some sort, and it is possible that you can make 
it work if that field is a vector of some AD type for which you can thread the 
remainder of the computations through AD evaluations.

I don't know whether that's going to work out any easier. but I did want to 
point it out as an option worth reading through.

Best
  W.

-- 
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                            www: http://www.math.colostate.edu/~bangerth/

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to