Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger1dCurve.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 25053554524b
children b91d23271481
comparison
equal deleted inserted replaced
492:6b95a894cbd7 493:6b8297f66c91
105 end 105 end
106 106
107 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 107 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
108 % u denotes the solution in the own domain 108 % u denotes the solution in the own domain
109 % v denotes the solution in the neighbour domain 109 % v denotes the solution in the neighbour domain
110 % [e_u,d_u,s_u] = obj.get_boundary_ops(boundary); 110 [e_u,d_u,s_u] = obj.get_boundary_ops(boundary);
111 % [e_v,d_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 111 [e_v,d_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
112 112
113 % a = -s_u* 1/2 * 1i ; 113 a = -s_u* 1/2 * 1i ;
114 % b = a'; 114 b = a';
115 115
116 % tau = b*d_u; 116 tau = b*d_u;
117 % sig = -a*e_u; 117 sig = -a*e_u;
118 118
119 % closure = obj.Hi * (tau*e_u' + sig*d_u'); 119 closure = obj.Hi * (tau*e_u' + sig*d_u');
120 % penalty = obj.Hi * (-tau*e_v' - sig*d_v'); 120 penalty = obj.Hi * (-tau*e_v' - sig*d_v');
121 end 121 end
122 122
123 % Ruturns the boundary ops and sign for the boundary specified by the string boundary. 123 % Ruturns the boundary ops and sign for the boundary specified by the string boundary.
124 % The right boundary is considered the positive boundary 124 % The right boundary is considered the positive boundary
125 function [e,d,s,p] = get_boundary_ops(obj,boundary) 125 function [e,d,s,p] = get_boundary_ops(obj,boundary)