Mercurial > repos > public > sbplib
comparison +scheme/Beam.m @ 223:96dedf892910 feature/beams
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 |
comparison
equal
deleted
inserted
replaced
222:e7e73173d44d | 223:96dedf892910 |
---|---|
86 tau = tau1*e+tau4*d3; | 86 tau = tau1*e+tau4*d3; |
87 sig = sig2*d1+sig3*d2; | 87 sig = sig2*d1+sig3*d2; |
88 | 88 |
89 closure = obj.Hi*(tau*e' + sig*d1'); | 89 closure = obj.Hi*(tau*e' + sig*d1'); |
90 | 90 |
91 penalty{1} = obj.Hi*tau; | 91 penalty{1} = -obj.Hi*tau; |
92 penalty{2} = obj.Hi*sig; | 92 penalty{2} = -obj.Hi*sig; |
93 otherwise % Unknown, boundary condition | 93 otherwise % Unknown, boundary condition |
94 error('No such boundary condition: type = %s',type); | 94 error('No such boundary condition: type = %s',type); |
95 end | 95 end |
96 end | 96 end |
97 | 97 |