コード例 #1
0
ファイル: SIR.py プロジェクト: oaarscorp/CME
OFSP_obj = OFSP_Solver(SIR_model, "SE1", 50, 1e-6)
OFSP_obj.step(0.01)

### Hybrid MRCE Computation
MRCE_obj = Hybrid_FSP_solver(SIR_model, stoc_vector, "MRCE", 1e-7)
MRCE_obj.set_initial_values(OFSP_obj.domain_states, OFSP_obj.p, t=0.01)

for t in T:
    MRCE_obj.step(t)
    MRCE_obj.print_stats
    #MRCE_obj.plot(inter=True)
    MRCE_obj.check_point()
    X = np.array([[150.0, 180.0], [1, 2]])
    MRCE_obj.probe_states(X)

MRCE_obj.plot_checked()

### Hybrid HL Computation
from pyme.Hybrid_FSP import Hybrid_FSP_solver
stoc_vector = np.array([True, False])

HL_obj = Hybrid_FSP_solver(SIR_model, stoc_vector, "HL", 1e-7)
HL_obj.set_initial_values(OFSP_obj.domain_states, OFSP_obj.p, t=0.01)
for t in T:
    HL_obj.step(t)
    HL_obj.print_stats
    HL_obj.check_point()
    X = np.array([[150.0, 180.0], [1, 2]])
    HL_obj.probe_states(X)

HL_obj.plot_checked()
コード例 #2
0
ファイル: SIR.py プロジェクト: SysSynBio/PyME
OFSP_obj = OFSP_Solver(SIR_model,"SE1",50,1e-6)
OFSP_obj.step(0.01)

### Hybrid MRCE Computation
MRCE_obj = Hybrid_FSP_solver(SIR_model,stoc_vector,"MRCE",1e-7)
MRCE_obj.set_initial_values(OFSP_obj.domain_states,OFSP_obj.p,t=0.01)

for t in T:
	MRCE_obj.step(t)
	MRCE_obj.print_stats
	#MRCE_obj.plot(inter=True)
	MRCE_obj.check_point()
	X = np.array([[150.0,180.0],[1,2]])
	MRCE_obj.probe_states(X)

MRCE_obj.plot_checked()


### Hybrid HL Computation
from pyme.Hybrid_FSP import Hybrid_FSP_solver
stoc_vector = np.array([True,False])

HL_obj = Hybrid_FSP_solver(SIR_model,stoc_vector,"HL",1e-7)
HL_obj.set_initial_values(OFSP_obj.domain_states,OFSP_obj.p,t=0.01)
for t in T:
	HL_obj.step(t)
	HL_obj.print_stats
	HL_obj.check_point()
	X = np.array([[150.0,180.0],[1,2]])
	HL_obj.probe_states(X)