--- a/+multiblock/Laplace.m	Sun Nov 04 12:36:30 2018 -0800
+++ b/+multiblock/Laplace.m	Wed Nov 21 18:29:29 2018 -0800
@@ -21,7 +21,7 @@
 
             obj.D = obj.mbDiffOp.D;
             obj.J = obj.jacobian();
-            obj.H = obj.mbDiffOp.H * obj.jacobian();
+            obj.H = obj.mbDiffOp.H;
         end
 
         function s = size(obj)
@@ -42,6 +42,10 @@
             op = getBoundaryOperator(obj.mbDiffOp, opName, boundary);
         end
 
+        function op = getBoundaryQuadrature(obj, boundary)
+            op = getBoundaryQuadrature(obj.mbDiffOp, boundary);
+        end
+
         function [closure, penalty] = boundary_condition(obj,boundary,type) % TODO: Change name to boundaryCondition
             [closure, penalty] = boundary_condition(obj.mbDiffOp, boundary, type);
         end