import numpy as np
from scipy.integrate import odeint
from DS_solver import bioDS
from cell_cycle import cell
test=cell(path='SinIR_delay/')
test.load()
test.set("Spo0A_p",0)
test.simulate(tmax=30)
test.update()
test.clear()
for i in range(200):
test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=50)
print(1)
test.parameterTurn(sigma=0.25)
test.display(['Ptap','Pmot'],save=True,fname='test')
#v_line=[2,6,8,12,14,18,20,24,26,30,32,36,38]
import numpy as np
from scipy.integrate import odeint
import pickle
from DS_solver import bioDS
from cell_cycle import cell
import matplotlib.pyplot as plt
for ppath in ['SinIR_delayf2_5/','SinIR_delayf2_4/','SinIR_delayf2_d/']:
print(ppath)
test=cell(path=ppath)
test.load()
test.set("Spo0A_p",0)
test.simulate(tmax=30)
test.update()
test.clear()
for i in range(20000):
test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=30)
test.parameterTurn(sigma=0.25)
if (i/1000)*1000==i:
print i
tt=test.trajectory[0]
yy=test.getTrajectory('Ptap')
zz=test.getTrajectory('Pmot')
pickle.dump([tt,yy,zz],open(ppath+"f.pickle",'wb'))
for threshold in [60,80,100,120,140,160,180,200,250,300,400,500]:
chain=[]
mot=[]
s=0
t=0
if yy[0]>threshold:
s=1
for i in range(1,len(tt)):
import numpy as np
from scipy.integrate import odeint
import pickle
from DS_solver import bioDS
from cell_cycle import cell
import matplotlib.pyplot as plt
test=cell(path='SinIR_feedback/')
test.load()
test.set("Spo0A_p",0)
test.simulate(tmax=30)
test.update()
test.clear()
for i in range(20000):
test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=30)
test.parameterTurn(sigma=0.25)
if (i/100)*100==i:
print i
tt=test.trajectory[0]
yy=test.getTrajectory('Ptap')
pickle.dump([tt,yy],open("SinIR_feedback/ttyy-fb_2.pickle",'wb'))
for threshold in [60,80,100,120,140,160,180,200,250,300,400,500]:
chain=[]
mot=[]
s=0
t=0
if yy[0]>threshold:
s=1
for i in range(1,len(tt)):
if s==0:
if yy[i]>threshold:
s=1
import numpy as np
from scipy.integrate import odeint
from DS_solver import bioDS
from cell_cycle import cell
test=cell(path='SinIR_delayf2_d/')
test.load()
test.set("Spo0A_p",0)
test.simulate(tmax=30)
test.update()
test.clear()
for i in range(200):
test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=50)
print(1)
test.parameterTurn(sigma=0.25)
test.display(['Ptap','Pmot'],save=True,fname='test')
#v_line=[2,6,8,12,14,18,20,24,26,30,32,36,38]