def gauss(d,**parms):
# gaussian, f(r) = exp(-(ep r)^2)
c = parms.get('centers')
#op = parms.get('operator','interp')
ep = parms.get('shapeparm',1)
DM = dmatrix(d, centers = c)
# eps_r = epsilon*r where epsilon may be an array or scalar
eps_r = dot(ep*eye(DM.shape[0]),DM)
return exp(-(eps_r)**2)
def mq(d,**parms):
# multiquadric, f(r) = sqrt(1 + (ep r)^2)
c = parms.get('centers')
#op = parms.get('operator','interp')
ep = parms.get('shapeparm',1)
DM = dmatrix(d, centers = c)
# eps_r = epsilon*r where epsilon may be an array or scalar
eps_r = dot(ep*eye(DM.shape[0]),DM)
return sqrt(1+(eps_r)**2)
