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

Protected Member Functions

void build_quad_mesh (unsigned int n_refinements=0)
 
void init ()
 
void clear ()
 
void setup_coupling_matrix (std::unique_ptr< CouplingMatrix > &coupling)
 

Protected Attributes

std::unique_ptr< MeshBase_mesh
 
std::unique_ptr< EquationSystems_es
 

Detailed Description

Definition at line 230 of file overlapping_coupling_test.C.

Member Function Documentation

◆ build_quad_mesh()

void OverlappingTestBase::build_quad_mesh ( unsigned int  n_refinements = 0)
inlineprotected

Definition at line 238 of file overlapping_coupling_test.C.

References libMesh::Elem::build_with_id(), libMesh::QUAD4, libMesh::Elem::set_node(), libMesh::Elem::subdomain_id(), and libMesh::MeshRefinement::uniformly_refine().

239  {
240  // We are making assumptions in various places about the presence
241  // of the elements on the current processor so we're restricting to
242  // ReplicatedMesh for now.
243  _mesh = std::make_unique<ReplicatedMesh>(*TestCommWorld);
244 
245  _mesh->set_mesh_dimension(2);
246 
247  _mesh->add_point( Point(0.0,0.0),0 );
248  _mesh->add_point( Point(1.0,0.0),1 );
249  _mesh->add_point( Point(1.0,1.0),2 );
250  _mesh->add_point( Point(0.0,1.0),3 );
251 
252  {
253  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 0));
254  elem->subdomain_id() = 1;
255 
256  for (unsigned int n=0; n<4; n++)
257  elem->set_node(n) = _mesh->node_ptr(n);
258  }
259 
260  _mesh->add_point( Point(1.0,2.0),4 );
261  _mesh->add_point( Point(0.0,2.0),5 );
262 
263  {
264  Elem * elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 1));
265  elem->subdomain_id() = 1;
266 
267  elem->set_node(0) = _mesh->node_ptr(3);
268  elem->set_node(1) = _mesh->node_ptr(2);
269  elem->set_node(2) = _mesh->node_ptr(4);
270  elem->set_node(3) = _mesh->node_ptr(5);
271  }
272 
273  _mesh->add_point( Point(0.0,0.0),6 );
274  _mesh->add_point( Point(1.0,0.0),7 );
275  _mesh->add_point( Point(1.0,2.0),8 );
276  _mesh->add_point( Point(0.0,2.0),9 );
277 
278  {
279  Elem* elem = _mesh->add_elem(Elem::build_with_id(QUAD4, 2));
280  elem->subdomain_id() = 2;
281 
282  elem->set_node(0) = _mesh->node_ptr(6);
283  elem->set_node(1) = _mesh->node_ptr(7);
284  elem->set_node(2) = _mesh->node_ptr(8);
285  elem->set_node(3) = _mesh->node_ptr(9);
286  }
287 
288  _mesh->partitioner() = std::make_unique<OverlappingTestPartitioner>();
289 
290  _mesh->prepare_for_use();
291 
292 #ifdef LIBMESH_ENABLE_AMR
293  if (n_refinements > 0)
294  {
295  MeshRefinement refine(*_mesh);
296  refine.uniformly_refine(n_refinements);
297  }
298 #endif // LIBMESH_ENABLE_AMR
299  }
virtual Node *& set_node(const unsigned int i)
Definition: elem.h:2381
This is the base class from which all geometric element types are derived.
Definition: elem.h:94
std::unique_ptr< MeshBase > _mesh
Implements (adaptive) mesh refinement algorithms for a MeshBase.
subdomain_id_type subdomain_id() const
Definition: elem.h:2391
A Point defines a location in LIBMESH_DIM dimensional Real space.
Definition: point.h:39

◆ clear()

void OverlappingTestBase::clear ( )
inlineprotected

Definition at line 321 of file overlapping_coupling_test.C.

322  {
323  _es.reset();
324  _mesh.reset();
325  }
std::unique_ptr< MeshBase > _mesh
std::unique_ptr< EquationSystems > _es

◆ init()

void OverlappingTestBase::init ( )
inlineprotected

Definition at line 301 of file overlapping_coupling_test.C.

References libMesh::System::add_variable(), libMesh::FIRST, and libMesh::LAGRANGE.

302  {
303  _es = std::make_unique<EquationSystems>(*_mesh);
304  LinearImplicitSystem & sys = _es->add_system<LinearImplicitSystem> ("SimpleSystem");
305 
306  std::set<subdomain_id_type> sub_one;
307  sub_one.insert(1);
308 
309  std::set<subdomain_id_type> sub_two;
310  sub_two.insert(2);
311 
312  sys.add_variable("U", FIRST, LAGRANGE, &sub_two);
313  sys.add_variable("L", FIRST, LAGRANGE, &sub_two);
314 
315  sys.add_variable("V", FIRST, LAGRANGE, &sub_one);
316  sys.add_variable("p", FIRST, LAGRANGE, &sub_one);
317 
318  _es->init();
319  }
Manages consistently variables, degrees of freedom, coefficient vectors, matrices and linear solvers ...
std::unique_ptr< EquationSystems > _es
unsigned int add_variable(std::string_view var, const FEType &type, const std::set< subdomain_id_type > *const active_subdomains=nullptr)
Adds the variable var to the list of variables for this system.
Definition: system.C:1305

◆ setup_coupling_matrix()

void OverlappingTestBase::setup_coupling_matrix ( std::unique_ptr< CouplingMatrix > &  coupling)
inlineprotected

Definition at line 327 of file overlapping_coupling_test.C.

References libMesh::System::n_vars(), and libMesh::System::variable_number().

328  {
329  System & system = _es->get_system("SimpleSystem");
330 
331  coupling = std::make_unique<CouplingMatrix>(system.n_vars());
332 
333  const unsigned int u_var = system.variable_number("U");
334  const unsigned int l_var = system.variable_number("L");
335  const unsigned int v_var = system.variable_number("V");
336  const unsigned int p_var = system.variable_number("p");
337 
338  // Only adding the overlapping couplings since the primary ones should
339  // be there by default.
340  (*coupling)(u_var,v_var) = true;
341  (*coupling)(l_var,v_var) = true;
342  (*coupling)(l_var,p_var) = true;
343  (*coupling)(v_var,u_var) = true;
344  (*coupling)(v_var,l_var) = true;
345  }
unsigned int variable_number(std::string_view var) const
Definition: system.C:1557
std::unique_ptr< EquationSystems > _es
Manages consistently variables, degrees of freedom, and coefficient vectors.
Definition: system.h:96
unsigned int n_vars() const
Definition: system.h:2349

Member Data Documentation

◆ _es

std::unique_ptr<EquationSystems> OverlappingTestBase::_es
protected

Definition at line 236 of file overlapping_coupling_test.C.

◆ _mesh

std::unique_ptr<MeshBase> OverlappingTestBase::_mesh
protected

Definition at line 234 of file overlapping_coupling_test.C.


The documentation for this class was generated from the following file: