66 const unsigned int n_elem_per_side = 5;
67 const std::unique_ptr<Elem> test_elem = Elem::build(elem_type);
68 const unsigned int ymax = test_elem->dim() > 1;
69 const unsigned int zmax = test_elem->dim() > 2;
70 const unsigned int ny = ymax * n_elem_per_side;
71 const unsigned int nz = zmax * n_elem_per_side;
85 std::set<dof_id_type> evaluable_elements;
87 for (
const auto & elem :
mesh.active_local_element_ptr_range())
90 evaluable_elements.insert(elem->id());
92 for (
unsigned int s1=0; s1 != elem->n_neighbors(); ++s1)
101 if (!evaluable_elements.count(n1->
id()))
104 evaluable_elements.insert(n1->
id());
107 for (
unsigned int s2=0; s2 != elem->n_neighbors(); ++s2)
114 if (!evaluable_elements.count(n2->
id()))
117 evaluable_elements.insert(n2->
id());
120 for (
unsigned int s3=0; s3 != elem->n_neighbors(); ++s3)
125 evaluable_elements.count(n3->
id()))
129 evaluable_elements.insert(n3->
id());
141 const std::size_t n_evaluable =
145 CPPUNIT_ASSERT_EQUAL(evaluable_elements.size(), n_evaluable);
This is the EquationSystems class.
libMesh::Parallel::Communicator * TestCommWorld
static constexpr Real TOLERANCE
This is the base class from which all geometric element types are derived.
DefaultCoupling & default_algebraic_ghosting()
Default algebraic ghosting functor.
Number point_value(unsigned int var, const Point &p, const bool insist_on_success=true, const NumericVector< Number > *sol=nullptr) const
Real distance(const Point &p)
void project_solution(FunctionBase< Number > *f, FunctionBase< Gradient > *g=nullptr) const
Projects arbitrary functions onto the current solution.
Number cubic_default_coupling_test(const Point &p, const Parameters &, const std::string &, const std::string &)
Manages consistently variables, degrees of freedom, and coefficient vectors.
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.
const Elem * neighbor_ptr(unsigned int i) const
void set_n_levels(unsigned int n_levels)
The Mesh class is a thin wrapper, around the ReplicatedMesh class by default.
processor_id_type processor_id() const
const DofMap & get_dof_map() const
processor_id_type processor_id() const
A Point defines a location in LIBMESH_DIM dimensional Real space.
Point vertex_average() const
bool is_evaluable(const DofObjectSubclass &obj, unsigned int var_num=libMesh::invalid_uint) const