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.
17 #include <GeomFill_SnglrFunc.ixx>
18 #include <Standard_NotImplemented.hxx>
19 #include <Precision.hxx>
21 GeomFill_SnglrFunc::GeomFill_SnglrFunc(const Handle(Adaptor3d_HCurve)& HC) :
22 myHCurve(HC), ratio(1)
26 void GeomFill_SnglrFunc::SetRatio(const Standard_Real Ratio)
31 Standard_Real GeomFill_SnglrFunc::FirstParameter() const
33 return myHCurve->FirstParameter();
36 Standard_Real GeomFill_SnglrFunc::LastParameter() const
38 return myHCurve->LastParameter();
41 Standard_Integer GeomFill_SnglrFunc::NbIntervals(const GeomAbs_Shape S)
43 GeomAbs_Shape HCS=GeomAbs_C0;
45 case GeomAbs_C0: HCS = GeomAbs_C2; break;
46 case GeomAbs_C1: HCS = GeomAbs_C3; break;
47 case GeomAbs_C2: HCS = GeomAbs_CN; break;
48 default: Standard_DomainError::Raise();
50 return myHCurve->NbIntervals(HCS);
53 void GeomFill_SnglrFunc::Intervals(TColStd_Array1OfReal& T,const GeomAbs_Shape S)
55 GeomAbs_Shape HCS=GeomAbs_C0;
57 case GeomAbs_C0: HCS = GeomAbs_C2; break;
58 case GeomAbs_C1: HCS = GeomAbs_C3; break;
59 case GeomAbs_C2: HCS = GeomAbs_CN; break;
60 default: Standard_DomainError::Raise();
62 myHCurve->Intervals(T, HCS);
65 Standard_Boolean GeomFill_SnglrFunc::IsPeriodic() const
67 return myHCurve->IsPeriodic();
70 Standard_Real GeomFill_SnglrFunc::Period() const
72 return myHCurve->Period();
76 gp_Pnt GeomFill_SnglrFunc::Value(const Standard_Real U) const
80 myHCurve->D2(U, C, DC, D2C);
82 return gp_Pnt(DC.Crossed(D2C).XYZ());
85 void GeomFill_SnglrFunc::D0(const Standard_Real U,gp_Pnt& P) const
89 myHCurve->D2(U, C, DC, D2C);
91 P = gp_Pnt(DC.Crossed(D2C).XYZ());
94 void GeomFill_SnglrFunc::D1(const Standard_Real U,gp_Pnt& P,gp_Vec& V) const
98 myHCurve->D3(U, C, DC, D2C, D3C);
100 P = gp_Pnt(DC.Crossed(D2C).XYZ());
104 void GeomFill_SnglrFunc::D2(const Standard_Real U,gp_Pnt& P,gp_Vec& V1,gp_Vec& V2) const
107 gp_Vec DC, D2C, D3C, D4C;
108 myHCurve->D3(U, C, DC, D2C, D3C);
109 P = gp_Pnt(DC.Crossed(D2C).XYZ());
110 V1 = DC.Crossed(D3C);
111 D4C = myHCurve->DN(U, 4);
112 V2 = D2C.Crossed(D3C) + DC.Crossed(D4C);
114 P.ChangeCoord() *= ratio;
119 void GeomFill_SnglrFunc::D3(const Standard_Real U,gp_Pnt& P,gp_Vec& V1,gp_Vec& V2,gp_Vec& V3) const
121 gp_Vec DC, D2C, D3C, D4C, D5C;
122 myHCurve->D3(U, P, DC, D2C, D3C);
123 D4C = myHCurve->DN(U, 4);
124 D5C = myHCurve->DN(U, 5);
125 P = gp_Pnt(DC.Crossed(D2C).XYZ()).ChangeCoord()*ratio;
126 V1 = DC.Crossed(D3C)*ratio;
127 V2 = (D2C.Crossed(D3C) + DC.Crossed(D4C))*ratio;
128 V3 = (DC.Crossed(D5C) + D2C.Crossed(D4C)*2)*ratio;
131 gp_Vec GeomFill_SnglrFunc::DN(const Standard_Real U,const Standard_Integer N) const
133 Standard_RangeError_Raise_if (N < 1, "Exception: Geom2d_OffsetCurve::DN(). N<1.");
135 gp_Vec D1C, D2C, D3C;
150 Standard_NotImplemented::Raise("Exception: Derivative order is greater than 3. "
151 "Cannot compute of derivative.");
158 Standard_Real GeomFill_SnglrFunc::Resolution(const Standard_Real R3D) const
160 return Precision::Parametric(R3D);
163 GeomAbs_CurveType GeomFill_SnglrFunc::GetType() const
165 return GeomAbs_OtherCurve;