Mercurial > repos > public > sbplib

comparison +scheme/Beam.m @ 968:a4ad90b37998 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 23 Dec 2018 14:39:31 +0100
parents a35ed1d124d3
children 2b1b944deae1 25d0efdb0f75
comparison
equal deleted inserted replaced
967:368a2773f78b 968:a4ad90b37998
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