On 08/24/2017 09:11 AM, Praveen C wrote:
I did this a while back, see here
https://bitbucket.org/cpraveen/deal_ii/src/865f1c963f1089a870054bebf7180c9712e93da8/step-15-sacado/?at=master
Ah, most excellent!
I suspect that this can be written even more concisely with the current
development version of deal.II. For example, I'm pretty sure you can write
std::vector<Tensor<1,dim,>Sacado::Fad::DFad<double>>
grad_u (n_q_points);
fe_values.get_function_gradients (solution, grad_u);
or something similar in place of what you do in lines 199, 207-217. Do
you want to see if something like this works now? I think it would be
excellent to have an example like this in step-15.
@jppelteret -- as the AD master, do you have any other suggestion?
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
