__init__(self, centers, potential_derivatives, modes, **opts): 

__call__(self, points, deriv_order=None, **kwargs):