Mercurial > repos > public > sbplib

view diags.m @ 1083:a22b23098021 feature/laplace1dUpdate

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Divide interface penalty strength by 4, to make it sharp. Stability verified by eigenvalue computations.
author Martin Almquist <malmquist@stanford.edu>
date Tue, 05 Mar 2019 12:41:50 -0800
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