Hello all, I am solving a saddle point problem like a nearly incompressible elasticity. The unknown field is u (displacement) and p (pressure). For Galerkin weak form, it is very convenient in deal.ii to use FEValuesExtractors::Vector or FEValuesExtractors::Scalar and use member function gradient/symmetric_gradient/divergence to quickly compute the 1st derivative related terms in the bilinear form.
However, I am trying to implement some stabilization technique that involves 2nd order derivatives of the shape function, like the divergence of the symmetric gradient of the displacement vector (actually, the shape function). For now, I can only think of selecting the components from the Hessian of the shape function to compute the div of the symmetric gradient. Is there any way in deal.ii to directly extract such values like divergence of the symmetric gradient of shape functions? If this is actually a quick question, a simple link to the reference or guide is also welcome. Thanks, Lixing -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/7a8fbef4-8ca0-4eb5-9cd2-8a6a9a29a1e9n%40googlegroups.com.
