--- a/+scheme/LaplaceCurvilinearNewCorner.m	Sun Mar 03 15:00:42 2019 -0800
+++ b/+scheme/LaplaceCurvilinearNewCorner.m	Sat May 04 14:44:34 2019 -0700
@@ -63,13 +63,13 @@
             % assert(isa(g, 'grid.Curvilinear'))
             if isa(a, 'function_handle')
                 a = grid.evalOn(g, a);
-                a = spdiag(a);
             end
+            a = spdiag(a);
 
             if isa(b, 'function_handle')
                 b = grid.evalOn(g, b);
-                b = spdiag(b);
             end
+            b = spdiag(b);
 
             % If b is scalar
             if length(b) == 1
@@ -305,7 +305,7 @@
                 case {'D','d','dirichlet'}
                     tuning = 1.0;
 
-                    sigma = 0;
+                    sigma = 0*b_b;
                     for i = 1:obj.dim
                         sigma = sigma + e'*J*K{i,m}*K{i,m}*e;
                     end
@@ -401,8 +401,8 @@
             b_b_v = e_v'*v.b*e_v;
 
             %--- Penalty strength tau -------------
-            sigma_u = 0;
-            sigma_v = 0;
+            sigma_u = 0*b_b_u;
+            sigma_v = 0*b_b_v;
             for i = 1:obj.dim
                 sigma_u = sigma_u + e_u'*J_u*K_u{i,m_u}*K_u{i,m_u}*e_u;
                 sigma_v = sigma_v + e_v'*J_v*K_v{i,m_v}*K_v{i,m_v}*e_v;