Salome HOME
Implemented [bos #35094] [EDF] (2023-T1) X,Y,Z to U,V.
[modules/geom.git] / src / GEOM_I / GEOM_IMeasureOperations_i.cc
index 0c09e09c4b1ebaccdf503e5e29c0762917fdee19..a740731b4b648c22badcb4e0f6115774649c1274 100644 (file)
@@ -21,6 +21,7 @@
 //
 
 #include <Standard_Stream.hxx>
+#include <TColStd_HArray1OfReal.hxx>
 
 #include "GEOM_IMeasureOperations_i.hh"
 
@@ -1243,3 +1244,383 @@ GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::SurfaceCurvatureByPointAndDirec
 
   return GetObject(anObject);
 }
+
+//=============================================================================
+/*!
+ *  XYZtoUV
+ */
+//=============================================================================
+GEOM::ListOfDouble* GEOM_IMeasureOperations_i::XYZtoUV
+                                        (GEOM::GEOM_Object_ptr     theSurf,
+                                         const GEOM::ListOfDouble& theXYZlist,
+                                         CORBA::Boolean            theIsNormalized)
+{
+  GEOM::ListOfDouble_var resUV = new GEOM::ListOfDouble;
+
+  //Set a not done flag
+  GetOperations()->SetNotDone();
+
+  //Get the reference shape
+  Handle(::GEOM_Object) aShape = GetObjectImpl(theSurf);
+  if (aShape.IsNull()) return resUV._retn();
+
+  //Get input XYZ list
+  Handle(TColStd_HArray1OfReal) aXYZlist =
+    new TColStd_HArray1OfReal (0, theXYZlist.length() - 1);
+  {
+    size_t nb = theXYZlist.length(); 
+    for (size_t i = 0; i < nb; ++i)
+      aXYZlist->SetValue(i, theXYZlist[i]);
+  }
+
+  //Call implementation
+  Handle(TColStd_HArray1OfReal) aUVlist =
+    GetOperations()->XYZtoUV(aShape, aXYZlist, theIsNormalized);
+
+  if (GetOperations()->IsDone()) {
+    resUV->length(aUVlist->Length());
+    int i0 = aUVlist->Lower();
+    int nb = aUVlist->Upper(); 
+    for (int i = i0; i <= nb; ++i)
+      resUV[ i-i0 ] = aUVlist->Value(i);
+  }
+
+  return resUV._retn();
+}
+
+//=============================================================================
+/*!
+ *  UVtoXYZ
+ */
+//=============================================================================
+GEOM::ListOfDouble* GEOM_IMeasureOperations_i::UVtoXYZ
+                                        (GEOM::GEOM_Object_ptr     theSurf,
+                                         const GEOM::ListOfDouble& theUVlist,
+                                         CORBA::Boolean            theIsNormalized)
+{
+  GEOM::ListOfDouble_var resXYZ = new GEOM::ListOfDouble;
+
+  //Set a not done flag
+  GetOperations()->SetNotDone();
+
+  //Get the reference shape
+  Handle(::GEOM_Object) aShape = GetObjectImpl(theSurf);
+  if (aShape.IsNull()) return resXYZ._retn();
+
+  //Get input UV list
+  Handle(TColStd_HArray1OfReal) aUVlist =
+    new TColStd_HArray1OfReal (0, theUVlist.length() - 1);
+  {
+    size_t nb = theUVlist.length(); 
+    for (size_t i = 0; i < nb; ++i)
+      aUVlist->SetValue(i, theUVlist[i]);
+  }
+
+  //Call implementation
+  Handle(TColStd_HArray1OfReal) aXYZlist =
+    GetOperations()->UVtoXYZ(aShape, aUVlist, theIsNormalized);
+
+  if (GetOperations()->IsDone()) {
+    resXYZ->length(aXYZlist->Length());
+    int i0 = aXYZlist->Lower();
+    int nb = aXYZlist->Upper(); 
+    for (int i = i0; i <= nb; ++i)
+      resXYZ[ i-i0 ] = aXYZlist->Value(i);
+  }
+
+  return resXYZ._retn();
+}
+
+//=============================================================================
+/*!
+ *  SelfIntersected2D
+ *  Find all self-intersected 2D curves.
+ *  \param theShape Shape for check.
+ */
+ //=============================================================================
+GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape* GEOM_IMeasureOperations_i::SelfIntersected2D(
+                                          const GEOM::GEOM_IMeasureOperations::CheckResults& theResuts)
+{
+  GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape_var aSeq =
+    new GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape;
+
+  // Perform patch face operation
+  std::list<GEOMImpl_IMeasureOperations::FailedChecks> aResults;
+  ConvertToList(theResuts, aResults);
+
+  std::list<GEOMImpl_IMeasureOperations::CoupleOfObjects> aSelfInters = GetOperations()->SelfIntersected2D(aResults);
+  if (!GetOperations()->IsDone() || aSelfInters.empty())
+    return aSeq._retn();
+
+  Standard_Integer aLength = aSelfInters.size();
+  aSeq->length(aLength);
+  std::list<GEOMImpl_IMeasureOperations::CoupleOfObjects>::iterator anIter(aSelfInters.begin());
+  for (Standard_Integer i = 0; i < aLength; i++, ++anIter)
+  {
+    aSeq[i].first = GetObject(Handle(::GEOM_Object)::DownCast((*anIter).first));
+    aSeq[i].second = GetObject(Handle(::GEOM_Object)::DownCast((*anIter).second));
+  }
+
+  return aSeq._retn();
+}
+
+//=============================================================================
+/*!
+ *  InterferingSubshapes
+ *  Find pairs of interfering sub-shapes, by default all pairs of interfering shapes are returned.
+ *  \param theShape Shape for check.
+ *  \param theShapeType1 Type of shape.
+ *  \param theShapeType2 Type of shape.
+ */
+ //=============================================================================
+GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape* GEOM_IMeasureOperations_i::InterferingSubshapes(
+    const GEOM::GEOM_IMeasureOperations::CheckResults& theResuts,
+    const CORBA::Long theShapeType1,
+    const CORBA::Long theShapeType2)
+{
+  GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape_var aSeq =
+    new GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape;
+
+  // Perform patch face operation
+  std::list<GEOMImpl_IMeasureOperations::FailedChecks> aResults;
+  ConvertToList(theResuts, aResults);
+
+  std::list<GEOMImpl_IMeasureOperations::CoupleOfObjects> aSelfInterf =
+    GetOperations()->InterferingSubshapes(aResults, theShapeType1, theShapeType2);
+  if (!GetOperations()->IsDone() || aSelfInterf.empty())
+    return aSeq._retn();
+
+  Standard_Integer aLength = aSelfInterf.size();
+  aSeq->length(aLength);
+  std::list<GEOMImpl_IMeasureOperations::CoupleOfObjects>::iterator anIter(aSelfInterf.begin());
+  for (Standard_Integer i = 0; i < aLength; i++, ++anIter)
+  {
+    aSeq[i].first = GetObject(Handle(::GEOM_Object)::DownCast((*anIter).first));
+    aSeq[i].second = GetObject(Handle(::GEOM_Object)::DownCast((*anIter).second));
+  }
+
+  return aSeq._retn();
+}
+
+//=============================================================================
+/*!
+ *  SmallEdges
+ *  Find edges, which are fully covered by tolerances of vertices.
+ *  \param theShape Shape for check.
+ */
+ //=============================================================================
+GEOM::ListOfGO* GEOM_IMeasureOperations_i::SmallEdges(const GEOM::GEOM_IMeasureOperations::CheckResults& theResuts)
+{
+  GEOM::ListOfGO_var aSeq = new GEOM::ListOfGO;
+
+  std::list<GEOMImpl_IMeasureOperations::FailedChecks> aResults;
+  ConvertToList(theResuts, aResults);
+
+  //Get the reference shape
+  Handle(TColStd_HSequenceOfTransient) aSmallEdges = GetOperations()->SmallEdges(aResults);
+  if (!GetOperations()->IsDone() || aSmallEdges.IsNull())
+    return aSeq._retn();
+
+  Standard_Integer aLength = aSmallEdges->Length();
+  aSeq->length(aLength);
+  for (Standard_Integer i = 1; i <= aLength; i++)
+    aSeq[i - 1] = GetObject(Handle(::GEOM_Object)::DownCast(aSmallEdges->Value(i)));
+
+  return aSeq._retn();
+}
+
+//=============================================================================
+/*!
+ *  DistantShapes
+ *  find remote objects (sub-shape on a shape).
+ *  \param theShape Shape for check.
+ *  \param theShapeType Type of shape.
+ *  \param theSubShapeType Type of sub-shape.
+ *  \param theTolerance tolerance.
+ */
+ //=============================================================================
+GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape* GEOM_IMeasureOperations_i::DistantShapes(
+    const GEOM::GEOM_IMeasureOperations::CheckResults& theResuts,
+    const CORBA::Long theShapeType,
+    const CORBA::Long theSubShapeType,
+    const CORBA::Double theTolerance)
+{
+  GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape_var aSeq =
+    new GEOM::GEOM_IMeasureOperations::SequenceOfPairOfShape;
+
+  std::list<GEOMImpl_IMeasureOperations::FailedChecks> aResults;
+  ConvertToList(theResuts, aResults);
+
+  // Perform patch face operation
+  std::list<GEOMImpl_IMeasureOperations::CoupleOfObjects> aDistantS =
+      GetOperations()->DistantShapes(aResults, theShapeType, theSubShapeType, theTolerance);
+  if (!GetOperations()->IsDone() || aDistantS.empty())
+    return aSeq._retn();
+
+  Standard_Integer aLength = aDistantS.size();
+  aSeq->length(aLength);
+  std::list<GEOMImpl_IMeasureOperations::CoupleOfObjects>::iterator anIter(aDistantS.begin());
+  for (Standard_Integer i = 0; i < aLength; i++, ++anIter)
+  {
+    aSeq[i].first = GetObject(Handle(::GEOM_Object)::DownCast((*anIter).first));
+    aSeq[i].second = GetObject(Handle(::GEOM_Object)::DownCast((*anIter).second));
+  }
+
+  return aSeq._retn();
+}
+
+//=============================================================================
+/*!
+ *  CheckConformityShape
+ *  Perform all required checks
+ */
+ //=============================================================================
+GEOM::GEOM_IMeasureOperations::CheckResults* GEOM_IMeasureOperations_i::CheckConformityShape(
+    GEOM::GEOM_Object_ptr theShape)
+{
+  GEOM::GEOM_IMeasureOperations::CheckResults_var aRes = new GEOM::GEOM_IMeasureOperations::CheckResults;
+
+  //Get the reference shape
+  Handle(::GEOM_Object) aShapeRef = GetObjectImpl(theShape);
+  if (aShapeRef.IsNull()) return nullptr;
+
+  std::list<GEOMImpl_IMeasureOperations::FailedChecks> aChecks;
+  GetOperations()->CheckConformityShape(aShapeRef, aChecks);
+
+  Standard_Integer aLength = aChecks.size();
+  aRes->length(aLength);
+  std::list<GEOMImpl_IMeasureOperations::FailedChecks>::const_iterator anIntIt = aChecks.cbegin();
+  for (Standard_Integer i = 0; i < aLength; i++, ++anIntIt)
+  {
+    aRes[i].type = (*anIntIt).TypeOfCheck;
+    aRes[i].failedShapes.first  = GetObject(Handle(::GEOM_Object)::DownCast((*anIntIt).FailedShapes.first));
+    aRes[i].failedShapes.second = GetObject(Handle(::GEOM_Object)::DownCast((*anIntIt).FailedShapes.second));
+  }
+
+  return aRes._retn();
+}
+
+//=============================================================================
+/*!
+ *  UpdateTolerance
+ *  Compute possible tolerance for the shape, minimize tolerance of shape as well
+ *  as tolerance of sub-shapes as much as possible
+ *  \param theShape Shape for check.
+ */
+ //=============================================================================
+CORBA::Double GEOM_IMeasureOperations_i::UpdateTolerance(GEOM::GEOM_Object_ptr theShape)
+{
+  //Get the reference shape
+  Handle(::GEOM_Object) aShapeRef = GetObjectImpl(theShape);
+  if (aShapeRef.IsNull())
+    return false;
+
+  return GetOperations()->UpdateTolerance(aShapeRef);
+}
+
+void GEOM_IMeasureOperations_i::ConvertToList(const GEOM::GEOM_IMeasureOperations::CheckResults& theResuts,
+                                              std::list<GEOMImpl_IMeasureOperations::FailedChecks>& theListOfResults)
+{
+  for (Standard_Integer i = 0; i < theResuts.length(); ++i)
+  {
+    GEOMImpl_IMeasureOperations::FailedChecks aCheck;
+    aCheck.TypeOfCheck = theResuts[i].type;
+    aCheck.FailedShapes.first = GetObjectImpl(theResuts[i].failedShapes.first);
+    aCheck.FailedShapes.second = GetObjectImpl(theResuts[i].failedShapes.second);
+    theListOfResults.push_back(aCheck);
+  }
+}
+
+//=============================================================================
+/*!
+ *  ShapeProximityCalculator
+ */
+//=============================================================================
+GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::ShapeProximityCalculator(GEOM::GEOM_Object_ptr theShape1,
+                                                                          GEOM::GEOM_Object_ptr theShape2)
+{
+  GEOM::GEOM_Object_var anEmptyCalc;
+
+  //Set a not done flag
+  GetOperations()->SetNotDone();
+
+  //Get the reference shape
+  Handle(::GEOM_Object) aShape1 = GetObjectImpl(theShape1);
+  Handle(::GEOM_Object) aShape2 = GetObjectImpl(theShape2);
+  if (aShape1.IsNull() || aShape2.IsNull())
+    return anEmptyCalc._retn();
+
+  Handle(::GEOM_Object) aCalculator = GetOperations()->ShapeProximityCalculator(aShape1, aShape2);
+  if (!GetOperations()->IsDone() || aCalculator.IsNull())
+    return anEmptyCalc._retn();
+
+  return GetObject(aCalculator);
+}
+
+//=============================================================================
+/*!
+ *  SetShapeSampling
+ */
+ //=============================================================================
+void GEOM_IMeasureOperations_i::SetShapeSampling(GEOM::GEOM_Object_ptr theCalculator,
+                                                 GEOM::GEOM_Object_ptr theShape,
+                                                 CORBA::Long theNbSamples)
+{
+  //Set a not done flag
+  GetOperations()->SetNotDone();
+
+  //Get the proximity calculator
+  Handle(::GEOM_Object) aCalc = GetObjectImpl(theCalculator);
+  if (aCalc.IsNull())
+    return ;
+  //Get the reference shape
+  Handle(::GEOM_Object) aShape = GetObjectImpl(theShape);
+  if (aShape.IsNull())
+    return ;
+
+  GetOperations()->SetShapeSampling(aCalc, aShape, theNbSamples);
+}
+
+//=============================================================================
+/*!
+ *  GetCoarseProximity
+ */
+ //=============================================================================
+CORBA::Double GEOM_IMeasureOperations_i::GetCoarseProximity(GEOM::GEOM_Object_ptr theCalculator)
+{
+  //Set a not done flag
+  GetOperations()->SetNotDone();
+
+  //Get the reference shape
+  Handle(::GEOM_Object) aCalc = GetObjectImpl(theCalculator);
+  if (aCalc.IsNull())
+    return -1.0;
+
+  Standard_Real aProximity = GetOperations()->GetCoarseProximity(aCalc);
+  if (!GetOperations()->IsDone())
+    return -1.0;
+
+  return aProximity;
+}
+
+//=============================================================================
+/*!
+ *  GetPreciseProximity
+ */
+ //=============================================================================
+CORBA::Double GEOM_IMeasureOperations_i::GetPreciseProximity(GEOM::GEOM_Object_ptr theCalculator)
+{
+  //Set a not done flag
+  GetOperations()->SetNotDone();
+
+  //Get the reference shape
+  Handle(::GEOM_Object) aCalc = GetObjectImpl(theCalculator);
+  if (aCalc.IsNull())
+    return -1.0;
+
+  Standard_Real aProximity = GetOperations()->GetPreciseProximity(aCalc);
+  if (!GetOperations()->IsDone())
+    return -1.0;
+
+  return aProximity;
+}