def pure_nugget(theta, d):
 theta = np.asarray(theta, dtype=np.float64)
 d = np.asarray(d, dtype=np.float64)
 n_eval = d.shape[0]
 r = np.zeros(n_eval)
 r[np.all((d == 0.0), axis=1)] = 1.0
 return r

def pure_nugget ( theta , d ) : theta = np . asarray ( theta , dtype = np . float64 ) d = np . asarray ( d , dtype = np . float64 ) n_eval = d . shape [ 0 ] r = np . zeros ( n_eval ) r [ np . all ( ( d == 0.0 ) , axis = 1 ) ] = 1.0 return r
spatial independence correlation model .