ksolve/SteadyState.cpp File Reference

#include "header.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 "../randnum/randnum.h"
#include "OdeSystem.h"
#include "VoxelPools.h"
#include "SteadyState.h"

Include dependency graph for SteadyState.cpp:

Classes
struct	reac_info
Functions
int	ss_func (const gsl_vector *x, void *params, gsl_vector *f)

---

Function Documentation

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.