Resolution of singularities of an idealistic filtration in dimension 3 after Benito-Villamayor

We establish an algorithm for resolution of singularities of an idealistic filtration in dimension 3 (at the local level) in positive characteristic, incorporating the method recently developed by Benito-Villamayor into our framework. Although (a global version of) our algorithm only implies embedded resolution of surfaces in the smooth ambient space of dimension 3, a classical result known before, we introduce some new invariant which effectively measures how much singularities are improved in the process of our algorithm and which strictly drops after each blow up. This is in contrast to the well-known Abhyankar-Moh pathology of the increase of the residual order under blow up and the phenomenon of the "Kangaroo" points observed by Hauser.