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)