Use scipy's Radau method to solve the odes for Zetas and Zetaw. This is an implicit method so we should be able to take larger step between each mechanics solves.
With the current implementation I am able to run the simple_demo and arrive at the same reference solution using mech_threshold=2.0 (default is 0.05) or dt_mech=0.1 (default is 1.0). This mean essentially that we solve the mechanics model every time step which is not any improvement.
Use scipy's Radau method to solve the odes for Zetas and Zetaw. This is an implicit method so we should be able to take larger step between each mechanics solves.
With the current implementation I am able to run the
simple_demoand arrive at the same reference solution usingmech_threshold=2.0(default is 0.05) ordt_mech=0.1(default is 1.0). This mean essentially that we solve the mechanics model every time step which is not any improvement.