def METHOD_NAME(dx, dy): 
    pp = np.mgrid[2*dx:1+dx:dx, 2*dy:2+dy:dy].reshape(2, -1).T
    x0 = np.arange(0, 4+dx, dx)
    x1 = np.arange(-dx/2, 4+dx+dx/2, dx) 
    y = np.arange(dy, 4+dy, dx)
    bp0 = np.column_stack((x0, np.zeros_like(x0)))
    bp1 = np.column_stack((x1, np.full_like(x1, dx)))
    bp = np.vstack((bp0,bp1))
    lp0 = np.column_stack((np.zeros_like(y), y))
    lp1 = np.column_stack((np.full_like(y, dy), y+dy/2))
    lp = np.vstack((lp0,lp1))
    rp0 = np.column_stack((np.full_like(y, 4), y)) 
    rp1 = np.column_stack((np.full_like(y, 4-dy), y+dy/2)) 
    rp = np.vstack((rp0,rp1))
    bpp = np.vstack((bp,lp,rp))
    return pp, bpp