Mercurial > repos > public > sbplib

comparison +scheme/Wave2dCurve.m @ 360:447ceb41fb65 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed penalty parameter in Wave2Dcurve.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 08 Dec 2016 11:11:28 +0100
parents 85c2fe06d551
children 32df00102268
comparison
equal deleted inserted replaced
359:0e198af1184b 360:447ceb41fb65
179 b1 = gamm*u.lambda./u.a11.^2; 179 b1 = gamm*u.lambda./u.a11.^2;
180 b2 = gamm*u.lambda./u.a22.^2; 180 b2 = gamm*u.lambda./u.a22.^2;
181 181
182 tau = -1./b1 - 1./b2; 182 tau = -1./b1 - 1./b2;
183 tau = tuning * spdiag(tau); 183 tau = tuning * spdiag(tau);
184 sig1 = 1/2; 184 sig1 = 1;
185 185
186 penalty_parameter_1 = halfnorm_inv_n*(tau + sig1*halfnorm_inv_t*F*e'*halfnorm_t)*e; 186 penalty_parameter_1 = halfnorm_inv_n*(tau + sig1*halfnorm_inv_t*F*e'*halfnorm_t)*e;
187 187
188 closure = obj.Ji*obj.c^2 * penalty_parameter_1*e'; 188 closure = obj.Ji*obj.c^2 * penalty_parameter_1*e';
189 penalty = -obj.Ji*obj.c^2 * penalty_parameter_1; 189 penalty = -obj.Ji*obj.c^2 * penalty_parameter_1;