Mercurial > repos > public > sbplib

view +sbp/diagInd.m @ 1275:8ff3a95ad7cc feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change metric computation in CurvilinearUpwind so that one can use odd orders for the diffOp
author Martin Almquist <malmquist@stanford.edu>
date Wed, 03 Jun 2020 10:29:53 -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