compute_coeffs_discrete¶
- gala.potential.scf.compute_coeffs_discrete(xyz, mass, nmax, lmax, r_s, skip_odd=False, skip_even=False, skip_m=False, compute_var=False)[source]¶
Compute the expansion coefficients for representing the density distribution of input points as a basis function expansion. The points,
xyz, are assumed to be samples from the density distribution.Computing the coefficients involves computing triple integrals which are computationally expensive. For an example of how to parallelize the computation of the coefficients, see
examples/parallel_compute_Anlm.py.- Parameters
- xyzarray_like
Samples from the density distribution. Should have shape
(n_samples, 3).- massarray_like
Mass of each sample. Should have shape
(n_samples,).- nmax
int Maximum value of
nfor the radial expansion.- lmax
int Maximum value of
lfor the spherical harmonics.- r_s
numeric Scale radius.
- skip_oddbool (optional)
Skip the odd terms in the angular portion of the expansion. For example, only take \(l=0, 2, 4, ...\)
- skip_evenbool (optional)
Skip the even terms in the angular portion of the expansion. For example, only take \(l=1, 3, 5, ...\)
- skip_mbool (optional)
Ignore terms with \(m > 0\).
- compute_varbool (optional)
Also compute the variances (and covariances) of the coefficients.
- Returns