> It happens when I call "distribute_dofs", because fe_system calls the
> following function:
>
> template <int dim>
> FiniteElementDomination::Domination
> FE_Nothing<dim> ::
> compare_for_face_domination (const FiniteElement<dim> & fe_other) const
> {
> if(dynamic_cast<const FE_Nothing<dim>*>(&fe_other) != 0)
> return FiniteElementDomination::either_element_can_dominate;
> else
> return FiniteElementDomination::this_element_dominates;
> }
>
>
> How should be the behavior of this function in the FE_Nothing case?
I tend to think that we should add another state to the
FiniteElementDomination enum, say no_requirements, that indicates that we
don't want to do anything on interfaces between FE_Nothing and anything else.
Then FESystem would also return no_requirements whenever each component of the
FESystem would return no_requirements.
Would that make sense?
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
