Mercurial > repos > public > sbplib

comparison +scheme/Scheme.m @ 1197:433c89bf19e0 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 07 Aug 2019 15:23:42 +0200
parents 5afc774fb7c4
children
comparison
equal deleted inserted replaced
1196:f6c571d8f22f 1197:433c89bf19e0
29 29
30 % type -- sets the type of interface, could be a string or a struct or something else 30 % type -- sets the type of interface, could be a string or a struct or something else
31 % depending on the particular scheme implementation 31 % depending on the particular scheme implementation
32 [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type) 32 [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
33 33
34 % TODO: op = getBoundaryOperator()?? 34 op = getBoundaryOperator(obj, opName, boundary)
35 % makes sense to have it available through a method instead of random properties 35 H_b= getBoundaryQuadrature(obj, boundary)
36 36
37 % Returns the number of degrees of freedom. 37 % Returns the number of degrees of freedom.
38 N = size(obj) 38 N = size(obj)
39 end 39 end
40
41 methods(Static)
42 % Calculates the matrcis need for the inteface coupling between
43 % boundary bound_u of scheme schm_u and bound_v of scheme schm_v.
44 % [uu, uv, vv, vu] = inteface_coupling(A,'r',B,'l')
45 function [uu, uv, vv, vu] = interface_coupling(schm_u,bound_u,schm_v,bound_v)
46 [uu,uv] = schm_u.interface(bound_u,schm_v,bound_v);
47 [vv,vu] = schm_v.interface(bound_v,schm_u,bound_u);
48 end
49 end
50 end 40 end