Reverted incorrect merge (missing files)
[blender.git] / extern / solid / src / convex / DT_Facet.h
1 /*
2  * SOLID - Software Library for Interference Detection
3  * 
4  * Copyright (C) 2001-2003  Dtecta.  All rights reserved.
5  *
6  * This library may be distributed under the terms of the Q Public License
7  * (QPL) as defined by Trolltech AS of Norway and appearing in the file
8  * LICENSE.QPL included in the packaging of this file.
9  *
10  * This library may be distributed and/or modified under the terms of the
11  * GNU General Public License (GPL) version 2 as published by the Free Software
12  * Foundation and appearing in the file LICENSE.GPL included in the
13  * packaging of this file.
14  *
15  * This library is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
16  * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
17  *
18  * Commercial use or any other use of this library not covered by either 
19  * the QPL or the GPL requires an additional license from Dtecta. 
20  * Please contact info@dtecta.com for enquiries about the terms of commercial
21  * use of this library.
22  */
23
24 #ifndef DT_FACET_H
25 #define DT_FACET_H
26
27 #include <string.h>
28 #include <vector>
29
30 #include <MT_Vector3.h>
31 #include <MT_Point3.h>
32
33 class DT_Facet;
34
35
36 class DT_Edge {
37 public:
38     DT_Edge() {}
39     DT_Edge(DT_Facet *facet, int index) : 
40         m_facet(facet), 
41         m_index(index) {}
42
43     DT_Facet *getFacet() const { return m_facet; }
44     int       getIndex() const { return m_index; }
45
46     int getSource() const;
47     int getTarget() const;
48
49 private:    
50     DT_Facet *m_facet;
51     int       m_index;
52 };
53
54 typedef std::vector<DT_Edge> DT_EdgeBuffer;
55
56
57 class DT_Facet {
58 public:
59     DT_Facet() {}
60     DT_Facet(int i0, int i1, int i2) 
61           :     m_obsolete(false) 
62     {
63                 m_indices[0] = i0; 
64                 m_indices[1] = i1; 
65                 m_indices[2] = i2;
66     }
67         
68     int operator[](int i) const { return m_indices[i]; } 
69
70     bool link(int edge0, DT_Facet *facet, int edge1);
71
72     
73     bool isObsolete() const { return m_obsolete; }
74     
75
76     bool computeClosest(const MT_Vector3 *verts);
77     
78     const MT_Vector3& getClosest() const { return m_closest; } 
79     
80     bool isClosestInternal() const
81         { 
82                 return m_lambda1 >= MT_Scalar(0.0) && 
83                         m_lambda2 >= MT_Scalar(0.0) && 
84                         m_lambda1 + m_lambda2 <= m_det;
85     } 
86
87     MT_Scalar getDist2() const { return m_dist2; }
88         
89     MT_Point3 getClosestPoint(const MT_Point3 *points) const 
90         {
91                 const MT_Point3& p0 = points[m_indices[0]];
92                 
93                 return p0 + (m_lambda1 * (points[m_indices[1]] - p0) + 
94                                          m_lambda2 * (points[m_indices[2]] - p0)) / m_det;
95     }
96     
97     void silhouette(const MT_Vector3& w, DT_EdgeBuffer& edgeBuffer) 
98         {
99                 edgeBuffer.clear();
100                 m_obsolete = true;
101                 m_adjFacets[0]->silhouette(m_adjEdges[0], w, edgeBuffer);
102                 m_adjFacets[1]->silhouette(m_adjEdges[1], w, edgeBuffer);
103                 m_adjFacets[2]->silhouette(m_adjEdges[2], w, edgeBuffer);
104     }
105         
106 private:
107     void silhouette(int index, const MT_Vector3& w, DT_EdgeBuffer& edgeBuffer);
108         
109     int         m_indices[3];
110     bool        m_obsolete;
111     DT_Facet   *m_adjFacets[3];
112     int         m_adjEdges[3];
113         
114     MT_Scalar   m_det;
115     MT_Scalar   m_lambda1;
116     MT_Scalar   m_lambda2;
117     MT_Vector3  m_closest;
118     MT_Scalar   m_dist2;
119 };
120
121
122 inline int incMod3(int i) { return ++i % 3; } 
123
124 inline int DT_Edge::getSource() const 
125 {
126     return (*m_facet)[m_index];
127 }
128
129 inline int DT_Edge::getTarget() const 
130 {
131     return (*m_facet)[incMod3(m_index)];
132 }
133
134 #endif