On 6/25/19 11:55 AM, Daniel Garcia-Sanchez wrote: > According to Alexander Grayver, deal.II interface does not permit dynamic > constraints (that is, constraints, which depend on a cell orientation). This > problem does not seem easy to solve.
That's only partly true. The DoFTools::make_hanging_node_constraints() class, and the interface it uses to the FiniteElement class, can not deal with such constraints. It assumes that the constraints have the same form on each cell with a hanging node. But one doesn't have to go through this function. One could just as well write a function that is specific to the FE_NedelecSZ element and that computes these kinds of constraints and puts them into AffineConstraints. I've done that for other elements before, though I'm not sure there is an example in deal.II itself. I'd be happy to provide my code as an example, though. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/990c3f87-454c-8fb2-18b6-03c256492da8%40colostate.edu. For more options, visit https://groups.google.com/d/optout.
