xi = np.asarray(xi, np.float_)