Mercurial > repos > public > sbplib

comparison +scheme/Beam.m @ 944:a35ed1d124d3 feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change from opts to type in a few schemes
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 05 Dec 2018 10:52:37 +0100
parents b9c98661ff5d
children 2b1b944deae1 25d0efdb0f75
comparison
equal deleted inserted replaced
943:21394c78c72e 944:a35ed1d124d3
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,opts) 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