Mercurial > repos > public > sbplib

view diags.m @ 536:6c6a3040a678 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add an abstrac class as help when defining transfinite multiblock grids
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 07 Aug 2017 09:57:03 +0200
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