Salome HOME
Merge remote branch 'origin/hydro/imps_2015'
[modules/geom.git] / src / GEOMUtils / GEOMUtils.hxx
1 // Copyright (C) 2007-2016  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
10 //
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 // Lesser General Public License for more details.
15 //
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 #ifndef _GEOMUtils_HXX_
24 #define _GEOMUtils_HXX_
25
26 #include <Standard_Macro.hxx>
27 #include <TopoDS_Shape.hxx>
28 #include <TopoDS_Vertex.hxx>
29
30 #include <TopTools_ListOfShape.hxx>
31
32 #include <TopAbs.hxx>
33
34 #include <gp_Ax3.hxx>
35 #include <gp_Vec.hxx>
36
37 #include <V3d_View.hxx>
38
39 #include <NCollection_DataMap.hxx>
40
41 #include <functional>
42
43 #include <map>
44 #include <vector>
45 #include <string>
46 #include <utility>
47
48 class Bnd_Box;
49
50 inline Standard_Boolean IsEqual (const TopoDS_Shape& S1, const TopoDS_Shape& S2)
51 {
52   return S1.IsSame(S2);
53 }
54
55 namespace GEOMUtils
56 {
57
58   /**
59    * This enumeration represents comparison conditions.
60    */
61   enum ComparisonCondition {
62     CC_GT, ///< Greater then
63     CC_GE, ///< Greater then or equal to
64     CC_LT, ///< Less then
65     CC_LE  ///< Less then or equal to
66   };
67
68   typedef std::vector<std::string> NodeLinks;
69   typedef std::map<std::string, NodeLinks> LevelInfo;
70   typedef std::vector<LevelInfo> LevelsList;
71   typedef std::map<std::string,std::pair<LevelsList,LevelsList> > TreeModel;
72
73   /*!
74    * \brief Compute numerical functor for the shape.
75    *
76    * Resulting value can be used to sort out shapes according to some parameter.
77    * 
78    * Returns a pair of two values (dist, functor) where
79    * - \a dist is a some value that is computed according to the center of mass of given shape;
80    * - \a functor is a numerical functor value
81    *
82    * The numerical functor is computed according to the shape's topological properties as follows:
83    * - orientation for vertices 
84    * - length for edges and wires
85    * - area for faces and shells
86    * - volume for solids, compounds, compsolids
87    *
88    * If \a isOldSorting parameter is set to \c true, for all cases linear properties of the shape
89    * are used (to support backward compatibility in some methods). By default, this parameter is
90    * set to \c false.
91    */
92   Standard_EXPORT std::pair<double, double> ShapeToDouble (const TopoDS_Shape& theShape,
93                                                            bool isOldSorting = false);
94
95   /*!
96    * \brief Get Local Coordinate System, corresponding to the given shape.
97    *
98    * Origin of the LCS is situated at the shape's center of mass.
99    * Axes of the LCS are obtained from shape's location or,
100    * if the shape is a planar face, from position of its plane.
101    */
102   Standard_EXPORT gp_Ax3 GetPosition (const TopoDS_Shape& theShape);
103
104   /*!
105    * \brief Get vector, defined by the given edge.
106    * \param theShape The edge.
107    * \param doConsiderOrientation If True, take into account the edge orientation.
108    * \note It is recommended to use doConsiderOrientation=Standard_False, because
109    *       the same edge can have different orientation depending on the way it was
110    *       extracted from a shape.
111    */
112   Standard_EXPORT gp_Vec GetVector (const TopoDS_Shape& theShape,
113                                     Standard_Boolean doConsiderOrientation);
114
115   /*!
116    * \brief Sort shapes in the list by their coordinates.
117    * \param SL The list of shapes to sort.
118    */
119   struct CompareShapes : public std::binary_function<TopoDS_Shape, TopoDS_Shape, bool>
120   {
121     CompareShapes (bool isOldSorting)
122       : myIsOldSorting(isOldSorting) {}
123
124     bool operator() (const TopoDS_Shape& lhs, const TopoDS_Shape& rhs);
125
126     typedef NCollection_DataMap<TopoDS_Shape, std::pair<double, double> > GEOMUtils_DataMapOfShapeDouble;
127     GEOMUtils_DataMapOfShapeDouble myMap;
128     bool myIsOldSorting;
129   };
130
131   /*!
132    * \brief Sort shapes by their centers of mass, using formula X*999 + Y*99 + Z*0.9
133    */
134   Standard_EXPORT void SortShapes (TopTools_ListOfShape& SL,
135                                    const Standard_Boolean isOldSorting = Standard_True);
136
137   /*!
138    * \brief Convert TopoDS_COMPSOLID to TopoDS_COMPOUND.
139    *
140    * If the argument shape is not of type TopoDS_COMPSOLID, this method returns it as is.
141    *
142    * \param theCompsolid The compsolid to be converted.
143    * \retval TopoDS_Shape Returns the resulting compound.
144    */
145   Standard_EXPORT TopoDS_Shape CompsolidToCompound (const TopoDS_Shape& theCompsolid);
146
147   /*!
148    * \brief Recursively extract all shapes from compounds and compsolids of the given shape into theList.
149    *
150    * If theShape is not compound or compsolid, theList will contain only theShape itself.
151    *
152    * \param theShape The shape to be exploded.
153    * \param theList Output parameter.
154    */
155   Standard_EXPORT void AddSimpleShapes (const TopoDS_Shape& theShape,
156                                         TopTools_ListOfShape& theList);
157
158   /*!
159    * \brief Build a triangulation on \a theShape if it is absent.
160    * \param theShape The shape to check/build triangulation on.
161    * \retval bool Returns false if the shape has no faces, i.e. impossible to build triangulation.
162    */
163   Standard_EXPORT bool CheckTriangulation (const TopoDS_Shape& theShape);
164   
165   /*!
166    * \brief Return type of shape for explode. In case of compound it will be a type of its first sub shape.
167    * \param theShape The shape to get type of.
168    * \retval TopAbs_ShapeEnum Return type of shape for explode.
169    */
170   Standard_EXPORT TopAbs_ShapeEnum GetTypeOfSimplePart (const TopoDS_Shape& theShape);
171
172   /*!
173    * \brief Find an edge of theShape, closest to thePoint.
174    *
175    * \param theShape The shape to explore.
176    * \param thePoint The point near the required edge.
177    * \retval TopoDS_Shape Returns the found edge or an empty shape if multiple edges found.
178    */
179   Standard_EXPORT TopoDS_Shape GetEdgeNearPoint (const TopoDS_Shape&  theShape,
180                                                  const TopoDS_Vertex& thePoint);
181
182   /*!
183    * \brief Compute precise bounding box of the shape based on the rough bounding box.
184    *
185    * \param theShape the shape.
186    * \param theBox rough bounding box on input; precise bounding box on output.
187    * \retval Standard_True in case of success; Standard_False otherwise.
188    */
189   Standard_EXPORT Standard_Boolean PreciseBoundingBox(const TopoDS_Shape &theShape, Bnd_Box &theBox);
190
191   /*!
192    * \brief Computes minumal distance between two shapes for singular cases
193    *        (workaround for bugs 19899, 19908 and 19910 from Mantis).
194    *
195    * \param aSh1 the first shape
196    * \param aSh2 the second shape
197    * \param Ptmp1 the output result point on the first shape
198    * \param Ptmp2 the output result point on the second shape
199    * \retval negative value if it is not a singular case; actual distance for singular case.
200    */
201   Standard_EXPORT Standard_Real GetMinDistanceSingular(const TopoDS_Shape& aSh1,
202                                                        const TopoDS_Shape& aSh2,
203                                                        gp_Pnt& Ptmp1, gp_Pnt& Ptmp2);
204   
205   /*!
206    * \brief Computes minumal distance between two shapes.
207    *
208    * \param theShape1 the first shape
209    * \param theShape2 the second shape
210    * \param thePnt1 the output result point on the first shape
211    * \param thePnt2 the output result point on the second shape
212    * \retval negative value in case of failure; otherwise the real distance.
213    */
214   Standard_EXPORT Standard_Real GetMinDistance(const TopoDS_Shape& theShape1,
215                                                const TopoDS_Shape& theShape2,
216                                                gp_Pnt& thePnt1, gp_Pnt& thePnt2);
217   
218   /*!
219    * \brief Returns the point clicked in 3D view.
220    *
221    * \param x The X coordinate in the view.
222    * \param y The Y coordinate in the view.
223    * \param theView View where the given point takes place.
224    * \retval gp_Pnt Returns the point clicked in 3D view
225    */
226   Standard_EXPORT gp_Pnt ConvertClickToPoint( int x, int y, Handle(V3d_View) theView );
227
228   /*!
229    * \brief Convert dependency tree data to the string representation
230    *
231    * \param tree dependency tree data
232    * \param dependencyStr output string
233    */
234   Standard_EXPORT void ConvertTreeToString( const TreeModel& tree,
235                                             std::string& dependencyStr );
236
237   /*!
238    * \brief Restore dependency tree data from the string representation
239    *
240    * \param dependencyStr string representation of tree data
241    * \param tree output dependency tree data
242    */
243   Standard_EXPORT void ConvertStringToTree( const std::string& dependencyStr,
244                                             TreeModel& tree );
245
246   /*!
247    * \brief Check shape
248    *
249    * \param shape input shape object
250    * \param checkGeometry when set to \c true, causes check of underlying geometry
251    *        in addition to the topology
252    * \return \c true if shape is valid or \c false otherwise
253    */
254   Standard_EXPORT bool CheckShape( TopoDS_Shape& shape, bool checkGeometry = false );
255
256   /*!
257    * \brief Check boolean and partition operations agruments
258    *
259    * \param theShape the agrument of an operation to be checked
260    * \return \c true if the agrument is valid for a boolean or partition
261    *         operation or \c false otherwise
262    */
263   Standard_EXPORT bool CheckBOPArguments(const TopoDS_Shape &theShape);
264   
265   /*!
266    * \brief Limit shape tolerance to the given value
267    *
268    * \param shape shape being fixed
269    * \param type topology type which tolerance is to be limited; TopAbs_SHAPE means
270    *             all types of topology
271    * \param tolerance expected tolerance value (1e-7 by default)
272    * \param checkGeometry check geometry validity of result
273    * \return \c true if resulting shape is valid
274    *
275    * \note Resulting tolerance of the shape is not mandatory equal to requested value
276    *       as it might be changed by fixshape operation in order to get valid shape where possible
277    * \note By default, result only checked for topology validity; check of geometry can be done by
278    *       passing \c true to \a checkGeometry parameter
279    */
280   Standard_EXPORT bool FixShapeTolerance( TopoDS_Shape& shape,
281                                           TopAbs_ShapeEnum type,
282                                           Standard_Real tolerance = Precision::Confusion(),
283                                           bool checkGeometry = false );
284
285   /*!
286    * \brief Limit shape tolerance to the given value
287    * This is overloaded function, it behaves exactly as previous one
288    */
289   Standard_EXPORT bool FixShapeTolerance( TopoDS_Shape& shape,
290                                           Standard_Real tolerance = Precision::Confusion(),
291                                           bool checkGeometry = false );
292
293   /*!
294    * \brief Limit shape tolerance to the given value
295    * This is overloaded function, it behaves exactly as previous one
296    */
297   Standard_EXPORT bool FixShapeTolerance( TopoDS_Shape& shape,
298                                           bool checkGeometry );
299
300   /*!
301    * \brief Fix curves of the given shape
302    * 
303    * The function checks each curve of the input shape in the following way:
304    * - compute deviation of the curve from the underlying surface in a set of points
305    *   computed with the certain discretization step value
306    * - find maximum tolerance between computed deviation values
307    * - limit tolerance of the curve with the computed maximum value
308    * 
309    * \param shape shape being fixed
310    * \return \c true if resulting shape is valid
311    */
312   Standard_EXPORT bool FixShapeCurves( TopoDS_Shape& shape );
313
314   /*!
315    * \brief Write shape to the BREP file
316    *
317    * \param source shape
318    * \return \c true if file was written or \c false otherwise
319    */
320   Standard_EXPORT bool Write( const TopoDS_Shape& shape,
321                               const char* fileName );
322   
323   /*!
324    * \brief Extract single SOLID from COMPSOLID or COMPOUND.
325    *
326    * If the argument shape is a COMPUND or COMPSOLID and there's
327    * only single simple-shape type inside, this sub-shape is returned as a result;
328    * otherwise, the shape is not changed.
329    *
330    * \param shape compound or compsolid being processed.
331    * \retval TopoDS_Shape resulting shape
332    */
333   Standard_EXPORT TopoDS_Shape ReduceCompound( const TopoDS_Shape& shape );
334
335   /*!
336    * \brief Generate triangulation for the shape.
337    *
338    * \param shape shape being meshed
339    * \param deflection deflection coefficient to be used
340    * \param forced if \c true, causes generation of mesh regardless it is already present in the shape
341    */
342   Standard_EXPORT void MeshShape( const TopoDS_Shape shape,
343                                   double deflection, bool forced = true );
344
345   /*!
346    * \brief Get default deflection coefficient used for triangulation
347    * \return default deflection value
348    */
349   Standard_EXPORT double DefaultDeflection();
350
351   /**
352    * \brief Check if the shape is not a closed wire or edge.
353    *
354    * This function is used for pipe creation algorithm to test if
355    * the pipe path is not closed. It returns false if theShape is a wire or
356    * an edge with coincident end vertices. Otherwise it returns true.
357    *
358    * \param theShape the shape to be tested.
359    * \return true if theShape is not a closed wire or edge.
360    */
361   Standard_EXPORT bool IsOpenPath(const TopoDS_Shape &theShape);
362
363   /**
364    * This function compares two tolerances. The shape tolerance (the first
365    * argument) is considered less than the reference tolerance (the second
366    * argument) if theTolShape < theTolRef - Tolerance(theTolRef). theTolShape is
367    * considered greater than theTolRef if theTolShape > theTolRef +
368    * Tolerance(theTolRef). Otherwise these tolerances are equal.
369    * Tolerance(theTolRef) = theTolRef*DEFAULT_TOLERANCE_TOLERANCE. But this value
370    * should not be greated than DEFAULT_MAX_TOLERANCE_TOLERANCE.
371    *
372    * \param theTolShape the shape tolerance
373    * \param theTolRef the reference tolerance
374    * \return -1 if theTolShape is less than theTolRef; 1 if theTolShape is greater
375    * than theTolRef; 0 if they are equal
376    */
377   Standard_EXPORT int CompareToleranceValues(const double theTolShape,
378                                              const double theTolRef);
379
380   /**
381    * Check if the comarison of tolerances fit the condition. The comparison of
382    * tolerances is performed using the function CompareToleranceValues.
383    *
384    * \param theCondition the condition
385    * \param theTolShape the shape tolerance
386    * \param theTolRef the reference tolerance
387    * \return true if the shape tolerance fits the condition; false otherwise.
388    */
389   Standard_EXPORT bool IsFitCondition(const ComparisonCondition theCondition,
390                                       const double              theTolShape,
391                                       const double              theTolRef);
392
393 };
394
395 #endif