Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger2dCurve.m @ 699:8f1eae1450b2 feature/optim

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
sign change on hamiltonian
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 16 Oct 2017 14:52:15 +0200
parents 0122cbe2e6d3
children de7f665e1d00
comparison
equal deleted inserted replaced
698:0122cbe2e6d3 699:8f1eae1450b2
129 129
130 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); 130 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);
131 end 131 end
132 132
133 function [Hamiltonian,n] = h_fun(obj,t) 133 function [Hamiltonian,n] = h_fun(obj,t)
134 Hamiltonian = obj.c^2*(obj.DUU + obj.DUV + obj.DVU + obj.DVV)*sqrt(obj.Ji); 134 Hamiltonian = -obj.c^2*(obj.DUU + obj.DUV + obj.DVU + obj.DVV)*sqrt(obj.Ji);
135 end 135 end
136 136
137 function [D ]= variable_update(obj,t) 137 function [D ]= variable_update(obj,t)
138 % Metric derivatives 138 % Metric derivatives
139 if(obj.t_up == t) 139 if(obj.t_up == t)