def METHOD_NAME():
    shape_v = [100, 100]
    cenv = shape_v[0] // 2
    r = 20
    a = np.zeros(shape_v)
    a[cenv - r : cenv + r, cenv - r : cenv + r] = 1.0
    sw_sigma = 0.50
    sw_threshold = 0.05
    new_sup_index = find_support(a, sw_sigma, sw_threshold)
    new_sup = np.zeros_like(a)
    new_sup[new_sup_index] = 1
    # the area of new support becomes larger
    assert np.sum(new_sup) == 1760