Mercurial > repos > public > sbplib

comparison +scheme/Wave2d.m @ 905:459eeb99130f feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Include type as (optional) input parameter in the interface method of all schemes.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 22 Nov 2018 22:03:44 -0800
parents a8ed986fcf57
children b9c98661ff5d
comparison
equal deleted inserted replaced
904:14b093a344eb 905:459eeb99130f
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