sigmalos2#

jeanspy.model.DSphModel.sigmalos2

Class and construction: DSphModel.

DSphModel.sigmalos2(R_pc, n=1024, n_kernel=128, ignore_RuntimeWarning=True, *, method='dequad')[source]#

Return the LOS velocity variance using the selected integration method.

Parameters:

  • R_pc (float or array_like) – Positive finite projected radius in pc; scalar or nonempty 1-D array. The model-dependent central limit at R=0 is not supported.

  • n (int, optional) – Outer quadrature order, default 1024.

  • n_kernel (int, optional) – Anisotropy-kernel quadrature order, default 128.

  • ignore_RuntimeWarning (bool, optional) – Suppress NumPy runtime warnings during integration, default true. Invalid integrands and nonfinite or negative variances still raise.

  • method ({"dequad"}, keyword-only, optional) – Integration method. The fixed double-exponential rule is currently the only supported choice and remains the default.

Returns:

scalar or ndarray – LOS velocity variance in (km/s)^2; scalar for scalar input, otherwise shape (N,).

Raises:

ValueError – Unsupported method, invalid projected radii, or nonphysical density, mass, integrand or variance.

Notes

This NumPy/SciPy calculation reads the stored model parameters. Refine n and n_kernel to check numerical convergence. The NumPy/SciPy inference model uses this entry point with its defaults; no extra callable or closure is needed when passing it to emcee.