Mercurial > repos > public > sbplib

view diags.m @ 1311:b2e84fe336d8 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Move common code for displacement and traction BC outside of switch statement
author Martin Almquist <malmquist@stanford.edu>
date Mon, 20 Jul 2020 22:21:57 -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