Mercurial > repos > public > sbplib

comparison +scheme/Beam2d.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 d095b5396103
children b9c98661ff5d
comparison
equal deleted inserted replaced
904:14b093a344eb 905:459eeb99130f
159 otherwise 159 otherwise
160 error('No such boundary condition: type = %s',type); 160 error('No such boundary condition: type = %s',type);
161 end 161 end
162 end 162 end
163 163
164 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 164 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
165 % u denotes the solution in the own domain 165 % u denotes the solution in the own domain
166 % v denotes the solution in the neighbour domain 166 % v denotes the solution in the neighbour domain
167 [e_u,d1_u,d2_u,d3_u,s_u,gamm_u,delt_u, halfnorm_inv] = obj.get_boundary_ops(boundary); 167 [e_u,d1_u,d2_u,d3_u,s_u,gamm_u,delt_u, halfnorm_inv] = obj.get_boundary_ops(boundary);
168 [e_v,d1_v,d2_v,d3_v,s_v,gamm_v,delt_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 168 [e_v,d1_v,d2_v,d3_v,s_v,gamm_v,delt_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
169 169