Mercurial > repos > public > sbplib

view diags.m @ 503:feebfca90080 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added the new abstract classes to hypsyst3D and Upwind1D to make them work with Jonatans changes
author Ylva Rydin <ylva.rydin@telia.com>
date Tue, 16 May 2017 09:18:10 +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