Mercurial > repos > public > sbplib
view find_elements.m @ 1270:9a7cac4202d0 feature/poroelastic
In CruvilinearAnisotropicUpwind: Make D1Upwind default opSet. Use D2Standard for metric derivatives. Bugfix global quadrature.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Sun, 31 May 2020 20:18:45 -0700 |
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