def METHOD_NAME(x, low, high, axis=-1):
    m1 = convert_angles(
        circmean(x, low=low, high=high, axis=axis, nan_policy="omit"), low, high
    )
    m2 = circ_mean(convert_angles(x, low, high), axis=axis)
    assert (np.round(m1, 4) == np.round(m2, 4)).all()