Found a small bug in integration. Refactored the basis functions to
prepare them for doing high order mappings.
Refactoring of the mesh interface to account for varying normal vector
and jacobian.
Higher order mappings are semi-functional. Still have some bugs!
Merge with isoparametric
Dislocations seem to function properly.
Updated some non-code.