Made BoundaryConditions functions templated in preparation for FEMContest API update.
Updated apply_neumann_normal. Need to think about the ones taking
a libMesh::Point.
Most of the conversion is done for FEMContext (I think). Need
to figure out how to deal with the cases where the flux is a
tensor.
Have BoundaryCondition methods take a RealGradient instead of a
Point. Still need to generalize, but this is the first step.
Generalized apply_neumann and apply_neumann_axisymmetric. Still
need to deal with the NeumannFuncObj cases.
Last fixes to get compilation with curent libMesh master. Still
need to deal with NeumannFuncObj. Going to solicit opinions.
Took @roystgnr's suggestion and added default typename for
BoundaryCondition methods. Went ahead and pulled out the template
part fo the function call where it was used.
Untabify.
Merge pull request #18 from pbauman/bound_cond_api_change

BoundaryConditions update to conform with FEMContext API update