Schwarz Contact: enable reconstruction of Schwarz contact Neumann and Dirichlet projection operators if the interface changes size and/or position during contact

[ikalash]

The N^src and N^dist shape functions used to define the projection operators will change in this case.

In looking at the code, we observed that the compute_normal function in interpolation.jl is limited to flat contact faces currently when using higher order elements. This will need to be modified to enable a more general contact scenario.