Mercurial > repos > public > sbplib
view +multiblock/local2globalPenalty.m @ 1286:4cb627c7fb90 feature/boundary_optimized_grids
Make D1Nonequidistant use the grid generation functions accurate/minimalBoundaryOptimizedGrid and remove grid generation from +implementations
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Wed, 01 Jul 2020 13:43:32 +0200 |
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