variables.cpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include "variables.hpp"

	#include "dfpm/meshes/mesh1d.hpp"

	namespace specmicp {
	namespace reactmicp {
	namespace systems {
	namespace satdiff {

	SaturatedVariables::SaturatedVariables(mesh::Mesh1DPtr the_mesh,
	RawDatabasePtr the_database):
	m_mesh(the_mesh),
	m_database(the_database)
	{}

	Vector SaturatedVariables::total_concentrations(index_t node) {
	index_t offset = offset_node(node);

	return porosity(node)*displacement().segment(offset+offset_aqueous_concentration(), nb_component())
	+ displacement().segment(offset+offset_solid_concentration(), nb_component());
	}


	} // end namespace satdiff
	} // end namespace systems
	} // end namespace reactmicp
	} // end namespace specmicp