unsaturated

boundary_conditions.hpp
View Options

			#ifndef SPECMICP_REACTMICP_SATURATED_DIFFUSION_EQUATIONNUMBERING_HPP
			#define SPECMICP_REACTMICP_SATURATED_DIFFUSION_EQUATIONNUMBERING_HPP

			#include <vector>
			#include "specmicp/equilibrium_data.hpp"

			namespace specmicp {
			namespace reactmicp {
			namespace systems {
			namespace siasaturated {

			//! \brief The different types of Boundary conditions accepted
			enum class BoundaryConditionType
			{
			Free, //!< No BC
			FixedComposition //!< Fixed composition for this node
			};

			//! \brief Store the boundary conditions to be parsed by the reactive transport solver
			struct SIABoundaryConditions
			{
			//! \brief list of the different initial state used
			//!
			//! The first one should be the one used everywhere
			std::vector<EquilibriumState> list_initial_states;
			//! \brief The BC type for each node
			std::vector<BoundaryConditionType> bs_types;
			//! \brief The initial state to use for each node
			std::vector<int> initial_states;

			SIABoundaryConditions(int nb_nodes):
			bs_types(nb_nodes, BoundaryConditionType::Free),
			initial_states(nb_nodes, 0)
			{}
			};


			} // end namespace siasaturated
			} // end namespace systems
			} // end namespace reactmicp
			} // end namespace specmicp


			#endif // SPECMICP_REACTMICP_SATURATED_DIFFUSION_EQUATIONNUMBERING_HPP