Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger2dCurve.m @ 709:f004b9e9d17a feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
changed a sign in penalty and closure, most probably everything blows up now...
author Ylva Rydin <ylva.rydin@telia.com>
date Tue, 21 Nov 2017 17:04:07 +0100
parents acb58769610e
children 023a95f63950
comparison
equal deleted inserted replaced
708:acb58769610e 709:f004b9e9d17a
248 248
249 a = @(t)spdiag(gamm_u(t)); 249 a = @(t)spdiag(gamm_u(t));
250 250
251 tau = s_u*1*1i/2; 251 tau = s_u*1*1i/2;
252 sig = -s_u*1*1i/2; 252 sig = -s_u*1*1i/2;
253 gamm = @(t) (-s_u*a(t))/2; 253 gamm = @(t) (s_u*a(t))/2;
254 254
255 penalty_parameter_1 = @(t) halfnorm_inv_u_n*(tau*halfnorm_inv_u_t*F_u(t)*e_u'*halfnorm_u_t*e_u); 255 penalty_parameter_1 = @(t) halfnorm_inv_u_n*(tau*halfnorm_inv_u_t*F_u(t)*e_u'*halfnorm_u_t*e_u);
256 penalty_parameter_2 = @(t) halfnorm_inv_u_n * e_u * (sig ); 256 penalty_parameter_2 = @(t) halfnorm_inv_u_n * e_u * (sig );
257 penalty_parameter_3 = @(t) halfnorm_inv_u_n * e_u * (gamm(t) ); 257 penalty_parameter_3 = @(t) halfnorm_inv_u_n * e_u * (gamm(t) );
258 258