Mercurial > repos > public > sbplib

comparison +scheme/Beam2d.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 78db023a7fe3
comparison
equal deleted inserted replaced
967:368a2773f78b 968:a4ad90b37998
159 otherwise 159 otherwise
160 error('No such boundary condition: type = %s',type); 160 error('No such boundary condition: type = %s',type);
161 end 161 end
162 end 162 end
163 163
164 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 164 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary, type)
165 % u denotes the solution in the own domain 165 % u denotes the solution in the own domain
166 % v denotes the solution in the neighbour domain 166 % v denotes the solution in the neighbour domain
167 [e_u,d1_u,d2_u,d3_u,s_u,gamm_u,delt_u, halfnorm_inv] = obj.get_boundary_ops(boundary); 167 [e_u,d1_u,d2_u,d3_u,s_u,gamm_u,delt_u, halfnorm_inv] = obj.get_boundary_ops(boundary);
168 [e_v,d1_v,d2_v,d3_v,s_v,gamm_v,delt_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 168 [e_v,d1_v,d2_v,d3_v,s_v,gamm_v,delt_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
169 169