Homec4science
Diffusion SpecMiCP / ReactMiCP (specmicp_bin) Browse

specmicp/src/specmicp/reduced_system_solver.hpp
99cf67ee7b3cspecmicp_bin

Last Change
Raw File

reduced_system_solver.hpp
View Options

/*-------------------------------------------------------
- Module : specmicp
- File : reduced_system_solver.hpp
- Author : Fabien Georget
Copyright (c) 2014, Fabien Georget, Princeton University
---------------------------------------------------------*/
#ifndef SPECMICP_SPECMICP_REDUCEDSYSTEMSOLVER_HPP
#define SPECMICP_SPECMICP_REDUCEDSYSTEMSOLVER_HPP
#include "database/data_container.hpp"
#include "reduced_system.hpp"
#include "micpsolver/micpsolver_structs.hpp"
//! \file reduced_system_solver.hpp Solve a reduced system
namespace specmicp {
struct ReducedSystemSolverOptions
{
micpsolver::NCPfunction ncp_function;
micpsolver::MiCPSolverOptions solver_options;
ReducedSystemSolverOptions():
ncp_function(micpsolver::NCPfunction::penalizedFB),
solver_options()
{}
};
//! \brief Solver a reduced system
//!
//! Take care of possibly non-existing component in the system
class ReducedSystemSolver
{
public:
//! Default constructor - you need to use another method to initialize it correctly
ReducedSystemSolver() {}
ReducedSystemSolver(std::shared_ptr<database::DataContainer> data,
Eigen::VectorXd& totconc);
//! \brief solve the problem using initial guess x
//!
//! @param[in,out] x, in -> initial guess, out -> solution
micpsolver::MiCPPerformance solve(Eigen::VectorXd& x);
//! Return the system used for the computation
std::shared_ptr<ReducedSystem> get_system() {return m_system;}
//! Set total aqueous concentrations
void total_aqueous_concentration(const Eigen::VectorXd& x, Eigen::VectorXd& totaq)
{
m_system->aqueous_tot_conc(x, totaq);
}
const ReducedSystemSolverOptions& get_options() const {return m_options;}
ReducedSystemSolverOptions& get_options() {return m_options;}
private:
//! \brief set up the true variable vector
void set_true_variable_vector(const Eigen::VectorXd& x);
//! \brief set up the true solution vector
//!
//! add zero components
void set_return_vector(Eigen::VectorXd& x);
//! \brief solve the problem
micpsolver::MiCPPerformance solve();
ReducedSystemSolverOptions m_options;
std::shared_ptr<database::DataContainer> m_data;
std::shared_ptr<ReducedSystem> m_system;
Eigen::VectorXd m_var;
};
inline micpsolver::MiCPPerformance solve_reduced_system(std::shared_ptr<database::DataContainer> data,
Eigen::VectorXd& totconc,
Eigen::VectorXd& x)
{
return ReducedSystemSolver(data, totconc).solve(x);
}
} // end namespace specmicp
#endif // SPECMICP_SPECMICP_REDUCEDSYSTEMSOLVER_HPP
c4science · Help