真实值为 x,测量值为 x0,
def rel_err(x, y): return np.max(np.abs(x-y)/np.maximum(1e-8, np.abs(x) + np.abs(y)))