In the context of hardware-oriented numerics, grid deformation, e.g. the relocation of grid points preserving the topology of the mesh, is a viable way for grid adaption and offers an alternative to the widespread h-adaptivity. We present a fast and accurate grid deformation method which -unlike the common MMPDE methods- works independently of the underlying problem and thus acts as black-box tool. It requires the solution of one global Poisson problem and one ODE for each grid point only. We show the convergence of our method. Moreover, exploiting the given grid hierarchy, this grid deformation method is of almost optimal complexity with respect to the number of vertices in the mesh.