Homec4science
Diffusion SpecMiCP / ReactMiCP (kinetics) Browse

specmicp/src/reactmicp/systems/saturated_diffusion/transport_neutrality_variables.hpp
ab2944bc8002kinetics

Last Change
Raw File

transport_neutrality_variables.hpp
View Options

#ifndef SPECMICP_REACTMICP_SATURATED_DIFFUSIONNEUTRALITY_VARIABLES_HPP
#define SPECMICP_REACTMICP_SATURATED_DIFFUSIONNEUTRALITY_VARIABLES_HPP
#include "reactmicp/systems/secondary_variables/secondary.hpp"
#include "specmicp/equilibrium_data.hpp"
namespace specmicp {
namespace reactmicp {
namespace systems {
namespace siasaturated {
//! Variables for the SIASaturatedReactiveTransportSolver
class TransportNeutralityVariables: public SecondaryVariables<TransportNeutralityVariables>
{
public:
TransportNeutralityVariables(
index_t nb_nodes,
RawDatabasePtr thedatabase
):
SecondaryVariables(nb_nodes, thedatabase),
m_data(thedatabase),
m_logb_component(nb_nodes*(m_data->nb_component+offset_component()))
{}
index_t get_ndf() const {return m_data->nb_component+offset_component();}
index_t get_dof_component(index_t node, index_t component) const {
specmicp_assert(component >= 1 and component < m_data->nb_component);
return (component+offset_component())+get_ndf()*node;
}
index_t get_edof_component(index_t enode, index_t component) const {
specmicp_assert(component >= 1 and component < m_data->nb_component);
return (component+offset_component())+get_ndf()*enode;
}
//! Return the concentration of 'component' at 'node'
scalar_t component_concentration(index_t node, index_t component) const {
return pow10(m_logb_component.coeff(get_dof_component(node, component)));
}
//! Return the concentration of 'component' at 'node'
scalar_t log_component_concentration(index_t node, index_t component) const {
specmicp_assert(component >= 1 and component < m_data->nb_component);
return m_logb_component(get_dof_component(node, component));
}
//! Return a reference to the concentration of 'component' at 'node'
scalar_t& log_component_concentration(index_t node, index_t component) {
specmicp_assert(component >= 1 and component < m_data->nb_component);
concentrations_have_changed(node);
return m_logb_component(get_dof_component(node, component));
}
//! Reference to the aqueous component concentrations
Vector& component_concentrations() {
return m_logb_component;
}
//! Eigen Expression for aqueous component concentrations of a node
Eigen::Block<Vector, Eigen::Dynamic, 1> component_concentrations(index_t node) {
return m_logb_component.segment(node*get_ndf(),get_ndf());
}
private:
// where the component informations are stored
index_t offset_component() const {return -1;}
RawDatabasePtr m_data;
Vector m_logb_component;
};
} // end namespace siasaturated
} // end namespace systems
} // end namespace reactmicp
} // end namespace specmicp
#endif // SPECMICP_REACTMICP_SATURATED_DIFFUSIONNEUTRALITY_VARIABLES_HPP
c4science · Help