Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 1047:6bc55a773e7c feature/getBoundaryOp

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix some mistakes from dc1bcbef2a86
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 22 Jan 2019 17:42:58 +0100
parents 19ed046aec52
children 84200bbae101
comparison
equal deleted inserted replaced
1046:19ed046aec52 1047:6bc55a773e7c
236 % neighbour_boundary is a string specifying which boundary to interface to. 236 % neighbour_boundary is a string specifying which boundary to interface to.
237 function [closure, penalty] = boundary_condition(obj, boundary, type, parameter) 237 function [closure, penalty] = boundary_condition(obj, boundary, type, parameter)
238 default_arg('type','neumann'); 238 default_arg('type','neumann');
239 default_arg('parameter', []); 239 default_arg('parameter', []);
240 240
241 [e, d] = obj.getBoundaryOperator({'e', 'd'}, boundary); 241 e = obj.getBoundaryOperator('e', boundary);
242 d = obj.getBoundaryOperator('d', boundary);
242 H_b = obj.getBoundaryQuadrature(boundary); 243 H_b = obj.getBoundaryQuadrature(boundary);
243 gamm = obj.getBoundaryBorrowing(boundary); 244 gamm = obj.getBoundaryBorrowing(boundary);
244 switch type 245 switch type
245 % Dirichlet boundary condition 246 % Dirichlet boundary condition
246 case {'D','d','dirichlet'} 247 case {'D','d','dirichlet'}
298 function [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type) 299 function [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type)
299 tuning = type.tuning; 300 tuning = type.tuning;
300 301
301 % u denotes the solution in the own domain 302 % u denotes the solution in the own domain
302 % v denotes the solution in the neighbour domain 303 % v denotes the solution in the neighbour domain
303 [e_u, d_u] = obj.getBoundaryOperator({'e', 'd'}, boundary); 304 e_u = obj.getBoundaryOperator('e', boundary);
305 d_u = obj.getBoundaryOperator('d', boundary);
304 H_b_u = obj.getBoundaryQuadrature(boundary); 306 H_b_u = obj.getBoundaryQuadrature(boundary);
305 I_u = obj.getBoundaryIndices(boundary); 307 I_u = obj.getBoundaryIndices(boundary);
306 gamm_u = obj.getBoundaryBorrowing(boundary); 308 gamm_u = obj.getBoundaryBorrowing(boundary);
307 309
308 [e_v, d_v] = neighbour_scheme.getBoundaryOperator({'e', 'd'}, neighbour_boundary); 310 e_v = neighbour_scheme.getBoundaryOperator('e', neighbour_boundary);
311 d_v = neighbour_scheme.getBoundaryOperator('d', neighbour_boundary);
309 H_b_v = neighbour_scheme.getBoundaryQuadrature(neighbour_boundary); 312 H_b_v = neighbour_scheme.getBoundaryQuadrature(neighbour_boundary);
310 I_v = neighbour_scheme.getBoundaryIndices(neighbour_boundary); 313 I_v = neighbour_scheme.getBoundaryIndices(neighbour_boundary);
311 gamm_v = neighbour_scheme.getBoundaryBorrowing(neighbour_boundary); 314 gamm_v = neighbour_scheme.getBoundaryBorrowing(neighbour_boundary);
312 315
313 u = obj; 316 u = obj;