Mercurial > repos > public > sbplib

comparison +scheme/Wave2d.m @ 955:92c3c170e90b

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merged in feature/utux2D (pull request #11) Feature/utux2D Approved-by: Jonatan Werpers <jonatan.werpers@it.uu.se> Approved-by: Martin Almquist <malmquist@stanford.edu> Approved-by: Vidar Stiernström <vidar.stiernstrom@it.uu.se>
author Jonatan Werpers <jonatan.werpers@it.uu.se>
date Mon, 10 Dec 2018 08:29:41 +0000
parents 706d1c2b4199
children 78db023a7fe3 c12b84fe9b00
comparison
equal deleted inserted replaced
951:a6c5e73ff44e 955:92c3c170e90b
156 otherwise 156 otherwise
157 error('No such boundary condition: type = %s',type); 157 error('No such boundary condition: type = %s',type);
158 end 158 end
159 end 159 end
160 160
161 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 161 function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type)
162 % u denotes the solution in the own domain 162 % u denotes the solution in the own domain
163 % v denotes the solution in the neighbour domain 163 % v denotes the solution in the neighbour domain
164 [e_u,d_u,s_u,gamm_u, halfnorm_inv] = obj.get_boundary_ops(boundary); 164 [e_u,d_u,s_u,gamm_u, halfnorm_inv] = obj.get_boundary_ops(boundary);
165 [e_v,d_v,s_v,gamm_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 165 [e_v,d_v,s_v,gamm_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
166 166