Mercurial > repos > public > sbplib

changeset 565:f4b0d0e84305 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add missing square
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Aug 2017 17:42:22 +0200
parents 9bf49338f8e6
children 9c98a0526afc
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:33:41 2017 +0200
+++ b/+scheme/LaplaceCurvilinear.m	Wed Aug 30 17:42:22 2017 +0200
@@ -185,10 +185,10 @@
             a22_n = spdiag(a22(I_n));
             a12_n = spdiag(a12(I_n));
 
-            s_w = sqrt((e_w'*x_v).^2 + (e_w'*y_v)');
-            s_e = sqrt((e_e'*x_v).^2 + (e_e'*y_v)');
-            s_s = sqrt((e_s'*x_u).^2 + (e_s'*y_u)');
-            s_n = sqrt((e_n'*x_u).^2 + (e_n'*y_u)');
+            s_w = sqrt((e_w'*x_v).^2 + (e_w'*y_v).^2);
+            s_e = sqrt((e_e'*x_v).^2 + (e_e'*y_v).^2);
+            s_s = sqrt((e_s'*x_u).^2 + (e_s'*y_u).^2);
+            s_n = sqrt((e_n'*x_u).^2 + (e_n'*y_u).^2);
 
             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')';