src/Geom/Geom_SurfaceOfRevolution.cxx

   1 // Created on: 1993-03-10
   2 // Created by: JCV
   3 // Copyright (c) 1993-1999 Matra Datavision
   4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   5 //
   6 // This file is part of Open CASCADE Technology software library.
   7 //
   8 // This library is free software; you can redistribute it and/or modify it under
   9 // the terms of the GNU Lesser General Public License version 2.1 as published
  10 // by the Free Software Foundation, with special exception defined in the file
  11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  12 // distribution for complete text of the license and disclaimer of any warranty.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17
  18 #include <BSplCLib.hxx>
  19 #include <BSplSLib.hxx>
  20 #include <Geom_BSplineCurve.hxx>
  21 #include <Geom_Circle.hxx>
  22 #include <Geom_Curve.hxx>
  23 #include <Geom_Geometry.hxx>
  24 #include <Geom_SurfaceOfRevolution.hxx>
  25 #include <Geom_UndefinedDerivative.hxx>
  26 #include <GeomEvaluator_SurfaceOfRevolution.hxx>
  27 #include <gp.hxx>
  28 #include <gp_Ax1.hxx>
  29 #include <gp_Ax2.hxx>
  30 #include <gp_Ax2d.hxx>
  31 #include <gp_Dir.hxx>
  32 #include <gp_GTrsf2d.hxx>
  33 #include <gp_Lin.hxx>
  34 #include <gp_Pnt.hxx>
  35 #include <gp_Trsf.hxx>
  36 #include <gp_Vec.hxx>
  37 #include <gp_XYZ.hxx>
  38 #include <Precision.hxx>
  39 #include <Standard_ConstructionError.hxx>
  40 #include <Standard_NotImplemented.hxx>
  41 #include <Standard_RangeError.hxx>
  42 #include <Standard_Type.hxx>
  43
  44 IMPLEMENT_STANDARD_RTTIEXT(Geom_SurfaceOfRevolution,Geom_SweptSurface)
  45
  46 #define  POLES    (poles->Array2())
  47 #define  WEIGHTS  (weights->Array2())
  48 #define  UKNOTS   (uknots->Array1())
  49 #define  VKNOTS   (vknots->Array1())
  50 #define  UFKNOTS  (ufknots->Array1())
  51 #define  VFKNOTS  (vfknots->Array1())
  52 #define  FMULTS   (BSplCLib::NoMults())
  53
  54 typedef Geom_SurfaceOfRevolution         SurfaceOfRevolution;
  55 typedef Geom_Curve                       Curve;
  56 typedef gp_Ax1  Ax1;
  57 typedef gp_Ax2  Ax2;
  58 typedef gp_Dir  Dir;
  59 typedef gp_Lin  Lin;
  60 typedef gp_Pnt  Pnt;
  61 typedef gp_Trsf Trsf;
  62 typedef gp_Vec  Vec;
  63 typedef gp_XYZ  XYZ;
  64
  65
  66
  67
  68
  69 //=======================================================================
  70 //function : Copy
  71 //purpose  :
  72 //=======================================================================
  73
  74 Handle(Geom_Geometry) Geom_SurfaceOfRevolution::Copy () const {
  75
  76   return new Geom_SurfaceOfRevolution (basisCurve, Axis());
  77 }
  78
  79
  80 //=======================================================================
  81 //function : Geom_SurfaceOfRevolution
  82 //purpose  :
  83 //=======================================================================
  84
  85 Geom_SurfaceOfRevolution::Geom_SurfaceOfRevolution
  86   (const Handle(Geom_Curve)& C ,
  87    const Ax1&           A1 ) : loc (A1.Location()) {
  88
  89   direction  = A1.Direction();
  90   SetBasisCurve(C);
  91 }
  92
  93
  94 //=======================================================================
  95 //function : UReverse
  96 //purpose  :
  97 //=======================================================================
  98
  99 void Geom_SurfaceOfRevolution::UReverse () {
 100
 101   direction.Reverse();
 102   myEvaluator->SetDirection(direction);
 103 }
 104
 105
 106 //=======================================================================
 107 //function : UReversedParameter
 108 //purpose  :
 109 //=======================================================================
 110
 111 Standard_Real Geom_SurfaceOfRevolution::UReversedParameter (const Standard_Real U) const {
 112
 113   return ( 2.*M_PI - U);
 114 }
 115
 116
 117 //=======================================================================
 118 //function : VReverse
 119 //purpose  :
 120 //=======================================================================
 121
 122 void Geom_SurfaceOfRevolution::VReverse () {
 123
 124   basisCurve->Reverse();
 125 }
 126
 127
 128 //=======================================================================
 129 //function : VReversedParameter
 130 //purpose  :
 131 //=======================================================================
 132
 133 Standard_Real Geom_SurfaceOfRevolution::VReversedParameter (const Standard_Real V) const {
 134
 135   return basisCurve->ReversedParameter(V);
 136 }
 137
 138
 139 //=======================================================================
 140 //function : Location
 141 //purpose  :
 142 //=======================================================================
 143
 144 const gp_Pnt& Geom_SurfaceOfRevolution::Location () const {
 145
 146   return loc;
 147 }
 148
 149 //=======================================================================
 150 //function : IsUPeriodic
 151 //purpose  :
 152 //=======================================================================
 153
 154 Standard_Boolean Geom_SurfaceOfRevolution::IsUPeriodic () const {
 155
 156   return Standard_True;
 157 }
 158
 159 //=======================================================================
 160 //function : IsCNu
 161 //purpose  :
 162 //=======================================================================
 163
 164 Standard_Boolean Geom_SurfaceOfRevolution::IsCNu (const Standard_Integer ) const  {
 165
 166   return Standard_True;
 167 }
 168
 169 //=======================================================================
 170 //function : Axis
 171 //purpose  :
 172 //=======================================================================
 173
 174 Ax1 Geom_SurfaceOfRevolution::Axis () const  {
 175
 176   return Ax1 (loc, direction);
 177 }
 178
 179 //=======================================================================
 180 //function : IsCNv
 181 //purpose  :
 182 //=======================================================================
 183
 184 Standard_Boolean Geom_SurfaceOfRevolution::IsCNv (const Standard_Integer N) const {
 185
 186   Standard_RangeError_Raise_if (N < 0, " ");
 187   return basisCurve->IsCN(N);
 188 }
 189
 190
 191 //=======================================================================
 192 //function : IsUClosed
 193 //purpose  :
 194 //=======================================================================
 195
 196 Standard_Boolean Geom_SurfaceOfRevolution::IsUClosed () const {
 197
 198   return Standard_True;
 199 }
 200
 201 //=======================================================================
 202 //function : IsVClosed
 203 //purpose  :
 204 //=======================================================================
 205
 206 Standard_Boolean Geom_SurfaceOfRevolution::IsVClosed () const
 207 {
 208   return basisCurve->IsClosed();
 209 }
 210
 211
 212 //=======================================================================
 213 //function : IsVPeriodic
 214 //purpose  :
 215 //=======================================================================
 216
 217 Standard_Boolean Geom_SurfaceOfRevolution::IsVPeriodic () const {
 218
 219   return basisCurve->IsPeriodic();
 220 }
 221
 222
 223 //=======================================================================
 224 //function : SetAxis
 225 //purpose  :
 226 //=======================================================================
 227
 228 void Geom_SurfaceOfRevolution::SetAxis (const Ax1& A1) {
 229
 230    direction = A1.Direction();
 231    loc = A1.Location();
 232    myEvaluator->SetAxis(A1);
 233 }
 234
 235
 236 //=======================================================================
 237 //function : SetDirection
 238 //purpose  :
 239 //=======================================================================
 240
 241 void Geom_SurfaceOfRevolution::SetDirection (const Dir& V) {
 242
 243    direction = V;
 244    myEvaluator->SetDirection(direction);
 245 }
 246
 247
 248 //=======================================================================
 249 //function : SetBasisCurve
 250 //purpose  :
 251 //=======================================================================
 252
 253 void Geom_SurfaceOfRevolution::SetBasisCurve (const Handle(Geom_Curve)& C) {
 254
 255    basisCurve = Handle(Geom_Curve)::DownCast(C->Copy());
 256    smooth     = C->Continuity();
 257    myEvaluator = new GeomEvaluator_SurfaceOfRevolution(basisCurve, direction, loc);
 258 }
 259
 260
 261 //=======================================================================
 262 //function : SetLocation
 263 //purpose  :
 264 //=======================================================================
 265
 266 void Geom_SurfaceOfRevolution::SetLocation (const Pnt& P) {
 267
 268    loc = P;
 269    myEvaluator->SetLocation(loc);
 270 }
 271
 272
 273 //=======================================================================
 274 //function : Bounds
 275 //purpose  :
 276 //=======================================================================
 277
 278 void Geom_SurfaceOfRevolution::Bounds ( Standard_Real& U1,
 279                                         Standard_Real& U2,
 280                                         Standard_Real& V1,
 281                                         Standard_Real& V2 ) const {
 282
 283   U1 = 0.0;
 284   U2 = 2.0 * M_PI;
 285   V1 = basisCurve->FirstParameter();
 286   V2 = basisCurve->LastParameter();
 287 }
 288
 289
 290 //=======================================================================
 291 //function : D0
 292 //purpose  :
 293 //=======================================================================
 294
 295 void Geom_SurfaceOfRevolution::D0
 296 (const Standard_Real U, const Standard_Real V, Pnt& P) const
 297 {
 298   myEvaluator->D0(U, V, P);
 299 }
 300
 301
 302 //=======================================================================
 303 //function : D1
 304 //purpose  :
 305 //=======================================================================
 306
 307 void Geom_SurfaceOfRevolution::D1
 308   (const Standard_Real U, const Standard_Real V,
 309          Pnt& P,
 310          Vec& D1U, Vec& D1V   ) const
 311 {
 312   myEvaluator->D1(U, V, P, D1U, D1V);
 313 }
 314
 315 //=======================================================================
 316 //function : D2
 317 //purpose  :
 318 //=======================================================================
 319
 320 void Geom_SurfaceOfRevolution::D2
 321   (const Standard_Real   U, const Standard_Real V,
 322          Pnt&   P,
 323          Vec& D1U, Vec& D1V,
 324          Vec& D2U, Vec& D2V, Vec& D2UV ) const
 325 {
 326   myEvaluator->D2(U, V, P, D1U, D1V, D2U, D2V, D2UV);
 327 }
 328
 329
 330
 331 //=======================================================================
 332 //function : D3
 333 //purpose  :
 334 //=======================================================================
 335
 336 void Geom_SurfaceOfRevolution::D3
 337   (const Standard_Real U, const Standard_Real V,
 338          Pnt& P,
 339          Vec& D1U, Vec& D1V,
 340          Vec& D2U, Vec& D2V, Vec& D2UV,
 341          Vec& D3U, Vec& D3V, Vec& D3UUV, Vec& D3UVV ) const
 342 {
 343   myEvaluator->D3(U, V, P, D1U, D1V, D2U, D2V, D2UV, D3U, D3V, D3UUV, D3UVV);
 344 }
 345
 346
 347 //=======================================================================
 348 //function : DN
 349 //purpose  :
 350 //=======================================================================
 351
 352 Vec Geom_SurfaceOfRevolution::DN (const Standard_Real    U , const Standard_Real    V,
 353                                   const Standard_Integer Nu, const Standard_Integer Nv) const
 354 {
 355   return myEvaluator->DN(U, V, Nu, Nv);
 356 }
 357
 358
 359 //=======================================================================
 360 //function : ReferencePlane
 361 //purpose  :
 362 //=======================================================================
 363
 364 Ax2 Geom_SurfaceOfRevolution::ReferencePlane() const {
 365
 366    throw Standard_NotImplemented();
 367 }
 368
 369
 370 //=======================================================================
 371 //function : UIso
 372 //purpose  :
 373 //=======================================================================
 374
 375 Handle(Geom_Curve) Geom_SurfaceOfRevolution::UIso (const Standard_Real U) const {
 376
 377    Handle(Geom_Curve) C = Handle(Geom_Curve)::DownCast(basisCurve->Copy());
 378    Ax1 RotAxis = Ax1 (loc, direction);
 379    C->Rotate (RotAxis, U);
 380    return C;
 381 }
 382
 383
 384 //=======================================================================
 385 //function : VIso
 386 //purpose  :
 387 //=======================================================================
 388
 389 Handle(Geom_Curve) Geom_SurfaceOfRevolution::VIso (const Standard_Real V) const {
 390
 391   Handle(Geom_Circle) Circ;
 392   Pnt Pc = basisCurve->Value (V);
 393   gp_Lin L1(loc,direction);
 394   Standard_Real Rad= L1.Distance(Pc);
 395
 396   Ax2 Rep ;
 397   if ( Rad > gp::Resolution()) {
 398     XYZ P  = Pc.XYZ();
 399     XYZ C;
 400     C.SetLinearForm((P-loc.XYZ()).Dot(direction.XYZ()),
 401                     direction.XYZ(), loc.XYZ() );
 402     P = P-C;
 403     if(P.Modulus() > gp::Resolution()) {
 404       gp_Dir D = P.Normalized();
 405       Rep = gp_Ax2(C, direction, D);
 406     }
 407     else
 408       Rep = gp_Ax2(C, direction);
 409   }
 410   else
 411     Rep = gp_Ax2(Pc, direction);
 412
 413   Circ   = new Geom_Circle (Rep, Rad);
 414   return Circ;
 415 }
 416
 417
 418 //=======================================================================
 419 //function : Transform
 420 //purpose  :
 421 //=======================================================================
 422
 423 void Geom_SurfaceOfRevolution::Transform (const Trsf& T) {
 424
 425   loc.Transform (T);
 426   direction.Transform (T);
 427   basisCurve->Transform (T);
 428   if(T.ScaleFactor()*T.HVectorialPart().Determinant() < 0.) UReverse();
 429   myEvaluator->SetDirection(direction);
 430   myEvaluator->SetLocation(loc);
 431 }
 432
 433 //=======================================================================
 434 //function : TransformParameters
 435 //purpose  :
 436 //=======================================================================
 437
 438 void Geom_SurfaceOfRevolution::TransformParameters(Standard_Real& ,
 439                                                    Standard_Real& V,
 440                                                    const gp_Trsf& T)
 441 const
 442 {
 443   V = basisCurve->TransformedParameter(V,T);
 444 }
 445
 446 //=======================================================================
 447 //function : ParametricTransformation
 448 //purpose  :
 449 //=======================================================================
 450
 451 gp_GTrsf2d Geom_SurfaceOfRevolution::ParametricTransformation
 452 (const gp_Trsf& T) const
 453 {
 454   gp_GTrsf2d T2;
 455   gp_Ax2d Axis(gp::Origin2d(),gp::DX2d());
 456   T2.SetAffinity(Axis, basisCurve->ParametricTransformation(T));
 457   return T2;
 458 }
 459