libMesh
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
GetBoundaryPointsSecondTest Class Reference
Inheritance diagram for GetBoundaryPointsSecondTest:
[legend]

Public Member Functions

 LIBMESH_CPPUNIT_TEST_SUITE (GetBoundaryPointsSecondTest)
 The goal of this test is the same as the previous, but now we use the second order mesh. More...
 
 CPPUNIT_TEST (testMesh)
 
 CPPUNIT_TEST_SUITE_END ()
 
void setUp ()
 
void testMesh ()
 
 LIBMESH_CPPUNIT_TEST_SUITE (GetBoundaryPointsTest)
 The goal of this test is to ensure that ReplicatedMesh::get_boundary_points returns a vector of boundary points enclosing the domain. More...
 
void tearDown ()
 

Protected Member Functions

void build_mesh ()
 

Protected Attributes

std::unique_ptr< ReplicatedMesh > _mesh
 

Detailed Description

Definition at line 233 of file boundary_points.C.

Member Function Documentation

◆ build_mesh()

void GetBoundaryPointsTest::build_mesh ( )
inlineprotectedinherited

Definition at line 30 of file boundary_points.C.

References libMesh::Elem::build_with_id(), libMesh::QUAD4, libMesh::Elem::set_node(), and TestCommWorld.

31  {
32  _mesh = std::make_unique<ReplicatedMesh>(*TestCommWorld, 2);
33 
34  // (0,1) (1,1) (2,1) (3,1) (4,1) (5,1) (6,1)
35  // x---------------x---------------x x-----------x-----------x-----------x
36  // | | | | | | |
37  // | | | | | | |
38  // | | | x-----------x-----------x-----------x
39  // | | | |(3,0.5) | void | |(6,0.5)
40  // | | | | | | |
41  // x---------------x---------------x x-----------x-----------x-----------x
42  // (0,0) (1,0) (2,0) | (3,0) | | |(6,0)
43  // | | | | | | |
44  // | | | | | | |
45  // | | | | | | |
46  // | | | | | | |
47  // x---------------x---------------x x-----------x-----------x-----------x
48  // (0,-1) (1,-1) (2,-1) (3,-1) (6,-1)
49 
50  _mesh->add_point( Point(0.0, 0.0), 0 );
51  _mesh->add_point( Point(1.0, 0.0), 1 );
52  _mesh->add_point( Point(1.0, 1.0), 2 );
53  _mesh->add_point( Point(0.0, 1.0), 3 );
54  _mesh->add_point( Point(0.0,-1.0), 4 );
55  _mesh->add_point( Point(1.0,-1.0), 5 );
56  _mesh->add_point( Point(2.0, 0.0), 6 );
57  _mesh->add_point( Point(2.0, 1.0), 7 );
58  _mesh->add_point( Point(2.0,-1.0), 8 );
59 
60  _mesh->add_point( Point(3.0,-1.0), 9 );
61  _mesh->add_point( Point(3.0, 0.0), 10 );
62  _mesh->add_point( Point(3.0, 0.5), 11 );
63  _mesh->add_point( Point(3.0, 1.0), 12 );
64  _mesh->add_point( Point(4.0,-1.0), 13 );
65  _mesh->add_point( Point(4.0, 0.0), 14 );
66  _mesh->add_point( Point(4.0, 0.5), 15 );
67  _mesh->add_point( Point(4.0, 1.0), 16 );
68  _mesh->add_point( Point(5.0,-1.0), 17 );
69  _mesh->add_point( Point(5.0, 0.0), 18 );
70  _mesh->add_point( Point(5.0, 0.5), 19 );
71  _mesh->add_point( Point(5.0, 1.0), 20 );
72  _mesh->add_point( Point(6.0,-1.0), 21 );
73  _mesh->add_point( Point(6.0, 0.0), 22 );
74  _mesh->add_point( Point(6.0, 0.5), 23 );
75  _mesh->add_point( Point(6.0, 1.0), 24 );
76 
77  {
78  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 0));
79  elem->set_node(0) = _mesh->node_ptr(0);
80  elem->set_node(1) = _mesh->node_ptr(1);
81  elem->set_node(2) = _mesh->node_ptr(2);
82  elem->set_node(3) = _mesh->node_ptr(3);
83  }
84  {
85  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 1));
86  elem->set_node(0) = _mesh->node_ptr(4);
87  elem->set_node(1) = _mesh->node_ptr(5);
88  elem->set_node(2) = _mesh->node_ptr(1);
89  elem->set_node(3) = _mesh->node_ptr(0);
90  }
91  {
92  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 2));
93  elem->set_node(0) = _mesh->node_ptr(7);
94  elem->set_node(1) = _mesh->node_ptr(2);
95  elem->set_node(2) = _mesh->node_ptr(1);
96  elem->set_node(3) = _mesh->node_ptr(6);
97  }
98  {
99  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 3));
100  elem->set_node(0) = _mesh->node_ptr(5);
101  elem->set_node(1) = _mesh->node_ptr(8);
102  elem->set_node(2) = _mesh->node_ptr(6);
103  elem->set_node(3) = _mesh->node_ptr(1);
104  }
105  {
106  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 4));
107  elem->set_node(0) = _mesh->node_ptr(9);
108  elem->set_node(1) = _mesh->node_ptr(13);
109  elem->set_node(2) = _mesh->node_ptr(14);
110  elem->set_node(3) = _mesh->node_ptr(10);
111  }
112  {
113  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 5));
114  elem->set_node(0) = _mesh->node_ptr(10);
115  elem->set_node(1) = _mesh->node_ptr(14);
116  elem->set_node(2) = _mesh->node_ptr(15);
117  elem->set_node(3) = _mesh->node_ptr(11);
118  }
119  {
120  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 6));
121  elem->set_node(0) = _mesh->node_ptr(11);
122  elem->set_node(1) = _mesh->node_ptr(15);
123  elem->set_node(2) = _mesh->node_ptr(16);
124  elem->set_node(3) = _mesh->node_ptr(12);
125  }
126  {
127  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 7));
128  elem->set_node(0) = _mesh->node_ptr(13);
129  elem->set_node(1) = _mesh->node_ptr(17);
130  elem->set_node(2) = _mesh->node_ptr(18);
131  elem->set_node(3) = _mesh->node_ptr(14);
132  }
133  // skip one element here
134  {
135  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 8));
136  elem->set_node(0) = _mesh->node_ptr(15);
137  elem->set_node(1) = _mesh->node_ptr(19);
138  elem->set_node(2) = _mesh->node_ptr(20);
139  elem->set_node(3) = _mesh->node_ptr(16);
140  }
141  {
142  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 9));
143  elem->set_node(0) = _mesh->node_ptr(17);
144  elem->set_node(1) = _mesh->node_ptr(21);
145  elem->set_node(2) = _mesh->node_ptr(22);
146  elem->set_node(3) = _mesh->node_ptr(18);
147  }
148  {
149  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 10));
150  elem->set_node(0) = _mesh->node_ptr(18);
151  elem->set_node(1) = _mesh->node_ptr(22);
152  elem->set_node(2) = _mesh->node_ptr(23);
153  elem->set_node(3) = _mesh->node_ptr(19);
154  }
155  {
156  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 11));
157  elem->set_node(0) = _mesh->node_ptr(19);
158  elem->set_node(1) = _mesh->node_ptr(23);
159  elem->set_node(2) = _mesh->node_ptr(24);
160  elem->set_node(3) = _mesh->node_ptr(20);
161  }
162 
163  // libMesh shouldn't renumber, or our based-on-initial-id
164  // assertions later may fail.
165  _mesh->allow_renumbering(false);
166 
167  _mesh->prepare_for_use();
168  }
libMesh::Elem::set_node
virtual Node *& set_node(const unsigned int i)
Definition: elem.h:2381
TestCommWorld
libMesh::Parallel::Communicator * TestCommWorld
Definition: driver.C:159
GetBoundaryPointsTest::_mesh
std::unique_ptr< ReplicatedMesh > _mesh
Definition: boundary_points.C:28
libMesh::Elem
This is the base class from which all geometric element types are derived.
Definition: elem.h:94
libMesh::Point
A Point defines a location in LIBMESH_DIM dimensional Real space.
Definition: point.h:39

◆ CPPUNIT_TEST()

GetBoundaryPointsSecondTest::CPPUNIT_TEST ( testMesh  )

◆ CPPUNIT_TEST_SUITE_END()

GetBoundaryPointsSecondTest::CPPUNIT_TEST_SUITE_END ( )

◆ LIBMESH_CPPUNIT_TEST_SUITE() [1/2]

GetBoundaryPointsTest::LIBMESH_CPPUNIT_TEST_SUITE ( GetBoundaryPointsTest  )
inherited

The goal of this test is to ensure that ReplicatedMesh::get_boundary_points returns a vector of boundary points enclosing the domain.

◆ LIBMESH_CPPUNIT_TEST_SUITE() [2/2]

GetBoundaryPointsSecondTest::LIBMESH_CPPUNIT_TEST_SUITE ( GetBoundaryPointsSecondTest  )

The goal of this test is the same as the previous, but now we use the second order mesh.

◆ setUp()

void GetBoundaryPointsSecondTest::setUp ( )
inline

Definition at line 250 of file boundary_points.C.

251  {
252 #if LIBMESH_DIM > 1
253  this->build_mesh();
254  _mesh->all_second_order(true);
255 #endif
256  }
GetBoundaryPointsTest::_mesh
std::unique_ptr< ReplicatedMesh > _mesh
Definition: boundary_points.C:28

◆ tearDown()

void GetBoundaryPointsTest::tearDown ( )
inlineinherited

Definition at line 178 of file boundary_points.C.

178 {}

◆ testMesh()

void GetBoundaryPointsSecondTest::testMesh ( )
inline

Definition at line 258 of file boundary_points.C.

259  {
260  LOG_UNIT_TEST;
261 
262  const auto points = _mesh->get_boundary_points();
263 
264  // There'd better be two disconnected subdomains
265  CPPUNIT_ASSERT_EQUAL( (std::size_t)2, points.size() );
266 
267  // The first key should better be element 0
268  auto it = points.find(0);
269  CPPUNIT_ASSERT(it != points.end());
270 
271  // There'd better be one boundary in the first subdomain
272  CPPUNIT_ASSERT_EQUAL( (std::size_t)1, it->second.size() );
273 
274  // There'd better be 16 points on the boundary of the first subdomain
275  CPPUNIT_ASSERT_EQUAL( (std::size_t)16, it->second[0].size() );
276 
277  // Boundary points should better start from (1,1)
278  CPPUNIT_ASSERT_EQUAL( Point(1,1), it->second[0][0] );
279 
280  // Boundary points should better include side middle points including (0.5,1)
281  CPPUNIT_ASSERT_EQUAL( Point(0.5,1), it->second[0][1] );
282 
283  // Boundary points should better end with (1.5,1)
284  CPPUNIT_ASSERT_EQUAL( Point(1.5,1), it->second[0][15] );
285 
286  // The second key should better be element 4
287  it = points.find(4);
288  CPPUNIT_ASSERT(it != points.end());
289 
290  // There'd better be two boundaries in the second subdomain due to the middle hole
291  CPPUNIT_ASSERT_EQUAL( (std::size_t)2, it->second.size() );
292 
293  // There'd better be 24 points on the first (outer) boundary of the second subdomain
294  CPPUNIT_ASSERT_EQUAL( (std::size_t)24, it->second[0].size() );
295 
296  // There'd better be 8 points on the first (inner) boundary of the second subdomain
297  CPPUNIT_ASSERT_EQUAL( (std::size_t)8, it->second[1].size() );
298 
299  // The first boundary points should better start from (3,-1)
300  CPPUNIT_ASSERT_EQUAL( Point(3,-1), it->second[0][0] );
301 
302  // The first boundary points should better include side middle points including (3.5,-1)
303  CPPUNIT_ASSERT_EQUAL( Point(3.5,-1), it->second[0][1] );
304 
305  // The first boundary points should better end with (3,-0.5)
306  CPPUNIT_ASSERT_EQUAL( Point(3,-0.5), it->second[0][23] );
307 
308  // The second boundary points should better start from (4,0)
309  CPPUNIT_ASSERT_EQUAL( Point(4,0), it->second[1][0] );
310 
311  // The second boundary points should better include side middle points including (4,0.25)
312  CPPUNIT_ASSERT_EQUAL( Point(4,0.25), it->second[1][1] );
313 
314  // The second boundary points should better end with (4.5,0)
315  CPPUNIT_ASSERT_EQUAL( Point(4.5,0), it->second[1][7] );
316  }
GetBoundaryPointsTest::_mesh
std::unique_ptr< ReplicatedMesh > _mesh
Definition: boundary_points.C:28
libMesh::Point
A Point defines a location in LIBMESH_DIM dimensional Real space.
Definition: point.h:39

Member Data Documentation

◆ _mesh

std::unique_ptr<ReplicatedMesh> GetBoundaryPointsTest::_mesh
protectedinherited

Definition at line 28 of file boundary_points.C.

The documentation for this class was generated from the following file:
Generated on Thu Apr 18 2024 15:18:01 for libMesh by   doxygen 1.8.14