Mercurial > repos > public > sbplib

changeset 564:9bf49338f8e6 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in creation of boundary inner products
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Aug 2017 17:33:41 +0200
parents 8f780429dab2
children f4b0d0e84305
files +scheme/LaplaceCurvilinear.m
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinear.m	Wed Aug 30 17:29:33 2017 +0200
+++ b/+scheme/LaplaceCurvilinear.m	Wed Aug 30 17:33:41 2017 +0200
@@ -199,10 +199,10 @@
             obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;
 
             %% Boundary inner products
-            H_w = H_v*s_w;
-            H_e = H_v*s_e;
-            H_s = H_u*s_s;
-            H_n = H_u*s_n;
+            obj.H_w = H_v*spdiag(s_w);
+            obj.H_e = H_v*spdiag(s_e);
+            obj.H_s = H_u*spdiag(s_s);
+            obj.H_n = H_u*spdiag(s_n);
 
             % Misc.
             obj.m = m;