Mercurial > repos > public > sbplib

comparison +multiblock/DiffOp.m @ 865:1cc5a0d26453 bcSetupExperiment

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use existing functions
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 25 Jul 2018 15:43:26 -0700
parents 517d6019d63c
children 57760d7088ad
comparison
equal deleted inserted replaced
864:c51b3f945c65 865:1cc5a0d26453
195 195
196 % Get the closure and penaly matrices 196 % Get the closure and penaly matrices
197 [blockClosure, blockPenalty] = obj.diffOps{I}.boundary_condition(name, type); 197 [blockClosure, blockPenalty] = obj.diffOps{I}.boundary_condition(name, type);
198 198
199 % Expand to matrix for full domain. 199 % Expand to matrix for full domain.
200 div = obj.blockmatrixDiv; 200 closure = multiblock.local2globalClosure(blockClosure, obj.blockmatrixDiv, I);
201 temp = blockmatrix.zero(div); 201 penalty = multiblock.local2globalPenalty(blockPenalty, obj.blockmatrixDiv, I);
202 temp{I,I} = blockClosure;
203 closure = blockmatrix.toMatrix(temp);
204
205 div{2} = size(blockPenalty, 2); % Penalty is a column vector
206 p = blockmatrix.zero(div);
207 p{I} = blockPenalty;
208 penalty = blockmatrix.toMatrix(p);
209 end 202 end
210 203
211 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 204 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
212 error('not implemented') 205 error('not implemented')
213 end 206 end