comparison isEquidistant.m @ 691:527a7b954f26 feature/quantumTriangles

Merge with default
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 14 Sep 2017 16:00:36 +0200
parents b43c4d841afe
children
comparison
equal deleted inserted replaced
520:f235284e2eb1 691:527a7b954f26
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