Hi All,
Is there a way to delete cells from a triangulation, such that no
degrees of freedom are distributed to them by the DoFHandler? I have
a selectively refined mesh (with hanging nodes), but only need to
solve a boundary value problem on a subregion of the whole domain.
One option would be to constrain the unnecessary degrees of freedom,
but I would end up constraining about 50% of the degrees of freedom.
I also thought of creating a new triangulation based on the vertex and
cell information for the desired cells, but this does not seem to be
feasible if the triangulation has hanging nodes, because neighborship
information is lost along the way.
Instead, is there a simple way I can mark certain cells as inactive?
It seems that the usual definition of an inactive cell is a cell that
has refined children, but this would not be the case for my
situation. Ideally, it would be nice if one could pass a
FilteredIterator to the DoFHandler, so that degrees of freedom would
be assigned to the requested subset of cells, rather than to the
active set. Would it be difficult to add this additional
functionality (in a modular way), or is the degree of freedom handling
too tightly linked to the definition of the active set?
Thanks for any ideas you may have,
Josh
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
