1 // Boost.Geometry (aka GGL, Generic Geometry Library)
2 // Unit Test
3
4 // Copyright (c) 2009-2012 Mateusz Loskot, London, UK.
5 // Copyright (c) 2017 Adam Wulkiewicz, Lodz, Poland.
6
7 // Use, modification and distribution is subject to the Boost Software License,
8 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
9 // http://www.boost.org/LICENSE_1_0.txt)
10
11
12 #include <geometry_test_common.hpp>
13
14 #include <boost/core/ignore_unused.hpp>
15
16 #include <boost/geometry/arithmetic/cross_product.hpp>
17
18 #include <boost/geometry/algorithms/assign.hpp>
19
20 #include <boost/geometry/geometries/point.hpp>
21 #include <boost/geometry/geometries/adapted/c_array.hpp>
22 #include <boost/geometry/geometries/adapted/boost_tuple.hpp>
23 #include <test_common/test_point.hpp>
24
25
26 template <typename P>
test_2d()27 void test_2d()
28 {
29 P p1;
30 bg::assign_values(p1, 20, 30);
31 P p2;
32 bg::assign_values(p2, 45, 70);
33 P c = bg::cross_product(p1, p2);
34
35 typedef typename bg::coordinate_type<P>::type scalar_type;
36 BOOST_CHECK_EQUAL(bg::get<0>(c), scalar_type(50));
37 }
38
39 template <typename P>
test_3d()40 void test_3d()
41 {
42 P p1;
43 bg::assign_values(p1, 20, 30, 10);
44 P p2;
45 bg::assign_values(p2, 45, 70, 20);
46 P c = bg::cross_product(p1, p2);
47
48 typedef typename bg::coordinate_type<P>::type scalar_type;
49 BOOST_CHECK_EQUAL(bg::get<0>(c), scalar_type(-100));
50 BOOST_CHECK_EQUAL(bg::get<1>(c), scalar_type(50));
51 BOOST_CHECK_EQUAL(bg::get<2>(c), scalar_type(50));
52 }
53
54 #ifdef TEST_FAIL_CROSS_PRODUCT
55 template <typename P>
test_4d()56 void test_4d()
57 {
58 P p1;
59 bg::assign_values(p1, 20, 30, 10);
60 bg::set<3>(p1, 15);
61 P p2;
62 bg::assign_values(p2, 45, 70, 20);
63 bg::set<3>(p2, 35);
64 P c = bg::cross_product(p1, p2);
65 boost::ignore_unused(c);
66 }
67 #endif
68
test_main(int,char * [])69 int test_main(int, char* [])
70 {
71 test_2d<bg::model::point<int, 2, bg::cs::cartesian> >();
72 test_2d<bg::model::point<float, 2, bg::cs::cartesian> >();
73 test_2d<bg::model::point<double, 2, bg::cs::cartesian> >();
74
75 test_3d<bg::model::point<int, 3, bg::cs::cartesian> >();
76 test_3d<bg::model::point<float, 3, bg::cs::cartesian> >();
77 test_3d<bg::model::point<double, 3, bg::cs::cartesian> >();
78
79 #ifdef TEST_FAIL_CROSS_PRODUCT
80 test_4d<bg::model::point<int, 4, bg::cs::cartesian> >();
81 test_4d<bg::model::point<float, 4, bg::cs::cartesian> >();
82 test_4d<bg::model::point<double, 4, bg::cs::cartesian> >();
83 #endif
84
85 return 0;
86 }
87
88