testDeltaNotch1.cpp
b7305f7b5a8bmaster

testDeltaNotch1.cpp
View Options

			#include "deltaNotch1.hpp"

			using namespace std;
			using namespace boost::numeric::odeint;

			//integrate<void (&)(const state_type&),state_type,double,double,double> mintegrate;
			void dummy(const state_type &x , state_type &dxdt , const double t){
			for(int i=0;i<5;i++){
			dxdt[i] = - 0.01*x[i];
			}
			}


			void DeltaNotch1::operator() (const state_type &y , state_type &dydt , const double t)
			{
			double ddi,dni,dnidi,dnidj,dsi;
			double di=y[0];
			double ni=y[1];
			double nidi=y[2];
			double nidj=y[3];
			double si = y[4];
			//// delta
			ddi= bd -gd*di; // production and degr.
			ddi = ddi+ktmnjdi-ktpdi*nj; // trans binding
			ddi = ddi+ks*njdi; //release after cleaving
			ddi = ddi-kcpnidi+kcm*nidi; //// cis binding
			//// notch
			dni = bn -gn*ni;
			dni = dni + ktmnidj -ktpdj*ni; // trans binding
			dni = dni - kcpnidi+kcm*nidi; // cis binding
			dni = dni + diffu(np+nn-2ni); //diffusion
			//// cis complex
			dnidi = kcpnidi - kcmnidi - kinidi;
			//// trans complex
			dnidj = -ktmnidj + ktpdjni - ksnidj;
			//// cleaved notch
			dsi = ksnidj -gssi;
			////
			dydt[0] = ddi;
			dydt[1] = dni;
			dydt[2] = dnidi;
			dydt[3] = dnidj;
			dydt[4] = dsi;
			}

			void writeState(const state_type y,const double t){
			cout<<t<<" ";
			for(uint i=0;i<y.size();i++){
			cout<<y[i]<<" ";
			}
			cout<<endl;

			}

			int main(int argc, char* arg[]){
			DeltaNotch1 dn;
			state_type y = {{1.0,1.0,0.0,0.0,0.0}};
			dn.setContext(1.02,0.98,0,0.3,0,0);
			// size_t nsteps=integrate(dummy,y,0.0 , 10.0 , 0.1);
			integrate(dn,y,0.0 , 10.0 , 0.1,writeState);
			// for(uint i=0;i<5;i++){
			// cout<<y[i]<<" ";
			// }
			// cout<<endl;
			}