Mercurial > repos > public > sbplib

diff +scheme/Schrodinger1dCurve.m @ 507:1341e8580d64 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Start over with the non-symmetric shrodnger to ceep track of changes.
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 29 May 2017 14:07:55 +0200
parents 2047ff9640b3
children 508b7493be94
line wrap: on
line diff
--- a/+scheme/Schrodinger1dCurve.m	Mon May 29 13:54:16 2017 +0200
+++ b/+scheme/Schrodinger1dCurve.m	Mon May 29 14:07:55 2017 +0200
@@ -17,7 +17,6 @@
         a
         a_xi
         Ji
-        J
         t_up
         x
         
@@ -86,8 +85,8 @@
         
         function [D] = d_fun(obj,t)
             obj.variable_update(t); % In driscretization?
-            D = sqrt(obj.Ji)*(-0.5*(obj.D1*obj.a*sqrt(obj.Ji) + obj.a*obj.D1*sqrt(obj.Ji)) + 1i*obj.D2(diag(obj.Ji))*obj.Ji + 1i*obj.V_mat);
-          %   D= obj.Ji*(-sqrt(obj.a)*obj.D1*sqrt(obj.a) + 0.5*obj.a_xi + 1i*obj.D2(diag(obj.Ji)) + 1i*obj.V_mat);
+            D= obj.Ji*(-0.5*(obj.D1*obj.a - obj.a_xi + obj.a*obj.D1) + 1i*obj.D2(diag(obj.Ji)) + 1i*obj.V_mat);
+            
         end
         
         
@@ -103,7 +102,6 @@
                 obj.a = sparse(diag((-ddt_x_l*( x_r - x_l) - (obj.x-x_l)*(ddt_x_r-ddt_x_l))/(x_r-x_l)));            
                 
                 obj.Ji = sparse(diag(1./(x_r - x_l + 0*obj.x)));
-                obj.J = sparse(x_r -x_l);
                 obj.a_xi = sparse(diag(-1*(ddt_x_r - ddt_x_l + 0*obj.x)));
                 obj.t_up = t;
             end
@@ -118,17 +116,17 @@
             switch type
                 % Dirichlet boundary condition
                 case {'D','d','dirichlet'}
-                    tau1 = @(t) s * 1i*obj.Ji(p,p)*d;
-                    tau2 = @(t) (1*s*obj.a(p,p))/2*e;
-                    closure = @(t)(obj.Ji*obj.Hi*tau1(t)*e' + obj.Hi*tau2(obj.a)*e');
+                    tau1 = @(t) s * 1i*obj.Ji(p,p)^2*d;
+                    tau2 = @(t) obj.Ji*(-1*s*obj.a(p,p) - abs(obj.a(p,p)))/4*e;
+                    closure = @(t) obj.Hi*tau1(t)*e' + obj.Hi*tau2(obj.a)*e';
                     
                     switch class(data)
                         case 'double'
-                            penalty = @(t) -(obj.Hi*tau1*data+obj.Hi*tau2(obj.a)*data);
+                            penalty = @(t) -obj.Ji*(obj.Hi*tau1*data+obj.Hi*tau2(obj.a)*data);
                             %                      case 'function_handle'
                             %                           penalty = @(t)-obj.Hi*tau*data(t);
                         otherwise
-                            error('Weird data argument!')
+                            error('Wierd data argument!')
                     end
                     
                     % Unknown, boundary condition