Mercurial > repos > public > sbplib
changeset 517:7a091a3527df feature/quantumTriangles
sign change in SAT-TERM
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Mon, 10 Jul 2017 09:27:58 +0200 |
| parents | afff85574ddb |
| children | 4709f2329372 |
| files | +scheme/Schrodinger2dCurve.m |
| diffstat | 1 files changed, 1 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/+scheme/Schrodinger2dCurve.m Fri Jul 07 14:52:47 2017 +0200 +++ b/+scheme/Schrodinger2dCurve.m Mon Jul 10 09:27:58 2017 +0200 @@ -136,7 +136,6 @@ else ti = parametrization.Ti.points(obj.p{1}(t),obj.p{2}(t),obj.p{3}(t),obj.p{4}(t)); ti_tau = parametrization.Ti.points(obj.p{5}(t),obj.p{6}(t),obj.p{7}(t),obj.p{8}(t)); - lcoords=points(obj.grid); [obj.xm,obj.ym]= ti.map(lcoords(1:obj.m_v:end,1),lcoords(1:obj.m_v,2)); [x_tau,y_tau]= ti_tau.map(lcoords(1:obj.m_v:end,1),lcoords(1:obj.m_v,2)); @@ -214,7 +213,7 @@ F = @(t)(s * a_n(t)*d_n' + s * a_t(t) *d_t')'; tau1 = 1; a = @(t)spdiag(g(t)); - tau2 = @(t) (-1*s*a(t))/2; + tau2 = @(t) (1*s*a(t))/2; penalty_parameter_1 = @(t) 1*1i*halfnorm_inv_n*halfnorm_inv_t*F(t)*e'*halfnorm_t*e; penalty_parameter_2 = @(t) halfnorm_inv_n*e*tau2(t);