test_uniform_mesh_1d.cpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#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);
	}
	}