changeset 1116:33c378e508d2

Bugfix minus sign in Laplace1d Dirichlet penalty.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 02 May 2019 13:25:14 -0700
parents 3230e4cbdbb4
children 433c89bf19e0 2924b3a9b921 dfb1b84b63f4 2a3bd78edb0e 0652b34f9f27
files +scheme/Laplace1d.m
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Laplace1d.m	Wed Apr 10 22:11:11 2019 +0200
+++ b/+scheme/Laplace1d.m	Thu May 02 13:25:14 2019 -0700
@@ -70,7 +70,7 @@
                     tau = tau1*e + tau2*d;
 
                     closure = obj.a*obj.Hi*tau*e';
-                    penalty = obj.a*obj.Hi*tau;
+                    penalty = -obj.a*obj.Hi*tau;
 
                 % Neumann boundary condition
                 case {'N','n','neumann'}