comparison +scheme/Wave2d.m @ 1018:eca4d9df9365 feature/advectionRV

Merge with default for latest changes
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Dec 2018 17:31:28 +0100
parents 706d1c2b4199
children 78db023a7fe3 c12b84fe9b00
comparison
equal deleted inserted replaced
1017:2d7c1333bd6c 1018:eca4d9df9365
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