1 // Created on: 1998-02-26
2 // Created by: Roman BORISOV
3 // Copyright (c) 1998-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
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.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
18 #include <Adaptor3d_Curve.hxx>
19 #include <GeomFill_SnglrFunc.hxx>
22 #include <Precision.hxx>
23 #include <Standard_DomainError.hxx>
24 #include <Standard_NotImplemented.hxx>
25 #include <Standard_OutOfRange.hxx>
27 GeomFill_SnglrFunc::GeomFill_SnglrFunc(const Handle(Adaptor3d_Curve)& HC) :
28 myHCurve(HC), ratio(1)
32 //=======================================================================
33 //function : ShallowCopy
35 //=======================================================================
37 Handle(Adaptor3d_Curve) GeomFill_SnglrFunc::ShallowCopy() const
39 Handle(GeomFill_SnglrFunc) aCopy = new GeomFill_SnglrFunc(myHCurve->ShallowCopy());
44 void GeomFill_SnglrFunc::SetRatio(const Standard_Real Ratio)
49 Standard_Real GeomFill_SnglrFunc::FirstParameter() const
51 return myHCurve->FirstParameter();
54 Standard_Real GeomFill_SnglrFunc::LastParameter() const
56 return myHCurve->LastParameter();
59 Standard_Integer GeomFill_SnglrFunc::NbIntervals(const GeomAbs_Shape S) const
61 GeomAbs_Shape HCS=GeomAbs_C0;
66 else if (S == GeomAbs_C1)
70 else if (S >= GeomAbs_C2)
74 return myHCurve->NbIntervals(HCS);
77 void GeomFill_SnglrFunc::Intervals(TColStd_Array1OfReal& T,const GeomAbs_Shape S) const
79 GeomAbs_Shape HCS=GeomAbs_C0;
84 else if (S == GeomAbs_C1)
88 else if (S >= GeomAbs_C2)
92 myHCurve->Intervals(T, HCS);
95 Standard_Boolean GeomFill_SnglrFunc::IsPeriodic() const
97 return myHCurve->IsPeriodic();
100 Standard_Real GeomFill_SnglrFunc::Period() const
102 return myHCurve->Period();
106 gp_Pnt GeomFill_SnglrFunc::Value(const Standard_Real U) const
110 myHCurve->D2(U, C, DC, D2C);
112 return gp_Pnt(DC.Crossed(D2C).XYZ());
115 void GeomFill_SnglrFunc::D0(const Standard_Real U,gp_Pnt& P) const
119 myHCurve->D2(U, C, DC, D2C);
121 P = gp_Pnt(DC.Crossed(D2C).XYZ());
124 void GeomFill_SnglrFunc::D1(const Standard_Real U,gp_Pnt& P,gp_Vec& V) const
128 myHCurve->D3(U, C, DC, D2C, D3C);
130 P = gp_Pnt(DC.Crossed(D2C).XYZ());
134 void GeomFill_SnglrFunc::D2(const Standard_Real U,gp_Pnt& P,gp_Vec& V1,gp_Vec& V2) const
137 gp_Vec DC, D2C, D3C, D4C;
138 myHCurve->D3(U, C, DC, D2C, D3C);
139 P = gp_Pnt(DC.Crossed(D2C).XYZ());
140 V1 = DC.Crossed(D3C);
141 D4C = myHCurve->DN(U, 4);
142 V2 = D2C.Crossed(D3C) + DC.Crossed(D4C);
144 P.ChangeCoord() *= ratio;
149 void GeomFill_SnglrFunc::D3(const Standard_Real U,gp_Pnt& P,gp_Vec& V1,gp_Vec& V2,gp_Vec& V3) const
151 gp_Vec DC, D2C, D3C, D4C, D5C;
152 myHCurve->D3(U, P, DC, D2C, D3C);
153 D4C = myHCurve->DN(U, 4);
154 D5C = myHCurve->DN(U, 5);
155 P = gp_Pnt(DC.Crossed(D2C).XYZ()).ChangeCoord()*ratio;
156 V1 = DC.Crossed(D3C)*ratio;
157 V2 = (D2C.Crossed(D3C) + DC.Crossed(D4C))*ratio;
158 V3 = (DC.Crossed(D5C) + D2C.Crossed(D4C)*2)*ratio;
161 gp_Vec GeomFill_SnglrFunc::DN(const Standard_Real U,const Standard_Integer N) const
163 Standard_RangeError_Raise_if (N < 1, "Exception: Geom2d_OffsetCurve::DN(). N<1.");
165 gp_Vec D1C, D2C, D3C;
180 throw Standard_NotImplemented("Exception: Derivative order is greater than 3. "
181 "Cannot compute of derivative.");
185 Standard_Real GeomFill_SnglrFunc::Resolution(const Standard_Real R3D) const
187 return Precision::Parametric(R3D);
190 GeomAbs_CurveType GeomFill_SnglrFunc::GetType() const
192 return GeomAbs_OtherCurve;