Mercurial > repos > public > sbplib

view +sbp/diagInd.m @ 921:19c05eefc8c6 feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add warning that interfaceNonConforming in LaplaceCurve only works for Cartesian grids.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 02 Dec 2018 16:27:49 -0800
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