Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger2dCurve.m @ 493:6b8297f66c91 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Commit before merge
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 23 Feb 2017 09:31:56 +0100
parents 6b95a894cbd7
children 4905446f165e
comparison
equal deleted inserted replaced
492:6b95a894cbd7 493:6b8297f66c91
139 J = x_u.*y_v - x_v.*y_u; 139 J = x_u.*y_v - x_v.*y_u;
140 a11 = 1./J.* (x_v.^2 + y_v.^2); 140 a11 = 1./J.* (x_v.^2 + y_v.^2);
141 a12 = -1./J .* (x_u.*x_v + y_u.*y_v); 141 a12 = -1./J .* (x_u.*x_v + y_u.*y_v);
142 a22 = 1./J .* (x_u.^2 + y_u.^2); 142 a22 = 1./J .* (x_u.^2 + y_u.^2);
143 143
144 obj.a11=a11; 144 obj.a11 = a11;
145 obj.a12=a12; 145 obj.a12 = a12;
146 obj.a22=a22; 146 obj.a22 = a22;
147 147
148 % Assemble full operators 148 % Assemble full operators
149 L_12 = spdiags(a12, 0, obj.m_tot, obj.m_tot); 149 L_12 = spdiags(a12, 0, obj.m_tot, obj.m_tot);
150 Duv = obj.Du*L_12*obj.Dv; 150 Duv = obj.Du*L_12*obj.Dv;
151 Dvu = obj.Dv*L_12*obj.Du; 151 Dvu = obj.Dv*L_12*obj.Du;