src/ShapeCustom/ShapeCustom_SweptToElementary.cxx

   1 // Copyright (c) 1999-2012 OPEN CASCADE SAS
   2 //
   3 // The content of this file is subject to the Open CASCADE Technology Public
   4 // License Version 6.5 (the "License"). You may not use the content of this file
   5 // except in compliance with the License. Please obtain a copy of the License
   6 // at http://www.opencascade.org and read it completely before using this file.
   7 //
   8 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
   9 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
  10 //
  11 // The Original Code and all software distributed under the License is
  12 // distributed on an "AS IS" basis, without warranty of any kind, and the
  13 // Initial Developer hereby disclaims all such warranties, including without
  14 // limitation, any warranties of merchantability, fitness for a particular
  15 // purpose or non-infringement. Please see the License for the specific terms
  16 // and conditions governing the rights and limitations under the License.
  17
  18 #include <ShapeCustom_SweptToElementary.ixx>
  19
  20 #include <BRep_Builder.hxx>
  21 #include <BRep_Tool.hxx>
  22 #include <BRep_TEdge.hxx>
  23 #include <BRep_ListIteratorOfListOfCurveRepresentation.hxx>
  24 #include <BRep_GCurve.hxx>
  25
  26 #include <Geom_SweptSurface.hxx>
  27 #include <Geom_SurfaceOfRevolution.hxx>
  28 #include <Geom_SurfaceOfLinearExtrusion.hxx>
  29 #include <Geom_RectangularTrimmedSurface.hxx>
  30 #include <Geom_OffsetSurface.hxx>
  31 #include <Geom_ToroidalSurface.hxx>
  32 #include <Geom_SphericalSurface.hxx>
  33 #include <Geom_CylindricalSurface.hxx>
  34 #include <Geom_ConicalSurface.hxx>
  35 #include <Geom_Plane.hxx>
  36
  37 #include <Adaptor3d_SurfaceOfRevolution.hxx>
  38 #include <Adaptor3d_SurfaceOfLinearExtrusion.hxx>
  39 #include <GeomAdaptor_HCurve.hxx>
  40
  41 #include <gp_Pln.hxx>
  42 #include <gp_Cylinder.hxx>
  43 #include <gp_Sphere.hxx>
  44 #include <gp_Cone.hxx>
  45 #include <gp_Torus.hxx>
  46
  47 #include <ShapeAnalysis_Surface.hxx>
  48
  49 #include <BRepTools.hxx>
  50
  51
  52 //=======================================================================
  53 //function : ShapeCustom_SweptToElementary
  54 //purpose  :
  55 //=======================================================================
  56
  57 ShapeCustom_SweptToElementary::ShapeCustom_SweptToElementary()
  58 {
  59 }
  60
  61
  62 //=======================================================================
  63 //function : IsToConvert
  64 //purpose  : auxilary (Analyze surface: is it to be converted?)
  65 //=======================================================================
  66 static Standard_Boolean IsToConvert (const Handle(Geom_Surface) &S,
  67                                      Handle(Geom_SweptSurface) &SS)
  68 {
  69   Handle(Geom_Surface) Stmp;
  70
  71   if(S->IsKind(STANDARD_TYPE(Geom_SweptSurface))) {
  72     SS = Handle(Geom_SweptSurface)::DownCast(S);
  73     return Standard_True;
  74   }
  75   if(S->IsKind(STANDARD_TYPE(Geom_RectangularTrimmedSurface))) {
  76     Handle(Geom_RectangularTrimmedSurface) RTS =
  77       Handle(Geom_RectangularTrimmedSurface)::DownCast(S);
  78     Stmp = RTS->BasisSurface();
  79   }
  80   else if(S->IsKind(STANDARD_TYPE(Geom_OffsetSurface))) {
  81     Handle(Geom_OffsetSurface) OS = Handle(Geom_OffsetSurface)::DownCast(S);
  82     Stmp = OS->BasisSurface();
  83   }
  84   if(Stmp.IsNull() ) return Standard_False;
  85   if(S->IsKind(STANDARD_TYPE(Geom_SweptSurface))) {
  86     SS = Handle(Geom_SweptSurface)::DownCast(Stmp);
  87     return Standard_True;
  88   }
  89   return Standard_False;
  90 }
  91
  92
  93 //=======================================================================
  94 //function : NewSurface
  95 //purpose  :
  96 //=======================================================================
  97
  98 Standard_Boolean ShapeCustom_SweptToElementary::NewSurface(const TopoDS_Face& F,
  99                                                            Handle(Geom_Surface)& S,
 100                                                            TopLoc_Location& L,
 101                                                            Standard_Real& Tol,
 102                                                            Standard_Boolean& RevWires,
 103                                                            Standard_Boolean& RevFace)
 104 {
 105   S = BRep_Tool::Surface(F,L);
 106   Handle(Geom_SweptSurface) SS;
 107   if(!IsToConvert(S,SS)) return Standard_False;
 108
 109   // case SurfaceOfRevolution
 110   if(SS->IsKind(STANDARD_TYPE(Geom_SurfaceOfRevolution))) {
 111     Handle(Geom_SurfaceOfRevolution) SR = Handle(Geom_SurfaceOfRevolution)::DownCast(SS);
 112     Handle(Geom_Curve) bc = SR->BasisCurve();
 113     gp_Ax1 ax1 = SR->Axis();
 114     Handle(GeomAdaptor_HCurve) HC = new GeomAdaptor_HCurve();
 115     HC->ChangeCurve().Load(bc,bc->FirstParameter(),bc->LastParameter());
 116     Adaptor3d_SurfaceOfRevolution AS(HC,ax1);
 117     switch(AS.GetType()){
 118     // skl 18.12.2003 - plane not used, problems in PRO14665.igs
 119     //case GeomAbs_Plane : {
 120     //  Handle(Geom_Plane) Pl = new Geom_Plane(AS.Plane());
 121     //  S = Pl;
 122     //} break;
 123     case GeomAbs_Cylinder : {
 124       Handle(Geom_CylindricalSurface) Cy =
 125         new Geom_CylindricalSurface(AS.Cylinder());
 126       S = Cy;
 127     } break;
 128     case GeomAbs_Sphere : {
 129       Handle(Geom_SphericalSurface) Sp =
 130         new Geom_SphericalSurface(AS.Sphere());
 131       S = Sp;
 132     } break;
 133     case GeomAbs_Cone : {
 134       Handle(Geom_ConicalSurface) Co =
 135         new Geom_ConicalSurface(AS.Cone());
 136       S = Co;
 137     } break;
 138     case GeomAbs_Torus : {
 139       Handle(Geom_ToroidalSurface) To =
 140         new Geom_ToroidalSurface(AS.Torus());
 141       S = To;
 142     } break;
 143     default : return Standard_False; break;
 144     }
 145   }
 146   // case SurfaceOfLinearExtrusion
 147   else if(SS->IsKind(STANDARD_TYPE(Geom_SurfaceOfLinearExtrusion))) {
 148     Handle(Geom_SurfaceOfLinearExtrusion) SLE =
 149       Handle(Geom_SurfaceOfLinearExtrusion)::DownCast(SS);
 150     Handle(Geom_Curve) bc = SLE->BasisCurve();
 151     gp_Dir dir = SLE->Direction();
 152     Handle(GeomAdaptor_HCurve) HC = new GeomAdaptor_HCurve();
 153     HC->ChangeCurve().Load(bc,bc->FirstParameter(),bc->LastParameter());
 154     Adaptor3d_SurfaceOfLinearExtrusion AS(HC,dir);
 155     switch(AS.GetType()){
 156     // skl 18.12.2003 - plane not used, problems in ims013.igs
 157     //case GeomAbs_Plane : {
 158     //  Handle(Geom_Plane) Pl = new Geom_Plane(AS.Plane());
 159     //  S = Pl;
 160     //} break;
 161     case GeomAbs_Cylinder : {
 162       Handle(Geom_CylindricalSurface) Cy =
 163         new Geom_CylindricalSurface(AS.Cylinder());
 164       S = Cy;
 165     } break;
 166     default : return Standard_False; break;
 167     }
 168   }
 169
 170   Tol = BRep_Tool::Tolerance(F);
 171   RevWires = Standard_False;
 172   RevFace = Standard_False;
 173   return Standard_True;
 174 }
 175
 176
 177 //=======================================================================
 178 //function : NewCurve
 179 //purpose  :
 180 //=======================================================================
 181
 182 Standard_Boolean ShapeCustom_SweptToElementary::NewCurve(const TopoDS_Edge& E,
 183                                                          Handle(Geom_Curve)& C,
 184                                                          TopLoc_Location& L,
 185                                                          Standard_Real& Tol)
 186 {
 187   //:p5 abv 26 Feb 99: force copying of edge if any its pcurve will be replaced
 188   Handle(BRep_TEdge)& TE = *((Handle(BRep_TEdge)*)&E.TShape());
 189
 190   // iterate on pcurves
 191   BRep_ListIteratorOfListOfCurveRepresentation itcr(TE->Curves());
 192   for ( ; itcr.More(); itcr.Next() ) {
 193     Handle(BRep_GCurve) GC = Handle(BRep_GCurve)::DownCast(itcr.Value());
 194     if ( GC.IsNull() || ! GC->IsCurveOnSurface() ) continue;
 195     Handle(Geom_Surface) S = GC->Surface();
 196     Handle(Geom_SweptSurface) SS;
 197     if(!IsToConvert(S,SS)) continue;
 198     Standard_Real f, l;
 199     C = BRep_Tool::Curve ( E, L, f, l );
 200     if ( ! C.IsNull() ) C = Handle(Geom_Curve)::DownCast ( C->Copy() );
 201     Tol = BRep_Tool::Tolerance ( E );
 202     return Standard_True;
 203   }
 204   return Standard_False;
 205 }
 206
 207
 208 //=======================================================================
 209 //function : NewPoint
 210 //purpose  :
 211 //=======================================================================
 212
 213 Standard_Boolean ShapeCustom_SweptToElementary::NewPoint(const TopoDS_Vertex& /*V*/,
 214                                                          gp_Pnt& /*P*/,Standard_Real& /*Tol*/)
 215 {
 216   // 3d points are never modified
 217   return Standard_False;
 218 }
 219
 220
 221 //=======================================================================
 222 //function : NewCurve2d
 223 //purpose  :
 224 //=======================================================================
 225
 226 Standard_Boolean ShapeCustom_SweptToElementary::NewCurve2d(const TopoDS_Edge& E,
 227                                                            const TopoDS_Face& F,
 228                                                            const TopoDS_Edge& NewE,
 229                                                            const TopoDS_Face& NewF,
 230                                                            Handle(Geom2d_Curve)& C,
 231                                                            Standard_Real& Tol)
 232 {
 233   TopLoc_Location L;
 234   Handle(Geom_Surface) S = BRep_Tool::Surface(F,L);
 235   Handle(Geom_SweptSurface) SS;
 236
 237   // just copy pcurve if either its surface is changing or edge was copied
 238   if ( !IsToConvert(S,SS) && E.IsSame(NewE) ) return Standard_False;
 239
 240   Standard_Real f, l;
 241   C = BRep_Tool::CurveOnSurface(E,F,f,l);
 242   if ( ! C.IsNull() ) {
 243     C = Handle(Geom2d_Curve)::DownCast ( C->Copy() );
 244     Handle(Geom_Surface) NS = BRep_Tool::Surface(NewF,L);
 245     if ( !NS.IsNull() && NS->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
 246       if(SS->IsKind(STANDARD_TYPE(Geom_SurfaceOfRevolution))) {
 247         Handle(Geom_SurfaceOfRevolution) SR = Handle(Geom_SurfaceOfRevolution)::DownCast(SS);
 248         Standard_Real U1,U2,V1,V2;
 249         SR->Bounds(U1,U2,V1,V2);
 250         gp_Pnt P0;
 251         SR->D0(U1,V1,P0);
 252         Handle(ShapeAnalysis_Surface) sas = new ShapeAnalysis_Surface(NS);
 253         gp_Pnt2d p2d = sas->ValueOfUV(P0,Precision::Confusion());
 254         gp_Vec2d shift(p2d.X()-U1,p2d.Y()-V1);
 255         C->Translate(shift);
 256       }
 257     }
 258     if ( !NS.IsNull() && NS->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ) {
 259       if(SS->IsKind(STANDARD_TYPE(Geom_SurfaceOfRevolution))) {
 260        Handle(Geom_SurfaceOfRevolution) SR = Handle(Geom_SurfaceOfRevolution)::DownCast(SS);
 261         gp_Pnt PR,PS;
 262         Handle(Geom_SphericalSurface) SPH = Handle(Geom_SphericalSurface)::DownCast(NS);
 263         Standard_Real US1,US2,VS1,VS2;
 264         SPH->Bounds(US1,US2,VS1,VS2);
 265         SPH->D0(US1,VS1,PS);
 266         Standard_Real UR1,UR2,VR1,VR2;
 267         SR->Bounds(UR1,UR2,VR1,VR2);
 268         SR->D0(UR1,VR1,PR);
 269         gp_Pnt P0 = SPH->Location();
 270         gp_Vec VS(P0,PS);
 271         gp_Vec VR(P0,PR);
 272         Standard_Real angle = VS.Angle(VR);
 273         gp_Vec2d shift(0,VS1-VR1+angle);
 274         C->Translate(shift);
 275       }
 276     }
 277   }
 278
 279   Tol = BRep_Tool::Tolerance ( E );
 280   return Standard_True;
 281 }
 282
 283
 284 //=======================================================================
 285 //function : NewParameter
 286 //purpose  :
 287 //=======================================================================
 288
 289 Standard_Boolean ShapeCustom_SweptToElementary::NewParameter(const TopoDS_Vertex& /*V*/,
 290                                                              const TopoDS_Edge& /*E*/,
 291                                                              Standard_Real& /*P*/,
 292                                                              Standard_Real& /*Tol*/)
 293 {
 294   return Standard_False;
 295 }
 296
 297
 298 //=======================================================================
 299 //function : Continuity
 300 //purpose  :
 301 //=======================================================================
 302
 303 GeomAbs_Shape ShapeCustom_SweptToElementary::Continuity(const TopoDS_Edge& E,
 304                                                         const TopoDS_Face& F1,
 305                                                         const TopoDS_Face& F2,
 306                                                         const TopoDS_Edge& /*NewE*/,
 307                                                         const TopoDS_Face& /*NewF1*/,
 308                                                         const TopoDS_Face& /*NewF2*/)
 309 {
 310   return BRep_Tool::Continuity(E,F1,F2);
 311 }
 312