galpy.df.quasiisothermaldf.vmomentdensity

quasiisothermaldf.vmomentdensity(R, z, n, m, o, nsigma=None, mc=False, nmc=10000, _returnmc=False, _vrs=None, _vts=None, _vzs=None, _rawgausssamples=False, gl=False, ngl=10, _returngl=False, _glqeval=None, _return_actions=False, _jr=None, _lz=None, _jz=None, _return_freqs=False, _rg=None, _kappa=None, _nu=None, _Omega=None, _sigmaR1=None, _sigmaz1=None, **kwargs)

NAME:

vmomentdensity

PURPOSE:

calculate the an arbitrary moment of the velocity distribution at R times the density

INPUT:

R - radius at which to calculate the moment(/ro)

n - vR^n

m - vT^m

o - vz^o

OPTIONAL INPUT:

nsigma - number of sigma to integrate the velocities over (when doing explicit numerical integral)

mc= if True, calculate using Monte Carlo integration

nmc= if mc, use nmc samples

gl= use Gauss-Legendre

_returngl= if True, return the evaluated DF

_return_actions= if True, return the evaluated actions (does not work with _returngl currently)

_return_freqs= if True, return the evaluated frequencies and rg (does not work with _returngl currently)

OUTPUT:

<vR^n vT^m x density> at R,z

HISTORY:

2012-08-06 - Written - Bovy (IAS@MPIA)