Mercurial > repos > public > sbplib

comparison +scheme/Beam2d.m @ 1033:037f203b9bf5 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with branch feature/advectioRV to utilize the +rv package
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 17 Jan 2019 10:44:12 +0100
parents a35ed1d124d3
children 78db023a7fe3
comparison
equal deleted inserted replaced
854:18162a0a5bb5 1033:037f203b9bf5
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