estimate_dt_n_steps¶
-
gala.dynamics.
estimate_dt_n_steps
(w0, potential, n_periods, n_steps_per_period, dE_threshold=1e-09, func=<function nanmax>)[source]¶ Estimate the timestep and number of steps to integrate an orbit for given its initial conditions and a potential object.
Parameters: w0 :
PhaseSpacePosition
, array_likeInitial conditions.
potential :
PotentialBase
The potential to integrate the orbit in.
n_periods : int
Number of (max) orbital periods to integrate for.
n_steps_per_period : int
Number of steps to take per (max) orbital period.
dE_threshold : numeric (optional)
Maximum fractional energy difference – used to determine initial timestep. Set to
None
to ignore this.func : callable (optional)
Determines which period to use. By default, this takes the maximum period using
nanmax()
. Other options could benanmin()
,nanmean()
,nanmedian()
.Returns: dt : float
The timestep.
n_steps : int
The number of timesteps to integrate for.