galpy.actionAngle.actionAngleTorus.xvJacobianFreqs

actionAngleTorus.xvJacobianFreqs(jr, jphi, jz, angler, anglephi, anglez, **kwargs)[source]

NAME:

xvJacobianFreqs

PURPOSE:

return [R,vR,vT,z,vz,phi], the Jacobian d [R,vR,vT,z,vz,phi] / d (J,angle), the Hessian dO/dJ, and frequencies Omega corresponding to a torus at multiple sets of angles

INPUT:

jr - radial action (scalar)

jphi - azimuthal action (scalar)

jz - vertical action (scalar)

angler - radial angle (array [N])

anglephi - azimuthal angle (array [N])

anglez - vertical angle (array [N])

tol= (object-wide value) goal for |dJ|/|J| along the torus

dJ= (object-wide value) action difference when computing derivatives (Hessian or Jacobian)

nosym= (False) if True, don’t explicitly symmetrize the Hessian (good to check errors)

OUTPUT:

([R,vR,vT,z,vz,phi], [N,6] array

d[R,vR,vT,z,vz,phi]/d[J,angle], –> (N,6,6) array

dO/dJ, –> (3,3) array

Omegar,Omegaphi,Omegaz, [N] arrays

Autofit error message)

HISTORY:

2016-07-19 - Written - Bovy (UofT)