unsaturated

transport_program.cpp
View Options

			#include "catch.hpp"

			#include "utils.hpp"
			#include "dfpm/meshes/uniform_mesh1d.hpp"
			#include "reactmicp/systems/saturated_diffusion/transport_program.hpp"
			#include "reactmicp/systems/saturated_diffusion/transport_parameters.hpp"

			using namespace specmicp;
			using namespace specmicp::reactmicp;
			using namespace specmicp::reactmicp::systems;
			using namespace specmicp::reactmicp::systems::siasaturated;

			TEST_CASE("Transport program", "[Transport, Diffusion, Finite Element, Solver]") {

			std::shared_ptr<database::DataContainer> database = get_test_carbo_database();
			int nb_nodes = 5;

			EquilibriumState initial_state = sample_carbo_composition(database);

			mesh::Mesh1DPtr themesh = mesh::uniform_mesh1d(nb_nodes, 1e-3, 0.001);
			auto parameters = std::make_shared<SaturatedDiffusionTransportParameters>(nb_nodes, 1e-6, 0.2);

			SECTION("Initialization") {
			SaturatedDiffusionProgram(themesh, database, parameters);
			}


			}