Mercurial > repos > public > sbplib

diff +scheme/Laplace1d.m @ 1116:33c378e508d2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix minus sign in Laplace1d Dirichlet penalty.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 02 May 2019 13:25:14 -0700
parents a22b23098021
children
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'}