Mercurial > repos > public > sbplib

view +sbp/diagInd.m @ 1131:ea225a4659fe feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add 2d tests for derivative of delta function, and confirm that they work.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 22 May 2019 15:31:26 -0700
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