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
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
// 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/DirectionH3.h $
// $Id: include/CGAL/Homogeneous/DirectionH3.h b26b07a1242 $
// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-Commercial
//
//
// Author(s) : Stefan Schirra
#ifndef CGAL_HOMOGENEOUS_DIRECTION_3_H
#define CGAL_HOMOGENEOUS_DIRECTION_3_H
#include <CGAL/kernel_config.h>
#include <CGAL/array.h>
#include <CGAL/Handle_for.h>
#include <CGAL/number_utils.h>
namespace CGAL {
template < class R_ >
class DirectionH3
{
typedef typename R_::RT RT;
typedef typename R_::FT FT;
typedef typename R_::Point_3 Point_3;
typedef typename R_::Vector_3 Vector_3;
typedef typename R_::Segment_3 Segment_3;
typedef typename R_::Line_3 Line_3;
typedef typename R_::Ray_3 Ray_3;
typedef std::array<RT, 4> Rep;
typedef typename R_::template Handle<Rep>::type Base;
Base base;
public:
typedef R_ R;
DirectionH3() {}
//DirectionH3(const Point_3 & p )
//: base(p) {}
DirectionH3(const Vector_3 & v )
{ *this = v.rep().direction(); }
DirectionH3(const Line_3 & l )
{ *this = l.rep().direction(); }
DirectionH3(const Ray_3 & r )
{ *this = r.rep().direction(); }
DirectionH3(const Segment_3 & s )
{ *this = s.rep().direction(); }
// the fourth argument is not documented. Should go away ?
DirectionH3(const RT& x, const RT& y,
const RT& z, const RT& w = RT(1) )
: base( w >= RT(0) ? CGAL::make_array(x, y, z, w)
: CGAL::make_array<RT>(-x, -y, -z, -w) ) {}
typename R::Boolean is_degenerate() const;
typename R::Boolean operator==( const DirectionH3<R>& d) const;
typename R::Boolean operator!=( const DirectionH3<R>& d) const;
Vector_3 to_vector() const;
Vector_3 vector() const { return to_vector(); }
const RT & dx() const { return get_pointee_or_identity(base)[0]; }
const RT & dy() const { return get_pointee_or_identity(base)[1]; }
const RT & dz() const { return get_pointee_or_identity(base)[2]; }
const RT & x() const { return get_pointee_or_identity(base)[0]; }
const RT & y() const { return get_pointee_or_identity(base)[1]; }
const RT & z() const { return get_pointee_or_identity(base)[2]; }
const RT & hx() const { return get_pointee_or_identity(base)[0]; }
const RT & hy() const { return get_pointee_or_identity(base)[1]; }
const RT & hz() const { return get_pointee_or_identity(base)[2]; }
};
template <class R >
CGAL_KERNEL_INLINE
typename R::Boolean
DirectionH3<R>::operator==( const DirectionH3<R>& d) const
{
return ( ( hx()*d.hy() == hy()*d.hx() )
&&( hx()*d.hz() == hz()*d.hx() )
&&( hy()*d.hz() == hz()*d.hy() )
&&( CGAL_NTS sign( hx() ) == CGAL_NTS sign( d.hx() ) )
&&( CGAL_NTS sign( hy() ) == CGAL_NTS sign( d.hy() ) )
&&( CGAL_NTS sign( hz() ) == CGAL_NTS sign( d.hz() ) ) );
}
template <class R >
inline
typename R::Boolean
DirectionH3<R>::operator!=( const DirectionH3<R>& d) const
{ return !operator==(d); }
template <class R >
CGAL_KERNEL_INLINE
typename R::Boolean
DirectionH3<R>::is_degenerate() const
{ return ((hx() == RT(0)) && (hy() == RT(0)) && (hz() == RT(0))); }
template <class R >
inline
typename DirectionH3<R>::Vector_3
DirectionH3<R>::to_vector() const
{ return Vector_3(dx(), dy(), dz(), RT(1)); }
template <class R>
CGAL_KERNEL_INLINE
DirectionH3<R>
cross_product( const DirectionH3<R>& d1,
const DirectionH3<R>& d2)
{ return cross_product(d1.to_vector(),d2.to_vector()).direction(); }
} //namespace CGAL
#endif // CGAL_HOMOGENEOUS_DIRECTION_3_H
|
