Calculate linear combinations of orthonormal cosine basis functions on two- or three-dimensional domains