An energy-minimizing interpolation for robust multigrid methods

We propose a robust interpolation for multigrid based on the concepts of en ergy minimization and approximation. It can handle PDE coefficients of vari ous types on structured or unstructured grids under one framework. The form ulation is general; it can be applied to any dimension. We demonstrate nume rically the effectiveness of the multigrid method in two dimensions by appl ying it to a discontinuous coefficient problem, an oscillatory coefficient problem, and an anisotropic problem. Empirically, the convergence rate is i ndependent of the coefficients of the underlying PDE, in addition to being independent of the mesh size. The proposed method is primarily designed for second-order elliptic PDEs, with possible extensions to other classes of p roblems such as integral equations.