11 #ifndef _STEADYSTATE_H
12 #define _STEADYSTATE_H
14 #include <boost/numeric/ublas/matrix.hpp>
16 #include "../randnum/RNG.h"
23 friend void testSteadyState();
44 double getTotal(
const unsigned int i )
const;
45 void setTotal(
const unsigned int i,
double val );
59 void settle(
bool forceSetup );
63 static void assignY(
double* S );
77 boost::numeric::ublas::matrix< value_type_ >& U
78 ,
const vector< value_type_ >& eliminatedTotal
79 , vector <value_type_> &yi
85 static void setMolN(
double y,
unsigned int i );
109 boost::numeric::ublas::matrix< value_type_ >
LU_;
110 boost::numeric::ublas::matrix< value_type_ >
Nr_;
111 boost::numeric::ublas::matrix< value_type_ >
gamma_;
136 #endif // _STEADYSTATE_H
void randomizeInitialCondition(const Eref &e)
unsigned int getMaxIter() const
unsigned int getNposEigenvalues() const
boost::numeric::ublas::matrix< value_type_ > gamma_
unsigned int getNiter() const
unsigned int getNumVarPools() const
double getConvergenceCriterion() const
unsigned int getRank() const
vector< double > eigenvalues_
double convergenceCriterion_
bool badStoichiometry() const
bool isInitialized() const
static void setMolN(double y, unsigned int i)
double getEigenvalue(const unsigned int i) const
unsigned int nPosEigenvalues_
void setEigenvalue(double val, const unsigned int i)
static const double DELTA
double getTotal(const unsigned int i) const
void setMaxIter(unsigned int value)
unsigned int nNegEigenvalues_
unsigned int getSolutionStatus() const
static const double EPSILON
unsigned int getStateType() const
void fitConservationRules(boost::numeric::ublas::matrix< value_type_ > &U, const vector< value_type_ > &eliminatedTotal, vector< value_type_ > &yi)
unsigned int getNnegEigenvalues() const
void setConvergenceCriterion(double value)
static const Cinfo * initCinfo()
boost::numeric::ublas::matrix< value_type_ > LU_
void classifyState(const double *T)
boost::numeric::ublas::matrix< value_type_ > Nr_
static void assignY(double *S)
unsigned int solutionStatus_
const Cinfo * initSteadyStateCinfo()
void settle(bool forceSetup)
void setTotal(const unsigned int i, double val)