Salome HOME
e8a8272a8e97fa697633f11d52e7db6b5e91625f
[modules/geom.git] / src / GEOM_I / GEOM_IMeasureOperations_i.cc
1 // Copyright (C) 2007-2014  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 #include <Standard_Stream.hxx>
24
25 #include "GEOM_IMeasureOperations_i.hh"
26
27 #include "utilities.h"
28 #include "OpUtil.hxx"
29
30 #include "GEOM_Engine.hxx"
31 #include "GEOM_Object.hxx"
32
33 /**
34  * This function converts shape errors from theErrorsFrom to theErrorsTo.
35  * Note that theErrorsTo is not cleared at first.
36  *
37  * \param theErrorsFrom errors to be converted.
38  * \param theErrorsTo result errors.
39  */
40 static void ConvertShapeError
41     (const GEOM::GEOM_IMeasureOperations::ShapeErrors         &theErrorsFrom,
42            std::list<GEOMImpl_IMeasureOperations::ShapeError> &theErrorsTo)
43 {
44   int aNbErr = theErrorsFrom.length();
45   int i = 0;
46
47   for (; i < aNbErr; i++) {
48     const GEOM::GEOM_IMeasureOperations::ShapeError anErr = theErrorsFrom[i];
49     const GEOM::GEOM_IMeasureOperations::ShapeErrorType aType = anErr.error;
50     const GEOM::ListOfLong anIncrims = anErr.incriminated;
51     GEOMImpl_IMeasureOperations::ShapeError anErrStruct;
52
53     switch (aType) {
54     case GEOM::GEOM_IMeasureOperations::InvalidPointOnCurve:
55       anErrStruct.error = BRepCheck_InvalidPointOnCurve;
56       break;
57     case GEOM::GEOM_IMeasureOperations::InvalidPointOnCurveOnSurface:
58       anErrStruct.error = BRepCheck_InvalidPointOnCurveOnSurface;
59       break;
60     case GEOM::GEOM_IMeasureOperations::InvalidPointOnSurface:
61       anErrStruct.error = BRepCheck_InvalidPointOnSurface;
62       break;
63     case GEOM::GEOM_IMeasureOperations::No3DCurve:
64       anErrStruct.error = BRepCheck_No3DCurve;
65       break;
66     case GEOM::GEOM_IMeasureOperations::Multiple3DCurve:
67       anErrStruct.error = BRepCheck_Multiple3DCurve;
68       break;
69     case GEOM::GEOM_IMeasureOperations::Invalid3DCurve:
70       anErrStruct.error = BRepCheck_Invalid3DCurve;
71       break;
72     case GEOM::GEOM_IMeasureOperations::NoCurveOnSurface:
73       anErrStruct.error = BRepCheck_NoCurveOnSurface;
74       break;
75     case GEOM::GEOM_IMeasureOperations::InvalidCurveOnSurface:
76       anErrStruct.error = BRepCheck_InvalidCurveOnSurface;
77       break;
78     case GEOM::GEOM_IMeasureOperations::InvalidCurveOnClosedSurface:
79       anErrStruct.error = BRepCheck_InvalidCurveOnClosedSurface;
80       break;
81     case GEOM::GEOM_IMeasureOperations::InvalidSameRangeFlag:
82       anErrStruct.error = BRepCheck_InvalidSameRangeFlag;
83       break;
84     case GEOM::GEOM_IMeasureOperations::InvalidSameParameterFlag:
85       anErrStruct.error = BRepCheck_InvalidSameParameterFlag;
86       break;
87     case GEOM::GEOM_IMeasureOperations::InvalidDegeneratedFlag:
88       anErrStruct.error = BRepCheck_InvalidDegeneratedFlag;
89       break;
90     case GEOM::GEOM_IMeasureOperations::FreeEdge:
91       anErrStruct.error = BRepCheck_FreeEdge;
92       break;
93     case GEOM::GEOM_IMeasureOperations::InvalidMultiConnexity:
94       anErrStruct.error = BRepCheck_InvalidMultiConnexity;
95       break;
96     case GEOM::GEOM_IMeasureOperations::InvalidRange:
97       anErrStruct.error = BRepCheck_InvalidRange;
98       break;
99     case GEOM::GEOM_IMeasureOperations::EmptyWire:
100       anErrStruct.error = BRepCheck_EmptyWire;
101       break;
102     case GEOM::GEOM_IMeasureOperations::RedundantEdge:
103       anErrStruct.error = BRepCheck_RedundantEdge;
104       break;
105     case GEOM::GEOM_IMeasureOperations::SelfIntersectingWire:
106       anErrStruct.error = BRepCheck_SelfIntersectingWire;
107       break;
108     case GEOM::GEOM_IMeasureOperations::NoSurface:
109       anErrStruct.error = BRepCheck_NoSurface;
110       break;
111     case GEOM::GEOM_IMeasureOperations::InvalidWire:
112       anErrStruct.error = BRepCheck_InvalidWire;
113       break;
114     case GEOM::GEOM_IMeasureOperations::RedundantWire:
115       anErrStruct.error = BRepCheck_RedundantWire;
116       break;
117     case GEOM::GEOM_IMeasureOperations::IntersectingWires:
118       anErrStruct.error = BRepCheck_IntersectingWires;
119       break;
120     case GEOM::GEOM_IMeasureOperations::InvalidImbricationOfWires:
121       anErrStruct.error = BRepCheck_InvalidImbricationOfWires;
122       break;
123     case GEOM::GEOM_IMeasureOperations::EmptyShell:
124       anErrStruct.error = BRepCheck_EmptyShell;
125       break;
126     case GEOM::GEOM_IMeasureOperations::RedundantFace:
127       anErrStruct.error = BRepCheck_RedundantFace;
128       break;
129     case GEOM::GEOM_IMeasureOperations::UnorientableShape:
130       anErrStruct.error = BRepCheck_UnorientableShape;
131       break;
132     case GEOM::GEOM_IMeasureOperations::NotClosed:
133       anErrStruct.error = BRepCheck_NotClosed;
134       break;
135     case GEOM::GEOM_IMeasureOperations::NotConnected:
136       anErrStruct.error = BRepCheck_NotConnected;
137       break;
138     case GEOM::GEOM_IMeasureOperations::SubshapeNotInShape:
139       anErrStruct.error = BRepCheck_SubshapeNotInShape;
140       break;
141     case GEOM::GEOM_IMeasureOperations::BadOrientation:
142       anErrStruct.error = BRepCheck_BadOrientation;
143       break;
144     case GEOM::GEOM_IMeasureOperations::BadOrientationOfSubshape:
145       anErrStruct.error = BRepCheck_BadOrientationOfSubshape;
146       break;
147     case GEOM::GEOM_IMeasureOperations::InvalidToleranceValue:
148       anErrStruct.error = BRepCheck_InvalidToleranceValue;
149       break;
150     case GEOM::GEOM_IMeasureOperations::CheckFail:
151       anErrStruct.error = BRepCheck_CheckFail;
152       break;
153     default:
154       break;
155     }
156
157     int ii = 0;
158     int aLen = anIncrims.length();
159
160     for (; ii < aLen; ii++) {
161       anErrStruct.incriminated.push_back(anIncrims[ii]);
162     }
163
164     theErrorsTo.push_back(anErrStruct);
165   }
166 }
167
168 /**
169  * This function converts shape errors from theErrorsFrom to theErrorsTo.
170  * Note that theErrorsTo is not cleared at first.
171  *
172  * \param theErrorsFrom errors to be converted.
173  * \param theErrorsTo result errors.
174  */
175 static void ConvertShapeError
176     (const std::list<GEOMImpl_IMeasureOperations::ShapeError> &theErrorsFrom,
177            GEOM::GEOM_IMeasureOperations::ShapeErrors_out     &theErrorsTo)
178 {
179   const int aNbErr = theErrorsFrom.size();
180   GEOM::GEOM_IMeasureOperations::ShapeErrors_var anErrArray =
181     new GEOM::GEOM_IMeasureOperations::ShapeErrors;
182
183   anErrArray->length(aNbErr);
184
185   // fill the local CORBA array with values from lists
186   std::list<GEOMImpl_IMeasureOperations::ShapeError>::const_iterator
187     anIt = theErrorsFrom.begin();
188   int i = 0;
189
190   for (; anIt != theErrorsFrom.end(); i++, anIt++) {
191     GEOM::GEOM_IMeasureOperations::ShapeError_var anErrStruct =
192       new GEOM::GEOM_IMeasureOperations::ShapeError;
193
194     switch (anIt->error) {
195     case BRepCheck_InvalidPointOnCurve:
196       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidPointOnCurve;
197       break;
198     case BRepCheck_InvalidPointOnCurveOnSurface:
199       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidPointOnCurveOnSurface;
200       break;
201     case BRepCheck_InvalidPointOnSurface:
202       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidPointOnSurface;
203       break;
204     case BRepCheck_No3DCurve:
205       anErrStruct->error = GEOM::GEOM_IMeasureOperations::No3DCurve;
206       break;
207     case BRepCheck_Multiple3DCurve:
208       anErrStruct->error = GEOM::GEOM_IMeasureOperations::Multiple3DCurve;
209       break;
210     case BRepCheck_Invalid3DCurve:
211       anErrStruct->error = GEOM::GEOM_IMeasureOperations::Invalid3DCurve;
212       break;
213     case BRepCheck_NoCurveOnSurface:
214       anErrStruct->error = GEOM::GEOM_IMeasureOperations::NoCurveOnSurface;
215       break;
216     case BRepCheck_InvalidCurveOnSurface:
217       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidCurveOnSurface;
218       break;
219     case BRepCheck_InvalidCurveOnClosedSurface:
220       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidCurveOnClosedSurface;
221       break;
222     case BRepCheck_InvalidSameRangeFlag:
223       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidSameRangeFlag;
224       break;
225     case BRepCheck_InvalidSameParameterFlag:
226       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidSameParameterFlag;
227       break;
228     case BRepCheck_InvalidDegeneratedFlag:
229       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidDegeneratedFlag;
230       break;
231     case BRepCheck_FreeEdge:
232       anErrStruct->error = GEOM::GEOM_IMeasureOperations::FreeEdge;
233       break;
234     case BRepCheck_InvalidMultiConnexity:
235       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidMultiConnexity;
236       break;
237     case BRepCheck_InvalidRange:
238       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidRange;
239       break;
240     case BRepCheck_EmptyWire:
241       anErrStruct->error = GEOM::GEOM_IMeasureOperations::EmptyWire;
242       break;
243     case BRepCheck_RedundantEdge:
244       anErrStruct->error = GEOM::GEOM_IMeasureOperations::RedundantEdge;
245       break;
246     case BRepCheck_SelfIntersectingWire:
247       anErrStruct->error = GEOM::GEOM_IMeasureOperations::SelfIntersectingWire;
248       break;
249     case BRepCheck_NoSurface:
250       anErrStruct->error = GEOM::GEOM_IMeasureOperations::NoSurface;
251       break;
252     case BRepCheck_InvalidWire:
253       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidWire;
254       break;
255     case BRepCheck_RedundantWire:
256       anErrStruct->error = GEOM::GEOM_IMeasureOperations::RedundantWire;
257       break;
258     case BRepCheck_IntersectingWires:
259       anErrStruct->error = GEOM::GEOM_IMeasureOperations::IntersectingWires;
260       break;
261     case BRepCheck_InvalidImbricationOfWires:
262       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidImbricationOfWires;
263       break;
264     case BRepCheck_EmptyShell:
265       anErrStruct->error = GEOM::GEOM_IMeasureOperations::EmptyShell;
266       break;
267     case BRepCheck_RedundantFace:
268       anErrStruct->error = GEOM::GEOM_IMeasureOperations::RedundantFace;
269       break;
270     case BRepCheck_UnorientableShape:
271       anErrStruct->error = GEOM::GEOM_IMeasureOperations::UnorientableShape;
272       break;
273     case BRepCheck_NotClosed:
274       anErrStruct->error = GEOM::GEOM_IMeasureOperations::NotClosed;
275       break;
276     case BRepCheck_NotConnected:
277       anErrStruct->error = GEOM::GEOM_IMeasureOperations::NotConnected;
278       break;
279     case BRepCheck_SubshapeNotInShape:
280       anErrStruct->error = GEOM::GEOM_IMeasureOperations::SubshapeNotInShape;
281       break;
282     case BRepCheck_BadOrientation:
283       anErrStruct->error = GEOM::GEOM_IMeasureOperations::BadOrientation;
284       break;
285     case BRepCheck_BadOrientationOfSubshape:
286       anErrStruct->error = GEOM::GEOM_IMeasureOperations::BadOrientationOfSubshape;
287       break;
288     case BRepCheck_InvalidToleranceValue:
289       anErrStruct->error = GEOM::GEOM_IMeasureOperations::InvalidToleranceValue;
290       break;
291     case BRepCheck_CheckFail:
292       anErrStruct->error = GEOM::GEOM_IMeasureOperations::CheckFail;
293       break;
294     default:
295       break;
296     }
297
298     std::list<int> aIndList = anIt->incriminated;
299     GEOM::ListOfLong_var anIncrims = new GEOM::ListOfLong();
300     anIncrims->length(aIndList.size());
301
302     std::list<int>::iterator anIndIt = aIndList.begin();
303     int jj = 0;
304
305     for (; anIndIt != aIndList.end(); jj++, anIndIt++) {
306       anIncrims[jj] = *anIndIt;
307     }
308
309     anErrStruct->incriminated = anIncrims;
310
311     anErrArray[i] = anErrStruct;
312   }
313
314   // initialize out-parameter with local array
315   theErrorsTo = anErrArray._retn();
316 }
317
318 //=============================================================================
319 /*!
320  *   constructor:
321  */
322 //=============================================================================
323 GEOM_IMeasureOperations_i::GEOM_IMeasureOperations_i (PortableServer::POA_ptr thePOA,
324                                                       GEOM::GEOM_Gen_ptr theEngine,
325                                                       ::GEOMImpl_IMeasureOperations* theImpl)
326 :GEOM_IOperations_i(thePOA, theEngine, theImpl)
327 {
328   MESSAGE("GEOM_IMeasureOperations_i::GEOM_IMeasureOperations_i");
329 }
330
331 //=============================================================================
332 /*!
333  *  destructor
334  */
335 //=============================================================================
336 GEOM_IMeasureOperations_i::~GEOM_IMeasureOperations_i()
337 {
338   MESSAGE("GEOM_IMeasureOperations_i::~GEOM_IMeasureOperations_i");
339 }
340
341 //=============================================================================
342 /*!
343  *  KindOfShape
344  */
345 //=============================================================================
346 GEOM::GEOM_IKindOfShape::shape_kind GEOM_IMeasureOperations_i::KindOfShape
347                                    (GEOM::GEOM_Object_ptr  theShape,
348                                     GEOM::ListOfLong_out   theIntegers,
349                                     GEOM::ListOfDouble_out theDoubles)
350 {
351   GEOMImpl_IMeasureOperations::ShapeKind aKind = GEOMImpl_IMeasureOperations::SK_NO_SHAPE;
352
353   // allocate the CORBA arrays
354   GEOM::ListOfLong_var anIntegersArray = new GEOM::ListOfLong();
355   GEOM::ListOfDouble_var aDoublesArray = new GEOM::ListOfDouble();
356
357   //Get the reference shape
358   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
359
360   if (!aShape.IsNull()) {
361     Handle(TColStd_HSequenceOfInteger) anIntegers = new TColStd_HSequenceOfInteger;
362     Handle(TColStd_HSequenceOfReal)    aDoubles   = new TColStd_HSequenceOfReal;
363
364     // Detect kind of shape and parameters
365     aKind = GetOperations()->KindOfShape(aShape, anIntegers, aDoubles);
366
367     int nbInts = anIntegers->Length();
368     int nbDbls = aDoubles->Length();
369
370     anIntegersArray->length(nbInts);
371     aDoublesArray->length(nbDbls);
372
373     for (int ii = 0; ii < nbInts; ii++) {
374       anIntegersArray[ii] = anIntegers->Value(ii + 1);
375     }
376     for (int id = 0; id < nbDbls; id++) {
377       aDoublesArray[id] = aDoubles->Value(id + 1);
378     }
379   }
380
381   // initialize out-parameters with local arrays
382   theIntegers = anIntegersArray._retn();
383   theDoubles  = aDoublesArray._retn();
384   return (GEOM::GEOM_IKindOfShape::shape_kind)aKind;
385 }
386
387 //=============================================================================
388 /*!
389  *  GetPosition
390  */
391 //=============================================================================
392 void GEOM_IMeasureOperations_i::GetPosition
393                  (GEOM::GEOM_Object_ptr theShape,
394                   CORBA::Double& Ox, CORBA::Double& Oy, CORBA::Double& Oz,
395                   CORBA::Double& Zx, CORBA::Double& Zy, CORBA::Double& Zz,
396                   CORBA::Double& Xx, CORBA::Double& Xy, CORBA::Double& Xz)
397 {
398   //Set a not done flag
399   GetOperations()->SetNotDone();
400
401   //Set default values: global CS
402   Ox = Oy = Oz = Zx = Zy = Xy = Xz = 0.;
403   Zz = Xx = 1.;
404
405   //Get the reference shape
406   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
407   if (aShape.IsNull()) return;
408
409   // Get shape parameters
410   GetOperations()->GetPosition(aShape, Ox,Oy,Oz, Zx,Zy,Zz, Xx,Xy,Xz);
411 }
412
413 //=============================================================================
414 /*!
415  *  GetCentreOfMass
416  */
417 //=============================================================================
418 GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::GetCentreOfMass
419                                               (GEOM::GEOM_Object_ptr theShape)
420 {
421   GEOM::GEOM_Object_var aGEOMObject;
422
423   //Set a not done flag
424   GetOperations()->SetNotDone();
425
426   //Get the reference shape
427   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
428   if (aShape.IsNull()) return aGEOMObject._retn();
429
430   // Make Point - centre of mass of theShape
431   Handle(GEOM_Object) anObject = GetOperations()->GetCentreOfMass(aShape);
432   if (!GetOperations()->IsDone() || anObject.IsNull())
433     return aGEOMObject._retn();
434
435   return GetObject(anObject);
436 }
437
438 //=============================================================================
439 /*!
440  *  GetVertexByIndex
441  */
442 //=============================================================================
443 GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::GetVertexByIndex
444   (GEOM::GEOM_Object_ptr theShape, CORBA::Long theIndex)
445 {
446   GEOM::GEOM_Object_var aGEOMObject;
447
448   //Set a not done flag
449   GetOperations()->SetNotDone();
450
451   //Get the reference shape
452   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
453   if ( aShape.IsNull() ) return aGEOMObject._retn();
454
455   // Get vertex by index
456   Handle(GEOM_Object) anObject = GetOperations()->GetVertexByIndex(aShape, theIndex);
457   if (!GetOperations()->IsDone() || anObject.IsNull())
458     return aGEOMObject._retn();
459
460   return GetObject(anObject);
461 }
462
463 //=============================================================================
464 /*!
465  *  GetNormal
466  */
467 //=============================================================================
468 GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::GetNormal
469                                        (GEOM::GEOM_Object_ptr theFace,
470                                         GEOM::GEOM_Object_ptr theOptionalPoint)
471 {
472   GEOM::GEOM_Object_var aGEOMObject;
473
474   //Set a not done flag
475   GetOperations()->SetNotDone();
476
477   //Get the reference shape
478   Handle(GEOM_Object) aFace = GetObjectImpl(theFace);
479   if (aFace.IsNull()) return aGEOMObject._retn();
480
481   // Get the OptionalPoint (can be not defined)
482   Handle(GEOM_Object) anOptionalPoint = GetObjectImpl(theOptionalPoint);
483
484   // Make Vector - normal to theFace (in point theOptionalPoint if the face is not planar)
485   Handle(GEOM_Object) anObject = GetOperations()->GetNormal(aFace, anOptionalPoint);
486   if (!GetOperations()->IsDone() || anObject.IsNull())
487     return aGEOMObject._retn();
488
489   return GetObject(anObject);
490 }
491
492 //=============================================================================
493 /*!
494  *  GetBasicProperties
495  */
496 //=============================================================================
497 void GEOM_IMeasureOperations_i::GetBasicProperties (GEOM::GEOM_Object_ptr theShape,
498                                                     CORBA::Double& theLength,
499                                                     CORBA::Double& theSurfArea,
500                                                     CORBA::Double& theVolume)
501 {
502   //Set a not done flag
503   GetOperations()->SetNotDone();
504
505   //Get the reference shape
506   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
507   if (aShape.IsNull()) return;
508
509   // Get shape parameters
510   GetOperations()->GetBasicProperties(aShape, theLength, theSurfArea, theVolume);
511 }
512
513 //=============================================================================
514 /*!
515  *  GetInertia
516  */
517 //=============================================================================
518 void GEOM_IMeasureOperations_i::GetInertia
519   (GEOM::GEOM_Object_ptr theShape,
520    CORBA::Double& I11, CORBA::Double& I12, CORBA::Double& I13,
521    CORBA::Double& I21, CORBA::Double& I22, CORBA::Double& I23,
522    CORBA::Double& I31, CORBA::Double& I32, CORBA::Double& I33,
523    CORBA::Double& Ix , CORBA::Double& Iy , CORBA::Double& Iz)
524 {
525   //Set a not done flag
526   GetOperations()->SetNotDone();
527
528   //Get the reference shape
529   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
530   if (aShape.IsNull()) return;
531
532   // Get shape parameters
533   GetOperations()->GetInertia(aShape,
534                               I11, I12, I13,
535                               I21, I22, I23,
536                               I31, I32, I33,
537                               Ix , Iy , Iz);
538 }
539
540 //=============================================================================
541 /*!
542  *  GetBoundingBox
543  */
544 //=============================================================================
545 void GEOM_IMeasureOperations_i::GetBoundingBox (GEOM::GEOM_Object_ptr theShape,
546                                                 CORBA::Boolean precise,
547                                                 CORBA::Double& Xmin, CORBA::Double& Xmax,
548                                                 CORBA::Double& Ymin, CORBA::Double& Ymax,
549                                                 CORBA::Double& Zmin, CORBA::Double& Zmax)
550 {
551   //Set a not done flag
552   GetOperations()->SetNotDone();
553
554   //Get the reference shape
555   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
556   if (aShape.IsNull()) return;
557
558   // Get shape parameters
559   GetOperations()->GetBoundingBox
560     (aShape, precise, Xmin, Xmax, Ymin, Ymax, Zmin, Zmax);
561 }
562
563 //=============================================================================
564 /*!
565  *  MakeBoundingBox
566  */
567 //=============================================================================
568 GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::MakeBoundingBox
569                                               (GEOM::GEOM_Object_ptr theShape,
570                                                CORBA::Boolean precise)
571 {
572   GEOM::GEOM_Object_var aGEOMObject;
573
574   //Set a not done flag
575   GetOperations()->SetNotDone();
576
577   //Get the reference shape
578   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
579   if (aShape.IsNull()) return aGEOMObject._retn();
580
581   // Make Box - bounding box of theShape
582   Handle(GEOM_Object) anObject =
583     GetOperations()->GetBoundingBox(aShape, precise);
584   if (!GetOperations()->IsDone() || anObject.IsNull())
585     return aGEOMObject._retn();
586
587   return GetObject(anObject);
588 }
589
590 //=============================================================================
591 /*!
592  *  GetTolerance
593  */
594 //=============================================================================
595 void GEOM_IMeasureOperations_i::GetTolerance
596                                 (GEOM::GEOM_Object_ptr theShape,
597                                  CORBA::Double& FaceMin, CORBA::Double& FaceMax,
598                                  CORBA::Double& EdgeMin, CORBA::Double& EdgeMax,
599                                  CORBA::Double& VertMin, CORBA::Double& VertMax)
600 {
601   //Set a not done flag
602   GetOperations()->SetNotDone();
603
604   //Get the reference shape
605   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
606   if (aShape.IsNull()) return;
607
608   // Get shape parameters
609   GetOperations()->GetTolerance(aShape,
610                                 FaceMin, FaceMax,
611                                 EdgeMin, EdgeMax,
612                                 VertMin, VertMax);
613 }
614
615 //=============================================================================
616 /*!
617  *  CheckShape
618  */
619 //=============================================================================
620 CORBA::Boolean GEOM_IMeasureOperations_i::CheckShape
621              (GEOM::GEOM_Object_ptr                          theShape,
622               GEOM::GEOM_IMeasureOperations::ShapeErrors_out theErrors)
623 {
624   //Set a not done flag
625   GetOperations()->SetNotDone();
626
627   if (CORBA::is_nil(theShape))
628   {
629     return 0;
630   }
631
632   //Get the reference shape
633   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
634
635   if (aShape.IsNull())
636   {
637     return 0;
638   }
639
640   std::list<GEOMImpl_IMeasureOperations::ShapeError> anErrList;
641   bool isOk = GetOperations()->CheckShape(aShape, false, anErrList);
642
643   ConvertShapeError(anErrList, theErrors);
644
645   return isOk;
646 }
647
648 CORBA::Boolean GEOM_IMeasureOperations_i::CheckShapeWithGeometry
649              (GEOM::GEOM_Object_ptr                          theShape,
650               GEOM::GEOM_IMeasureOperations::ShapeErrors_out theErrors)
651 {
652   //Set a not done flag
653   GetOperations()->SetNotDone();
654
655   if (CORBA::is_nil(theShape))
656   {
657     return 0;
658   }
659
660   //Get the reference shape
661   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
662
663   if (aShape.IsNull())
664   {
665     return 0;
666   }
667
668   std::list<GEOMImpl_IMeasureOperations::ShapeError> anErrList;
669   bool isOk = GetOperations()->CheckShape(aShape, true, anErrList);
670
671   ConvertShapeError(anErrList, theErrors);
672
673   return isOk;
674 }
675
676 //=============================================================================
677 /*!
678  *  PrintShapeErrors
679  */
680 //=============================================================================
681 char* GEOM_IMeasureOperations_i::PrintShapeErrors
682              (      GEOM::GEOM_Object_ptr                       theShape,
683               const GEOM::GEOM_IMeasureOperations::ShapeErrors &theErrors)
684 {
685   //Get the reference shape
686   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
687
688   if (aShape.IsNull()) {
689     return NULL;
690   }
691
692   // Convert the errors sequence
693   std::list<GEOMImpl_IMeasureOperations::ShapeError> anErrList;
694
695   ConvertShapeError(theErrors, anErrList);
696
697   TCollection_AsciiString aDescr =
698     GetOperations()->PrintShapeErrors(aShape, anErrList);
699
700   return CORBA::string_dup(aDescr.ToCString());
701 }
702
703 //=============================================================================
704 /*!
705  *  CheckSelfIntersections
706  */
707 //=============================================================================
708 CORBA::Boolean GEOM_IMeasureOperations_i::CheckSelfIntersections (GEOM::GEOM_Object_ptr theShape,
709                                                                   CORBA::Long           theCheckLevel,
710                                                                   GEOM::ListOfLong_out  theIntersections)
711 {
712   // Set a not done flag
713   GetOperations()->SetNotDone();
714
715   bool isGood = false;
716
717   // Allocate the CORBA arrays
718   GEOM::ListOfLong_var anIntegersArray = new GEOM::ListOfLong();
719
720   // Get the reference shape
721   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
722
723   if (!aShape.IsNull()) {
724     GEOMImpl_IMeasureOperations::SICheckLevel aCheckLevel;
725
726     switch(theCheckLevel) {
727     case GEOM::SI_V_V:
728       aCheckLevel = GEOMImpl_IMeasureOperations::SI_V_V;
729       break;
730     case GEOM::SI_V_E:
731       aCheckLevel = GEOMImpl_IMeasureOperations::SI_V_E;
732       break;
733     case GEOM::SI_E_E:
734       aCheckLevel = GEOMImpl_IMeasureOperations::SI_E_E;
735       break;
736     case GEOM::SI_V_F:
737       aCheckLevel = GEOMImpl_IMeasureOperations::SI_V_F;
738       break;
739     case GEOM::SI_E_F:
740       aCheckLevel = GEOMImpl_IMeasureOperations::SI_E_F;
741       break;
742     case GEOM::SI_ALL:
743     default:
744       aCheckLevel = GEOMImpl_IMeasureOperations::SI_ALL;
745       break;
746     }
747
748     Handle(TColStd_HSequenceOfInteger) anIntegers = new TColStd_HSequenceOfInteger;
749
750     // Detect self-intersections
751     isGood = GetOperations()->CheckSelfIntersections
752       (aShape, aCheckLevel, anIntegers);
753
754     int nbInts = anIntegers->Length();
755
756     anIntegersArray->length(nbInts);
757
758     for (int ii = 0; ii < nbInts; ii++) {
759       anIntegersArray[ii] = anIntegers->Value(ii + 1);
760     }
761   }
762
763   // Initialize out-parameters with local arrays
764   theIntersections = anIntegersArray._retn();
765   return isGood;
766 }
767
768 //=============================================================================
769 /*!
770  *  FastIntersect
771  */
772 //=============================================================================
773 CORBA::Boolean GEOM_IMeasureOperations_i::FastIntersect (GEOM::GEOM_Object_ptr theShape1,
774                                                          GEOM::GEOM_Object_ptr theShape2,
775                                                          CORBA::Double         theTolerance,
776                                                          CORBA::Float          theDeflection,
777                                                          GEOM::ListOfLong_out  theIntersections1,
778                                                          GEOM::ListOfLong_out  theIntersections2)
779 {
780   // Set a not done flag
781   GetOperations()->SetNotDone();
782
783   bool isGood = false;
784
785   // Allocate the CORBA arrays
786   GEOM::ListOfLong_var anIntegersArray1 = new GEOM::ListOfLong();
787   GEOM::ListOfLong_var anIntegersArray2 = new GEOM::ListOfLong();
788
789   // Get the reference shape
790   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
791   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
792   
793   if (!aShape1.IsNull() && !aShape2.IsNull()) {
794     Handle(TColStd_HSequenceOfInteger) anIntegers1 = new TColStd_HSequenceOfInteger;
795     Handle(TColStd_HSequenceOfInteger) anIntegers2 = new TColStd_HSequenceOfInteger;
796
797     // Detect intersections
798     isGood = GetOperations()->FastIntersect
799       (aShape1, aShape2, theTolerance, theDeflection, anIntegers1, anIntegers2);
800
801     int nbInts1 = anIntegers1->Length();
802     int nbInts2 = anIntegers2->Length();
803
804     anIntegersArray1->length(nbInts1);
805     anIntegersArray2->length(nbInts2);
806
807     for (int ii = 0; ii < nbInts1; ii++) {
808       anIntegersArray1[ii] = anIntegers1->Value(ii + 1);
809     }
810     for (int ii = 0; ii < nbInts2; ii++) {
811       anIntegersArray2[ii] = anIntegers2->Value(ii + 1);
812     }
813   }
814
815   // Initialize out-parameters with local arrays
816   theIntersections1 = anIntegersArray1._retn();
817   theIntersections2 = anIntegersArray2._retn();
818   return isGood;
819 }
820
821 //=============================================================================
822 /*!
823  *  IsGoodForSolid
824  */
825 //=============================================================================
826 char* GEOM_IMeasureOperations_i::IsGoodForSolid (GEOM::GEOM_Object_ptr theShape)
827 {
828   //Set a not done flag
829   GetOperations()->SetNotDone();
830
831   //Get the reference shape
832   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
833   if (aShape.IsNull()) return CORBA::string_dup("WRN_NULL_OBJECT_OR_SHAPE");
834
835   // Get shape parameters
836   TCollection_AsciiString aDescription = GetOperations()->IsGoodForSolid(aShape);
837   return CORBA::string_dup(aDescription.ToCString());
838 }
839
840 //=============================================================================
841 /*!
842  *  WhatIs
843  */
844 //=============================================================================
845 char* GEOM_IMeasureOperations_i::WhatIs (GEOM::GEOM_Object_ptr theShape)
846 {
847   //Set a not done flag
848   GetOperations()->SetNotDone();
849
850   //Get the reference shape
851   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
852   if (aShape.IsNull()) return NULL;
853
854   // Get shape parameters
855   TCollection_AsciiString aDescription = GetOperations()->WhatIs(aShape);
856   return CORBA::string_dup(aDescription.ToCString());
857 }
858
859 //=============================================================================
860 /*!
861  *  AreCoordsInside
862  */
863 //=============================================================================
864 GEOM::ListOfBool* GEOM_IMeasureOperations_i::AreCoordsInside (GEOM::GEOM_Object_ptr theShape,
865                                                               const GEOM::ListOfDouble& theCoords,
866                                                               CORBA::Double tolerance)
867 {
868   //Set a not done flag
869   GetOperations()->SetNotDone();
870   
871   unsigned int nb_points = theCoords.length()/3;
872
873   GEOM::ListOfBool_var aResults = new GEOM::ListOfBool;
874   aResults->length(nb_points);
875   
876   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
877   
878   std::vector<double> tmp(3*nb_points);
879   for (int i = 0; i < 3*nb_points; i++)
880     tmp[i] = theCoords[i];
881   std::vector<bool> res = GetOperations()->AreCoordsInside(aShape, tmp, tolerance);
882   for (int i = 0; i < nb_points; i++)
883     aResults[i] = i < res.size() ? res[i] : false;
884   return aResults._retn();
885 }
886
887 //=============================================================================
888 /*!
889  *  GetMinDistance
890  */
891 //=============================================================================
892 CORBA::Double GEOM_IMeasureOperations_i::GetMinDistance
893   (GEOM::GEOM_Object_ptr theShape1, GEOM::GEOM_Object_ptr theShape2,
894    CORBA::Double& X1, CORBA::Double& Y1, CORBA::Double& Z1,
895    CORBA::Double& X2, CORBA::Double& Y2, CORBA::Double& Z2)
896 {
897   //Set a not done flag
898   GetOperations()->SetNotDone();
899
900   //Get the reference shape
901   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
902   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
903   if (aShape1.IsNull() || aShape2.IsNull()) return -1.0;
904
905   // Get shape parameters
906   return GetOperations()->GetMinDistance(aShape1, aShape2, X1, Y1, Z1, X2, Y2, Z2);
907 }
908
909 //=============================================================================
910 /*!
911  *  ClosestPoints
912  */
913 //=============================================================================
914 CORBA::Long GEOM_IMeasureOperations_i::ClosestPoints
915   (GEOM::GEOM_Object_ptr theShape1, GEOM::GEOM_Object_ptr theShape2,
916    GEOM::ListOfDouble_out theCoords)
917 {
918   //Set a not done flag
919   GetOperations()->SetNotDone();
920
921   // allocate the CORBA array
922   int nbSols = 0;
923   GEOM::ListOfDouble_var aDoublesArray = new GEOM::ListOfDouble();
924
925   //Get the reference shape
926   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
927   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
928
929   if (!aShape1.IsNull() && !aShape2.IsNull()) {
930     Handle(TColStd_HSequenceOfReal) aDoubles = new TColStd_HSequenceOfReal;
931     // Get shape parameters
932     nbSols = GetOperations()->ClosestPoints(aShape1, aShape2, aDoubles);
933     int nbDbls = aDoubles->Length();
934     aDoublesArray->length(nbDbls);
935     for (int id = 0; id < nbDbls; id++) {
936       aDoublesArray[id] = aDoubles->Value(id + 1);
937     }
938   }
939
940   theCoords = aDoublesArray._retn();
941   return nbSols;
942 }
943
944 //=============================================================================
945 /*!
946  *  PointCoordinates
947  */
948 //=============================================================================
949 void GEOM_IMeasureOperations_i::PointCoordinates (GEOM::GEOM_Object_ptr theShape,
950                                                   CORBA::Double& X, CORBA::Double& Y, CORBA::Double& Z)
951
952 {
953   //Set a not done flag
954   GetOperations()->SetNotDone();
955
956   //Get the reference shape
957   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
958   if (aShape.IsNull())
959     return;
960
961   // Get shape parameters
962   GetOperations()->PointCoordinates( aShape, X, Y, Z );
963 }
964
965 //=============================================================================
966 /*!
967  *  GetAngle
968  */
969 //=============================================================================
970 CORBA::Double GEOM_IMeasureOperations_i::GetAngle (GEOM::GEOM_Object_ptr theShape1,
971                                                    GEOM::GEOM_Object_ptr theShape2)
972 {
973   //Set a not done flag
974   GetOperations()->SetNotDone();
975
976   //Get the reference shapes
977   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
978   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
979   if (aShape1.IsNull() || aShape2.IsNull()) return -1.0;
980
981   // Get the angle
982   return GetOperations()->GetAngle(aShape1, aShape2);
983 }
984
985 //=============================================================================
986 /*!
987  *  GetAngle
988  */
989 //=============================================================================
990 CORBA::Double GEOM_IMeasureOperations_i::GetAngleBtwVectors (GEOM::GEOM_Object_ptr theShape1,
991                                                              GEOM::GEOM_Object_ptr theShape2)
992 {
993   //Set a not done flag
994   GetOperations()->SetNotDone();
995
996   //Get the reference shapes
997   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
998   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
999   if (aShape1.IsNull() || aShape2.IsNull()) return -1.0;
1000
1001   // Get the angle
1002   return GetOperations()->GetAngleBtwVectors(aShape1, aShape2);
1003 }
1004
1005
1006 //=============================================================================
1007 /*!
1008  *  CurveCurvatureByParam
1009  */
1010 //=============================================================================
1011 CORBA::Double GEOM_IMeasureOperations_i::CurveCurvatureByParam
1012                        (GEOM::GEOM_Object_ptr theCurve, CORBA::Double theParam)
1013 {
1014   //Set a not done flag
1015   GetOperations()->SetNotDone();
1016
1017   //Get the reference shape
1018   Handle(GEOM_Object) aShape = GetObjectImpl(theCurve);
1019   if(aShape.IsNull()) return -1.0;
1020
1021   return GetOperations()->CurveCurvatureByParam(aShape,theParam);
1022 }
1023
1024 //=============================================================================
1025 /*!
1026  *  CurveCurvatureByPoint
1027  */
1028 //=============================================================================
1029 CORBA::Double GEOM_IMeasureOperations_i::CurveCurvatureByPoint
1030                (GEOM::GEOM_Object_ptr theCurve, GEOM::GEOM_Object_ptr thePoint)
1031 {
1032   //Set a not done flag
1033   GetOperations()->SetNotDone();
1034
1035   //Get the reference shape
1036   Handle(GEOM_Object) aShape = GetObjectImpl(theCurve);
1037   Handle(GEOM_Object) aPoint = GetObjectImpl(thePoint);
1038   if( aShape.IsNull() || aPoint.IsNull() ) return -1.0;
1039
1040   return GetOperations()->CurveCurvatureByPoint(aShape,aPoint);
1041 }
1042
1043
1044 //=============================================================================
1045 /*!
1046  *  MaxSurfaceCurvatureByParam
1047  */
1048 //=============================================================================
1049 CORBA::Double GEOM_IMeasureOperations_i::MaxSurfaceCurvatureByParam
1050                                                 (GEOM::GEOM_Object_ptr theSurf,
1051                                                  CORBA::Double theUParam,
1052                                                  CORBA::Double theVParam)
1053 {
1054   //Set a not done flag
1055   GetOperations()->SetNotDone();
1056
1057   //Get the reference shape
1058   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
1059   if(aShape.IsNull()) return -1.0;
1060
1061   return GetOperations()->MaxSurfaceCurvatureByParam(aShape,theUParam,theVParam);
1062 }
1063
1064 //=============================================================================
1065 /*!
1066  *  MaxSurfaceCurvatureByPoint
1067  */
1068 //=============================================================================
1069 CORBA::Double GEOM_IMeasureOperations_i::MaxSurfaceCurvatureByPoint
1070                 (GEOM::GEOM_Object_ptr theSurf, GEOM::GEOM_Object_ptr thePoint)
1071 {
1072   //Set a not done flag
1073   GetOperations()->SetNotDone();
1074
1075   //Get the reference shape
1076   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
1077   Handle(GEOM_Object) aPoint = GetObjectImpl(thePoint);
1078   if( aShape.IsNull() || aPoint.IsNull() ) return -1.0;
1079
1080   return GetOperations()->MaxSurfaceCurvatureByPoint(aShape,aPoint);
1081 }
1082
1083 //=============================================================================
1084 /*!
1085  *  MinSurfaceCurvatureByParam
1086  */
1087 //=============================================================================
1088 CORBA::Double GEOM_IMeasureOperations_i::MinSurfaceCurvatureByParam
1089                                                 (GEOM::GEOM_Object_ptr theSurf,
1090                                                  CORBA::Double theUParam,
1091                                                  CORBA::Double theVParam)
1092 {
1093   //Set a not done flag
1094   GetOperations()->SetNotDone();
1095
1096   //Get the reference shape
1097   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
1098   if (aShape.IsNull()) return -1.0;
1099
1100   return GetOperations()->MinSurfaceCurvatureByParam(aShape,theUParam,theVParam);
1101 }
1102
1103 //=============================================================================
1104 /*!
1105  *  MinSurfaceCurvatureByPoint
1106  */
1107 //=============================================================================
1108 CORBA::Double GEOM_IMeasureOperations_i::MinSurfaceCurvatureByPoint
1109                 (GEOM::GEOM_Object_ptr theSurf, GEOM::GEOM_Object_ptr thePoint)
1110 {
1111   //Set a not done flag
1112   GetOperations()->SetNotDone();
1113
1114   //Get the reference shape
1115   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
1116   Handle(GEOM_Object) aPoint = GetObjectImpl(thePoint);
1117   if (aShape.IsNull() || aPoint.IsNull()) return -1.0;
1118
1119   return GetOperations()->MinSurfaceCurvatureByPoint(aShape,aPoint);
1120 }