Mercurial > repos > public > sbplib

view diags.m @ 1131:ea225a4659fe feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add 2d tests for derivative of delta function, and confirm that they work.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 22 May 2019 15:31:26 -0700
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