Mercurial > repos > public > sbplib

view diags.m @ 1166:17bb3c46472c feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add function for constructing differential operator on coarser mesh. Not sure if it will be usefull tho.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 27 Jun 2019 17:15:11 +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