ActiveElementFractionTest.C File Reference

Go to the source code of this file.

Functions
	TEST (ActiveElementFractionTest, UnitSquare)

Function Documentation

TEST()

TEST	(	ActiveElementFractionTest	,
		UnitSquare
	)

Definition at line 22 of file ActiveElementFractionTest.C.

{
  Parallel::Communicator comm(MPI_COMM_SELF);
  auto mesh = std::make_unique<libMesh::SerialMesh>(comm);
 
  Node * n0 = mesh->add_point(Point(0.0, 0.0, 0.0), 0);
  Node * n1 = mesh->add_point(Point(1.0, 0.0, 0.0), 1);
  Node * n2 = mesh->add_point(Point(1.0, 1.0, 0.0), 2);
  Node * n3 = mesh->add_point(Point(0.0, 1.0, 0.0), 3);
 
  auto quad = new Quad4();
  quad->set_node(0, n0);
  quad->set_node(1, n1);
  quad->set_node(2, n2);
  quad->set_node(3, n3);
  const Elem * elem = mesh->add_elem(quad);
 
  mesh->prepare_for_use();
 
  const auto all_active = [](const Point &) { return true; };
  const auto none_active = [](const Point &) { return false; };
  const auto left_half = [](const Point & p) { return p(0) < 0.5; };
 
  // SECOND-order Gauss has 2 points per dimension, symmetric about the centre with none on it.
  EXPECT_NEAR(SBMUtils::activeElementFraction(*elem, SECOND, all_active), 1.0, 1e-12);
  EXPECT_NEAR(SBMUtils::activeElementFraction(*elem, SECOND, none_active), 0.0, 1e-12);
  EXPECT_NEAR(SBMUtils::activeElementFraction(*elem, SECOND, left_half), 0.5, 1e-12);
}