constants#

Module containing functions for computing the constants of motion for an orbit in Kerr Spacetime. Constants of motion are computed using the method described in Appendix B of Schmidt.

Functions

angular_momentum(a, p, e, x[, E])

Computes the dimensionless angular momentum of a bound orbit with the given parameters

apex_from_constants(a, E, L, Q)

Computes the orbital parameters \((a,p,e,x)\) for a stable bound orbit with the given constants of motion

carter_constant(a, p, e, x[, E, L])

Computes the dimensionless carter constant of a bound orbit with the given parameters

constants_of_motion(a, p, e, x)

Computes the dimensionless energy, angular momentum, and Carter constant of a bound orbit with the given parameters

energy(a, p, e, x)

Computes the dimensionless energy of a bound orbit with the given parameters

fast_separatrix(a[, grid_spacing])

Constructs a faster separatrix function for a given value of \(a\) by interpolating over a grid of \(e\) and \(x\) values.

is_stable(a, p, e, x)

Tests whether or not the given orbital parameters define a stable bound orbit

plunging_radial_roots(a, E, L, Q)

Computes the radial roots for a plunging orbit.

scale_constants(constants, M, mu)

Scales the dimensionless constants of motion to the given mass parameters

separatrix(a, e, x)

Returns the value of p at the separatrix for the given orbital parameters computed using the bracked root finding method described in Stein and Warburton

stable_polar_roots(a, p, e, x[, constants])

Computes the polar roots for a stable bound orbit as defined in equation 10 of Fujita and Hikida.

stable_radial_roots(a, p, e, x[, constants])

Computes the radial roots for a stable bound orbit as defined in equation 10 of Fujita and Hikida.

valid_params(a, e, x)

Tests whether the given parameters fall into the allowed ranges