Be able to construct global dof index to derivative map for given dual number

[lindsayad]

Reason

@joshuahansel is coupling physics on disconnected different-dimension meshes (I think) where the only way he can get the Jacobians right in an AD setting is to be able to map derivatives using global dof indices (there is nothing local about his dependence)

Design

Create a utility that takes a dual number, DofMap, and element, and generates the map from global dof indices to derivatives. This is admittedly a slow way to do this. If we want to optimize this later, we could at a minimum generate a local to global index mapping for the current element one time and then use that map for all dual numbers supplied as arguments to it.

Impact

Help users (@joshuahansel ) do non-local derivative coupling