SteadyStateGsl.cpp File Reference

#include "../basecode/header.h"
#include "../basecode/global.h"
#include "SparseMatrix.h"
#include "KinSparseMatrix.h"
#include "RateTerm.h"
#include "FuncTerm.h"
#include "VoxelPoolsBase.h"
#include "../mesh/VoxelJunction.h"
#include "XferInfo.h"
#include "ZombiePoolInterface.h"
#include "Stoich.h"
#include "OdeSystem.h"
#include "VoxelPools.h"
#include "SteadyStateGsl.h"

Include dependency graph for SteadyStateGsl.cpp:

Go to the source code of this file.

Classes
struct	reac_info

Functions
static bool	checkAboveZero (const vector< double > &y)
static double	invop (double x)
static bool	isSolutionPositive (const vector< double > &x)
static double	op (double x)
int	ss_func (const gsl_vector *x, void *params, gsl_vector *f)

Variables
static const Cinfo *	steadyStateCinfo = SteadyState::initCinfo()

Function Documentation

static bool checkAboveZero ( const vector< double > &  y )

static

Definition at line 1086 of file SteadyStateGsl.cpp.

{
    for ( vector< double >::const_iterator
            i = y.begin(); i != y.end(); ++i )
    {
        if ( *i < 0.0 )
            return false;
    }
    return true;
}

static double invop ( double  x )

static

Definition at line 674 of file SteadyStateGsl.cpp.

{
    if ( x > 0.0 )
        return sqrt( x );
    return 0.0;
}

static bool isSolutionPositive ( const vector< double > &  x )

static

Definition at line 825 of file SteadyStateGsl.cpp.

{
    for ( vector< double >::const_iterator
            i = x.begin(); i != x.end(); ++i )
    {
        if ( *i < 0.0 )
        {
            cout << "Warning: SteadyState iteration gave negative concs\n";
            return false;
        }
    }
    return true;
}

static double op ( double  x )

static

Definition at line 669 of file SteadyStateGsl.cpp.

Referenced by SetGet0::checkOpClass(), SetGet1< A >::checkOpClass(), PostMaster::clearPendingSetGet(), localGet(), OpFunc::ops(), SetGet0::set(), SetGet1< A >::set(), SetGet2< L, A >::set(), SetGet3< A1, A2, A3 >::set(), SetGet4< A1, A2, A3, A4 >::set(), SetGet5< A1, A2, A3, A4, A5 >::set(), SetGet6< A1, A2, A3, A4, A5, A6 >::set(), SetGet1< A >::setVec(), SetGet2< L, A >::setVec(), and wildcardFieldComparison().

{
    return x * x;
}

Here is the caller graph for this function:

int ss_func	(	const gsl_vector *	x,
		void *	params,
		gsl_vector *	f
	)

This program works out a steady-state value for a reaction system. It uses GSL heavily, and isn't even compiled if the flag isn't set. It finds the ss value closest to the initial conditions.

If you want to find multiple stable states, it is best to do this in Python as it gives a lot of flexibility in working out how to find steady states. Likewise, if you want to carry out a dose-response calculation.

Variable Documentation

const Cinfo* steadyStateCinfo = SteadyState::initCinfo()

static

Definition at line 314 of file SteadyStateGsl.cpp.