We haven't taken a look at this code in a long time, but when doing multilake calibrations or different types of calibrations, the compare_to_field workflow is pretty slow. This is probably a low priority for the time being, but probably some easy wins using more modern packages to do this type of matching/interpolating within the output matrix.
Before starting on this it would be good to benchmark with an example.