where N l is the number of positive roots in the Gauss-Legendre quadrature, M t is the number of points in the upper hemisphere of the spherical t-design grid, and ε is the error. For computational efficiency, we express the complex exponentials as the real and imaginary trigonometric terms and for a given atom, we define a basis function at a specific cubature node (α, n, m) as: