Mercurial > repos > public > sbplib
comparison isEquidistant.m @ 538:95e41f7a6f1a feature/grids
Merge default
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 08 Aug 2017 08:26:49 +0200 |
| parents | b43c4d841afe |
| children |
comparison
equal
deleted
inserted
replaced
| 537:a70d5387d2ca | 538:95e41f7a6f1a |
|---|---|
| 1 % Tests if consecutive elements of vector v are euidistant | |
| 2 function b = isEquidistant(v) | |
| 3 if length(v) < 2 | |
| 4 error('sbplib:isEquidistant:inputTooShort', 'Input vector is too short'); | |
| 5 end | |
| 6 | |
| 7 tol = 1e-8; | |
| 8 | |
| 9 d = v(2:end) - v(1:end-1); | |
| 10 err = abs(d - d(1)); | |
| 11 | |
| 12 relErr = err./abs(d); | |
| 13 | |
| 14 I_zero = find(d < tol); | |
| 15 | |
| 16 relErr(I_zero) = err(I_zero); | |
| 17 | |
| 18 b = all(relErr < tol); | |
| 19 end |