Mercurial > repos > public > sbplib

comparison +multiblock/local2globalPenalty.m @ 1018:eca4d9df9365 feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default for latest changes
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Dec 2018 17:31:28 +0100
parents ba10f24bf476
children
comparison
equal deleted inserted replaced
1017:2d7c1333bd6c 1018:eca4d9df9365
1 % Takes the block-local penalty and turns it into a global penalty
2 % local -- The local penalty
3 % div -- block matrix division for the diffOp
4 % I -- Index of blockmatrix block
5 function penalty = local2globalPenalty(local, div, I)
6 penaltyDiv = {div{1}, size(local,2)};
7 penalty_bm = blockmatrix.zero(penaltyDiv);
8 penalty_bm{I,1} = local;
9
10 penalty = blockmatrix.toMatrix(penalty_bm);
11 end