Closed RudolfWeeber closed 6 months ago
Furhter discussions: np.any() rather than np.all() also, termination can just be controlled in system.integator.set_steepest_descent(..., f_max=1/syste.time_step) system.integrator.run(big_integer) if np.any(np.abs(system.part.all().f)>1/system_time_Step): raise Exception("Overlap removal did not converge")
Different criteria are used in different tutorials. I think
while np.all(np.abs(system.part.all().f)>1/system.time_Step)
is preferrable