On 5/16/19 2:41 PM, 'Maxi Miller' via deal.II User Group wrote: > If I understand that correctly, the correct order for doing that would be: > -> Interpolate boundary condition without any component mask, giving me a > vector > -> Loop over all elements in the vector, and cross-check with the mask (using > extract_dofs()) if they should be interpolated in that way > -> If yes, copy them to a second vector > -> Repeat with the next boundary condition and new component mask > -> After finishing: Replace first vector with second vector (or use a > temporary vector as first vector) > Is that correct? > In addition: How do I treat different boundary markings, i.e. > dirichlet-conditions on 0, neumann-conditions on 1, and similar stuff?
Something of this sort. Of course, if you had sorted your DoFs by component (for example if you wanted to put things into a block vector), then the check whether a particular DoF belongs to a particular component becomes pretty trivial. Best W> -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/6b4003a6-04cf-9d51-6383-72034f9818ed%40colostate.edu. For more options, visit https://groups.google.com/d/optout.