Mercurial > repos > public > sbplib

comparison +scheme/Scheme.m @ 1072:6468a5f6ec79 feature/grids/LaplaceSquared

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 12 Feb 2019 17:12:42 +0100
parents 5afc774fb7c4
children
comparison
equal deleted inserted replaced
1071:92cb03e64ca4 1072:6468a5f6ec79
24 % interfaced to. 24 % interfaced to.
25 % neighbour_boundary is a string specifying which boundary to 25 % neighbour_boundary is a string specifying which boundary to
26 % interface to. 26 % interface to.
27 % penalty may be a cell array if there are several penalties with different weights 27 % penalty may be a cell array if there are several penalties with different weights
28 [closure, penalty] = boundary_condition(obj,boundary,type) % TODO: Change name to boundaryCondition 28 [closure, penalty] = boundary_condition(obj,boundary,type) % TODO: Change name to boundaryCondition
29 [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
30 29
31 % TODO: op = getBoundaryOperator()?? 30 % type -- sets the type of interface, could be a string or a struct or something else
32 % makes sense to have it available through a method instead of random properties 31 % depending on the particular scheme implementation
32 [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
33
34 op = getBoundaryOperator(obj, opName, boundary)
35 H_b= getBoundaryQuadrature(obj, boundary)
33 36
34 % Returns the number of degrees of freedom. 37 % Returns the number of degrees of freedom.
35 N = size(obj) 38 N = size(obj)
36 end 39 end
37
38 methods(Static)
39 % Calculates the matrcis need for the inteface coupling between
40 % boundary bound_u of scheme schm_u and bound_v of scheme schm_v.
41 % [uu, uv, vv, vu] = inteface_coupling(A,'r',B,'l')
42 function [uu, uv, vv, vu] = interface_coupling(schm_u,bound_u,schm_v,bound_v)
43 [uu,uv] = schm_u.interface(bound_u,schm_v,bound_v);
44 [vv,vu] = schm_v.interface(bound_v,schm_u,bound_u);
45 end
46 end
47 end 40 end