Mercurial > repos > public > sbplib

comparison +multiblock/local2globalPenalty.m @ 968:a4ad90b37998 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 23 Dec 2018 14:39:31 +0100
parents ba10f24bf476
children
comparison
equal deleted inserted replaced
967:368a2773f78b 968:a4ad90b37998
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