eps = np.finfo(float).eps