Daniel,
Thanks again. To summarize, I want to implement the gradients and the
laplacian functions for the SymmetricTensor class.
The laplacian (trace of hessian) for the fe_values block extractor is
straightforward, but the gradients (of the solution field) for the
SymmetricTensor has some iss
Daniel,
thank you for the explanation. I understand that the components are in
space (x/y/z), and (*shape_gradient_ptr++)[d] would be derivative of phi in
each direction; is that correct?
(compared to vectors) the implementation of divergence function for the
tensors is as following;
Hi Wolfgang,
In the fe_values.cc, the function
"FEValuesViews::internal::do_function_derivatives" (for vectors), and in
the fe_values.h, the type functions "vector::gradient" and
"vector::symmetric_gradient"; are not clear to me because I'm not familiar
with the data structure.
Inside the fe_