Mercurial > repos > public > sbplib

comparison +scheme/Wave2dCurve.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 706d1c2b4199
children
comparison
equal deleted inserted replaced
967:368a2773f78b 968:a4ad90b37998
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);