Mercurial > repos > public > sbplib
diff +scheme/Schrodinger2dCurve.m @ 708:acb58769610e feature/quantumTriangles
fixed error in diffOp
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Tue, 21 Nov 2017 16:51:51 +0100 |
| parents | e89715fe6a6e |
| children | f004b9e9d17a |
line wrap: on
line diff
--- a/+scheme/Schrodinger2dCurve.m Fri Nov 10 14:22:56 2017 +0100 +++ b/+scheme/Schrodinger2dCurve.m Tue Nov 21 16:51:51 2017 +0100 @@ -211,10 +211,9 @@ % neighbour_boundary is a string specifying which boundary to interface to. function [closure, penalty,closureHamiltonian,penaltyHamiltonian] = boundary_condition(obj, boundary,~) [e, d_n, d_t, coeff_t, coeff_n s, halfnorm_inv_n, halfnorm_inv_t, halfnorm_t,g] = obj.get_boundary_ops(boundary); - a_t = @(t) spdiag(coeff_t(t)); a_n = @(t) spdiag(coeff_n(t)); - + boundary F = @(t)(s * a_n(t)*d_n' + s * a_t(t) *d_t')'; tau1 = 1; @@ -258,8 +257,10 @@ penalty_parameter_3 = @(t) halfnorm_inv_u_n * e_u * (gamm(t) ); - closure =@(t) sqrt(Ji_u)*obj.c^2 * ( penalty_parameter_1(t)*e_u' + penalty_parameter_2(t)*F_u(t)' + penalty_parameter_3(t)*e_u')*sqrt(Ji_u); - penalty =@(t) sqrt(Ji_u)*obj.c^2 * ( -penalty_parameter_1(t)*e_v' - penalty_parameter_2(t)*F_v(t)' - penalty_parameter_3(t)*e_v')*sqrt(Ji_v); + closure =@(t) sqrt(Ji_u)*obj.c^2 * ( penalty_parameter_1(t)*e_u'... + + penalty_parameter_2(t)*F_u(t)' + penalty_parameter_3(t)*e_u')*sqrt(Ji_u); + penalty =@(t) sqrt(Ji_u)*obj.c^2 * ( -penalty_parameter_1(t)*e_v'... + - penalty_parameter_2(t)*F_v(t)' - penalty_parameter_3(t)*e_v')*sqrt(Ji_v); end