David,
the code snippet you posted is not going to work. You can't use a FEValues
without initializing it with a cell_iterator.
If all the information you need for your postprocessing is contained in a
data vector based on a single DoFHandler, then DataPostprocessor is what
you want to use.
If
Hi,
Well, I found that the FiniteElement base class itself has a method
Tensor< 1, dim > FiniteElement< dim, spacedim >::shape_grad (const unsigned
int i, const Point ) const
which can compute shape gradient values at a given point, but the result it
get is on the reference cell, not the