Mercurial > repos > public > sbplib

changeset 67:446d67a49cd8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed some errors in scheme.Schrodinger.m
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 20 Nov 2015 10:28:27 +0100
parents dbc50fa58ca6
children f9a65e62c5e5
files +scheme/Schrodinger.m
diffstat 1 files changed, 3 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
diff -r dbc50fa58ca6 -r 446d67a49cd8 +scheme/Schrodinger.m
--- a/+scheme/Schrodinger.m	Fri Nov 20 10:24:21 2015 +0100
+++ b/+scheme/Schrodinger.m	Fri Nov 20 10:28:27 2015 +0100
@@ -46,7 +46,7 @@
 
             V_mat = spdiags(V_vec,0,m,m);
 
-            obj.D = 1i * obj.D2 - 1i * V;
+            obj.D = 1i * obj.D2 - 1i * V_mat;
 
             obj.m = m;
             obj.h = h;
@@ -96,11 +96,11 @@
             [e_u,d_u,s_u] = obj.get_boundary_ops(boundary);
             [e_v,d_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
 
-            a =  s* 1/2 * 1i ;
+            a =  -s_u* 1/2 * 1i ;
             b =  a';
 
             tau = b*d_u;
-            sig = a*e_u;
+            sig = -a*e_u;
 
             closure = obj.Hi * (tau*e_u' + sig*d_u');
             penalty = obj.Hi * (-tau*e_v' - sig*d_v');