Mercurial > repos > public > sbplib

comparison +scheme/Laplace1d.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 0c504a21432d
comparison
equal deleted inserted replaced
1046:19ed046aec52 1047:6bc55a773e7c
86 end 86 end
87 87
88 function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) 88 function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type)
89 % u denotes the solution in the own domain 89 % u denotes the solution in the own domain
90 % v denotes the solution in the neighbour domain 90 % v denotes the solution in the neighbour domain
91 [e_u, d_u] = obj.getBoundaryOperator({'e', 'd'}, boundary); 91 e_u = obj.getBoundaryOperator('e', boundary);
92 d_u = obj.getBoundaryOperator('d', boundary);
92 s_u = obj.getBoundarySign(boundary); 93 s_u = obj.getBoundarySign(boundary);
93 94
94 [e_v, d_v] = neighbour_scheme.getBoundaryOperator({'e', 'd'}, neighbour_boundary); 95 e_v = neighbour_scheme.getBoundaryOperator('e', neighbour_boundary);
96 d_v = neighbour_scheme.getBoundaryOperator('d', neighbour_boundary);
95 s_v = neighbour_scheme.getBoundarySign(neighbour_boundary); 97 s_v = neighbour_scheme.getBoundarySign(neighbour_boundary);
96 98
97 a_u = obj.a; 99 a_u = obj.a;
98 a_v = neighbour_scheme.a; 100 a_v = neighbour_scheme.a;
99 101