[deal.II] Re: Specifying internal "boundary" conditions

Hi Melanie, we also have a few examples in ASPECT that probably do what you are aiming for. E.g. the cookbook prescribed_velocity (https://github.com/geodynamics/aspect/blob/master/cookbooks/prescribed_velocity/prescribed_velocity.cc#L295) uses the `add_line` and `set_inhomogeneity` functions t

[deal.II] Re: Specifying internal "boundary" conditions

Dear Bruno, Thank you, that helps a lot. Melanie On Wednesday, September 18, 2019 at 1:58:39 PM UTC-4, Bruno Turcksin wrote: > > Mélanie, > > I am not sure I understand exactly what you want to do but often people > who want internal boundary conditions, really just want to impose some > const

[deal.II] Re: Specifying internal "boundary" conditions

Mélanie, I am not sure I understand exactly what you want to do but often people who want internal boundary conditions, really just want to impose some constraints to a few dofs. This can be done using AffineConstraints . I