galpy.df.quasiisothermaldf.__call__

quasiisothermaldf.__call__(*args, **kwargs)[source]

NAME:

__call__

PURPOSE:

return the DF

INPUT:

Either:

a)(jr,lz,jz) tuple; each can be a Quantity
where:

jr - radial action lz - z-component of angular momentum jz - vertical action

  1. R,vR,vT,z,vz

  2. Orbit instance: initial condition used if that’s it, orbit(t) if there is a time given as well

log= if True, return the natural log

+scipy.integrate.quadrature kwargs

func= function of (jr,lz,jz) to multiply f with (useful for moments)

OUTPUT:

value of DF

HISTORY:

2012-07-25 - Written - Bovy (IAS@MPIA)

NOTE:

For Miyamoto-Nagai/adiabatic approximation this seems to take about 30 ms / evaluation in the extended Solar neighborhood For a MWPotential/adiabatic approximation this takes about 50 ms / evaluation in the extended Solar neighborhood

For adiabatic-approximation grid this seems to take about 0.67 to 0.75 ms / evaluation in the extended Solar neighborhood (includes some out of the grid)

up to 200x faster when called with vector R,vR,vT,z,vz