Mercurial > repos > public > sbplib

comparison +scheme/Divergence.m @ 1124:c2d281633e14 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Divergence and Gradient do not work for multi-block. Not sure how to fix, left error message for now.
author Martin Almquist <malmquist@stanford.edu>
date Sat, 11 May 2019 17:52:08 -0700
parents 8984b12feba6
children
comparison
equal deleted inserted replaced
1123:ff39d6692489 1124:c2d281633e14
96 % type Struct that specifies the interface coupling. 96 % type Struct that specifies the interface coupling.
97 % Fields: 97 % Fields:
98 % -- tuning: penalty strength, defaults to 1.2 98 % -- tuning: penalty strength, defaults to 1.2
99 % -- interpolation: type of interpolation, default 'none' 99 % -- interpolation: type of interpolation, default 'none'
100 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type) 100 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
101 101 error('Not implemented')
102 [m, n] = size(obj.D);
103 closure = sparse(m, n);
104
105 [m, n] = size(neighbour_scheme.D);
106 penalty = sparse(m, n);
107 end 102 end
108 103
109 function N = size(obj) 104 function N = size(obj)
110 N = obj.dim*prod(obj.m); 105 N = obj.dim*prod(obj.m);
111 end 106 end