Mercurial > repos > public > sbplib

view diags.m @ 1269:5fa2f62a58a1 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add class for upwind D1 based on the standard norm matrices
author Martin Almquist <malmquist@stanford.edu>
date Sun, 31 May 2020 20:16:54 -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