def err_metric(a, b, atol=1e-290):
    m = abs(b - a) / (atol + abs(b))
    m[np.isinf(b) & (a == b)] = 0
    return m
dataset/ETHPy150Open scipy/scipy/scipy/special/c_misc/struve_convergence.py/err_metric/swapped operands of `-`