Mercurial > repos > public > sbplib

comparison +scheme/Wave2dCurve.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 706d1c2b4199
children
comparison
equal deleted inserted replaced
854:18162a0a5bb5 1033:037f203b9bf5
241 otherwise 241 otherwise
242 error('No such boundary condition: type = %s',type); 242 error('No such boundary condition: type = %s',type);
243 end 243 end
244 end 244 end
245 245
246 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 246 function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type)
247 % u denotes the solution in the own domain 247 % u denotes the solution in the own domain
248 % v denotes the solution in the neighbour domain 248 % v denotes the solution in the neighbour domain
249 tuning = 1.2; 249 tuning = 1.2;
250 % tuning = 20.2; 250 % tuning = 20.2;
251 [e_u, d_n_u, d_t_u, coeff_n_u, coeff_t_u, s_u, gamm_u, halfnorm_inv_u_n, halfnorm_inv_u_t, halfnorm_u_t, I_u] = obj.get_boundary_ops(boundary); 251 [e_u, d_n_u, d_t_u, coeff_n_u, coeff_t_u, s_u, gamm_u, halfnorm_inv_u_n, halfnorm_inv_u_t, halfnorm_u_t, I_u] = obj.get_boundary_ops(boundary);