A collocation code for computing gaussians on a grid of the form: out_Lp = x^l y^m z^n \sum_i coeff_i e^(exponent_i * (|center - p|)^2) , where the returned matrix dimension are the angular momentum (L) by number of requested points (p).