Mercurial > repos > public > sbplib
diff +multiblock/local2globalPenalty.m @ 902:306f5b3cd7bc
Merged in bcSetupExperiment (pull request #10)
Improved bulk boundary condition handling
Approved-by: Martin Almquist <malmquist@stanford.edu>
| author | Jonatan Werpers <jonatan.werpers@it.uu.se> |
|---|---|
| date | Sat, 24 Nov 2018 14:55:48 +0000 |
| parents | ba10f24bf476 |
| children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+multiblock/local2globalPenalty.m Sat Nov 24 14:55:48 2018 +0000 @@ -0,0 +1,11 @@ +% Takes the block-local penalty and turns it into a global penalty +% local -- The local penalty +% div -- block matrix division for the diffOp +% I -- Index of blockmatrix block +function penalty = local2globalPenalty(local, div, I) + penaltyDiv = {div{1}, size(local,2)}; + penalty_bm = blockmatrix.zero(penaltyDiv); + penalty_bm{I,1} = local; + + penalty = blockmatrix.toMatrix(penalty_bm); +end