Mercurial > repos > public > sbplib

comparison +scheme/Beam.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 2b1b944deae1 25d0efdb0f75
comparison
equal deleted inserted replaced
854:18162a0a5bb5 1033:037f203b9bf5
17 gamm 17 gamm
18 delt 18 delt
19 alphaII 19 alphaII
20 alphaIII 20 alphaIII
21 21
22 opt 22 opt % TODO: Get rid of this and use the interface type instead
23 end 23 end
24 24
25 methods 25 methods
26 function obj = Beam(grid, order, alpha, opsGen, opt) 26 function obj = Beam(grid, order, alpha, opsGen, opt)
27 default_arg('alpha', -1); 27 default_arg('alpha', -1);
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