unsaturated

transport_program.cpp
View Options

			#include "catch.hpp"

			#include "utils.hpp"
			#include "reactmicp/systems/saturated_diffusion/transport_program.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);

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

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


			}