Mercurial > repos > public > sbplib

changeset 560:6132c52bf923 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 29 Aug 2017 13:04:07 +0200
parents e3f587bfd7f8 (diff) 54c775c3348a (current diff)
children 3a13916f8ff0
files +scheme/LaplaceCurvilinear.m
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinear.m	Tue Aug 29 12:48:30 2017 +0200
+++ b/+scheme/LaplaceCurvilinear.m	Tue Aug 29 13:04:07 2017 +0200
@@ -185,10 +185,10 @@
             a22_n = spdiag(a22(I_n));
             a12_n = spdiag(a12(I_n));
 
-            d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')';
-            d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')';
-            d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')';
-            d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')';
+            obj.d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')';
+            obj.d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')';
+            obj.d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')';
+            obj.d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')';
 
             obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;
             obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;