1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
|
// Copyright (c) 1999
// Utrecht University (The Netherlands),
// ETH Zurich (Switzerland),
// INRIA Sophia-Antipolis (France),
// Max-Planck-Institute Saarbruecken (Germany),
// and Tel-Aviv University (Israel). All rights reserved.
//
// This file is part of CGAL (www.cgal.org)
//
// $URL: https://github.com/CGAL/cgal/blob/v6.1/Homogeneous_kernel/include/CGAL/Homogeneous/DirectionH2.h $
// $Id: include/CGAL/Homogeneous/DirectionH2.h b26b07a1242 $
// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-Commercial
//
//
// Author(s) : Stefan Schirra
#ifndef CGAL_HOMOGENEOUS_DIRECTION_2_H
#define CGAL_HOMOGENEOUS_DIRECTION_2_H
#include <CGAL/array.h>
#include <CGAL/Handle_for.h>
#include <CGAL/kernel_config.h>
#include <CGAL/number_type_basic.h>
namespace CGAL {
template < class R_ >
class DirectionH2
{
typedef DirectionH2<R_> Self;
typedef typename R_::FT FT;
typedef typename R_::RT RT;
typedef typename R_::Point_2 Point_2;
typedef typename R_::Vector_2 Vector_2;
typedef typename R_::Line_2 Line_2;
typedef typename R_::Ray_2 Ray_2;
typedef typename R_::Segment_2 Segment_2;
typedef std::array<RT, 3> Rep;
typedef typename R_::template Handle<Rep>::type Base;
Base base;
public:
typedef R_ R;
typedef const RT& Homogeneous_coordinate_type;
const Self&
rep() const
{
return static_cast<const Self& >(*this);
}
DirectionH2() {}
DirectionH2(const RT& x, const RT& y)
: base{x, y, RT(1)} {}
// TODO Not documented : should not exist, not used.
// we should also change array<RT, 3> -> array<RT, 2>
DirectionH2(const RT& x, const RT& y, const RT& w )
: base( w > RT(0) ? CGAL::make_array(x, y, w)
: CGAL::make_array<RT>(-x, -y, -w) ) {}
typename R_::Boolean operator==( const DirectionH2<R>& d) const;
typename R_::Boolean operator!=( const DirectionH2<R>& d) const;
Vector_2 to_vector() const;
const RT & x() const { return CGAL::get_pointee_or_identity(base)[0]; }
const RT & y() const { return CGAL::get_pointee_or_identity(base)[1]; }
const RT & delta(int i) const;
const RT & dx() const { return CGAL::get_pointee_or_identity(base)[0]; }
const RT & dy() const { return CGAL::get_pointee_or_identity(base)[1]; }
};
template <class R >
CGAL_KERNEL_INLINE
typename R::Boolean
DirectionH2<R>::operator==( const DirectionH2<R>& d) const
{
return ( ( x() * d.y() == y() * d.x() )
&&( CGAL_NTS sign( x() ) == CGAL_NTS sign( d.x() ) )
&&( CGAL_NTS sign( y() ) == CGAL_NTS sign( d.y() ) ) );
}
template <class R >
inline
typename R::Boolean
DirectionH2<R>::operator!=( const DirectionH2<R>& d) const
{ return !(*this == d); }
} //namespace CGAL
#include <CGAL/Homogeneous/predicates_on_directionsH2.h>
namespace CGAL {
template <class R >
CGAL_KERNEL_INLINE
typename DirectionH2<R>::Vector_2
DirectionH2<R>::to_vector() const
{ return Vector_2(dx(), dy()); }
} //namespace CGAL
#endif // CGAL_HOMOGENEOUS_DIRECTION_2_H
|
