src/BlockFix/BlockFix_BlockFixAPI.cxx

   1 // Copyright (C) 2007-2024  CEA, EDF, 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 //  File:    BlockFix_BlockFixAPI.cxx
  24 //  Created: Tue Dec  7 11:59:05 2004
  25 //  Author:  Pavel DURANDIN
  26
  27 #include <BlockFix_BlockFixAPI.hxx>
  28
  29 #include <BlockFix.hxx>
  30 #include <BlockFix_UnionFaces.hxx>
  31 #include <BlockFix_UnionEdges.hxx>
  32
  33 #include <ShapeUpgrade_RemoveLocations.hxx>
  34 #include <TopoDS_Edge.hxx>
  35 #include <ShapeUpgrade_UnifySameDomain.hxx>
  36
  37 #include <Precision.hxx>
  38
  39 #include <BRepTools.hxx>
  40 #include <BRepBuilderAPI_Copy.hxx>
  41 #include <BRepCheck_Analyzer.hxx>
  42
  43 #include <Basics_OCCTVersion.hxx>
  44
  45 IMPLEMENT_STANDARD_RTTIEXT(BlockFix_BlockFixAPI, Standard_Transient)
  46
  47 //=======================================================================
  48 //function : BlockFix_BlockFixAPI
  49 //purpose  :
  50 //=======================================================================
  51 BlockFix_BlockFixAPI::BlockFix_BlockFixAPI()
  52 {
  53   myTolerance = Precision::Confusion();
  54   myOptimumNbFaces = 6;
  55 }
  56
  57 //=======================================================================
  58 //function : ~BlockFix_BlockFixAPI
  59 //purpose  :
  60 //=======================================================================
  61 BlockFix_BlockFixAPI::~BlockFix_BlockFixAPI() {}
  62
  63 //=======================================================================
  64 //function : Perform
  65 //purpose  :
  66 //=======================================================================
  67 void BlockFix_BlockFixAPI::Perform()
  68 {
  69   // processing spheres with degenerativities
  70   TopoDS_Shape aShape = Shape();
  71   // Copy the shape to avoid modification of initial shape
  72   // while trying the approach with small rotation
  73   BRepBuilderAPI_Copy aMC (aShape);
  74   if (!aMC.IsDone()) return;
  75   TopoDS_Shape aSCopy = aMC.Shape();
  76   TopoDS_Shape aNewShape = BlockFix::RotateSphereSpace(aSCopy, myTolerance, Standard_True);
  77   BRepCheck_Analyzer ana (aNewShape, false);
  78   if (ana.IsValid()) {
  79     if (aNewShape == aSCopy)
  80       myShape = aShape;
  81     else
  82       myShape = aNewShape;
  83   }
  84   else {
  85     myShape = BlockFix::RotateSphereSpace(aShape, myTolerance, Standard_False);
  86   }
  87
  88   // try to approximate non-canonic surfaces
  89   // with singularities on boundaries by filling
  90   myShape = BlockFix::RefillProblemFaces(myShape);
  91
  92   // faces unification
  93   TopoDS_Shape aResult = myShape;
  94 #if OCC_VERSION_LARGE < 0x07050301
  95   BlockFix_UnionFaces aFaceUnifier;
  96   aFaceUnifier.GetTolerance() = myTolerance;
  97   aFaceUnifier.GetOptimumNbFaces() = myOptimumNbFaces;
  98   aResult = aFaceUnifier.Perform(aResult);
  99 #else
 100   if (myOptimumNbFaces > 1) {
 101     // use old algo BlockFix_UnionFaces for exactly given resulting number of faces
 102     BlockFix_UnionFaces aFaceUnifier;
 103     aFaceUnifier.GetTolerance() = myTolerance;
 104     aFaceUnifier.GetOptimumNbFaces() = myOptimumNbFaces;
 105     aResult = aFaceUnifier.Perform(aResult);
 106   }
 107   else if (myOptimumNbFaces != -1) {
 108     // use OCCT algo ShapeUpgrade_UnifySameDomain
 109     ShapeUpgrade_UnifySameDomain Unifier;
 110     //only faces
 111     Standard_Boolean isUnifyEdges = Standard_False;
 112     Standard_Boolean isUnifyFaces = Standard_True;
 113     Standard_Boolean isConcatBSplines = Standard_True;
 114     Unifier.Initialize(myShape, isUnifyEdges, isUnifyFaces, isConcatBSplines);
 115     //Unifier.SetLinearTolerance(myTolerance);
 116     Unifier.SetLinearTolerance(Precision::Confusion());
 117     Unifier.SetAngularTolerance(Precision::Confusion());
 118     //Unifier.SetAngularTolerance(1e-5);
 119     //Unifier.SetAngularTolerance(0.1);
 120     Unifier.Build();
 121     aResult = Unifier.Shape();
 122   }
 123   else {
 124     // myOptimumNbFaces == -1 means do not union faces
 125   }
 126 #endif
 127
 128   // avoid problem with degenerated edges appearance
 129   // due to shape quality regress
 130   ShapeUpgrade_RemoveLocations RemLoc;
 131   RemLoc.Remove(aResult);
 132   aResult = RemLoc.GetResult();
 133
 134   // edges unification
 135 #if OCC_VERSION_LARGE < 0x07050301
 136   BlockFix_UnionEdges anEdgeUnifier;
 137   myShape = anEdgeUnifier.Perform(aResult,myTolerance);
 138 #else
 139   ShapeUpgrade_UnifySameDomain Unifier;
 140   Standard_Boolean isUnifyEdges = Standard_True;
 141   Standard_Boolean isUnifyFaces = Standard_False; //only edges
 142   Standard_Boolean isConcatBSplines = Standard_True;
 143   Unifier.Initialize(aResult, isUnifyEdges, isUnifyFaces, isConcatBSplines);
 144   Unifier.SetLinearTolerance(myTolerance);
 145   Unifier.SetAngularTolerance(1e-5);
 146   Unifier.Build();
 147   myShape = Unifier.Shape();
 148 #endif
 149
 150   TopoDS_Shape aRes = BlockFix::FixRanges(myShape,myTolerance);
 151   myShape = aRes;
 152 }