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
|
// Copyright (c) 2014
// INRIA Saclay-Ile de France (France)
//
// This file is part of CGAL (www.cgal.org); you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 3 of the License,
// or (at your option) any later version.
//
// Licensees holding a valid commercial license may use this file in
// accordance with the commercial license agreement provided with the software.
//
// This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
// WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
//
// $URL$
// $Id$
//
// Author(s) : Marc Glisse
#ifndef CGAL_DETVEC_H
#define CGAL_DETVEC_H
#include <CGAL/determinant.h>
#include <CGAL/predicates/sign_of_determinant.h>
namespace CGAL {
// TODO: determine whether it is better to pass them by lines or columns.
template <class NT, class Vector> inline
NT determinant_of_vectors(Vector const&a, Vector const&b){
return determinant<NT>(a[0],a[1],b[0],b[1]);
}
template <class NT, class Vector> inline
typename Sgn<NT>::result_type
sign_of_determinant_of_vectors(Vector const&a, Vector const&b){
return sign_of_determinant<NT>(a[0],a[1],b[0],b[1]);
}
template <class NT, class Vector>
NT determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c){
return determinant<NT>(a[0],a[1],a[2],b[0],b[1],b[2],c[0],c[1],c[2]);
}
template <class NT, class Vector>
typename Sgn<NT>::result_type
sign_of_determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c){
return sign_of_determinant<NT>(a[0],a[1],a[2],b[0],b[1],b[2],c[0],c[1],c[2]);
}
template <class NT, class Vector>
NT determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c, Vector const&d){
return determinant<NT>(
a[0],a[1],a[2],a[3],
b[0],b[1],b[2],b[3],
c[0],c[1],c[2],c[3],
d[0],d[1],d[2],d[3]);
}
template <class NT, class Vector>
typename Sgn<NT>::result_type
sign_of_determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c, Vector const&d){
return sign_of_determinant<NT>(
a[0],a[1],a[2],a[3],
b[0],b[1],b[2],b[3],
c[0],c[1],c[2],c[3],
d[0],d[1],d[2],d[3]);
}
template <class NT, class Vector>
NT determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c, Vector const&d, Vector const&e){
return determinant<NT>(
a[0],a[1],a[2],a[3],a[4],
b[0],b[1],b[2],b[3],b[4],
c[0],c[1],c[2],c[3],c[4],
d[0],d[1],d[2],d[3],d[4],
e[0],e[1],e[2],e[3],e[4]);
}
template <class NT, class Vector>
typename Sgn<NT>::result_type
sign_of_determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c, Vector const&d, Vector const&e){
return sign_of_determinant<NT>(
a[0],a[1],a[2],a[3],a[4],
b[0],b[1],b[2],b[3],b[4],
c[0],c[1],c[2],c[3],c[4],
d[0],d[1],d[2],d[3],d[4],
e[0],e[1],e[2],e[3],e[4]);
}
template <class NT, class Vector>
NT determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c, Vector const&d, Vector const&e, Vector const&f){
return determinant<NT>(
a[0],a[1],a[2],a[3],a[4],a[5],
b[0],b[1],b[2],b[3],b[4],b[5],
c[0],c[1],c[2],c[3],c[4],c[5],
d[0],d[1],d[2],d[3],d[4],d[5],
e[0],e[1],e[2],e[3],e[4],e[5],
f[0],f[1],f[2],f[3],f[4],f[5]);
}
template <class NT, class Vector>
typename Sgn<NT>::result_type
sign_of_determinant_of_vectors(Vector const&a, Vector const&b,
Vector const&c, Vector const&d, Vector const&e, Vector const&f){
return sign_of_determinant<NT>(
a[0],a[1],a[2],a[3],a[4],a[5],
b[0],b[1],b[2],b[3],b[4],b[5],
c[0],c[1],c[2],c[3],c[4],c[5],
d[0],d[1],d[2],d[3],d[4],d[5],
e[0],e[1],e[2],e[3],e[4],e[5],
f[0],f[1],f[2],f[3],f[4],f[5]);
}
}
#endif
|
