annotate find_elements.m @ 835:008496ca38f3
feature/burgers1d
Compute the residual in between each runge-kutta stage.
Note: It is not clear whether the correct residual is used when computing the stages. Must investigate further.
| author |
Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
| date |
Thu, 13 Sep 2018 18:14:54 +0200 |
| parents |
48b6fb693025 |
| children |
|
| rev |
line source |
|
0
|
1 % I = find_elements(a,b)
|
|
|
2 % Finds the index of elements a in b.
|
|
|
3 % a and b have to be in the same order.
|
|
|
4 function I = find_elements(a,b)
|
|
|
5 I = [];
|
|
|
6
|
|
|
7 j = 1;
|
|
|
8 for i = 1:length(a)
|
|
|
9 while b(j) ~= a(i)
|
|
|
10 j = j + 1;
|
|
|
11 end
|
|
|
12 I(end+1) = j;
|
|
|
13 j = j+1;
|
|
|
14 end
|
|
|
15
|
|
|
16 assert(length(I) == length(a),'Expected %d but got %d elements',length(a),length(I))
|
|
|
17 end |
