DGtal 1.4.0
Loading...
Searching...
No Matches
testPointVector.cpp
Go to the documentation of this file.
1
32#include <cstdio>
33#include <cmath>
34#include <iostream>
35#include <fstream>
36#include <vector>
37#include <type_traits>
38#include <functional>
39
40#include "DGtal/base/Common.h"
41#include "DGtal/kernel/PointVector.h"
42
43#include "DGtalCatch.h"
44
45using namespace DGtal;
46using namespace std;
47
48// Compare the value of the two parameters and also their component type.
49#define COMPARE_VALUE_AND_TYPE(expr, check) \
50 REQUIRE( (expr) == (check) ); \
51 REQUIRE( ( std::is_same<decltype(expr)::Component, decltype(check)::Component>::value ) );
52
53TEST_CASE( "2D Point Vector Unit tests" )
54{
55 using Real = double;
56 using Integer = DGtal::int32_t;
57
58 typedef PointVector<2, Integer> Point2D;
59 typedef PointVector<2, Real> RealPoint2D;
60 typedef PointVector<3, Integer> Point3D;
61 typedef PointVector<3, Real> RealPoint3D;
62
63 Integer t1[] = {1,2};
64 Integer t2[] = {5,4};
65 Real t3[] = {1.5,2.5};
66 Real t4[] = {5.5,4.5};
67
68 Point2D p1( t1 );
69 Point2D p2( t2 );
70 RealPoint2D p3(t3);
71 RealPoint2D p4(t4);
72
73 Point3D p1_3d( p1[0], p1[1] );
74 Point3D p2_3d( p2[0], p2[1] );
75 RealPoint3D p3_3d( p3[0], p3[1] );
76 RealPoint3D p4_3d( p4[0], p4[1] );
77
78 SECTION("Cross products with integers")
79 {
80 COMPARE_VALUE_AND_TYPE( p1.crossProduct(p2), p1_3d.crossProduct(p2_3d) );
81 COMPARE_VALUE_AND_TYPE( crossProduct(p1, p2), crossProduct(p1_3d, p2_3d) );
82 COMPARE_VALUE_AND_TYPE( p2.crossProduct(p1), p2_3d.crossProduct(p1_3d) );
83 COMPARE_VALUE_AND_TYPE( crossProduct(p2, p1), crossProduct(p2_3d, p1_3d) );
84 }
85
86 SECTION("Cross products with reals")
87 {
88 COMPARE_VALUE_AND_TYPE( p3.crossProduct(p4), p3_3d.crossProduct(p4_3d) );
89 COMPARE_VALUE_AND_TYPE( crossProduct(p3, p4), crossProduct(p3_3d, p4_3d) );
90 COMPARE_VALUE_AND_TYPE( p4.crossProduct(p3), p4_3d.crossProduct(p3_3d) );
91 COMPARE_VALUE_AND_TYPE( crossProduct(p4, p3), crossProduct(p4_3d, p3_3d) );
92 }
93
94 SECTION("Cross products with mixed integers/reals")
95 {
96 COMPARE_VALUE_AND_TYPE( p1.crossProduct(p3), p1_3d.crossProduct(p3_3d) );
97 COMPARE_VALUE_AND_TYPE( crossProduct(p1, p3), crossProduct(p1_3d, p3_3d) );
98 COMPARE_VALUE_AND_TYPE( p3.crossProduct(p1), p3_3d.crossProduct(p1_3d) );
99 COMPARE_VALUE_AND_TYPE( crossProduct(p3, p1), crossProduct(p3_3d, p1_3d) );
100 }
101 SECTION("Access data() of internal container")
102 {
103 const auto d = p1_3d.data();
104 CHECK(d[0] == p1[0]);
105 CHECK(d[1] == p1[1]);
106 }
107}
108
109TEST_CASE( "3D Point Vector Unit tests" )
110{
111 using Real = double;
112 using Integer = DGtal::int32_t;
113
114 typedef PointVector<3, Integer> Point;
115 typedef PointVector<3, Real> RealPoint;
116
117 Integer t1[] = {1,2,3};
118 Integer t2[] = {5,4,3};
119 Real t3[] = {1.5,2.5,3.5};
120 Real t4[] = {5.5,4.5,3.5};
121
122 Point p1( t1 );
123 Point p2( t2 );
124 RealPoint p3(t3);
125 RealPoint p4(t4);
126
127 SECTION("Cross products with integers")
128 {
129 COMPARE_VALUE_AND_TYPE( p1.crossProduct(p2), Point(-6, 12, -6) );
130 COMPARE_VALUE_AND_TYPE( crossProduct(p1, p2), Point(-6, 12, -6) );
131 COMPARE_VALUE_AND_TYPE( p2.crossProduct(p1), Point(6, -12, 6) );
132 COMPARE_VALUE_AND_TYPE( crossProduct(p2, p1), Point(6, -12, 6) );
133 }
134
135 SECTION("Cross products with reals")
136 {
137 COMPARE_VALUE_AND_TYPE( p3.crossProduct(p4), RealPoint(-7., 14., -7.) );
138 COMPARE_VALUE_AND_TYPE( crossProduct(p3, p4), RealPoint(-7., 14., -7.) );
139 COMPARE_VALUE_AND_TYPE( p4.crossProduct(p3), RealPoint(7., -14., 7.) );
140 COMPARE_VALUE_AND_TYPE( crossProduct(p4, p3), RealPoint(7., -14., 7.) );
141 }
142
143 SECTION("Cross products with mixed integers/reals")
144 {
145 COMPARE_VALUE_AND_TYPE( p1.crossProduct(p3), RealPoint(-0.5, 1., -0.5) );
146 COMPARE_VALUE_AND_TYPE( crossProduct(p1, p3), RealPoint(-0.5, 1., -0.5) );
147 COMPARE_VALUE_AND_TYPE( p3.crossProduct(p1), RealPoint(0.5, -1., 0.5) );
148 COMPARE_VALUE_AND_TYPE( crossProduct(p3, p1), RealPoint(0.5, -1., 0.5) );
149 }
150}
151
152TEST_CASE( "4D Point Vector Unit tests" )
153{
154 using Real = double;
155 using Integer = DGtal::int32_t;
156
157 typedef PointVector<4, Integer> Point;
158 typedef PointVector<4, Real> RealPoint;
159
160 const Real pi = std::acos(Real(-1));
161
162 Integer t1[] = {1,2,3,4};
163 Integer t2[] = {5,4,3,2};
164 Real t3[] = {1.0,-1.0,2.0,-2.0};
165 Real t4[] = {5.5,-4.5,3.5,2.5};
166
167 Point p1( t1 );
168 Point p1bis( t1 );
169 Point p2( t2 );
170 RealPoint p3(t3);
171 RealPoint p4(t4);
172 RealPoint p1r(p1);
173 RealPoint p2r(p2);
174
175 SECTION("Construction")
176 {
177 REQUIRE( p1 == Point( {1, 2, 3, 4} ) );
178 REQUIRE( p1r == RealPoint( {1., 2., 3., 4.} ) );
179
180 REQUIRE( p1 == Point( p1r ) );
181 REQUIRE( p1 == Point( RealPoint( {0.5, 2.1, 3.1, 3.9} ), DGtal::functors::Round<>() ) );
182 REQUIRE( p1 == Point( Point(0, 0, 1, 3), Point(1, 2, 2, 1), std::plus<Integer>() ) );
183 }
184
185 SECTION("Assignments")
186 {
187 Point dummy1;
188 dummy1 = p1;
189 REQUIRE( p1 == dummy1 );
190
191 Point dummy2(1, 3, 3, 5);
192 dummy2.partialCopy( Point(0, 2, 0, 4), {1, 3} );
193 REQUIRE( p1 == dummy2 );
194
195 Point dummy3(2, 2, 1, 4);
196 dummy3.partialCopyInv( Point(1, 0, 3, 0), {1, 3} );
197 REQUIRE( p1 == dummy3 );
198
199 RealPoint dummy1r;
200 dummy1r = p1;
201 REQUIRE( p1r == dummy1r );
202
203 RealPoint dummy2r(1, 3, 3, 5);
204 dummy2r.partialCopy( Point(0, 2, 0, 4), {1, 3} );
205 REQUIRE( p1r == dummy2r );
206
207 RealPoint dummy3r(2, 2, 1, 4);
208 dummy3r.partialCopyInv( Point(1, 0, 3, 0), {1, 3} );
209 REQUIRE( p1r == dummy3r );
210
211 Point dummy4(1, 3, 3, 5);
212 dummy4.partialCopy( RealPoint(0, 1.5, 0, 4.1), {1, 3}, DGtal::functors::Round<>() );
213 REQUIRE( p1 == dummy4 );
214
215 Point dummy5(2, 2, 1, 4);
216 dummy5.partialCopyInv( RealPoint(1.1, 0, 2.5, 0), {1, 3}, DGtal::functors::Round<>() );
217 REQUIRE( p1 == dummy5 );
218 }
219
220 SECTION("Comparisons")
221 {
222 // Partial equality
223 REQUIRE( p1.partialEqual( RealPoint(0, 2, 0, 4), {1, 3} ) );
224 REQUIRE( ! p1.partialEqual( RealPoint(0, 1, 0, 4), {1, 3} ) );
225 REQUIRE( p1.partialEqualInv( RealPoint(1, 0, 3, 0), {1, 3} ) );
226 REQUIRE( ! p1.partialEqualInv( RealPoint(1, 0, 2, 0), {1, 3} ) );
227
228 // Two equal points of same type
229 REQUIRE( p1 == p1bis );
230 REQUIRE( p1bis == p1 );
231 REQUIRE( ! (p1 != p1bis) );
232 REQUIRE( ! (p1bis != p1) );
233 REQUIRE( p1 <= p1bis );
234 REQUIRE( p1 >= p1bis );
235 REQUIRE( p1bis >= p1 );
236 REQUIRE( p1bis <= p1 );
237 REQUIRE( ! (p1 < p1bis) );
238 REQUIRE( ! (p1 > p1bis) );
239 REQUIRE( ! (p1bis > p1) );
240 REQUIRE( ! (p1bis < p1) );
241
242 // Two equal points of different types
243 REQUIRE( p1 == p1r );
244 REQUIRE( p1r == p1 );
245 REQUIRE( ! (p1 != p1r) );
246 REQUIRE( ! (p1r != p1) );
247 REQUIRE( p1 <= p1r );
248 REQUIRE( p1 >= p1r );
249 REQUIRE( p1r >= p1 );
250 REQUIRE( p1r <= p1 );
251 REQUIRE( ! (p1 < p1r) );
252 REQUIRE( ! (p1 > p1r) );
253 REQUIRE( ! (p1r > p1) );
254 REQUIRE( ! (p1r < p1) );
255
256 // Two ordered points of same type
257 REQUIRE( ! (p1 == p2) );
258 REQUIRE( ! (p2 == p1) );
259 REQUIRE( p1 != p2 );
260 REQUIRE( p2 != p1 );
261 REQUIRE( p1 <= p2 );
262 REQUIRE( ! (p1 >= p2) );
263 REQUIRE( p2 >= p1 );
264 REQUIRE( ! (p2 <= p1) );
265 REQUIRE( p1 < p2 );
266 REQUIRE( ! (p1 > p2) );
267 REQUIRE( p2 > p1 );
268 REQUIRE( ! (p2 < p1) );
269
270 // Two ordered points of different types
271 REQUIRE( ! (p1 == p2r) );
272 REQUIRE( ! (p2r == p1) );
273 REQUIRE( p1 != p2r );
274 REQUIRE( p2r != p1 );
275 REQUIRE( p1 <= p2r );
276 REQUIRE( ! (p1 >= p2r) );
277 REQUIRE( p2r >= p1 );
278 REQUIRE( ! (p2r <= p1) );
279 REQUIRE( p1 < p2r );
280 REQUIRE( ! (p1 > p2r) );
281 REQUIRE( p2r > p1 );
282 REQUIRE( ! (p2r < p1) );
283 }
284
285 SECTION("Min/Max of vector components")
286 {
287 REQUIRE( p3.max() == 2.0 );
288 REQUIRE( p3.min() == -2.0 );
289 REQUIRE( *p3.maxElement() == 2.0 );
290 REQUIRE( *p3.minElement() == -2.0 );
291 }
292
293 Point aPoint;
294 aPoint[ 3 ] = 0;
295 aPoint[ 2 ] = 2;
296 aPoint[ 1 ] = -1;
297 aPoint[ 0 ] = 3;
298
299 SECTION("Testing norms")
300 {
301 RealPoint normalized = aPoint.getNormalized();
302 CAPTURE( normalized );
303 REQUIRE( aPoint.norm ( Point::L_1 ) == 6 );
304 REQUIRE( aPoint.norm ( Point::L_infty ) == 3 );
305 REQUIRE( aPoint.squaredNorm() == Approx(aPoint.norm()*aPoint.norm()) );
306 REQUIRE( normalized[0] == Approx( 0.801784) );
307 REQUIRE( normalized[1] == Approx( -0.267261) );
308 REQUIRE( normalized[2] == Approx( 0.534522) );
309 REQUIRE( normalized[3] == Approx( 0.0) );
310 }
311
312 SECTION("PointVector Iterator")
313 {
314 PointVector<25,int> aPoint25;
315 for (unsigned int i=0;i<25;++i)
316 aPoint25[i] = i;
317
318 int sum = 0;
319 for (PointVector<25,int>::ConstIterator it = aPoint25.begin() ; it != aPoint25.end(); ++it)
320 sum += (*it);
321
322 CAPTURE(aPoint25);
323 CAPTURE(sum);
324 REQUIRE( sum == 300 );
325 }
326
327 SECTION("Arithmetical operators with integers")
328 {
329 COMPARE_VALUE_AND_TYPE( p1 + p2, Point(6,6,6,6) );
330 COMPARE_VALUE_AND_TYPE( p1 - p2, Point(-4,-2,0,2) );
331 COMPARE_VALUE_AND_TYPE( p1 * p2, Point(5,8,9,8) );
332 COMPARE_VALUE_AND_TYPE( p2 / p1, Point(5,2,1,0) );
333
334 COMPARE_VALUE_AND_TYPE( p1 + 2, Point(3,4,5,6) );
335 COMPARE_VALUE_AND_TYPE( 2 + p1, Point(3,4,5,6) );
336 COMPARE_VALUE_AND_TYPE( p1 - 2, Point(-1,0,1,2) );
337 COMPARE_VALUE_AND_TYPE( 2 - p1, Point(1,0,-1,-2) );
338 COMPARE_VALUE_AND_TYPE( p1 * 2, Point(2,4,6,8) );
339 COMPARE_VALUE_AND_TYPE( 2 * p1, Point(2,4,6,8) );
340 COMPARE_VALUE_AND_TYPE( p1 / 2, Point(0,1,1,2) );
341 COMPARE_VALUE_AND_TYPE( 2 / p1, Point(2,1,0,0) );
342
343 COMPARE_VALUE_AND_TYPE( -p1, Point(-1,-2,-3,-4) );
344
345 p1 *= 2; COMPARE_VALUE_AND_TYPE( p1, Point(2,4,6,8) );
346 p1 += 2; COMPARE_VALUE_AND_TYPE( p1, Point(4,6,8,10) );
347 p1 -= 2; COMPARE_VALUE_AND_TYPE( p1, Point(2,4,6,8) );
348 p1 /= 2; COMPARE_VALUE_AND_TYPE( p1, Point(1,2,3,4) );
349
350 p1 *= p2; COMPARE_VALUE_AND_TYPE( p1, Point(5,8,9,8) );
351 p1 += p2; COMPARE_VALUE_AND_TYPE( p1, Point(10,12,12,10) );
352 p1 -= p2; COMPARE_VALUE_AND_TYPE( p1, Point(5,8,9,8) );
353 p1 /= p2; COMPARE_VALUE_AND_TYPE( p1, Point(1,2,3,4) );
354 }
355
356 SECTION("Other operators with integers")
357 {
358 COMPARE_VALUE_AND_TYPE( p1.inf(p2), Point(1,2,3,2) );
359 COMPARE_VALUE_AND_TYPE( p1.sup(p2), Point(5,4,3,4) );
360 COMPARE_VALUE_AND_TYPE( inf(p1, p2), Point(1,2,3,2) );
361 COMPARE_VALUE_AND_TYPE( sup(p1, p2), Point(5,4,3,4) );
362
363 REQUIRE( p1.dot(p2) == 30 );
364 REQUIRE( dotProduct(p1, p2) == 30 );
365
366 REQUIRE( p1.cosineSimilarity(p1) == Approx(0.).margin(0.000001));
367 REQUIRE( p1.cosineSimilarity(-p1) == Approx(pi).margin(0.000001));
368 REQUIRE( p1.cosineSimilarity( Point(-2,1,-4,3) ) == Approx(pi/2).margin(0.000001) );
369 REQUIRE( cosineSimilarity(p1, p1) == Approx(0.).margin(0.000001) );
370 REQUIRE( cosineSimilarity(p1, -p1) == Approx(pi).margin(0.000001) );
371 REQUIRE( cosineSimilarity(p1, Point(-2,1,-4,3)) == Approx(pi/2).margin(0.000001) );
372
373 REQUIRE( p1.isLower(p2) == false );
374 REQUIRE( isLower(p1, p2) == false );
375 REQUIRE( p2.isUpper(p1) == false );
376 REQUIRE( isUpper(p2, p1) == false );
377 p1[3] = 2;
378 REQUIRE( p1.isLower(p2) == true );
379 REQUIRE( isLower(p1, p2) == true );
380 REQUIRE( p2.isUpper(p1) == true );
381 REQUIRE( isUpper(p2, p1) == true );
382 }
383
384 SECTION("Arithmetical Operators with reals")
385 {
386 COMPARE_VALUE_AND_TYPE( p3 + p4, RealPoint(6.5,-5.5,5.5,0.5) );
387 COMPARE_VALUE_AND_TYPE( p3 - p4, RealPoint(-4.5,3.5,-1.5,-4.5) );
388 COMPARE_VALUE_AND_TYPE( p3 * p4, RealPoint(5.5,4.5,7.0,-5.0) );
389 COMPARE_VALUE_AND_TYPE( p4 / p3, RealPoint(5.5,4.5,1.75,-1.25) );
390
391 COMPARE_VALUE_AND_TYPE( p3 + 2., RealPoint(3.,1.,4.,0.) );
392 COMPARE_VALUE_AND_TYPE( 2. + p3, RealPoint(3.,1.,4.,0.) );
393 COMPARE_VALUE_AND_TYPE( p3 - 2., RealPoint(-1.,-3.,0.,-4.) );
394 COMPARE_VALUE_AND_TYPE( 2. - p3, RealPoint(1.,3.,0.,4.) );
395 COMPARE_VALUE_AND_TYPE( p3 * 2., RealPoint(2.,-2.,4.,-4.) );
396 COMPARE_VALUE_AND_TYPE( 2. * p3, RealPoint(2.,-2.,4.,-4.) );
397 COMPARE_VALUE_AND_TYPE( p3 / 2., RealPoint(0.5,-0.5,1.,-1.) );
398 COMPARE_VALUE_AND_TYPE( 2. / p3, RealPoint(2.,-2.,1.,-1.) );
399
400 COMPARE_VALUE_AND_TYPE( -p3, RealPoint(-1.,1.,-2.,2.) );
401
402 p3 *= 2.5; COMPARE_VALUE_AND_TYPE( p3, RealPoint(2.5, -2.5, 5., -5.) );
403 p3 += 2.5; COMPARE_VALUE_AND_TYPE( p3, RealPoint(5., 0., 7.5, -2.5) );
404 p3 -= 2.5; COMPARE_VALUE_AND_TYPE( p3, RealPoint(2.5, -2.5, 5., -5.) );
405 p3 /= 2.5; COMPARE_VALUE_AND_TYPE( p3, RealPoint(1., -1., 2., -2.) );
406
407 p3 *= p4; COMPARE_VALUE_AND_TYPE( p3, RealPoint(5.5, 4.5, 7., -5.) );
408 p3 += p4; COMPARE_VALUE_AND_TYPE( p3, RealPoint(11, 0., 10.5, -2.5) );
409 p3 -= p4; COMPARE_VALUE_AND_TYPE( p3, RealPoint(5.5, 4.5, 7., -5.) );
410 p3 /= p4; COMPARE_VALUE_AND_TYPE( p3, RealPoint(1., -1., 2., -2.) );
411 }
412
413 SECTION("Other operators with reals")
414 {
415 COMPARE_VALUE_AND_TYPE( p3.inf(p4), RealPoint(1.,-4.5,2.,-2.) );
416 COMPARE_VALUE_AND_TYPE( p3.sup(p4), RealPoint(5.5,-1.,3.5,2.5) );
417 COMPARE_VALUE_AND_TYPE( inf(p3, p4), RealPoint(1.,-4.5,2.,-2.) );
418 COMPARE_VALUE_AND_TYPE( sup(p3, p4), RealPoint(5.5,-1.,3.5,2.5) );
419
420 REQUIRE( p3.dot(p4) == 12. );
421 REQUIRE( dotProduct(p3, p4) == 12. );
422
423 REQUIRE( p3.cosineSimilarity(p3) == Approx(0.).margin(0.000001) );
424 REQUIRE( p3.cosineSimilarity(-p3) == Approx(pi).margin(0.000001) );
425 REQUIRE( p3.cosineSimilarity( RealPoint(1.0,1.0,2.0,2.0) ) == Approx(pi/2).margin(0.000001) );
426 REQUIRE( cosineSimilarity(p3, p3) == Approx(0.).margin(0.000001) );
427 REQUIRE( cosineSimilarity(p3, -p3) == Approx(pi).margin(0.000001) );
428 REQUIRE( cosineSimilarity(p3, RealPoint(1.0,1.0,2.0,2.0)) == Approx(pi/2).margin(0.000001) );
429
430 REQUIRE( p3.isLower(p4) == false );
431 REQUIRE( isLower(p3, p4) == false );
432 REQUIRE( p4.isUpper(p3) == false );
433 REQUIRE( isUpper(p4, p3) == false );
434 p4[1] = -p4[1];
435 REQUIRE( p3.isLower(p4) == true );
436 REQUIRE( isLower(p3, p4) == true );
437 REQUIRE( p4.isUpper(p3) == true );
438 REQUIRE( isUpper(p4, p3) == true );
439 }
440
441 SECTION("Arithmetical Operators with mixed integers/reals")
442 {
443 COMPARE_VALUE_AND_TYPE( p1 + p4, RealPoint(6.5,-2.5,6.5,6.5) );
444 COMPARE_VALUE_AND_TYPE( p4 + p1, RealPoint(6.5,-2.5,6.5,6.5) );
445 COMPARE_VALUE_AND_TYPE( p1 - p4, RealPoint(-4.5,6.5,-0.5,1.5) );
446 COMPARE_VALUE_AND_TYPE( p4 - p1, RealPoint(4.5,-6.5,0.5,-1.5) );
447 COMPARE_VALUE_AND_TYPE( p1 * p4, RealPoint(5.5,-9.0,10.5,10.0) );
448 COMPARE_VALUE_AND_TYPE( p4 * p1, RealPoint(5.5,-9.0,10.5,10.0) );
449 COMPARE_VALUE_AND_TYPE( p1 / p3, RealPoint(1.,-2.,1.5,-2.0) );
450 COMPARE_VALUE_AND_TYPE( p3 / p1, RealPoint(1.,-0.5,2./3.,-0.5) );
451
452 COMPARE_VALUE_AND_TYPE( p3 + 2, RealPoint(3.,1.,4.,0.) );
453 COMPARE_VALUE_AND_TYPE( 2 + p3, RealPoint(3.,1.,4.,0.) );
454 COMPARE_VALUE_AND_TYPE( p3 - 2, RealPoint(-1.,-3.,0.,-4.) );
455 COMPARE_VALUE_AND_TYPE( 2 - p3, RealPoint(1.,3.,0.,4.) );
456 COMPARE_VALUE_AND_TYPE( p3 * 2, RealPoint(2.,-2.,4.,-4.) );
457 COMPARE_VALUE_AND_TYPE( 2 * p3, RealPoint(2.,-2.,4.,-4.) );
458 COMPARE_VALUE_AND_TYPE( p3 / 2, RealPoint(0.5,-0.5,1.,-1.) );
459 COMPARE_VALUE_AND_TYPE( 2 / p3, RealPoint(2.,-2.,1.,-1.) );
460
461 COMPARE_VALUE_AND_TYPE( p1 + 2.5, RealPoint(3.5,4.5,5.5,6.5) );
462 COMPARE_VALUE_AND_TYPE( 2.5 + p1, RealPoint(3.5,4.5,5.5,6.5) );
463 COMPARE_VALUE_AND_TYPE( p1 - 2.5, RealPoint(-1.5,-0.5,0.5,1.5) );
464 COMPARE_VALUE_AND_TYPE( 2.5 - p1, RealPoint(1.5,0.5,-0.5,-1.5) );
465 COMPARE_VALUE_AND_TYPE( p1 * 2.5, RealPoint(2.5,5.,7.5,10.) );
466 COMPARE_VALUE_AND_TYPE( 2.5 * p1, RealPoint(2.5,5.,7.5,10.) );
467 COMPARE_VALUE_AND_TYPE( p1 / 0.5, RealPoint(2.,4.,6.,8.) );
468 COMPARE_VALUE_AND_TYPE( 2. / p1, RealPoint(2.,1.,2./3.,0.5) );
469
470 p3 *= 2; COMPARE_VALUE_AND_TYPE( p3, RealPoint(2, -2, 4., -4.) );
471 p3 += 2; COMPARE_VALUE_AND_TYPE( p3, RealPoint(4., 0., 6., -2.) );
472 p3 -= 2; COMPARE_VALUE_AND_TYPE( p3, RealPoint(2, -2, 4., -4.) );
473 p3 /= 2; COMPARE_VALUE_AND_TYPE( p3, RealPoint(1., -1., 2., -2.) );
474
475 p3 *= p1; COMPARE_VALUE_AND_TYPE( p3, RealPoint(1., -2., 6., -8.) );
476 p3 += p1; COMPARE_VALUE_AND_TYPE( p3, RealPoint(2., 0., 9., -4.) );
477 p3 -= p1; COMPARE_VALUE_AND_TYPE( p3, RealPoint(1., -2., 6., -8.) );
478 p3 /= p1; COMPARE_VALUE_AND_TYPE( p3, RealPoint(1., -1., 2., -2.) );
479 }
480
481 SECTION("Other operators with mixed integers/reals")
482 {
483 COMPARE_VALUE_AND_TYPE( p1.inf(p3), RealPoint(1.,-1.,2.,-2.) );
484 COMPARE_VALUE_AND_TYPE( p3.inf(p1), RealPoint(1.,-1.,2.,-2.) );
485 COMPARE_VALUE_AND_TYPE( p1.sup(p3), RealPoint(1.,2.,3.,4.) );
486 COMPARE_VALUE_AND_TYPE( p3.sup(p1), RealPoint(1.,2.,3.,4.) );
487 COMPARE_VALUE_AND_TYPE( inf(p1, p3), RealPoint(1.,-1.,2.,-2.) );
488 COMPARE_VALUE_AND_TYPE( inf(p3, p1), RealPoint(1.,-1.,2.,-2.) );
489 COMPARE_VALUE_AND_TYPE( sup(p1, p3), RealPoint(1.,2.,3.,4.) );
490 COMPARE_VALUE_AND_TYPE( sup(p3, p1), RealPoint(1.,2.,3.,4.) );
491
492 REQUIRE( p4.dot(p1) == 17.0 );
493 REQUIRE( dotProduct(p4, p1) == 17.0 );
494 REQUIRE( dotProduct(p1, p4) == 17.0 );
495
496 REQUIRE( p1.cosineSimilarity(RealPoint(p1)) == Approx(0.).margin(0.000001) );
497 REQUIRE( p1.cosineSimilarity(-RealPoint(p1)) == Approx(pi).margin(0.000001) );
498 REQUIRE( p1.cosineSimilarity( RealPoint(-2,1,-4,3) ) == Approx(pi/2).margin(0.000001) );
499 REQUIRE( cosineSimilarity(p1, RealPoint(p1)) == Approx(0.).margin(0.000001) );
500 REQUIRE( cosineSimilarity(p1, -RealPoint(p1)) == Approx(pi).margin(0.000001) );
501 REQUIRE( cosineSimilarity(p1, RealPoint(-2,1,-4,3)) == Approx(pi/2).margin(0.000001) );
502
503 REQUIRE( p3.cosineSimilarity(Point(1,-1,2,-2)) == Approx(0.).margin(0.000001) );
504 REQUIRE( p3.cosineSimilarity(-Point(1,-1,2,-2)) == Approx(pi).margin(0.000001) );
505 REQUIRE( p3.cosineSimilarity( Point(1,1,2,2) ) == Approx(pi/2).margin(0.000001) );
506 REQUIRE( cosineSimilarity(p3, Point(1,-1,2,-2)) == Approx(0.).margin(0.000001) );
507 REQUIRE( cosineSimilarity(p3, -Point(1,-1,2,-2)) == Approx(pi).margin(0.000001) );
508 REQUIRE( cosineSimilarity(p3, Point(1,1,2,2)) == Approx(pi/2).margin(0.000001) );
509
510 REQUIRE( p2.isLower(p4) == false );
511 REQUIRE( isLower(p2, p4) == false );
512 REQUIRE( p4.isUpper(p2) == false );
513 REQUIRE( isUpper(p4, p2) == false );
514 p4[1] = -p4[1];
515 REQUIRE( p2.isLower(p4) == true );
516 REQUIRE( isLower(p2, p4) == true );
517 REQUIRE( p4.isUpper(p2) == true );
518 REQUIRE( isUpper(p4, p2) == true );
519 }
520
521}
522
523
524TEST_CASE("Benchmarking","[.benchmark]")
525{
526 using Integer = DGtal::int32_t;
527 typedef PointVector<3, Integer> Point;
528 Point p1 = {1,2,3,4};
529 Point p2 = {3,4,5,6};
530
531 using Real = double;
532 typedef PointVector<3, Real> RPoint;
533 RPoint rp1 = {1,2,3,4};
534 RPoint rp2 = {3,4,5,6};
535
536 CHECK(p1.dot(p2) == 26);
537 CHECK(rp1.dot(rp2) == Approx(26));
538
539 BENCHMARK("Dot product int")
540 {
541 return p1.dot(p2);
542 };
543
544 BENCHMARK("Dot product double (with int->double cast)")
545 {
546 return rp1.dot(p2);
547 };
548
549 BENCHMARK("Dot product double")
550 {
551 return rp1.dot(rp2);
552 };
553
554}
555
DGtal::PointVector
Aim: Implements basic operations that will be used in Point and Vector classes.
Definition PointVector.h:593
DGtal::PointVector::max
Component max() const
DGtal::PointVector::isLower
bool isLower(const PointVector< dim, OtherComponent, OtherStorage > &p) const
Return true if this point is below a given point.
DGtal::PointVector::dot
auto dot(const PointVector< dim, OtherComponent, OtherStorage > &v) const -> decltype(DGtal::dotProduct(*this, v))
Dot product with a PointVector.
DGtal::PointVector::partialCopyInv
Self & partialCopyInv(const PointVector< dim, OtherComponent, OtherContainer > &pv, const std::vector< Dimension > &dimensions)
Partial copy of a given PointVector.
DGtal::PointVector::maxElement
Iterator maxElement()
DGtal::PointVector::inf
auto inf(const PointVector< dim, OtherComponent, OtherStorage > &aPoint) const -> decltype(DGtal::inf(*this, aPoint))
Implements the infimum (or greatest lower bound).
DGtal::PointVector::sup
auto sup(const PointVector< dim, OtherComponent, OtherStorage > &aPoint) const -> decltype(DGtal::sup(*this, aPoint))
Implements the supremum (or least upper bound).
DGtal::PointVector::crossProduct
auto crossProduct(const PointVector< dim, OtherComponent, OtherStorage > &v) const -> decltype(DGtal::crossProduct(*this, v))
Cross product with a PointVector.
DGtal::PointVector::isUpper
bool isUpper(const PointVector< dim, OtherComponent, OtherStorage > &p) const
Return true if this point is upper a given point.
DGtal::PointVector::min
Component min() const
DGtal::PointVector::ConstIterator
Container::const_iterator ConstIterator
Constant iterator type.
Definition PointVector.h:641
DGtal::PointVector::cosineSimilarity
double cosineSimilarity(const PointVector< dim, OtherComponent, OtherStorage > &v) const
Positive angle between two vectors, deduced from their scalar product.
DGtal::PointVector::minElement
Iterator minElement()
DGtal::PointVector::partialCopy
Self & partialCopy(const PointVector< dim, OtherComponent, OtherContainer > &pv, const std::vector< Dimension > &dimensions)
Partial copy of a given PointVector.
DGtal
DGtal is the top-level namespace which contains all DGtal functions and types.
Definition ClosedIntegerHalfPlane.h:49
DGtal::crossProduct
auto crossProduct(PointVector< 3, LeftEuclideanRing, LeftContainer > const &lhs, PointVector< 3, RightEuclideanRing, RightContainer > const &rhs) -> decltype(DGtal::constructFromArithmeticConversion(lhs, rhs))
Cross product of two 3D Points/Vectors.
DGtal::dotProduct
DGtal::ArithmeticConversionType< LeftEuclideanRing, RightEuclideanRing > dotProduct(PointVector< ptDim, LeftEuclideanRing, LeftContainer > const &lhs, PointVector< ptDim, RightEuclideanRing, RightContainer > const &rhs)
Dot product between two points/vectors.
DGtal::isUpper
bool isUpper(PointVector< ptDim, LeftEuclideanRing, LeftContainer > const &lhs, PointVector< ptDim, RightEuclideanRing, RightContainer > const &rhs)
Return true if the first point is upper the second point.
DGtal::cosineSimilarity
double cosineSimilarity(PointVector< ptDim, LeftEuclideanRing, LeftContainer > const &lhs, PointVector< ptDim, RightEuclideanRing, RightContainer > const &rhs)
Positive angle between two vectors, deduced from their scalar product.
DGtal::inf
auto inf(PointVector< ptDim, LeftEuclideanRing, LeftContainer > const &lhs, PointVector< ptDim, RightEuclideanRing, RightContainer > const &rhs) -> decltype(DGtal::constructFromArithmeticConversion(lhs, rhs))
Implements the infimum (or greatest lower bound).
DGtal::int32_t
boost::int32_t int32_t
signed 32-bit integer.
Definition BasicTypes.h:72
DGtal::sup
auto sup(PointVector< ptDim, LeftEuclideanRing, LeftContainer > const &lhs, PointVector< ptDim, RightEuclideanRing, RightContainer > const &rhs) -> decltype(DGtal::constructFromArithmeticConversion(lhs, rhs))
Implements the supremum (or least upper bound).
DGtal::isLower
bool isLower(PointVector< ptDim, LeftEuclideanRing, LeftContainer > const &lhs, PointVector< ptDim, RightEuclideanRing, RightContainer > const &rhs)
Return true if the first point is below the second point.
std
STL namespace.
RealPoint3D
Z3i::RealPoint RealPoint3D
Definition sphereCotangentLaplaceOperator.cpp:70
DGtal::functors::Round
Functor that rounds to the nearest integer.
Definition BasicFunctors.h:152
BENCHMARK
BENCHMARK(BM_StringCreation)
Point
MyPointD Point
Definition testClone2.cpp:383
CAPTURE
CAPTURE(thicknessHV)
TEST_CASE
TEST_CASE("2D Point Vector Unit tests")
Definition testPointVector.cpp:53
COMPARE_VALUE_AND_TYPE
#define COMPARE_VALUE_AND_TYPE(expr, check)
Definition testPointVector.cpp:49
aPoint
const Point aPoint(3, 4)
SECTION
SECTION("Testing constant forward iterators")
Definition testSimpleRandomAccessRangeFromPoint.cpp:66
REQUIRE
REQUIRE(domain.isInside(aPoint))
RealPoint
PointVector< 3, double > RealPoint
Definition testTriangulatedSurface.cpp:51