Hello Markus,
Thank you. FEValuesExtractors/FEValuesViews look a bit complicated for
at the moment. I will try to discover them later though. Currently I do
this:
typedef std::pair< std::pair< unsigned int, unsigned int >, unsigned int
> BaseIndex;
BaseIndex base_i = fe.system_to_base_index(i);
BaseIndex base_j = fe.system_to_base_index(j);
if (base_i.first.second == base_j.first.second)
{
// Diagonal
}
else
{
// Off-diagonal
}
Is this correct?
On 27.04.2012 10:39, Markus Bürg wrote:
Hello Alexander,
FiniteElement::system_to_base_index is the equivalent for
non-primitive shape functions. However, another option would be to use
FEValuesExtractors/FEValuesViews. Then handling with shape functions
becomes much more natural again (as in the scalar case).
Best Regards,
Markus
--
Regards,
Alexander
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
