Finite element discretizations of problems in computational physics often
rely on adaptive mesh refinement (AMR) to preferentially resolve regions
containing important features during simulation. However, these spatial
refinement strategies are often heuristic and rely on domain-specif