Mercurial > repos > public > sbplib
view +multiblock/local2globalPenalty.m @ 1252:8fc2f9a4c882 feature/poroelastic
Add (temporary) fix for e-s , w-n, and x-x couplings
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Tue, 07 Jan 2020 12:49:09 -0800 |
parents | ba10f24bf476 |
children |
line wrap: on
line source
% 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