test_uniform_mesh_1d.cpp
114f05dd6425specmicp_bin

test_uniform_mesh_1d.cpp
View Options

			#include "catch.hpp"

			#include "dfpm/meshes/uniform_mesh1d.hpp"

			using namespace specmicp::mesh;

			TEST_CASE("Uniform mesh 1d", "[ReactMiCP],[Mesh],[uniform],[1D]") {

			SECTION("test uniform mesh 1d") {

			double nb_node = 5;
			double dx = 1;
			double cross = 2;

			UniformMesh1D mesh(nb_node, dx, cross);

			REQUIRE(mesh.nb_nodes() == nb_node);
			REQUIRE(mesh.nb_elements() == nb_node - 1);

			REQUIRE(mesh.get_dx(0) == dx);
			REQUIRE(mesh.get_face_area(2) == cross);
			REQUIRE(mesh.get_volume_element(2) == dx*cross);

			int check_nb_node = 0;
			for (int node: mesh.range_nodes()) {++check_nb_node;}
			REQUIRE(check_nb_node == nb_node);

			int check_nb_element = 0;
			for (int element: mesh.range_elements()) {++check_nb_element;}
			REQUIRE(check_nb_element == mesh.nb_elements());

			int check_nen = 0;
			for (int node: mesh.range_nen()) {++check_nen;}
			REQUIRE(check_nen == 2);
			}
			}