added scipy solver back into risk_simulator

tapios / risk-networks

Code for risk networks: a blend of compartmental models, graphs, data assimilation and semi-supervised learning

Other

2 stars 2 forks source link

added scipy solver back into risk_simulator #162

Closed odunbar closed 4 years ago

odunbar commented 4 years ago

Replace forward and backward solver in risk_simulator with RK45 scipy solver for larger timestepping Coauthored with @jinlong83 @dburov190

dburov190 commented 4 years ago

I slightly tweaked the formatting and namings (like, removed self. from a couple variables because they are only used locally; renamed do_step -> compute_rhs; etc.).

There are two potential issues:

the one I mentioned above, whether we need to evaluate the closure every little time-step (which is easy to implement, just stick self.eval_closure inside compute_rhs, but I'm willing to give this crude way a try)
no clipping

Anyway, I'm gonna merge this and see if it works (: