changeset 563:8f780429dab2 feature/grids/laplace_refactor

Fix sign error in Neumann BC
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Aug 2017 17:29:33 +0200
parents 11d8d6ccbcd7
children 9bf49338f8e6
files +scheme/LaplaceCurvilinear.m
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinear.m	Wed Aug 30 15:43:51 2017 +0200
+++ b/+scheme/LaplaceCurvilinear.m	Wed Aug 30 17:29:33 2017 +0200
@@ -260,7 +260,7 @@
                 case {'N','n','neumann'}
                     tau1 = -1;
                     tau2 = 0;
-                    tau = s*obj.a*obj.Ji*(tau1*e + tau2*d);
+                    tau = obj.a*obj.Ji*(tau1*e + tau2*d);
 
                     closure = halfnorm_inv*tau*d';
                     penalty = -halfnorm_inv*tau;