Mercurial > repos > public > sbplib

view find_elements.m @ 1065:c2bd7f15da48 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add scheme where only corner points are multiplied by dim. Does not seem to decrease spectral radius much, but is stable.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 17 Jan 2019 18:57:54 -0800
parents 48b6fb693025
children
line wrap: on
line source

% I = find_elements(a,b)
% Finds the index of elements a in b.
% a and b have to be in the same order.
function I = find_elements(a,b)
    I = [];

    j = 1;
    for i = 1:length(a)
        while b(j) ~= a(i)
            j = j + 1;
        end
        I(end+1) = j;
        j = j+1;
    end

    assert(length(I) == length(a),'Expected %d but got %d elements',length(a),length(I))
end