// Copyright (c) 1995-1999 Matra Datavision // Copyright (c) 1999-2014 OPEN CASCADE SAS // // This file is part of Open CASCADE Technology software library. // // This library is free software; you can redistribute it and/or modify it under // the terms of the GNU Lesser General Public License version 2.1 as published // by the Free Software Foundation, with special exception defined in the file // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT // distribution for complete text of the license and disclaimer of any warranty. // // Alternatively, this file may be used under the terms of Open CASCADE // commercial license or contractual agreement. //JCV 16/10/91 #include #include #include #include #include static const Standard_Integer TheUDegree = 2; static const Standard_Integer TheVDegree = 2; static const Standard_Integer MaxNbUKnots = 4; static const Standard_Integer MaxNbVKnots = 4; static const Standard_Integer MaxNbUPoles = 7; static const Standard_Integer MaxNbVPoles = 7; static void ComputePoles ( const Standard_Real R, const Standard_Real r, const Standard_Real U1, const Standard_Real U2, const Standard_Real V1, const Standard_Real V2, TColgp_Array2OfPnt& Poles) { Standard_Real deltaU = U2 - U1; Standard_Real deltaV = V2 - V1; Standard_Integer i, j; // Number of spans : maximum opening = 150 degrees ( = PI / 1.2 rds) Standard_Integer nbUSpans = (Standard_Integer)IntegerPart( 1.2 * deltaU / M_PI) + 1; Standard_Integer nbVSpans = (Standard_Integer)IntegerPart( 1.2 * deltaV / M_PI) + 1; Standard_Real AlfaU = deltaU / ( nbUSpans * 2); Standard_Real AlfaV = deltaV / ( nbVSpans * 2); Standard_Integer nbVP = 2 * nbVSpans + 1; Standard_Real x[MaxNbVPoles]; Standard_Real z[MaxNbVPoles]; x[0] = R + r * Cos( V1); z[0] = r * Sin( V1); Standard_Real VStart = V1; for ( i = 1; i <= nbVSpans; i++) { x[2*i-1] = R + r * Cos( VStart + AlfaV) / Cos( AlfaV); z[2*i-1] = r * Sin( VStart + AlfaV) / Cos( AlfaV); x[2*i] = R + r * Cos( VStart + 2 * AlfaV); z[2*i] = r * Sin( VStart + 2 * AlfaV); VStart += 2*AlfaV; } Standard_Real UStart = U1; for ( j = 0; j <= nbVP-1; j++) { Poles( 1, j+1) = gp_Pnt(x[j]*Cos(UStart),x[j]*Sin(UStart),z[j]); } for ( i = 1; i <= nbUSpans; i++) { for ( j = 0; j<= nbVP-1; j++) { Poles( 2*i, j+1) = gp_Pnt( x[j] * Cos(UStart+AlfaU) / Cos(AlfaU), x[j] * Sin(UStart+AlfaU) / Cos(AlfaU), z[j] ); Poles(2*i+1,j+1) = gp_Pnt( x[j] * Cos(UStart+2*AlfaU), x[j] * Sin(UStart+2*AlfaU), z[j] ); } UStart += 2*AlfaU; } } //======================================================================= //function : Convert_TorusToBSplineSurface //purpose : //======================================================================= Convert_TorusToBSplineSurface::Convert_TorusToBSplineSurface (const gp_Torus& T, const Standard_Real U1, const Standard_Real U2, const Standard_Real V1, const Standard_Real V2) : Convert_ElementarySurfaceToBSplineSurface (MaxNbUPoles, MaxNbVPoles, MaxNbUKnots, MaxNbVKnots, TheUDegree , TheVDegree) { Standard_Real deltaU = U2 - U1; Standard_Real deltaV = V2 - V1; Standard_DomainError_Raise_if( (deltaU>2*M_PI) || (deltaU<0.) || (deltaV>2*M_PI) || (deltaV<0.), "Convert_TorusToBSplineSurface"); isuperiodic = Standard_False; isvperiodic = Standard_False; Standard_Integer i,j; // construction of the torus in the reference mark xOy. // Number of spans : maximum opening = 150 degrees ( = PI / 1.2 rds) Standard_Integer nbUSpans = (Standard_Integer)IntegerPart( 1.2 * deltaU / M_PI) + 1; Standard_Integer nbVSpans = (Standard_Integer)IntegerPart( 1.2 * deltaV / M_PI) + 1; Standard_Real AlfaU = deltaU / ( nbUSpans * 2); Standard_Real AlfaV = deltaV / ( nbVSpans * 2); nbUPoles = 2 * nbUSpans + 1; nbVPoles = 2 * nbVSpans + 1; nbUKnots = nbUSpans + 1; nbVKnots = nbVSpans + 1; Standard_Real R = T.MajorRadius(); Standard_Real r = T.MinorRadius(); ComputePoles( R, r, U1, U2, V1, V2, poles); for ( i = 1; i<= nbUKnots; i++) { uknots(i) = U1 + (i-1) * 2 * AlfaU; umults(i) = 2; } umults(1)++; umults(nbUKnots)++; for ( i = 1; i<= nbVKnots; i++) { vknots(i) = V1 + (i-1) * 2 * AlfaV; vmults(i) = 2; } vmults(1)++; vmults(nbVKnots)++; // Replace the bspline in the reference of the torus. // and calculate the weight of the bspline. Standard_Real W1, W2; gp_Trsf Trsf; Trsf.SetTransformation( T.Position(), gp::XOY()); for ( i = 1; i <= nbUPoles; i++) { if ( i % 2 == 0) W1 = Cos(AlfaU); else W1 = 1.; for ( j = 1; j <= nbVPoles; j++) { if ( j % 2 == 0) W2 = Cos(AlfaV); else W2 = 1.; weights( i, j) = W1 * W2; poles( i, j).Transform( Trsf); } } } //======================================================================= //function : Convert_TorusToBSplineSurface //purpose : //======================================================================= Convert_TorusToBSplineSurface::Convert_TorusToBSplineSurface (const gp_Torus& T, const Standard_Real Param1, const Standard_Real Param2, const Standard_Boolean UTrim ) : Convert_ElementarySurfaceToBSplineSurface (MaxNbUPoles, MaxNbVPoles, MaxNbUKnots, MaxNbVKnots, TheUDegree , TheVDegree) { #ifndef No_Exception Standard_Real delta = Param2 - Param1; #endif Standard_DomainError_Raise_if( (delta>2*M_PI) || (delta<0.), "Convert_TorusToBSplineSurface"); Standard_Integer i, j; Standard_Real deltaU, deltaV; isuperiodic = !UTrim; isvperiodic = UTrim; Standard_Real R = T.MajorRadius(); Standard_Real r = T.MinorRadius(); Standard_Real W1, W2, CosU, CosV; if ( isuperiodic) { ComputePoles(R, r, 0, 2.*M_PI, Param1, Param2, poles); nbUPoles = 6; nbUKnots = 4; deltaV = Param2 - Param1; Standard_Integer nbVSpans = (Standard_Integer)IntegerPart( 1.2 * deltaV / M_PI) + 1; Standard_Real AlfaV = deltaV / ( nbVSpans * 2); nbVPoles = 2 * nbVSpans + 1; nbVKnots = nbVSpans + 1; for ( i = 1; i <= nbUKnots; i++) { uknots(i) = ( i-1) * 2. * M_PI /3.; umults(i) = 2; } for ( i = 1; i <= nbVKnots; i++) { vknots(i) = Param1 + (i-1) * 2 * AlfaV; vmults(i) = 2; } vmults(1)++; vmults(nbVKnots)++; CosU = 0.5; // = Cos(pi /3) CosV = Cos(AlfaV); } else { ComputePoles(R, r, Param1, Param2, 0., 2.*M_PI, poles); nbVPoles = 6; nbVKnots = 4; deltaU = Param2 - Param1; Standard_Integer nbUSpans = (Standard_Integer)IntegerPart( 1.2 * deltaU / M_PI) + 1; Standard_Real AlfaU = deltaU / ( nbUSpans * 2); nbUPoles = 2 * nbUSpans + 1; nbUKnots = nbUSpans + 1; for ( i = 1; i <= nbVKnots; i++) { vknots(i) = ( i-1) * 2. * M_PI /3.; vmults(i) = 2; } for ( i = 1; i <= nbUKnots; i++) { uknots(i) = Param1 + (i-1) * 2 * AlfaU; umults(i) = 2; } umults(1)++; umults(nbUKnots)++; CosV = 0.5; // = Cos(pi /3) CosU = Cos(AlfaU); } // Replace the bspline in the reference of the torus. // and calculate the weight of the bspline. gp_Trsf Trsf; Trsf.SetTransformation( T.Position(), gp::XOY()); for ( i = 1; i <= nbUPoles; i++) { if ( i % 2 == 0) W1 = CosU; else W1 = 1.; for ( j = 1; j <= nbVPoles; j++) { if ( j % 2 == 0) W2 = CosV; else W2 = 1.; weights( i, j) = W1 * W2; poles( i, j).Transform( Trsf); } } } //======================================================================= //function : Convert_TorusToBSplineSurface //purpose : //======================================================================= Convert_TorusToBSplineSurface::Convert_TorusToBSplineSurface (const gp_Torus& T ) : Convert_ElementarySurfaceToBSplineSurface (MaxNbUPoles, MaxNbVPoles, MaxNbUKnots, MaxNbVKnots, TheUDegree , TheVDegree ) { isuperiodic = Standard_True; isvperiodic = Standard_True; Standard_Real W1, W2; Standard_Integer i, j; nbUPoles = 6; nbVPoles = 6; nbUKnots = 4; nbVKnots = 4; // Construction of the Torus in the reference mark xOy. Standard_Real R = T.MajorRadius(); Standard_Real r = T.MinorRadius(); ComputePoles( R, r, 0., 2.*M_PI, 0., 2.*M_PI, poles); uknots( 1) = vknots( 1) = 0.; uknots( 2) = vknots( 2) = 2. * M_PI / 3.; uknots( 3) = vknots( 3) = 4. * M_PI / 3.; uknots( 4) = vknots( 4) = 2. * M_PI; for ( i = 1; i <= 4; i++) { umults( i) = vmults( i) = 2; } // Replace the bspline in the mark of the torus. // and calculate the weight of the bspline. gp_Trsf Trsf; Trsf.SetTransformation( T.Position(), gp::XOY()); for ( i = 1; i <= nbUPoles; i++) { if ( i % 2 == 0) W1 = 0.5; else W1 = 1.; for ( j = 1; j <= nbVPoles; j++) { if ( j % 2 == 0) W2 = 0.5; else W2 = 1.; weights( i, j) = W1 * W2; poles( i, j).Transform( Trsf); } } }