Mercurial > repos > public > sbplib

view diags.m @ 711:397d1b22cc37 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Changed closure and penalty ao that sparse matrices are created in DiffOpTimeDep
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 30 Nov 2017 18:01:21 +0100
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