For a one-dimensional unit domain (0,1), if one needs to set non-homogenous Neumann BC at a particular end of the domain (say the left end, x=0),
Looking at step-7 and adapting it a little bit, does the following code look reasonable to mark the boundary ID of the left edge as '1' for later setting its flux? for (const auto &cell : triangulation.cell_iterators <https://dealii.org/developer/doxygen/deal.II/group__CPP11.html#ga7295c669f32e78c45446ddc236f19136> ()) for (const auto &face : cell->face_iterators()) { const auto center = face->center(); if (std::fabs(center(0)) < 1e-12) face->set_boundary_id(1); } This is based on my assumption that center(0) method of a face iterator shall return the actual co-ordinate itself (since in the 1D case, the face degenerates to a single point and its center is the same). On a related note, when we use the arrow operator, under which page should I look up the documentation. Searching for center() returns a whole lot of results in doxygen. -- 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/415017fc-c2ec-4463-b51d-feee4d7493b2%40googlegroups.com.
