Mercurial > repos > public > sbplib

changeset 520:f235284e2eb1 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
changed sign in penalty parameter
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 07 Aug 2017 13:20:48 +0200
parents 0de024556427
children 527a7b954f26
files +scheme/Schrodinger2dCurve.m
diffstat 1 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Schrodinger2dCurve.m	Mon Aug 07 11:56:53 2017 +0200
+++ b/+scheme/Schrodinger2dCurve.m	Mon Aug 07 13:20:48 2017 +0200
@@ -125,7 +125,7 @@
             %  1i*obj.c^2*(obj.DUU + obj.DUV + obj.DVU + obj.DVV)); (ols
             %  not skew sym disc
             
-            D = sqrt(obj.Ji)*(-1/2*(obj.b1*obj.Du + obj.Du*obj.b1) - 1/2*(obj.b2*obj.Dv + obj.Dv*obj.b2) + 1i*obj.c^2*(obj.DUU + 0*obj.DUV + 0*obj.DVU + 0*obj.DVV))*sqrt(obj.Ji); 
+            D = sqrt(obj.Ji)*(-1/2*(obj.b1*obj.Du + obj.Du*obj.b1) - 1/2*(obj.b2*obj.Dv + obj.Dv*obj.b2) + 1i*obj.c^2*(obj.DUU + obj.DUV + obj.DVU + obj.DVV))*sqrt(obj.Ji); 
         end
         
         
@@ -220,7 +220,7 @@
                     penalty_parameter_2 = @(t) halfnorm_inv_n*e*tau2(t);
 
                     closure = @(t)  sqrt(obj.Ji)*(obj.c^2 * penalty_parameter_1(t)*e' +  penalty_parameter_2(t)*e')*sqrt(obj.Ji);
-                    penalty = @(t) -sqrt(obj.Ji)*(obj.c^2 * penalty_parameter_1(t)*e' -  penalty_parameter_2(t)*e')*sqrt(obj.Ji);
+                    penalty = @(t) -sqrt(obj.Ji)*(obj.c^2 * penalty_parameter_1(t)*e' +  penalty_parameter_2(t)*e')*sqrt(obj.Ji);
                 
         end