Added Function::evaluate() variant with 3D support (for testing).
Fixed a bug in the BC imposition, and added more debugging output.  Not yet working for non-constant exact solutions.