Mercurial > repos > public > sbplib

comparison +scheme/Beam.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 4ced7d47bd1f
children b9c98661ff5d
comparison
equal deleted inserted replaced
904:14b093a344eb 905:459eeb99130f
168 otherwise % Unknown, boundary condition 168 otherwise % Unknown, boundary condition
169 error('No such boundary condition: type = %s',type); 169 error('No such boundary condition: type = %s',type);
170 end 170 end
171 end 171 end
172 172
173 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 173 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
174 % u denotes the solution in the own domain 174 % u denotes the solution in the own domain
175 % v denotes the solution in the neighbour domain 175 % v denotes the solution in the neighbour domain
176 [e_u,d1_u,d2_u,d3_u,s_u] = obj.get_boundary_ops(boundary); 176 [e_u,d1_u,d2_u,d3_u,s_u] = obj.get_boundary_ops(boundary);
177 [e_v,d1_v,d2_v,d3_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 177 [e_v,d1_v,d2_v,d3_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
178 178