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()
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)