I am using libMesh to develop a multilevel solver. I understand the adaptive mesh in libMesh is hierarchical and embeded, which is good for multilevel/multigrid type of solvers. But I am wondering if there is an easy way to construct linear systems on meshes at different levels and represent vectors on those levels. That is essential for multilevel/multigrid methods. If the current implementation of libMesh is not easy for this, could you shed some light on what extension I can make to keep things easy?
      Thank you very much!