Checkpointing work on BoundaryConditions conforming to updated NeumannFuncObj.
Might end up pulling out the flux vector API so we don't have to deal with
3rd order tensors...