This code is a modification of MFEM, a C++ library for finite element methods. The modification includes: 4D meshes and finite elements, multigrid and time-slabbing solvers and a natural setting for considering (constrained) first-order system least-squares formulations.