Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger2dCurve.m @ 519:0de024556427 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
removed y-direction for debugging
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 07 Aug 2017 11:56:53 +0200
parents 4709f2329372
children f235284e2eb1
comparison
equal deleted inserted replaced
518:4709f2329372 519:0de024556427
123 % D = obj.Ji*(-1/2*(obj.b1*obj.Du-obj.b1_u+obj.Du*obj.b1) - 123 % D = obj.Ji*(-1/2*(obj.b1*obj.Du-obj.b1_u+obj.Du*obj.b1) -
124 % 1/2*(obj.b2*obj.Dv - obj.b2_v +obj.Dv*obj.b2) + 124 % 1/2*(obj.b2*obj.Dv - obj.b2_v +obj.Dv*obj.b2) +
125 % 1i*obj.c^2*(obj.DUU + obj.DUV + obj.DVU + obj.DVV)); (ols 125 % 1i*obj.c^2*(obj.DUU + obj.DUV + obj.DVU + obj.DVV)); (ols
126 % not skew sym disc 126 % not skew sym disc
127 127
128 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); 128 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);
129 end 129 end
130 130
131 131
132 function [D ]= variable_update(obj,t) 132 function [D ]= variable_update(obj,t)
133 % Metric derivatives 133 % Metric derivatives