void operator()(const Device *ctx, const int &ng, const int &lmax, const FPTYPE &SQRT2, const FPTYPE &PI, const FPTYPE &PI_HALF, const FPTYPE &FOUR_PI, const FPTYPE &SQRT_INVERSE_FOUR_PI, const FPTYPE *g, FPTYPE *p, FPTYPE *ylm)
YLM_REAL::Real spherical harmonics ylm(G) up to l=lmax Use Numerical recursive algorithm as given in ...