Mercurial > repos > public > sbplib

view diags.m @ 680:cd1a76c38565 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove schemes elasticDilation and elasticShear because they are subsets of elastic.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 05 Feb 2018 14:58:48 -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