Mercurial > repos > public > sbplib

diff +scheme/LaplaceCurvilinear.m @ 1110:9ddd2e89d018 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix in interfaceNonconforming(). Remove unused properties gamm_u, gamm_v.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 10 Apr 2019 13:27:50 -0700
parents 01d28cfafe7c
children 0e5dc5dfa559
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinear.m	Wed Apr 10 11:28:23 2019 -0700
+++ b/+scheme/LaplaceCurvilinear.m	Wed Apr 10 13:27:50 2019 -0700
@@ -47,8 +47,6 @@
 
         % Temporary, only used for nonconforming interfaces but should be removed.
         lambda
-        gamm_u, gamm_v
-
     end
 
     methods
@@ -262,8 +260,6 @@
 
             % Temporary
             obj.lambda = lambda;
-            obj.gamm_u = h_u*ops_u.borrowing.M.d1;
-            obj.gamm_v = h_v*ops_v.borrowing.M.d1;
         end
 
 
@@ -429,17 +425,17 @@
 
             % u denotes the solution in the own domain
             % v denotes the solution in the neighbour domain
-            e_u    = obj.getBoundaryOperator('e', boundary);
-            d_u    = obj.getBoundaryOperator('d', boundary);
-            H_b_u  = obj.getBoundaryQuadrature(boundary);
-            I_u    = obj.getBoundaryIndices(boundary);
-            gamm_u = obj.getBoundaryBorrowing(boundary);
+            e_u         = obj.getBoundaryOperator('e', boundary);
+            d_u         = obj.getBoundaryOperator('d', boundary);
+            H_b_u       = obj.getBoundaryQuadrature(boundary);
+            I_u         = obj.getBoundaryIndices(boundary);
+            [~, gamm_u] = obj.getBoundaryBorrowing(boundary);
 
-            e_v    = neighbour_scheme.getBoundaryOperator('e', neighbour_boundary);
-            d_v    = neighbour_scheme.getBoundaryOperator('d', neighbour_boundary);
-            H_b_v  = neighbour_scheme.getBoundaryQuadrature(neighbour_boundary);
-            I_v    = neighbour_scheme.getBoundaryIndices(neighbour_boundary);
-            gamm_v = neighbour_scheme.getBoundaryBorrowing(neighbour_boundary);
+            e_v         = neighbour_scheme.getBoundaryOperator('e', neighbour_boundary);
+            d_v         = neighbour_scheme.getBoundaryOperator('d', neighbour_boundary);
+            H_b_v       = neighbour_scheme.getBoundaryQuadrature(neighbour_boundary);
+            I_v         = neighbour_scheme.getBoundaryIndices(neighbour_boundary);
+            [~, gamm_v] = neighbour_scheme.getBoundaryBorrowing(neighbour_boundary);
 
 
             % Find the number of grid points along the interface