This leads (for example) to us doing unnecessary matvecs (`resample_by_mat`) instead of indirect memory access (`resample_by_picking`) in nodal DG. Solution started in #105 cc @nchristensen