Mercurial > repos > public > sbplib

view +multiblock/local2globalClosure.m @ 1168:af3c4eb0cbbd feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add method for constructing diffops to multistage approach of computing the residual viscosity. Differentiate between constructDiffOpsMultiGrid and constructDiffOpsMultiStage
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 28 Jun 2019 13:12:29 +0200
parents ba10f24bf476
children
line wrap: on
line source

% Takes the block-local closures and turns it into a global closure
%   local -- The local closure
%   div   -- block matrix division for the diffOp
%   I     -- Index of blockmatrix block
function closure = local2globalClosure(local, div, I)
    closure_bm = blockmatrix.zero(div);
    closure_bm{I,I} = local;

    closure = blockmatrix.toMatrix(closure_bm);
end