Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add new tests from SubDiv branch into main. #22

Merged
merged 3 commits into from
Oct 25, 2024
Merged

Add new tests from SubDiv branch into main. #22

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
5ae02ce
Updated BoundaryMeshingTests to include tests to cover for genRect()
TheBEllis Oct 25, 2024
b94420f
Updated BoundaryMeshingTests to include tests to cover for genRect()
TheBEllis Oct 25, 2024
afb0594
Merge pull request #21 from aurora-multiphysics/main
TheBEllis Oct 25, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
132 changes: 116 additions & 16 deletions test/src/BoundaryMeshingTests.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,109 @@ TEST_CASE("BoundaryGenerator", "[cube]") {
mesh.read(std::string(USER_DIR) +
"/test/testingMeshes/unit_cube/unit_cube.e");

SECTION("genSquare") {
SECTION("checkBoundary") {

// Boundary Parameters
// Boundary should be too small resulting in exception being thrown
double length = 0.1;
int subdivisions = 10;

//
BoundaryGenerator bound_generator(mesh, surface_mesh, boundary_mesh);
Eigen::MatrixXd boundary_vertices;
Eigen::MatrixXi boundary_elems;

// Generate square made of edge elements
REQUIRE_THROWS_AS(bound_generator.checkBoundary(length),
std::invalid_argument);
}

SECTION("genRect") {

// Boundary Parameters
// Boundary should be too small resulting in exception being thrown
double x_dim = 1;
double y_dim = 1.5;
int x_subdiv = 2;
int y_subdiv = 3;


//
BoundaryGenerator bound_generator(mesh, surface_mesh, boundary_mesh);
Eigen::MatrixXd boundary_verts;
Eigen::MatrixXi boundary_elems;

Eigen::MatrixXd correct_boundary_verts(10, 2);
correct_boundary_verts(0, 0) = -0.5;
correct_boundary_verts(0, 1) = -0.75;

correct_boundary_verts(1, 0) = 0.0;
correct_boundary_verts(1, 1) = -0.75;

correct_boundary_verts(2, 0) = 0.5;
correct_boundary_verts(2, 1) = -0.75;

correct_boundary_verts(3, 0) = 0.5;
correct_boundary_verts(3, 1) = -0.25;

correct_boundary_verts(4, 0) = 0.5;
correct_boundary_verts(4, 1) = 0.25;

correct_boundary_verts(5, 0) = 0.5;
correct_boundary_verts(5, 1) = 0.75;

correct_boundary_verts(6, 0) = 0.0;
correct_boundary_verts(6, 1) = 0.75;

correct_boundary_verts(7, 0) = -0.5;
correct_boundary_verts(7, 1) = 0.75;

correct_boundary_verts(8, 0) = -0.5;
correct_boundary_verts(8, 1) = 0.25;

correct_boundary_verts(9, 0) = -0.5;
correct_boundary_verts(9, 1) = -0.25;


Eigen::MatrixXi correct_boundary_elems(10, 2);
correct_boundary_elems(0, 0) = 0;
correct_boundary_elems(0, 1) = 1;

correct_boundary_elems(1, 0) = 1;
correct_boundary_elems(1, 1) = 2;

correct_boundary_elems(2, 0) = 2;
correct_boundary_elems(2, 1) = 3;

correct_boundary_elems(3, 0) = 3;
correct_boundary_elems(3, 1) = 4;

correct_boundary_elems(4, 0) = 4;
correct_boundary_elems(4, 1) = 5;

correct_boundary_elems(5, 0) = 5;
correct_boundary_elems(5, 1) = 6;

correct_boundary_elems(6, 0) = 6;
correct_boundary_elems(6, 1) = 7;

correct_boundary_elems(7, 0) = 7;
correct_boundary_elems(7, 1) = 8;

correct_boundary_elems(8, 0) = 8;
correct_boundary_elems(8, 1) = 9;

correct_boundary_elems(9, 0) = 0;
correct_boundary_elems(9, 1) = 9;

// Generate boundary elements of rectangle
bound_generator.genRect(boundary_verts, boundary_elems, x_dim, y_dim, x_subdiv, y_subdiv);

REQUIRE(boundary_verts == correct_boundary_verts);
REQUIRE(boundary_elems == correct_boundary_elems);
}

SECTION("genSquare") {

// Boundary Parameters
double length = 10;
Expand All @@ -34,22 +136,20 @@ TEST_CASE("BoundaryGenerator", "[cube]") {
REQUIRE(square_vertices.rows() == (subdivisions * 4));
}

SECTION("getSeedingPoints") {
// SECTION("genTriangulatedRect") {

// Boundary Parameters
// Boundary should be too small resulting in exception being thrown
double length = 0.1;
int subdivisions = 10;

//
BoundaryGenerator bound_generator(mesh, surface_mesh, boundary_mesh);
Eigen::MatrixXd boundary_vertices;
Eigen::MatrixXi boundary_elems;
// // Boundary Parameters
// // Boundary should be too small resulting in exception being thrown
// double length = 0.1;
// int subdivisions = 10;

// Generate square made of edge elements
REQUIRE_THROWS_AS(bound_generator.checkBoundary(length),
std::invalid_argument);
}
// //
// BoundaryGenerator bound_generator(mesh, surface_mesh, boundary_mesh);
// Eigen::MatrixXd boundary_vertices;
// Eigen::MatrixXi boundary_elems;

SECTION("throw") {}
// // Generate square made of edge elements
// REQUIRE_THROWS_AS(bound_generator.checkBoundary(length),
// std::invalid_argument);
// }
}