Thomas, > The idea was to use a hp_collector like > > fe_1 ( FE_Q<dim>(1), 1, FE_Nothing<dim>(), 1 ), > fe_2 ( FE_Nothing<dim>(), 1, FE_Q<dim>(1), 1 ), > > fe_collection.push_back (fe_1); > fe_collection.push_back (fe_2);
This is exactly what I always thought the FENothing should be used for! > in this way we would have double degrees of freedom only on the interface, > BUT we have the following problem: > > The violated condition was: > domination != FiniteElementDomination::neither_element_dominates > > that makes sense in general, but not in my case, as I will never try to > restrict my FE_Q to FE_Nothing. Right. I think that the check in this place is the right thing to do in general, and that we should special case FE_Nothing somehow. Can you say where exactly this exception is produced? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
