// Create the objects that describe the boundaries and attach them
// to the triangulation as the ones to use on boundaries marked
- // with boundary indicators 8 and 9
+ // with boundary indicators 8 and 9
const double inner_radius = 1.5;
const double outer_radius = 2.5;
// Then loop over all faces of the domain and, if for the position
// of the center of a face the following holds then set boundary
// indicators:
- // - if y>3 and z<=2.5 or z>=5 then use boundary indicator 8
- // - if y>3 and 2.5<=z<=5 then use boundary indicator 9
+ // - if y>3 and z<=2.5 or z>=5 then use boundary indicator 8
+ // - if y>3 and 2.5<=z<=5 then use boundary indicator 9
typename Triangulation<3>::active_cell_iterator
cell = triangulation.begin_active(),
endc = triangulation.end();
if (cell->face(f)->at_boundary())
{
- if ((face_center[2] <= 2.5 || face_center[2] >= 5) &&
+ if ((face_center[2] <= 2.5 || face_center[2] >= 5) &&
face_center[1] >= 3)
cell->face(f)->set_boundary_indicator(8);
- if (face_center[2] >= 2.5 &&
+ if (face_center[2] >= 2.5 &&
face_center[2] <= 5
&& face_center[1] >= 3)
cell->face(f)->set_boundary_indicator(9);
Point<3>(0,0,1),
Point<3>(0,3,0));
triangulation.set_boundary (9, outer_cylinder);
-@encode
+@endcode
This yields an improvement, though it is still not quite correct:
src="http://www.dealii.org/images/steps/developer/step-49.yuhan.5.png"
alt="" height="200px">
-This is already better. However, something is still going wrong on the
+This is already better. However, something is still going wrong on the
front left face. On second look, we can also see that the faces where
the geometry widens have been refined at the bottom, that there is one
transition face that looks wrong because it has a triangle rather than
<img
src="http://www.dealii.org/images/steps/developer/step-49.yuhan.7.png"
alt="" height="200px">
-
+
So, yes!, this is finally what we were looking for!
-
\ No newline at end of file
+