Salome HOME
Issue #1922: After partition, some edges are not displayed below the connected shape...
[modules/shaper.git] / src / GeomAlgoAPI / GeomAlgoAPI_Placement.cpp
index 6e28fc603c5231023467fcbed35e795957c7fbbc..58128437f08beed518f8d23b05305ca277580a6f 100644 (file)
@@ -4,10 +4,13 @@
 // Created:     2 Dec 2014
 // Author:      Artem ZHIDKOV
 
-#include <GeomAlgoAPI_Placement.h>
+#include "GeomAlgoAPI_Placement.h"
+
 #include <GeomAlgoAPI_DFLoader.h>
 
+#include <GeomAPI_Dir.h>
 #include <GeomAPI_Edge.h>
+#include <GeomAPI_Face.h>
 #include <GeomAPI_Lin.h>
 #include <GeomAPI_Pnt.h>
 #include <GeomAPI_Pln.h>
 #include <GeomAPI_XYZ.h>
 
 #include <BRepBuilderAPI_Transform.hxx>
+#include <BRepClass3d_SolidClassifier.hxx>
+#include <BRepGProp.hxx>
 #include <gp_Trsf.hxx>
 #include <gp_Quaternion.hxx>
-#include <TopExp_Explorer.hxx>
-#include <BRepCheck_Analyzer.hxx>
-#include <BRepClass3d_SolidClassifier.hxx>
 #include <GProp_GProps.hxx>
-#include <BRepGProp.hxx>
 #include <Precision.hxx>
 
-GeomAlgoAPI_Placement::GeomAlgoAPI_Placement(
-  std::shared_ptr<GeomAPI_Shape> theSourceSolid,
-  std::shared_ptr<GeomAPI_Shape> theDestSolid,
-  std::shared_ptr<GeomAPI_Shape> theSourceShape,
-  std::shared_ptr<GeomAPI_Shape> theDestShape,
-  bool theIsReverse,
-  bool theIsCentering, 
-  bool theSimpleTransform)
-  : myDone(false),
-  myShape(new GeomAPI_Shape())
+GeomAlgoAPI_Placement::GeomAlgoAPI_Placement(const std::shared_ptr<GeomAPI_Shape> theSourceSolid,
+                                             const std::shared_ptr<GeomAPI_Shape> theDestSolid,
+                                             const std::shared_ptr<GeomAPI_Shape> theSourceShape,
+                                             const std::shared_ptr<GeomAPI_Shape> theDestShape,
+                                             const bool theIsReverse,
+                                             const bool theIsCentering,
+                                             const bool theSimpleTransform)
 {
-  build(theSourceSolid, theDestSolid, theSourceShape, theDestShape, 
+  build(theSourceSolid, theDestSolid, theSourceShape, theDestShape,
     theIsReverse, theIsCentering, theSimpleTransform);
 }
 
-void GeomAlgoAPI_Placement::build(
-  const std::shared_ptr<GeomAPI_Shape>& theSourceSolid,
-  const std::shared_ptr<GeomAPI_Shape>& theDestSolid,
-  const std::shared_ptr<GeomAPI_Shape>& theSourceShape,
-  const std::shared_ptr<GeomAPI_Shape>& theDestShape,
-  bool theIsReverse,
-  bool theIsCentering,
-  bool theSimpleTransform)
+void GeomAlgoAPI_Placement::build(const std::shared_ptr<GeomAPI_Shape>& theSourceSolid,
+                                  const std::shared_ptr<GeomAPI_Shape>& theDestSolid,
+                                  const std::shared_ptr<GeomAPI_Shape>& theSourceShape,
+                                  const std::shared_ptr<GeomAPI_Shape>& theDestShape,
+                                  const bool theIsReverse,
+                                  const bool theIsCentering,
+                                  const bool theSimpleTransform)
 {
   // Filling the parameters of the objects
   static const int aNbObjects = 2;
@@ -57,15 +54,20 @@ void GeomAlgoAPI_Placement::build(
   bool hasDirection[aNbObjects];
   std::shared_ptr<GeomAPI_Shape> aShapes[aNbObjects] = {theSourceShape, theDestShape};
 
+  GProp_GProps aProps;
+  static const double aPropEps = 1.e-4;
   for (int i = 0; i < aNbObjects; i++) {
     if (aShapes[i]->isFace()) {
       std::shared_ptr<GeomAPI_Face> aFace(new GeomAPI_Face(aShapes[i]));
       std::shared_ptr<GeomAPI_Pln> aPlane = aFace->getPlane();
       std::shared_ptr<GeomAPI_Dir> aDir = aPlane->direction();
-      std::shared_ptr<GeomAPI_Pnt> aLoc = aPlane->location();
-      aSrcDstPoints[i].SetCoord(aLoc->x(), aLoc->y(), aLoc->z());
       aSrcDstNormals[i].SetCoord(aDir->x(), aDir->y(), aDir->z());
-    } else if (aShapes[i]->isEdge()) {
+
+      BRepGProp::SurfaceProperties(aFace->impl<TopoDS_Face>(), aProps, aPropEps);
+      gp_Pnt aLoc = aProps.CentreOfMass();
+      aSrcDstPoints[i].SetCoord(aLoc.X(), aLoc.Y(), aLoc.Z());
+    }
+    else if (aShapes[i]->isEdge()) {
       std::shared_ptr<GeomAPI_Edge> anEdge(new GeomAPI_Edge(aShapes[i]));
       std::shared_ptr<GeomAPI_Lin> aLine = anEdge->line();
       std::shared_ptr<GeomAPI_Dir> aDir = aLine->direction();
@@ -74,7 +76,8 @@ void GeomAlgoAPI_Placement::build(
       std::shared_ptr<GeomAPI_XYZ> aLoc = aFirstPnt->xyz()->added(aLastPnt->xyz())->multiplied(0.5);
       aSrcDstPoints[i].SetCoord(aLoc->x(), aLoc->y(), aLoc->z());
       aSrcDstDirections[i].SetCoord(aDir->x(), aDir->y(), aDir->z());
-    } else if (aShapes[i]->isVertex()) {
+    }
+    else if (aShapes[i]->isVertex()) {
       std::shared_ptr<GeomAPI_Vertex> aVertex(new GeomAPI_Vertex(aShapes[i]));
       std::shared_ptr<GeomAPI_Pnt> aPnt = aVertex->point();
       aSrcDstPoints[i].SetCoord(aPnt->x(), aPnt->y(), aPnt->z());
@@ -114,11 +117,13 @@ void GeomAlgoAPI_Placement::build(
       if (hasDirection[0] || hasDirection[1]) { // plane - line
         int anInd = hasDirection[0] ? 0 : 1;
         gp_Vec aVec = aSrcDstNormals[1 - anInd].Crossed(aSrcDstDirections[anInd]);
-        if (aVec.SquareMagnitude() < Precision::SquareConfusion()) { // normal and direction are collinear
+        if (aVec.SquareMagnitude() < Precision::SquareConfusion()) {
+          // normal and direction are collinear
           aVec = aSrcDstNormals[1 - anInd].Crossed(
             gp_Vec(aSrcDstPoints[1 - anInd], aSrcDstPoints[anInd]));
-          if (aVec.SquareMagnitude() < Precision::SquareConfusion()) { // normal and points direction are collinear
-            if (Abs(aSrcDstNormals[1 - anInd].Y()) >= Precision::Confusion() || 
+          if (aVec.SquareMagnitude() < Precision::SquareConfusion()) {
+            // normal and points direction are collinear
+            if (Abs(aSrcDstNormals[1 - anInd].Y()) >= Precision::Confusion() ||
               Abs(aSrcDstNormals[1 - anInd].Z()) >= Precision::Confusion())
               aVec = gp::DX();
             else
@@ -158,7 +163,7 @@ void GeomAlgoAPI_Placement::build(
         gp_Vec aVec(aSrcDstPoints[anInd], aSrcDstPoints[1 - anInd]);
         aVec.Cross(aSrcDstDirections[anInd]);
         if (aVec.SquareMagnitude() < Precision::SquareConfusion()) { // point is on line
-          if (Abs(aSrcDstDirections[1 - anInd].Y()) >= Precision::Confusion() || 
+          if (Abs(aSrcDstDirections[1 - anInd].Y()) >= Precision::Confusion() ||
             Abs(aSrcDstDirections[1 - anInd].Z()) >= Precision::Confusion())
             aVec = gp::DX();
           else
@@ -192,82 +197,35 @@ void GeomAlgoAPI_Placement::build(
 
   if (theSimpleTransform) { // just add transformation
     TopLoc_Location aDelta(aTrsf);
-    TopoDS_Shape aResult = aSourceShape.Moved(aDelta);
-    myShape->setImpl(new TopoDS_Shape(aResult));
     // store the accumulated information about the result and this delta
-    myTrsf = std::shared_ptr<GeomAPI_Trsf>(new GeomAPI_Trsf(new gp_Trsf(aTrsf * aSourceShape.Location().Transformation())));
-    myDone = true; // it is allways true for simple transformation generation
+    myTrsf.reset(new GeomAPI_Trsf(new gp_Trsf(aTrsf)));
+    TopoDS_Shape aResult = aSourceShape.Moved(aDelta);
+    std::shared_ptr<GeomAPI_Shape> aShape(new GeomAPI_Shape());
+    aShape->setImpl(new TopoDS_Shape(aResult));
+    this->setShape(aShape);
+    this->setDone(true); // it is allways true for simple transformation generation
   } else { // internal rebuild of the shape
     // Transform the shape with copying it
     BRepBuilderAPI_Transform* aBuilder = new BRepBuilderAPI_Transform(aSourceShape, aTrsf, true);
-    if(aBuilder) {
-      setImpl(aBuilder);
-      myDone = aBuilder->IsDone() == Standard_True;
-      if (myDone) {
-        TopoDS_Shape aResult = aBuilder->Shape();
-        // fill data map to keep correct orientation of sub-shapes 
-        for (TopExp_Explorer Exp(aResult,TopAbs_FACE); Exp.More(); Exp.Next()) {
-          std::shared_ptr<GeomAPI_Shape> aCurrentShape(new GeomAPI_Shape());
-          aCurrentShape->setImpl(new TopoDS_Shape(Exp.Current()));
-          myMap.bind(aCurrentShape, aCurrentShape);
-        }
-        myShape->setImpl(new TopoDS_Shape(aResult));
-        myMkShape = new GeomAlgoAPI_MakeShape (aBuilder);
-      }
+    if(!aBuilder) {
+      return;
     }
-  }
-}
-
-//============================================================================
-const bool GeomAlgoAPI_Placement::isValid() const
-{
-  if (myShape.get()) { // only for not-simple transform
-    BRepCheck_Analyzer aChecker(myShape->impl<TopoDS_Shape>());
-    return (aChecker.IsValid() == Standard_True);
-  }
-  return true;
-}
+    this->setImpl(aBuilder);
+    this->setBuilderType(OCCT_BRepBuilderAPI_MakeShape);
+    if(aBuilder->IsDone() != Standard_True) {
+      return;
+    }
+    TopoDS_Shape aResult = aBuilder->Shape();
 
-//============================================================================
-const bool GeomAlgoAPI_Placement::hasVolume() const
-{
-  bool hasVolume(false);
-  if(isValid()) {
-    const TopoDS_Shape& aRShape = myShape->impl<TopoDS_Shape>();
-    GProp_GProps aGProp;
-    BRepGProp::VolumeProperties(aRShape, aGProp);
-    if(aGProp.Mass() > Precision::Confusion()) 
-      hasVolume = true;        
+    std::shared_ptr<GeomAPI_Shape> aShape(new GeomAPI_Shape());
+    aShape->setImpl(new TopoDS_Shape(aResult));
+    this->setShape(aShape);
+    this->setDone(true);
   }
-  return hasVolume;
-}
-
-//============================================================================
-const std::shared_ptr<GeomAPI_Shape>& GeomAlgoAPI_Placement::shape () const 
-{
-  return myShape;
-}
-
-//============================================================================
-void GeomAlgoAPI_Placement::mapOfShapes (GeomAPI_DataMapOfShapeShape& theMap) const
-{
-  theMap = myMap;
-}
-
-//============================================================================
-GeomAlgoAPI_MakeShape * GeomAlgoAPI_Placement::makeShape() const
-{
-  return myMkShape;
 }
 
+//=================================================================================================
 std::shared_ptr<GeomAPI_Trsf> GeomAlgoAPI_Placement::transformation() const
 {
   return myTrsf;
 }
-
-//============================================================================
-GeomAlgoAPI_Placement::~GeomAlgoAPI_Placement()
-{
-  if (!empty())
-    myMap.clear();
-}