//
#include <Standard_Stream.hxx>
+#include <TColStd_HArray1OfReal.hxx>
#include "GEOM_IMeasureOperations_i.hh"
return GetOperations()->GetAngleBtwVectors(aShape1, aShape2);
}
+//=============================================================================
+/*!
+* PatchFace
+*/
+//=============================================================================
+GEOM::ListOfGO* GEOM_IMeasureOperations_i::PatchFace(GEOM::GEOM_Object_ptr theShape)
+{
+ GEOM::ListOfGO_var aSeq = new GEOM::ListOfGO;
+
+ //Get the reference shape
+ Handle(::GEOM_Object) aShapeRef = GetObjectImpl(theShape);
+ if (aShapeRef.IsNull()) return aSeq._retn();
+
+ // Perform patch face operation
+ Handle(TColStd_HSequenceOfTransient) aHSeq =
+ GetOperations()->PatchFace(aShapeRef);
+ if (!GetOperations()->IsDone() || aHSeq.IsNull())
+ return aSeq._retn();
+
+ Standard_Integer aLength = aHSeq->Length();
+ aSeq->length(aLength);
+ for (Standard_Integer i = 1; i <= aLength; i++)
+ aSeq[i - 1] = GetObject(Handle(::GEOM_Object)::DownCast(aHSeq->Value(i)));
+
+ return aSeq._retn();
+}
//=============================================================================
/*!
return GetOperations()->MinSurfaceCurvatureByPoint(aShape,aPoint);
}
+
+//=============================================================================
+/*!
+ * SurfaceCurvatureByPointAndDirection
+ */
+//=============================================================================
+GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::SurfaceCurvatureByPointAndDirection
+ (GEOM::GEOM_Object_ptr theSurf,
+ GEOM::GEOM_Object_ptr thePoint,
+ GEOM::GEOM_Object_ptr theDirection)
+{
+ GEOM::GEOM_Object_var aGEOMObject;
+
+ //Set a not done flag
+ GetOperations()->SetNotDone();
+
+ //Get the reference shape
+ Handle(::GEOM_Object) aShape = GetObjectImpl(theSurf);
+ Handle(::GEOM_Object) aPoint = GetObjectImpl(thePoint);
+ Handle(::GEOM_Object) aDirection = GetObjectImpl(theDirection);
+ if (aShape.IsNull() || aPoint.IsNull() || aDirection.IsNull()) return aGEOMObject._retn();
+
+ Handle(::GEOM_Object) anObject =
+ GetOperations()->SurfaceCurvatureByPointAndDirection(aShape,aPoint,aDirection);
+ if (!GetOperations()->IsDone() || anObject.IsNull())
+ return aGEOMObject._retn();
+
+ 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;
+}