/src/contrib/geom-5.1.2.7/src/BlockFix/BlockFix_PeriodicSurfaceModifier.cpp

  1//  Copyright (C) 2007-2008  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.
 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// File:      BlockFix_PeriodicSurfaceModifier.cxx
 23// Created:   15.12.04 10:08:50
 24// Author:    Sergey KUUL
 25//
 26#include <BlockFix_PeriodicSurfaceModifier.ixx>
 27
 28#include <BRep_Builder.hxx>
 29#include <BRep_Tool.hxx>
 30#include <BRepTools.hxx>
 31#include <Geom_CylindricalSurface.hxx>
 32#include <Geom_SphericalSurface.hxx>
 33#include <ShapeFix_Edge.hxx>
 34#include <TopExp.hxx>
 35
 36
 37//=======================================================================
 38//function : BlockFix_PeriodicSurfaceModifier()
 39//purpose  : Constructor
 40//=======================================================================
 41
 42BlockFix_PeriodicSurfaceModifier::BlockFix_PeriodicSurfaceModifier (  )
 43{
 44  myMapOfFaces.Clear();
 45  myMapOfSurfaces.Clear();
 46}
 47
 48
 49//=======================================================================
 50//function : SetTolerance
 51//purpose  : 
 52//=======================================================================
 53
 54void BlockFix_PeriodicSurfaceModifier::SetTolerance(const Standard_Real Tol)
 55{
 56  myTolerance = Tol;
 57}
 58
 59
 60//=======================================================================
 61//function : ModifySurface
 62//purpose  : auxilary
 63//=======================================================================
 64
 65static Standard_Boolean ModifySurface(const TopoDS_Face& aFace,
 66                                      const Handle(Geom_Surface)& aSurface,
 67                                      Handle(Geom_Surface)& aNewSurface)
 68{
 69  Handle(Geom_Surface) S = aSurface;
 70
 71  if(S->IsKind(STANDARD_TYPE(Geom_CylindricalSurface))) {
 72    Handle(Geom_CylindricalSurface) aCyl = 
 73      Handle(Geom_CylindricalSurface)::DownCast(S);
 74    Standard_Real Umin, Umax, Vmin, Vmax;
 75    BRepTools::UVBounds(aFace, Umin, Umax, Vmin, Vmax);
 76    if( Umin<-Precision::PConfusion() || Umax>2*PI+Precision::PConfusion() ) {
 77      gp_Ax3 ax3 = aCyl->Position();
 78      gp_Ax1 NDir = ax3.Axis();
 79      gp_Ax3 newax3 = ax3.Rotated(NDir,Umin-Precision::PConfusion());
 80      Handle(Geom_CylindricalSurface) aNewCyl =
 81        new Geom_CylindricalSurface(newax3,aCyl->Radius());
 82      aNewSurface = aNewCyl;
 83      return Standard_True;
 84    }
 85  }
 86  
 87  if(S->IsKind(STANDARD_TYPE(Geom_SphericalSurface))) {
 88    Handle(Geom_SphericalSurface) aSphere = Handle(Geom_SphericalSurface)::DownCast(S);
 89    Standard_Real Umin, Umax, Vmin, Vmax;
 90    BRepTools::UVBounds(aFace, Umin, Umax, Vmin, Vmax);
 91    if( Umin<-Precision::PConfusion() || Umax>2*PI+Precision::PConfusion() ) {
 92      gp_Ax3 ax3 = aSphere->Position();
 93      gp_Ax1 NDir = ax3.Axis();
 94      gp_Ax3 newax3 = ax3.Rotated(NDir,Umin-Precision::PConfusion());
 95      Handle(Geom_SphericalSurface) aNewSphere = new Geom_SphericalSurface(newax3,aSphere->Radius());
 96      aNewSurface = aNewSphere;
 97      return Standard_True;
 98    }
 99  }
100
101  return Standard_False;
102}
103                                      
104
105//=======================================================================
106//function : NewSurface
107//purpose  : 
108//=======================================================================
109
110Standard_Boolean BlockFix_PeriodicSurfaceModifier::NewSurface(const TopoDS_Face& F,
111                                                              Handle(Geom_Surface)& S,
112                                                              TopLoc_Location& L,Standard_Real& Tol,
113                                                              Standard_Boolean& RevWires,
114                                                              Standard_Boolean& RevFace) 
115{
116  TopLoc_Location LS;
117  Handle(Geom_Surface) SIni = BRep_Tool::Surface(F, LS);
118  
119  if(ModifySurface(F, SIni, S)) {
120  
121    RevWires = Standard_False;
122    RevFace = Standard_False;
123    
124    L = LS;
125    Tol = BRep_Tool::Tolerance(F);
126     
127    Standard_Integer anIndex = myMapOfSurfaces.Add(S);
128    myMapOfFaces.Bind(F,anIndex);
129    return Standard_True;
130  }
131
132  return Standard_False;
133}
134
135
136//=======================================================================
137//function : NewCurve
138//purpose  : 
139//=======================================================================
140
141Standard_Boolean BlockFix_PeriodicSurfaceModifier::NewCurve(const TopoDS_Edge& /*E*/,
142                                                            Handle(Geom_Curve)& /*C*/,
143                                                            TopLoc_Location& /*L*/,
144                                                            Standard_Real& /*Tol*/) 
145{
146  return Standard_False;
147}
148
149
150//=======================================================================
151//function : NewPoint
152//purpose  : 
153//=======================================================================
154
155Standard_Boolean BlockFix_PeriodicSurfaceModifier::NewPoint(const TopoDS_Vertex& /*V*/,
156                                                            gp_Pnt& /*P*/,
157                                                            Standard_Real& /*Tol*/) 
158{
159  return Standard_False;
160}
161
162
163//=======================================================================
164//function : NewCurve2d
165//purpose  : 
166//=======================================================================
167
168Standard_Boolean BlockFix_PeriodicSurfaceModifier::NewCurve2d(const TopoDS_Edge& E,
169                                                              const TopoDS_Face& F,
170                                                              const TopoDS_Edge& /*NewE*/,
171                                                              const TopoDS_Face& /*NewF*/,
172                                                              Handle(Geom2d_Curve)& C,
173                                                              Standard_Real& Tol) 
174{
175  //check if undelying surface of the face was modified
176  if(myMapOfFaces.IsBound(F)) {
177    Standard_Integer anIndex = myMapOfFaces.Find(F);
178    
179    Handle(Geom_Surface) aNewSurf = Handle(Geom_Surface)::DownCast(myMapOfSurfaces.FindKey(anIndex));
180    
181    Standard_Real f,l;
182    TopLoc_Location LC, LS;
183    Handle(Geom_Curve) C3d = BRep_Tool::Curve ( E, LC, f, l );
184    Handle(Geom_Surface) S = BRep_Tool::Surface(F, LS);
185  
186    //taking into accound the orientation of the seam
187    C = BRep_Tool::CurveOnSurface(E,F,f,l);
188    Tol = BRep_Tool::Tolerance(E);
189     
190    BRep_Builder B;
191    TopoDS_Edge TempE;
192    B.MakeEdge(TempE);
193    B.Add(TempE, TopExp::FirstVertex(E));
194    B.Add(TempE, TopExp::LastVertex(E));
195
196    if(!C3d.IsNull()) 
197      B.UpdateEdge(TempE, Handle(Geom_Curve)::DownCast(C3d->Transformed(LC.Transformation())), Precision::Confusion());
198    B.Range(TempE, f, l);
199    
200    Handle(ShapeFix_Edge) sfe = new ShapeFix_Edge;
201    Handle(Geom_Surface) STemp = Handle(Geom_Surface)::DownCast(aNewSurf->Transformed(LS.Transformation()));
202    TopLoc_Location LTemp;
203    LTemp.Identity();
204    
205    Standard_Boolean isClosed = BRep_Tool::IsClosed (E, F);
206    Standard_Real aWorkTol = 2*myTolerance+Tol;
207    sfe->FixAddPCurve(TempE, STemp, LTemp, isClosed, Max(Precision::Confusion(), aWorkTol));
208    sfe->FixSameParameter(TempE);
209      
210    //keep the orientation of original edge
211    TempE.Orientation(E.Orientation());
212    C = BRep_Tool::CurveOnSurface(TempE, STemp, LTemp, f, l);
213    
214    //surface was modified
215    return Standard_True;
216  }
217
218  return Standard_False;
219}
220
221
222//=======================================================================
223//function : NewParameter
224//purpose  : 
225//=======================================================================
226
227Standard_Boolean BlockFix_PeriodicSurfaceModifier::NewParameter(const TopoDS_Vertex& /*V*/,
228                                                                const TopoDS_Edge& /*E*/,
229                                                                Standard_Real& /*P*/,
230                                                                Standard_Real& /*Tol*/) 
231{
232  return Standard_False;
233}
234
235
236//=======================================================================
237//function : Continuity
238//purpose  : 
239//=======================================================================
240
241GeomAbs_Shape BlockFix_PeriodicSurfaceModifier::Continuity(const TopoDS_Edge& E,
242                                                           const TopoDS_Face& F1,
243                                                           const TopoDS_Face& F2,
244                                                           const TopoDS_Edge& /*NewE*/,
245                                                           const TopoDS_Face& /*NewF1*/,
246                                                           const TopoDS_Face& /*NewF2*/) 
247{
248  return BRep_Tool::Continuity(E,F1,F2);
249}
250