Homec4science
Diffusion tisue modeling (master) Browse

tismod/include/deltaMechanicalForces.hpp
b7305f7b5a8bmaster

Last Change
Raw File

deltaMechanicalForces.hpp
View Options

#ifndef __DELTA_MECHANICAL_FORCES_H__
#define __DELTA_MECHANICAL_FORCES_H__
#include <iostream>
#include <boost/array.hpp>
#include "mechanicalForces.hpp"
#include "cellInteraction.hpp"
//#include "lattice.hpp"
//#include "observer.hpp"
class DeltaMechanicalForces;
//template <class Observer>
class DeltaMechanicalForces:public MechanicalForces, public CellInteraction{
protected:
Lattice *lattice;
// Observer *observer;
double gamma_d;
double h_dc;
double dprod;
double hfdp;
double del_sig;
double tau_sig;
double tau_mech;
protected:
double deltacisfunc(double delta) {double a=pow(h_dc,hill_co);double b=pow(delta,hill_co);return a/(a+b);}
virtual bool isHC(Cell *c){return m_signal[c->id()]<0.05;}
virtual bool isSC(Cell *c){return m_signal[c->id()]>=2;}
// double weightfunc(double s)
// {double a=pow(h_we,hill_co);double b=pow(s,hill_co);return we_min+we_r*a/(a+b);}
public:
DeltaMechanicalForces(){initParameters();}
void initParameters();
void operator()(const state_type &y, state_type &dydt, const double t);
virtual float getWeight(Boundary *b);
// void setState(state_type& s){MechanicalForces::setState(s);CellInteraction::setState(s);}
virtual void initState(state_type &s);
virtual void initState(state_type &s, double *delta);
virtual double getDeltaProduction(double si){return dprod/(1+pow(del_sig*si,hfdp));}
virtual uint cellStateDim(){return 1;}
virtual uint BoundaryStateDim(){return 0;}
virtual uint setStateVectorIndex();
virtual void setLattice(Lattice *lat);
void fillDeltaVector(double *v, state_type& s);
// void fillSignalVector(double *v); // in Mechanical forces
void updateDelta(const state_type &y, state_type &dydt);
virtual void updateSignal(const state_type &y);
// void findT1Transitions();
// virtual void checkT1Transition(Boundary *b);
// bool fromState(const state_type &y);
};
// void myWrapper(Boundary *b, void *args){ checkT1TransitionWrapper<DeltaMechanicalForces>(b, args);}
//checkT1TransitionWrapper<DeltaMechanicalForces>
#endif
c4science · Help