Homec4science
Diffusion tisue modeling (master) Browse

tismod/src/debug/deb_coupling.cpp
b7305f7b5a8bmaster

Last Change
Raw File

deb_coupling.cpp
View Options

#include <math.h>
#include <iostream>
#define PIf 3.14159265358979323846264338327950288419716939937510582097494459230781640628621
using namespace std;
class CollierCoupling{
protected:
struct fparams {
double amp;
double hdeg;//a^deg
double deg;
public:
void init_hill(double amp, double half,double deg){
this->amp=amp;this->deg=deg;this->hdeg=pow(half,deg);}
void init_cauchy(double amp, double half,double deg){
this->amp=amp;this->deg=deg;this->hdeg=half;}
} fhp,ghp;
public:
virtual void init(){}
virtual double f_func(double x){}
virtual double g_func(double x){}
double f_hill(double x)
{double xh = pow(x,fhp.deg);return fhp.amp*xh/(xh+fhp.hdeg);}
double g_hill(double x)
{double xh = pow(x,ghp.deg);return ghp.amp*ghp.hdeg/(xh+ghp.hdeg);}
double f_cauchy(double x){return fhp.amp*(0.5+1/PI*atan((x-fhp.hdeg)/fhp.deg));};
double g_cauchy(double x){return ghp.amp*(0.5-1/PI*atan((x-ghp.hdeg)/ghp.deg));};
void set_hill_params(double amp, double half,double deg)
{fhp.amp=amp;fhp.deg=deg;fhp.hdeg=pow(half,deg);}
};
class HillCoupling : public CollierCoupling{
public:
virtual void init(){fhp.init_hill(1.0,0.2,3.0);ghp = fhp;}
virtual double f_func(double x){return f_hill(x);}
virtual double g_func(double x){return g_hill(x);}
};
class HillCauchyCoupling : public CollierCoupling{
public:
virtual void init(){
fhp.init_hill(pow(10,0.6139),pow(10,-3.8524),pow(10,-0.2898));
ghp.init_cauchy(pow(10,-1.0172),pow(10,-0.3977),0.0001);}
virtual double f_func(double x){return f_hill(x);}
virtual double g_func(double x){return g_cauchy(x);}
};
int main(int argc, char *argv[]){
HillCauchyCoupling sig;
sig.init();
for (double x=0.0;x<4.5;x+=0.01){
double y [4];
y[0] = sig.g_func(x);
y[1] =sig.f_func(y[0]);
y[2] =sig.g_func(y[1]);
y[3] =sig.f_func(y[2]);
cout<<x<<":"<<y[0]<<" "<<y[1]<<" "<<y[2]<<" "<<y[3]<<" "<< ((y[3]>x)?"+":"-") <<endl;
}
return 0;
}
c4science · Help