The current functions to find neighborhoods for a given cell (both vertex and face neighbors) are hideously slow when considering the AMR mesh. This might be easiest to circumvent by using a VTK backend, with a couple of options:
Other suggestions welcome.