Mercurial > repos > public > sbplib

diff +scheme/Beam.m @ 223:96dedf892910 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed sign error in beam BC penalties.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 28 Jun 2016 14:05:17 +0200
parents 5df8d20281fe
children 05304d7792dc
line wrap: on
line diff
--- a/+scheme/Beam.m	Tue Jun 28 13:21:02 2016 +0200
+++ b/+scheme/Beam.m	Tue Jun 28 14:05:17 2016 +0200
@@ -88,8 +88,8 @@
 
                     closure = obj.Hi*(tau*e' + sig*d1');
 
-                    penalty{1} = obj.Hi*tau;
-                    penalty{2} = obj.Hi*sig;
+                    penalty{1} = -obj.Hi*tau;
+                    penalty{2} = -obj.Hi*sig;
                 otherwise % Unknown, boundary condition
                     error('No such boundary condition: type = %s',type);
             end