Mercurial > repos > public > sbplib

view diags.m @ 926:60a3bc79835a feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make multiblock.setAllInterfaceOptions set opts only for those entries where grid.connections indicates interface couplings.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 03 Dec 2018 12:02:27 -0800
parents 0b87863cb29c
children
line wrap: on
line source

function A = diags(B,d,m,n)
    assert(size(B,1) == m);

    A = repmat(B(:,1)*0, [1, n]);

    for i = 1:size(B,2)
        A(:,d(i)+ (1:m)) = A(:,d(i)+ (1:m)) + diag(B(:,i));
    end
end