Mercurial > repos > public > sbplib
diff +scheme/Schrodinger2dCurve.m @ 695:ba0d31ce4121 feature/quantumTriangles
Changed jacobian at right side of the penalty term in 2D, seems like it is working
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Tue, 19 Sep 2017 09:51:08 +0200 |
| parents | 1157375c678a |
| children | 0122cbe2e6d3 |
line wrap: on
line diff
--- a/+scheme/Schrodinger2dCurve.m Tue Sep 19 09:47:14 2017 +0200 +++ b/+scheme/Schrodinger2dCurve.m Tue Sep 19 09:51:08 2017 +0200 @@ -251,11 +251,11 @@ 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_v)*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); + 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 - function [e, d_n, d_t, coeff_t,coeff_n, s, halfnorm_inv_n, halfnorm_inv_t, halfnorm_t,g] = get_boundary_ops(obj, boundary) + function [e, d_n, d_t, coeff_t,coeff_n, s, halfnorm_inv_n, halfnorm_inv_t, halfnorm_t,g,Ji] = get_boundary_ops(obj, boundary) % gridMatrix = zeros(obj.m(2),obj.m(1)); % gridMatrix(:) = 1:numel(gridMatrix);