boundary_conditions.hpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#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