Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger2dCurve.m @ 518:4709f2329372 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added J to propertied
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 07 Aug 2017 09:02:49 +0200
parents 7a091a3527df
children 0de024556427
comparison
equal deleted inserted replaced
517:7a091a3527df 518:4709f2329372
35 t_up 35 t_up
36 36
37 a11, a12, a22 37 a11, a12, a22
38 m_tot, m_u, m_v 38 m_tot, m_u, m_v
39 p 39 p
40 Ji 40 Ji, J
41 end 41 end
42 42
43 methods 43 methods
44 function obj = Schrodinger2dCurve(g ,order, opSet,p) 44 function obj = Schrodinger2dCurve(g ,order, opSet,p)
45 default_arg('opSet',@sbp.D2Variable); 45 default_arg('opSet',@sbp.D2Variable);
151 x_v = obj.Dv*x; 151 x_v = obj.Dv*x;
152 y_u = obj.Du*y; 152 y_u = obj.Du*y;
153 y_v = obj.Dv*y; 153 y_v = obj.Dv*y;
154 154
155 J = (x_u.*y_v - x_v.*y_u); 155 J = (x_u.*y_v - x_v.*y_u);
156 obj.J = spdiags(J, 0, obj.m_tot, obj.m_tot);
156 157
157 Ji = spdiags(1./J, 0, obj.m_tot, obj.m_tot); 158 Ji = spdiags(1./J, 0, obj.m_tot, obj.m_tot);
158 obj.Ji = Ji; 159 obj.Ji = Ji;
159 160
160 a11 = Ji* (x_v.^2 + y_v.^2); 161 a11 = Ji* (x_v.^2 + y_v.^2);