is_stable#

kerrgeopy.initial_conditions.is_stable(a, initial_position, initial_velocity, constants=None, tol=0.0001)[source]#

Tests whether the orbit with the given initial conditions is stable

Parameters:

  • a (double) – spin parameter

  • initial_position (tuple(double,double,double,double)) – initial position \((t_0,r_0,\theta_0,\phi_0)\)

  • initial_velocity (tuple(double,double,double,double)) – initial four-velocity \((u^t_0,u^r_0,u^\theta_0,u^\phi_0)\)

  • constants (tuple(double,double,double), optional) – tuple of constants of motion \((\mathcal{E},\mathcal{L},\mathcal{Q})\) can be passed in to avoid recomputing them

  • tol (double, optional) – numerical tolerance used when comparing \(r_0\) to radial roots, defaults to 1e-4

Returns:

True if the orbit is stable, False otherwise

Return type:

bool