Mercurial > repos > public > sbplib

view diags.m @ 736:d991f8a0b61c feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add grid.blockEvalOn for functions with different functional expressions for different grid blocks
author Martin Almquist <malmquist@stanford.edu>
date Wed, 25 Apr 2018 14:56:17 -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