Hi, I am soon going to be turning from the idealized problem I have been working with for the last 2 years to one with more realistic geometry and parameters, and it might call parts of the domain to be "masked out" - based on some criterion certain dofs (possibly large regions of them) will be taken out of the FE solve and the relevant quantities will be solved for by more locally-defined relations.
If I still want to use a rectangular grid (and I think I do, rather than creating a custom mesh), one of the ways I can think of to handle this is to use a ConstraintMatrix - basically set the DoF values to zero in those masked regions. (I haven't completely thought out how I will do it yet.) But this seems like an awkward way to do this, if the number of entries in the ConstraintMatrix is on the order of the number of DoFs, right? Is there some cell-wise property in the triangulation that can be used so that the above step isn't necessary - so that, when distribute_dofs() is called, it knows that certain dofs on certain cells will not be part of the calculation, and so the distribution, number, and ordering of dofs is calculated accordingly? Thanks very much Dan _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
