Mercurial > repos > public > sbplib

view +sbp/diagInd.m @ 527:763d50caddf3 feature/boundaryGroup

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Move the description of boundary identifiers to the grid package and update documentation for grid.Grid.getBoundaryNames
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 03 Aug 2017 14:04:36 +0200
parents 32b39dc44474
children
line wrap: on
line source

function A = diagInd(d,n,m)
    A = zeros(n,length(d));
    for i = 1:length(d)
        i0 = 1;
        j0 = d(i)+1;

        I = i0 + (0:(n-1))';
        J = j0 + (0:(n-1))';

        A(:,i) = matInd2VecInd(I,J,n);

    end
end

function I = matInd2VecInd(i,j,n)
    I = i + (j-1)*n;
end