Add in ApplyBCPatch macro system

[KoenigKrote]

In progress to #9 #14 #15

Note: Primary work and changes live in bc_faces.h and in richards_jacobian_eval.c, rest of macros are setting up loop changes to work with this system.

Current progress:

• Set up ForBCStructNumPatches macro
• Design a macro system for applying BC contributions (ApplyBCPatch)
• Set up BCStructPatchLoopX, GrGeomPatchLoopX, and GrGeomOctreeFaceLoopX macro to handle ApplyBCPatch macro
• Convert RichardsJacobian Symmetric contribution loop to use new macro system
• Convert RichardsJacobian Boundary Condition Contribution loops to use new macro system
• Convert RichardsJacobian JC Submatrix loop to use new macro system
• DirichletBC Pressure Contribution Macro #14 DirichletBC Pressure macro written, needs to be removed from this PR and into its own
• [?] Boundary Condition Macro System - NL Function Eval #15 Apply macro system to NlFunctionEval
  • Precondition / Postcondition macros
• Compiles correctly
• Document system
• Decide how to handle outer iteration spaces, specify them directly or use Do_ macro verbage?
• Decide on specific naming schemes
• Decide on specific header file organization
• Determine which macro systems can be used in both Richards and NL (i.e. where can we remove "Richards" from the Do_* macros and also use them in NlFunctionEval)
• Verify output
• Check performance changes (in serial)

[KoenigKrote]

It would make the most sense the finalize and merge adding boundary condition macros before continuing to review this procedure, since we want to add/merge additional OverlandFlow cases anyways