Feature request: Refiner should return old to new element mapping

It would be useful to know how elements in the coarse mesh map to elements in the fine mesh. For instance, when we are doing iterative refinement for QBX, we can avoid a lot of work by not considering elements that already pass the refinement criteria.