b311480e |
1 | // Created on: 1993-05-14 |
2 | // Created by: Joelle CHAUVET |
3 | // Copyright (c) 1993-1999 Matra Datavision |
973c2be1 |
4 | // Copyright (c) 1999-2014 OPEN CASCADE SAS |
b311480e |
5 | // |
973c2be1 |
6 | // This file is part of Open CASCADE Technology software library. |
b311480e |
7 | // |
d5f74e42 |
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 |
973c2be1 |
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. |
b311480e |
13 | // |
973c2be1 |
14 | // Alternatively, this file may be used under the terms of Open CASCADE |
15 | // commercial license or contractual agreement. |
7fd59977 |
16 | |
17 | // Modified: Thu Nov 26 16:37:18 1998 |
7fd59977 |
18 | // correction in NbUIntervals for SurfaceOfLinearExtrusion |
19 | // (PRO16346) |
20 | |
21 | #define No_Standard_RangeError |
22 | #define No_Standard_OutOfRange |
23 | #define PosTol (Precision::PConfusion()*0.5) |
24 | |
7fd59977 |
25 | |
42cf5bc1 |
26 | #include <Adaptor3d_HCurve.hxx> |
7fd59977 |
27 | #include <Adaptor3d_HSurface.hxx> |
42cf5bc1 |
28 | #include <BSplCLib.hxx> |
29 | #include <BSplSLib_Cache.hxx> |
30 | #include <Geom_BezierSurface.hxx> |
31 | #include <Geom_BSplineSurface.hxx> |
32 | #include <Geom_Circle.hxx> |
33 | #include <Geom_ConicalSurface.hxx> |
34 | #include <Geom_Curve.hxx> |
35 | #include <Geom_CylindricalSurface.hxx> |
7fd59977 |
36 | #include <Geom_OffsetSurface.hxx> |
7fd59977 |
37 | #include <Geom_Plane.hxx> |
42cf5bc1 |
38 | #include <Geom_RectangularTrimmedSurface.hxx> |
7fd59977 |
39 | #include <Geom_SphericalSurface.hxx> |
42cf5bc1 |
40 | #include <Geom_Surface.hxx> |
7fd59977 |
41 | #include <Geom_SurfaceOfLinearExtrusion.hxx> |
42cf5bc1 |
42 | #include <Geom_SurfaceOfRevolution.hxx> |
43 | #include <Geom_ToroidalSurface.hxx> |
44 | #include <GeomAdaptor_Curve.hxx> |
45 | #include <GeomAdaptor_HCurve.hxx> |
46 | #include <GeomAdaptor_HSurface.hxx> |
47 | #include <GeomAdaptor_Surface.hxx> |
6b84c3f7 |
48 | #include <GeomEvaluator_OffsetSurface.hxx> |
49 | #include <GeomEvaluator_SurfaceOfExtrusion.hxx> |
50 | #include <GeomEvaluator_SurfaceOfRevolution.hxx> |
42cf5bc1 |
51 | #include <gp_Ax1.hxx> |
7fd59977 |
52 | #include <gp_Circ.hxx> |
42cf5bc1 |
53 | #include <gp_Cone.hxx> |
54 | #include <gp_Cylinder.hxx> |
55 | #include <gp_Dir.hxx> |
7fd59977 |
56 | #include <gp_Lin.hxx> |
42cf5bc1 |
57 | #include <gp_Pln.hxx> |
58 | #include <gp_Pnt.hxx> |
59 | #include <gp_Sphere.hxx> |
60 | #include <gp_Torus.hxx> |
7fd59977 |
61 | #include <gp_Trsf.hxx> |
42cf5bc1 |
62 | #include <gp_Vec.hxx> |
7fd59977 |
63 | #include <Precision.hxx> |
42cf5bc1 |
64 | #include <Standard_ConstructionError.hxx> |
65 | #include <Standard_DomainError.hxx> |
7fd59977 |
66 | #include <Standard_NoSuchObject.hxx> |
041bfce9 |
67 | #include <Standard_NullObject.hxx> |
42cf5bc1 |
68 | #include <Standard_OutOfRange.hxx> |
69 | #include <TColStd_Array1OfInteger.hxx> |
70 | #include <TColStd_Array1OfReal.hxx> |
71 | #include <TColStd_HArray1OfInteger.hxx> |
7fd59977 |
72 | |
7fd59977 |
73 | //======================================================================= |
74 | //function : LocalContinuity |
75 | //purpose : |
76 | //======================================================================= |
77 | |
78 | GeomAbs_Shape LocalContinuity(Standard_Integer Degree, |
79 | Standard_Integer Nb, |
80 | TColStd_Array1OfReal& TK, |
81 | TColStd_Array1OfInteger& TM, |
82 | Standard_Real PFirst, |
83 | Standard_Real PLast, |
84 | Standard_Boolean IsPeriodic) |
85 | { |
86 | Standard_DomainError_Raise_if( (TK.Length()!=Nb || TM.Length()!=Nb )," "); |
87 | Standard_Integer Index1 = 0; |
88 | Standard_Integer Index2 = 0; |
89 | Standard_Real newFirst, newLast; |
90 | BSplCLib::LocateParameter(Degree,TK,TM,PFirst,IsPeriodic,1,Nb,Index1,newFirst); |
91 | BSplCLib::LocateParameter(Degree,TK,TM,PLast, IsPeriodic,1,Nb,Index2,newLast ); |
92 | const Standard_Real EpsKnot = Precision::PConfusion(); |
93 | if (Abs(newFirst-TK(Index1+1))< EpsKnot) Index1++; |
94 | if (Abs(newLast -TK(Index2 ))< EpsKnot) Index2--; |
95 | // attention aux courbes peridiques. |
96 | if ( (IsPeriodic) && (Index1 == Nb) ) |
97 | Index1 = 1; |
98 | |
99 | if (Index2!=Index1) |
100 | { |
101 | Standard_Integer i, Multmax = TM(Index1+1); |
102 | for (i = Index1+1; i<=Index2; i++) { |
103 | if (TM(i)>Multmax) Multmax=TM(i); |
104 | } |
105 | Multmax = Degree - Multmax; |
106 | if (Multmax <= 0) return GeomAbs_C0; |
107 | switch (Multmax) { |
108 | case 1: return GeomAbs_C1; |
109 | case 2: return GeomAbs_C2; |
110 | case 3: return GeomAbs_C3; |
111 | } |
112 | } |
113 | return GeomAbs_CN; |
114 | } |
115 | |
116 | //======================================================================= |
117 | //function : Load |
118 | //purpose : |
119 | //======================================================================= |
120 | |
041bfce9 |
121 | void GeomAdaptor_Surface::load(const Handle(Geom_Surface)& S, |
7fd59977 |
122 | const Standard_Real UFirst, |
123 | const Standard_Real ULast, |
124 | const Standard_Real VFirst, |
125 | const Standard_Real VLast, |
126 | const Standard_Real TolU, |
127 | const Standard_Real TolV) |
128 | { |
7fd59977 |
129 | myTolU = TolU; |
130 | myTolV = TolV; |
131 | myUFirst = UFirst; |
132 | myULast = ULast; |
133 | myVFirst = VFirst; |
134 | myVLast = VLast; |
135 | |
136 | if ( mySurface != S) { |
137 | mySurface = S; |
6b84c3f7 |
138 | mySurfaceCache = Handle(BSplSLib_Cache)(); |
139 | myNestedEvaluator = Handle(GeomEvaluator_Surface)(); |
140 | |
7fd59977 |
141 | const Handle(Standard_Type)& TheType = S->DynamicType(); |
142 | if ( TheType == STANDARD_TYPE(Geom_BezierSurface)) |
143 | mySurfaceType = GeomAbs_BezierSurface; |
144 | else if (TheType == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) { |
c5f3a425 |
145 | Load(Handle(Geom_RectangularTrimmedSurface)::DownCast (S)->BasisSurface(), |
94f71cad |
146 | UFirst,ULast,VFirst,VLast); |
7fd59977 |
147 | } |
148 | else if ( TheType == STANDARD_TYPE(Geom_Plane)) |
149 | mySurfaceType = GeomAbs_Plane; |
150 | else if ( TheType == STANDARD_TYPE(Geom_CylindricalSurface)) |
151 | mySurfaceType = GeomAbs_Cylinder; |
152 | else if ( TheType == STANDARD_TYPE(Geom_ConicalSurface)) |
153 | mySurfaceType = GeomAbs_Cone; |
154 | else if ( TheType == STANDARD_TYPE(Geom_SphericalSurface)) |
155 | mySurfaceType = GeomAbs_Sphere; |
156 | else if ( TheType == STANDARD_TYPE(Geom_ToroidalSurface)) |
157 | mySurfaceType = GeomAbs_Torus; |
6b84c3f7 |
158 | else if (TheType == STANDARD_TYPE(Geom_SurfaceOfRevolution)) |
159 | { |
7fd59977 |
160 | mySurfaceType = GeomAbs_SurfaceOfRevolution; |
6b84c3f7 |
161 | Handle(Geom_SurfaceOfRevolution) myRevSurf = |
162 | Handle(Geom_SurfaceOfRevolution)::DownCast(mySurface); |
163 | // Create nested adaptor for base curve |
164 | Handle(Geom_Curve) aBaseCurve = myRevSurf->BasisCurve(); |
165 | Handle(GeomAdaptor_HCurve) aBaseAdaptor = new GeomAdaptor_HCurve(aBaseCurve); |
166 | // Create corresponding evaluator |
167 | myNestedEvaluator = new GeomEvaluator_SurfaceOfRevolution( |
168 | Handle(Adaptor3d_HCurve)::DownCast(aBaseAdaptor), |
169 | myRevSurf->Direction(), myRevSurf->Location()); |
170 | } |
171 | else if (TheType == STANDARD_TYPE(Geom_SurfaceOfLinearExtrusion)) |
172 | { |
7fd59977 |
173 | mySurfaceType = GeomAbs_SurfaceOfExtrusion; |
6b84c3f7 |
174 | Handle(Geom_SurfaceOfLinearExtrusion) myExtSurf = |
175 | Handle(Geom_SurfaceOfLinearExtrusion)::DownCast(mySurface); |
176 | // Create nested adaptor for base curve |
177 | Handle(Geom_Curve) aBaseCurve = myExtSurf->BasisCurve(); |
178 | Handle(GeomAdaptor_HCurve) aBaseAdaptor = new GeomAdaptor_HCurve(aBaseCurve); |
179 | // Create corresponding evaluator |
180 | myNestedEvaluator = new GeomEvaluator_SurfaceOfExtrusion( |
181 | Handle(Adaptor3d_HCurve)::DownCast(aBaseAdaptor), myExtSurf->Direction()); |
182 | } |
7fd59977 |
183 | else if ( TheType == STANDARD_TYPE(Geom_BSplineSurface)) { |
184 | mySurfaceType = GeomAbs_BSplineSurface; |
6b84c3f7 |
185 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
94f71cad |
186 | // Create cache for B-spline |
6b84c3f7 |
187 | mySurfaceCache = new BSplSLib_Cache( |
94f71cad |
188 | myBspl->UDegree(), myBspl->IsUPeriodic(), myBspl->UKnotSequence(), |
189 | myBspl->VDegree(), myBspl->IsVPeriodic(), myBspl->VKnotSequence(), |
190 | myBspl->Poles(), myBspl->Weights()); |
7fd59977 |
191 | } |
6b84c3f7 |
192 | else if (TheType == STANDARD_TYPE(Geom_OffsetSurface)) |
193 | { |
7fd59977 |
194 | mySurfaceType = GeomAbs_OffsetSurface; |
6b84c3f7 |
195 | Handle(Geom_OffsetSurface) myOffSurf = Handle(Geom_OffsetSurface)::DownCast(mySurface); |
196 | // Create nested adaptor for base surface |
197 | Handle(Geom_Surface) aBaseSurf = myOffSurf->BasisSurface(); |
198 | Handle(GeomAdaptor_HSurface) aBaseAdaptor = |
199 | new GeomAdaptor_HSurface(aBaseSurf, myUFirst, myULast, myVFirst, myVLast, myTolU, myTolV); |
200 | myNestedEvaluator = new GeomEvaluator_OffsetSurface( |
201 | aBaseAdaptor, myOffSurf->Offset(), myOffSurf->OsculatingSurface()); |
202 | } |
7fd59977 |
203 | else |
204 | mySurfaceType = GeomAbs_OtherSurface; |
205 | } |
206 | } |
207 | |
208 | // -- |
209 | // -- Global methods - Apply to the whole Surface. |
210 | // -- |
211 | |
212 | |
213 | //======================================================================= |
214 | //function : UContinuity |
215 | //purpose : |
216 | //======================================================================= |
217 | |
218 | GeomAbs_Shape GeomAdaptor_Surface::UContinuity() const |
219 | { |
220 | switch (mySurfaceType) |
221 | { |
222 | case GeomAbs_BSplineSurface: |
223 | { |
6b84c3f7 |
224 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
7fd59977 |
225 | const Standard_Integer N = myBspl->NbUKnots(); |
226 | TColStd_Array1OfReal TK(1,N); |
227 | TColStd_Array1OfInteger TM(1,N); |
228 | myBspl->UKnots(TK); |
229 | myBspl->UMultiplicities(TM); |
230 | return LocalContinuity(myBspl->UDegree(), myBspl->NbUKnots(), TK, TM, |
231 | myUFirst, myULast, IsUPeriodic()); |
232 | } |
566f8441 |
233 | case GeomAbs_OffsetSurface: |
7fd59977 |
234 | { |
235 | switch(BasisSurface()->UContinuity()) |
236 | { |
566f8441 |
237 | case GeomAbs_CN : |
238 | case GeomAbs_C3 : return GeomAbs_CN; |
239 | case GeomAbs_G2 : |
240 | case GeomAbs_C2 : return GeomAbs_C1; |
241 | case GeomAbs_G1 : |
242 | case GeomAbs_C1 : return GeomAbs_C0; |
243 | case GeomAbs_C0 : break; |
7fd59977 |
244 | } |
245 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::UContinuity"); |
566f8441 |
246 | break; |
7fd59977 |
247 | } |
566f8441 |
248 | case GeomAbs_SurfaceOfExtrusion: |
7fd59977 |
249 | { |
6b84c3f7 |
250 | Handle(Geom_SurfaceOfLinearExtrusion) myExtSurf = |
251 | Handle(Geom_SurfaceOfLinearExtrusion)::DownCast(mySurface); |
252 | GeomAdaptor_Curve GC(myExtSurf->BasisCurve(), myUFirst, myULast); |
7fd59977 |
253 | return GC.Continuity(); |
254 | } |
566f8441 |
255 | case GeomAbs_OtherSurface: |
256 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::UContinuity"); |
257 | case GeomAbs_Plane: |
258 | case GeomAbs_Cylinder: |
259 | case GeomAbs_Cone: |
260 | case GeomAbs_Sphere: |
261 | case GeomAbs_Torus: |
262 | case GeomAbs_BezierSurface: |
263 | case GeomAbs_SurfaceOfRevolution: break; |
7fd59977 |
264 | } |
265 | return GeomAbs_CN; |
266 | } |
267 | |
268 | //======================================================================= |
269 | //function : VContinuity |
270 | //purpose : |
271 | //======================================================================= |
272 | |
273 | GeomAbs_Shape GeomAdaptor_Surface::VContinuity() const |
274 | { |
275 | switch (mySurfaceType) |
276 | { |
277 | case GeomAbs_BSplineSurface: |
278 | { |
6b84c3f7 |
279 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
7fd59977 |
280 | const Standard_Integer N = myBspl->NbVKnots(); |
281 | TColStd_Array1OfReal TK(1,N); |
282 | TColStd_Array1OfInteger TM(1,N); |
283 | myBspl->VKnots(TK); |
284 | myBspl->VMultiplicities(TM); |
285 | return LocalContinuity(myBspl->VDegree(), myBspl->NbVKnots(), TK, TM, |
286 | myVFirst, myVLast, IsVPeriodic()); |
287 | } |
566f8441 |
288 | case GeomAbs_OffsetSurface: |
7fd59977 |
289 | { |
290 | switch(BasisSurface()->VContinuity()) |
291 | { |
566f8441 |
292 | case GeomAbs_CN : |
293 | case GeomAbs_C3 : return GeomAbs_CN; |
294 | case GeomAbs_G2 : |
295 | case GeomAbs_C2 : return GeomAbs_C1; |
296 | case GeomAbs_G1 : |
297 | case GeomAbs_C1 : return GeomAbs_C0; |
298 | case GeomAbs_C0 : break; |
7fd59977 |
299 | } |
300 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::VContinuity"); |
301 | break; |
302 | } |
566f8441 |
303 | case GeomAbs_SurfaceOfRevolution: |
7fd59977 |
304 | { |
6b84c3f7 |
305 | Handle(Geom_SurfaceOfRevolution) myRevSurf = |
306 | Handle(Geom_SurfaceOfRevolution)::DownCast(mySurface); |
307 | GeomAdaptor_Curve GC(myRevSurf->BasisCurve(), myVFirst, myVLast); |
7fd59977 |
308 | return GC.Continuity(); |
309 | } |
566f8441 |
310 | case GeomAbs_OtherSurface: |
311 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::VContinuity"); |
312 | case GeomAbs_Plane: |
313 | case GeomAbs_Cylinder: |
314 | case GeomAbs_Cone: |
315 | case GeomAbs_Sphere: |
316 | case GeomAbs_Torus: |
317 | case GeomAbs_BezierSurface: |
318 | case GeomAbs_SurfaceOfExtrusion: break; |
7fd59977 |
319 | } |
320 | return GeomAbs_CN; |
321 | } |
322 | |
323 | //======================================================================= |
324 | //function : NbUIntervals |
325 | //purpose : |
326 | //======================================================================= |
327 | |
328 | Standard_Integer GeomAdaptor_Surface::NbUIntervals(const GeomAbs_Shape S) const |
329 | { |
330 | switch (mySurfaceType) |
331 | { |
332 | case GeomAbs_BSplineSurface: |
333 | { |
6b84c3f7 |
334 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
7fd59977 |
335 | GeomAdaptor_Curve myBasisCurve |
336 | (myBspl->VIso(myBspl->VKnot(myBspl->FirstVKnotIndex())),myUFirst,myULast); |
337 | return myBasisCurve.NbIntervals(S); |
338 | } |
566f8441 |
339 | case GeomAbs_SurfaceOfExtrusion: |
7fd59977 |
340 | { |
6b84c3f7 |
341 | Handle(Geom_SurfaceOfLinearExtrusion) myExtSurf = |
342 | Handle(Geom_SurfaceOfLinearExtrusion)::DownCast(mySurface); |
343 | GeomAdaptor_Curve myBasisCurve(myExtSurf->BasisCurve(), myUFirst, myULast); |
344 | if (myBasisCurve.GetType() == GeomAbs_BSplineCurve) |
7fd59977 |
345 | return myBasisCurve.NbIntervals(S); |
346 | break; |
347 | } |
566f8441 |
348 | case GeomAbs_OffsetSurface: |
7fd59977 |
349 | { |
350 | GeomAbs_Shape BaseS = GeomAbs_CN; |
351 | switch(S) |
352 | { |
353 | case GeomAbs_G1: |
354 | case GeomAbs_G2: Standard_DomainError::Raise("GeomAdaptor_Curve::NbUIntervals"); |
355 | case GeomAbs_C0: BaseS = GeomAbs_C1; break; |
356 | case GeomAbs_C1: BaseS = GeomAbs_C2; break; |
357 | case GeomAbs_C2: BaseS = GeomAbs_C3; break; |
566f8441 |
358 | case GeomAbs_C3: |
359 | case GeomAbs_CN: break; |
7fd59977 |
360 | } |
6b84c3f7 |
361 | Handle(Geom_OffsetSurface) myOffSurf = Handle(Geom_OffsetSurface)::DownCast(mySurface); |
362 | GeomAdaptor_Surface Sur(myOffSurf->BasisSurface()); |
7fd59977 |
363 | return Sur.NbUIntervals(BaseS); |
364 | } |
566f8441 |
365 | case GeomAbs_Plane: |
366 | case GeomAbs_Cylinder: |
367 | case GeomAbs_Cone: |
368 | case GeomAbs_Sphere: |
369 | case GeomAbs_Torus: |
370 | case GeomAbs_BezierSurface: |
371 | case GeomAbs_OtherSurface: |
372 | case GeomAbs_SurfaceOfRevolution: break; |
7fd59977 |
373 | } |
374 | return 1; |
375 | } |
376 | |
377 | //======================================================================= |
378 | //function : NbVIntervals |
379 | //purpose : |
380 | //======================================================================= |
381 | |
382 | Standard_Integer GeomAdaptor_Surface::NbVIntervals(const GeomAbs_Shape S) const |
383 | { |
384 | switch (mySurfaceType) |
385 | { |
386 | case GeomAbs_BSplineSurface: |
387 | { |
6b84c3f7 |
388 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
7fd59977 |
389 | GeomAdaptor_Curve myBasisCurve |
390 | (myBspl->UIso(myBspl->UKnot(myBspl->FirstUKnotIndex())),myVFirst,myVLast); |
391 | return myBasisCurve.NbIntervals(S); |
392 | } |
6b84c3f7 |
393 | case GeomAbs_SurfaceOfRevolution: |
7fd59977 |
394 | { |
6b84c3f7 |
395 | Handle(Geom_SurfaceOfRevolution) myRevSurf = |
396 | Handle(Geom_SurfaceOfRevolution)::DownCast(mySurface); |
397 | GeomAdaptor_Curve myBasisCurve(myRevSurf->BasisCurve(), myVFirst, myVLast); |
398 | if (myBasisCurve.GetType() == GeomAbs_BSplineCurve) |
7fd59977 |
399 | return myBasisCurve.NbIntervals(S); |
400 | break; |
401 | } |
6b84c3f7 |
402 | case GeomAbs_OffsetSurface: |
7fd59977 |
403 | { |
404 | GeomAbs_Shape BaseS = GeomAbs_CN; |
405 | switch(S) |
406 | { |
407 | case GeomAbs_G1: |
408 | case GeomAbs_G2: Standard_DomainError::Raise("GeomAdaptor_Curve::NbVIntervals"); |
409 | case GeomAbs_C0: BaseS = GeomAbs_C1; break; |
410 | case GeomAbs_C1: BaseS = GeomAbs_C2; break; |
411 | case GeomAbs_C2: BaseS = GeomAbs_C3; break; |
566f8441 |
412 | case GeomAbs_C3: |
413 | case GeomAbs_CN: break; |
7fd59977 |
414 | } |
6b84c3f7 |
415 | Handle(Geom_OffsetSurface) myOffSurf = Handle(Geom_OffsetSurface)::DownCast(mySurface); |
416 | GeomAdaptor_Surface Sur(myOffSurf->BasisSurface()); |
7fd59977 |
417 | return Sur.NbVIntervals(BaseS); |
6b84c3f7 |
418 | } |
566f8441 |
419 | case GeomAbs_Plane: |
420 | case GeomAbs_Cylinder: |
421 | case GeomAbs_Cone: |
422 | case GeomAbs_Sphere: |
423 | case GeomAbs_Torus: |
424 | case GeomAbs_BezierSurface: |
425 | case GeomAbs_OtherSurface: |
426 | case GeomAbs_SurfaceOfExtrusion: break; |
7fd59977 |
427 | } |
428 | return 1; |
429 | } |
430 | |
431 | //======================================================================= |
432 | //function : UIntervals |
433 | //purpose : |
434 | //======================================================================= |
435 | |
436 | void GeomAdaptor_Surface::UIntervals(TColStd_Array1OfReal& T, const GeomAbs_Shape S) const |
437 | { |
438 | Standard_Integer myNbUIntervals = 1; |
439 | |
440 | switch (mySurfaceType) |
441 | { |
442 | case GeomAbs_BSplineSurface: |
443 | { |
6b84c3f7 |
444 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
7fd59977 |
445 | GeomAdaptor_Curve myBasisCurve |
446 | (myBspl->VIso(myBspl->VKnot(myBspl->FirstVKnotIndex())),myUFirst,myULast); |
447 | myNbUIntervals = myBasisCurve.NbIntervals(S); |
448 | myBasisCurve.Intervals(T,S); |
449 | break; |
450 | } |
6b84c3f7 |
451 | case GeomAbs_SurfaceOfExtrusion: |
7fd59977 |
452 | { |
453 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
454 | (Handle(Geom_SurfaceOfLinearExtrusion)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
455 | if (myBasisCurve.GetType() == GeomAbs_BSplineCurve) |
456 | { |
457 | myNbUIntervals = myBasisCurve.NbIntervals(S); |
458 | myBasisCurve.Intervals(T,S); |
459 | } |
460 | break; |
461 | } |
6b84c3f7 |
462 | case GeomAbs_OffsetSurface: |
7fd59977 |
463 | { |
464 | GeomAbs_Shape BaseS = GeomAbs_CN; |
465 | switch(S) |
466 | { |
467 | case GeomAbs_G1: |
468 | case GeomAbs_G2: Standard_DomainError::Raise("GeomAdaptor_Curve::UIntervals"); |
469 | case GeomAbs_C0: BaseS = GeomAbs_C1; break; |
470 | case GeomAbs_C1: BaseS = GeomAbs_C2; break; |
471 | case GeomAbs_C2: BaseS = GeomAbs_C3; break; |
566f8441 |
472 | case GeomAbs_C3: |
473 | case GeomAbs_CN: break; |
7fd59977 |
474 | } |
6b84c3f7 |
475 | Handle(Geom_OffsetSurface) myOffSurf = Handle(Geom_OffsetSurface)::DownCast(mySurface); |
476 | GeomAdaptor_Surface Sur(myOffSurf->BasisSurface()); |
7fd59977 |
477 | myNbUIntervals = Sur.NbUIntervals(BaseS); |
478 | Sur.UIntervals(T, BaseS); |
479 | } |
566f8441 |
480 | case GeomAbs_Plane: |
481 | case GeomAbs_Cylinder: |
482 | case GeomAbs_Cone: |
483 | case GeomAbs_Sphere: |
484 | case GeomAbs_Torus: |
485 | case GeomAbs_BezierSurface: |
486 | case GeomAbs_OtherSurface: |
487 | case GeomAbs_SurfaceOfRevolution: break; |
7fd59977 |
488 | } |
489 | |
490 | T(T.Lower()) = myUFirst; |
491 | T(T.Lower() + myNbUIntervals) = myULast; |
492 | } |
493 | |
494 | //======================================================================= |
495 | //function : VIntervals |
496 | //purpose : |
497 | //======================================================================= |
498 | |
499 | void GeomAdaptor_Surface::VIntervals(TColStd_Array1OfReal& T, const GeomAbs_Shape S) const |
500 | { |
501 | Standard_Integer myNbVIntervals = 1; |
502 | |
503 | switch (mySurfaceType) |
504 | { |
505 | case GeomAbs_BSplineSurface: |
506 | { |
6b84c3f7 |
507 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
7fd59977 |
508 | GeomAdaptor_Curve myBasisCurve |
509 | (myBspl->UIso(myBspl->UKnot(myBspl->FirstUKnotIndex())),myVFirst,myVLast); |
510 | myNbVIntervals = myBasisCurve.NbIntervals(S); |
511 | myBasisCurve.Intervals(T,S); |
512 | break; |
513 | } |
6b84c3f7 |
514 | case GeomAbs_SurfaceOfRevolution: |
7fd59977 |
515 | { |
6b84c3f7 |
516 | Handle(Geom_SurfaceOfRevolution) myRevSurf = |
517 | Handle(Geom_SurfaceOfRevolution)::DownCast(mySurface); |
518 | GeomAdaptor_Curve myBasisCurve(myRevSurf->BasisCurve(), myVFirst, myVLast); |
7fd59977 |
519 | if (myBasisCurve.GetType() == GeomAbs_BSplineCurve) |
520 | { |
521 | myNbVIntervals = myBasisCurve.NbIntervals(S); |
522 | myBasisCurve.Intervals(T,S); |
523 | } |
524 | break; |
525 | } |
6b84c3f7 |
526 | case GeomAbs_OffsetSurface: |
7fd59977 |
527 | { |
528 | GeomAbs_Shape BaseS = GeomAbs_CN; |
529 | switch(S) |
530 | { |
531 | case GeomAbs_G1: |
532 | case GeomAbs_G2: Standard_DomainError::Raise("GeomAdaptor_Curve::VIntervals"); |
533 | case GeomAbs_C0: BaseS = GeomAbs_C1; break; |
534 | case GeomAbs_C1: BaseS = GeomAbs_C2; break; |
535 | case GeomAbs_C2: BaseS = GeomAbs_C3; break; |
566f8441 |
536 | case GeomAbs_C3: |
537 | case GeomAbs_CN: break; |
7fd59977 |
538 | } |
6b84c3f7 |
539 | Handle(Geom_OffsetSurface) myOffSurf = Handle(Geom_OffsetSurface)::DownCast(mySurface); |
540 | GeomAdaptor_Surface Sur(myOffSurf->BasisSurface()); |
7fd59977 |
541 | myNbVIntervals = Sur.NbVIntervals(BaseS); |
542 | Sur.VIntervals(T, BaseS); |
543 | } |
566f8441 |
544 | case GeomAbs_Plane: |
545 | case GeomAbs_Cylinder: |
546 | case GeomAbs_Cone: |
547 | case GeomAbs_Sphere: |
548 | case GeomAbs_Torus: |
549 | case GeomAbs_BezierSurface: |
550 | case GeomAbs_OtherSurface: |
551 | case GeomAbs_SurfaceOfExtrusion: break; |
7fd59977 |
552 | } |
553 | |
554 | T(T.Lower()) = myVFirst; |
555 | T(T.Lower() + myNbVIntervals) = myVLast; |
556 | } |
557 | |
558 | //======================================================================= |
559 | //function : UTrim |
560 | //purpose : |
561 | //======================================================================= |
562 | |
563 | Handle(Adaptor3d_HSurface) GeomAdaptor_Surface::UTrim(const Standard_Real First, |
564 | const Standard_Real Last , |
565 | const Standard_Real Tol ) const |
566 | { |
567 | return Handle(GeomAdaptor_HSurface) |
568 | (new GeomAdaptor_HSurface(mySurface,First,Last,myVFirst,myVLast,Tol,myTolV)); |
569 | } |
570 | |
571 | //======================================================================= |
572 | //function : VTrim |
573 | //purpose : |
574 | //======================================================================= |
575 | |
576 | Handle(Adaptor3d_HSurface) GeomAdaptor_Surface::VTrim(const Standard_Real First, |
577 | const Standard_Real Last , |
578 | const Standard_Real Tol ) const |
579 | { |
580 | return Handle(GeomAdaptor_HSurface) |
581 | (new GeomAdaptor_HSurface(mySurface,myUFirst,myULast,First,Last,myTolU,Tol)); |
582 | } |
583 | |
584 | //======================================================================= |
585 | //function : IsUClosed |
586 | //purpose : |
587 | //======================================================================= |
588 | |
589 | Standard_Boolean GeomAdaptor_Surface::IsUClosed() const |
590 | { |
591 | if (!mySurface->IsUClosed()) |
592 | return Standard_False; |
593 | |
594 | Standard_Real U1,U2,V1,V2; |
595 | mySurface->Bounds(U1,U2,V1,V2); |
596 | if (mySurface->IsUPeriodic()) |
597 | return (Abs(Abs(U1-U2)-Abs(myUFirst-myULast))<Precision::PConfusion()); |
598 | |
599 | return ( Abs(U1-myUFirst)<Precision::PConfusion() |
600 | && Abs(U2-myULast )<Precision::PConfusion() ); |
601 | } |
602 | |
603 | //======================================================================= |
604 | //function : IsVClosed |
605 | //purpose : |
606 | //======================================================================= |
607 | |
608 | Standard_Boolean GeomAdaptor_Surface::IsVClosed() const |
609 | { |
610 | if (!mySurface->IsVClosed()) |
611 | return Standard_False; |
612 | |
613 | Standard_Real U1,U2,V1,V2; |
614 | mySurface->Bounds(U1,U2,V1,V2); |
615 | if (mySurface->IsVPeriodic()) |
616 | return (Abs(Abs(V1-V2)-Abs(myVFirst-myVLast))<Precision::PConfusion()); |
617 | |
618 | return ( Abs(V1-myVFirst)<Precision::PConfusion() |
619 | && Abs(V2-myVLast )<Precision::PConfusion() ); |
620 | } |
621 | |
622 | //======================================================================= |
623 | //function : IsUPeriodic |
624 | //purpose : |
625 | //======================================================================= |
626 | |
627 | Standard_Boolean GeomAdaptor_Surface::IsUPeriodic() const |
628 | { |
629 | return (mySurface->IsUPeriodic()); |
630 | } |
631 | |
632 | //======================================================================= |
633 | //function : UPeriod |
634 | //purpose : |
635 | //======================================================================= |
636 | |
637 | Standard_Real GeomAdaptor_Surface::UPeriod() const |
638 | { |
639 | Standard_NoSuchObject_Raise_if(!IsUPeriodic()," "); |
640 | return mySurface->UPeriod(); |
641 | } |
642 | |
643 | //======================================================================= |
644 | //function : IsVPeriodic |
645 | //purpose : |
646 | //======================================================================= |
647 | |
648 | Standard_Boolean GeomAdaptor_Surface::IsVPeriodic() const |
649 | { |
650 | return (mySurface->IsVPeriodic()); |
651 | } |
652 | |
653 | //======================================================================= |
654 | //function : VPeriod |
655 | //purpose : |
656 | //======================================================================= |
657 | |
658 | Standard_Real GeomAdaptor_Surface::VPeriod() const |
659 | { |
660 | Standard_NoSuchObject_Raise_if(!IsVPeriodic()," "); |
661 | return mySurface->VPeriod(); |
662 | } |
663 | |
94f71cad |
664 | //======================================================================= |
665 | //function : RebuildCache |
666 | //purpose : |
667 | //======================================================================= |
668 | void GeomAdaptor_Surface::RebuildCache(const Standard_Real theU, |
669 | const Standard_Real theV) const |
670 | { |
6b84c3f7 |
671 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
672 | mySurfaceCache->BuildCache(theU, theV, |
94f71cad |
673 | myBspl->UDegree(), myBspl->IsUPeriodic(), myBspl->UKnotSequence(), |
674 | myBspl->VDegree(), myBspl->IsVPeriodic(), myBspl->VKnotSequence(), |
675 | myBspl->Poles(), myBspl->Weights()); |
676 | } |
677 | |
7fd59977 |
678 | //======================================================================= |
679 | //function : Value |
680 | //purpose : |
681 | //======================================================================= |
682 | |
683 | gp_Pnt GeomAdaptor_Surface::Value(const Standard_Real U, |
684 | const Standard_Real V) const |
685 | { |
6b84c3f7 |
686 | gp_Pnt aValue; |
687 | D0(U, V, aValue); |
688 | return aValue; |
7fd59977 |
689 | } |
690 | |
691 | //======================================================================= |
692 | //function : D0 |
693 | //purpose : |
694 | //======================================================================= |
695 | |
696 | void GeomAdaptor_Surface::D0(const Standard_Real U, |
697 | const Standard_Real V, gp_Pnt& P) const |
698 | { |
6b84c3f7 |
699 | switch (mySurfaceType) |
94f71cad |
700 | { |
6b84c3f7 |
701 | case GeomAbs_BSplineSurface: |
702 | if (!mySurfaceCache.IsNull()) |
703 | { |
704 | if (!mySurfaceCache->IsCacheValid(U, V)) |
705 | RebuildCache(U, V); |
706 | mySurfaceCache->D0(U, V, P); |
707 | } |
708 | else |
709 | mySurface->D0(U, V, P); |
710 | break; |
94f71cad |
711 | |
6b84c3f7 |
712 | case GeomAbs_OffsetSurface: |
713 | case GeomAbs_SurfaceOfExtrusion: |
714 | case GeomAbs_SurfaceOfRevolution: |
715 | Standard_NoSuchObject_Raise_if(myNestedEvaluator.IsNull(), |
716 | "GeomAdaptor_Surface::D0: evaluator is not initialized"); |
717 | myNestedEvaluator->D0(U, V, P); |
718 | break; |
719 | |
720 | default: |
721 | mySurface->D0(U, V, P); |
722 | } |
7fd59977 |
723 | } |
724 | |
725 | |
726 | //======================================================================= |
727 | //function : D1 |
728 | //purpose : |
729 | //======================================================================= |
730 | |
731 | void GeomAdaptor_Surface::D1(const Standard_Real U, |
732 | const Standard_Real V, |
94f71cad |
733 | gp_Pnt& P, |
734 | gp_Vec& D1U, |
735 | gp_Vec& D1V ) const |
7fd59977 |
736 | { |
94f71cad |
737 | Standard_Integer Ideb, Ifin, IVdeb, IVfin, USide=0, VSide=0; |
7fd59977 |
738 | Standard_Real u = U, v = V; |
739 | if (Abs(U-myUFirst) <= myTolU) {USide= 1; u = myUFirst;} |
740 | else if (Abs(U-myULast) <= myTolU) {USide= -1; u = myULast;} |
741 | if (Abs(V-myVFirst) <= myTolV) {VSide= 1; v = myVFirst;} |
742 | else if (Abs(V-myVLast) <= myTolV) {VSide= -1; v = myVLast;} |
743 | |
744 | switch(mySurfaceType) { |
6b84c3f7 |
745 | case GeomAbs_BSplineSurface: { |
746 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
747 | if ((USide != 0 || VSide != 0) && |
94f71cad |
748 | IfUVBound(u, v, Ideb, Ifin, IVdeb, IVfin, USide, VSide)) |
749 | myBspl->LocalD1(u, v, Ideb, Ifin, IVdeb, IVfin, P, D1U, D1V); |
750 | else if (!mySurfaceCache.IsNull()) |
751 | { |
752 | if (!mySurfaceCache->IsCacheValid(U, V)) |
753 | RebuildCache(U, V); |
754 | mySurfaceCache->D1(U, V, P, D1U, D1V); |
7fd59977 |
755 | } |
94f71cad |
756 | else |
757 | myBspl->D1(u, v, P, D1U, D1V); |
758 | break; |
6b84c3f7 |
759 | } |
94f71cad |
760 | |
761 | case GeomAbs_SurfaceOfExtrusion: |
94f71cad |
762 | case GeomAbs_SurfaceOfRevolution: |
94f71cad |
763 | case GeomAbs_OffsetSurface: |
6b84c3f7 |
764 | Standard_NoSuchObject_Raise_if(myNestedEvaluator.IsNull(), |
765 | "GeomAdaptor_Surface::D1: evaluator is not initialized"); |
766 | myNestedEvaluator->D1(u, v, P, D1U, D1V); |
94f71cad |
767 | break; |
6b84c3f7 |
768 | |
94f71cad |
769 | default: |
770 | mySurface->D1(u, v, P, D1U, D1V); |
7fd59977 |
771 | } |
772 | } |
773 | |
774 | //======================================================================= |
775 | //function : D2 |
776 | //purpose : |
777 | //======================================================================= |
778 | |
779 | void GeomAdaptor_Surface::D2(const Standard_Real U, |
94f71cad |
780 | const Standard_Real V, |
781 | gp_Pnt& P, |
782 | gp_Vec& D1U, |
783 | gp_Vec& D1V, |
784 | gp_Vec& D2U, |
785 | gp_Vec& D2V, |
786 | gp_Vec& D2UV) const |
7fd59977 |
787 | { |
94f71cad |
788 | Standard_Integer Ideb, Ifin, IVdeb, IVfin, USide=0, VSide=0; |
7fd59977 |
789 | Standard_Real u = U, v = V; |
790 | if (Abs(U-myUFirst) <= myTolU) {USide= 1; u = myUFirst;} |
791 | else if (Abs(U-myULast) <= myTolU) {USide= -1; u = myULast;} |
792 | if (Abs(V-myVFirst) <= myTolV) {VSide= 1; v = myVFirst;} |
793 | else if (Abs(V-myVLast) <= myTolV) {VSide= -1; v = myVLast;} |
794 | |
94f71cad |
795 | switch(mySurfaceType) { |
6b84c3f7 |
796 | case GeomAbs_BSplineSurface: { |
797 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
798 | if ((USide != 0 || VSide != 0) && |
799 | IfUVBound(u, v, Ideb, Ifin, IVdeb, IVfin, USide, VSide)) |
94f71cad |
800 | myBspl->LocalD2(u, v, Ideb, Ifin, IVdeb, IVfin, P, D1U, D1V, D2U, D2V, D2UV); |
801 | else if (!mySurfaceCache.IsNull()) |
7fd59977 |
802 | { |
94f71cad |
803 | if (!mySurfaceCache->IsCacheValid(U, V)) |
804 | RebuildCache(U, V); |
805 | mySurfaceCache->D2(U, V, P, D1U, D1V, D2U, D2V, D2UV); |
806 | } |
6b84c3f7 |
807 | else myBspl->D2(u, v, P, D1U, D1V, D2U, D2V, D2UV); |
94f71cad |
808 | break; |
6b84c3f7 |
809 | } |
94f71cad |
810 | |
6b84c3f7 |
811 | case GeomAbs_SurfaceOfExtrusion: |
812 | case GeomAbs_SurfaceOfRevolution: |
813 | case GeomAbs_OffsetSurface : |
814 | Standard_NoSuchObject_Raise_if(myNestedEvaluator.IsNull(), |
815 | "GeomAdaptor_Surface::D2: evaluator is not initialized"); |
816 | myNestedEvaluator->D2(u, v, P, D1U, D1V, D2U, D2V, D2UV); |
94f71cad |
817 | break; |
818 | |
6b84c3f7 |
819 | default: { mySurface->D2(u, v, P, D1U, D1V, D2U, D2V, D2UV); |
94f71cad |
820 | break;} |
821 | } |
7fd59977 |
822 | } |
823 | |
824 | |
825 | //======================================================================= |
826 | //function : D3 |
827 | //purpose : |
828 | //======================================================================= |
829 | |
830 | void GeomAdaptor_Surface::D3(const Standard_Real U, const Standard_Real V, |
831 | gp_Pnt& P, gp_Vec& D1U, gp_Vec& D1V, |
832 | gp_Vec& D2U, gp_Vec& D2V, gp_Vec& D2UV, |
833 | gp_Vec& D3U, gp_Vec& D3V, gp_Vec& D3UUV, |
834 | gp_Vec& D3UVV) const |
835 | { |
836 | Standard_Integer Ideb,Ifin,IVdeb,IVfin,USide=0,VSide=0; |
837 | Standard_Real u = U, v = V; |
838 | if (Abs(U-myUFirst) <= myTolU) {USide= 1; u = myUFirst;} |
839 | else if (Abs(U-myULast) <= myTolU) {USide= -1; u = myULast;} |
840 | if (Abs(V-myVFirst) <= myTolV) {VSide= 1; v = myVFirst;} |
841 | else if (Abs(V-myVLast) <= myTolV) {VSide= -1; v = myVLast;} |
842 | |
843 | switch(mySurfaceType) { |
6b84c3f7 |
844 | case GeomAbs_BSplineSurface: { |
845 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
846 | if ((USide == 0) && (VSide == 0)) |
847 | myBspl->D3(u, v, P, D1U, D1V, D2U, D2V, D2UV, D3U, D3V, D3UUV, D3UVV); |
848 | else { |
849 | if (IfUVBound(u, v, Ideb, Ifin, IVdeb, IVfin, USide, VSide)) |
850 | myBspl->LocalD3(u, v, Ideb, Ifin, IVdeb, IVfin, |
851 | P, D1U, D1V, D2U, D2V, D2UV, D3U, D3V, D3UUV, D3UVV); |
852 | else |
853 | myBspl->D3(u, v, P, D1U, D1V, D2U, D2V, D2UV, D3U, D3V, D3UUV, D3UVV); |
7fd59977 |
854 | } |
6b84c3f7 |
855 | break; |
856 | } |
857 | |
858 | case GeomAbs_SurfaceOfExtrusion: |
859 | case GeomAbs_SurfaceOfRevolution: |
860 | case GeomAbs_OffsetSurface: |
861 | Standard_NoSuchObject_Raise_if(myNestedEvaluator.IsNull(), |
862 | "GeomAdaptor_Surface::D3: evaluator is not initialized"); |
863 | myNestedEvaluator->D3(u, v, P, D1U, D1V, D2U, D2V, D2UV, D3U, D3V, D3UUV, D3UVV); |
864 | break; |
865 | |
866 | default : { mySurface->D3(u,v,P,D1U,D1V,D2U,D2V,D2UV,D3U,D3V,D3UUV,D3UVV); |
7fd59977 |
867 | break;} |
868 | } |
869 | } |
870 | |
871 | //======================================================================= |
872 | //function : DN |
873 | //purpose : |
874 | //======================================================================= |
875 | |
876 | gp_Vec GeomAdaptor_Surface::DN(const Standard_Real U, |
877 | const Standard_Real V, |
878 | const Standard_Integer Nu, |
879 | const Standard_Integer Nv) const |
880 | { |
881 | Standard_Integer Ideb,Ifin,IVdeb,IVfin,USide=0,VSide=0; |
882 | Standard_Real u = U, v = V; |
883 | if (Abs(U-myUFirst) <= myTolU) {USide= 1; u = myUFirst;} |
884 | else if (Abs(U-myULast) <= myTolU) {USide= -1; u = myULast;} |
885 | if (Abs(V-myVFirst) <= myTolV) {VSide= 1; v = myVFirst;} |
886 | else if (Abs(V-myVLast) <= myTolV) {VSide= -1; v = myVLast;} |
887 | |
888 | switch(mySurfaceType) |
889 | { |
6b84c3f7 |
890 | case GeomAbs_BSplineSurface: { |
891 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
892 | if ((USide == 0) && (VSide == 0)) |
893 | return myBspl->DN(u, v, Nu, Nv); |
7fd59977 |
894 | else { |
6b84c3f7 |
895 | if (IfUVBound(u, v, Ideb, Ifin, IVdeb, IVfin, USide, VSide)) |
896 | return myBspl->LocalDN(u, v, Ideb, Ifin, IVdeb, IVfin, Nu, Nv); |
7fd59977 |
897 | else |
6b84c3f7 |
898 | return myBspl->DN(u, v, Nu, Nv); |
7fd59977 |
899 | } |
6b84c3f7 |
900 | } |
901 | |
7fd59977 |
902 | case GeomAbs_SurfaceOfExtrusion: |
7fd59977 |
903 | case GeomAbs_SurfaceOfRevolution: |
7fd59977 |
904 | case GeomAbs_OffsetSurface: |
6b84c3f7 |
905 | Standard_NoSuchObject_Raise_if(myNestedEvaluator.IsNull(), |
906 | "GeomAdaptor_Surface::DN: evaluator is not initialized"); |
907 | return myNestedEvaluator->DN(u, v, Nu, Nv); |
7fd59977 |
908 | |
909 | case GeomAbs_Plane: |
910 | case GeomAbs_Cylinder: |
911 | case GeomAbs_Cone: |
912 | case GeomAbs_Sphere: |
913 | case GeomAbs_Torus: |
914 | case GeomAbs_BezierSurface: |
915 | case GeomAbs_OtherSurface: |
916 | default: |
917 | break; |
918 | } |
919 | |
920 | return mySurface->DN(u,v, Nu, Nv); |
921 | } |
922 | |
923 | |
924 | //======================================================================= |
925 | //function : UResolution |
926 | //purpose : |
927 | //======================================================================= |
928 | |
929 | Standard_Real GeomAdaptor_Surface::UResolution(const Standard_Real R3d) const |
930 | { |
931 | Standard_Real Res = 0.; |
932 | |
933 | switch (mySurfaceType) |
934 | { |
935 | case GeomAbs_SurfaceOfExtrusion: |
936 | { |
937 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
938 | (Handle(Geom_SurfaceOfLinearExtrusion)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
939 | return myBasisCurve.Resolution(R3d); |
940 | } |
941 | case GeomAbs_Torus: |
942 | { |
a1eb3afd |
943 | Handle(Geom_ToroidalSurface) S (Handle(Geom_ToroidalSurface)::DownCast (mySurface)); |
7fd59977 |
944 | const Standard_Real R = S->MajorRadius() + S->MinorRadius(); |
945 | if(R>Precision::Confusion()) |
946 | Res = R3d/(2.*R); |
947 | break; |
948 | } |
949 | case GeomAbs_Sphere: |
950 | { |
a1eb3afd |
951 | Handle(Geom_SphericalSurface) S (Handle(Geom_SphericalSurface)::DownCast (mySurface)); |
7fd59977 |
952 | const Standard_Real R = S->Radius(); |
953 | if(R>Precision::Confusion()) |
954 | Res = R3d/(2.*R); |
955 | break; |
956 | } |
957 | case GeomAbs_Cylinder: |
958 | { |
a1eb3afd |
959 | Handle(Geom_CylindricalSurface) S (Handle(Geom_CylindricalSurface)::DownCast (mySurface)); |
7fd59977 |
960 | const Standard_Real R = S->Radius(); |
961 | if(R>Precision::Confusion()) |
962 | Res = R3d/(2.*R); |
963 | break; |
964 | } |
965 | case GeomAbs_Cone: |
966 | { |
967 | if (myVLast - myVFirst > 1.e10) { |
968 | // Pas vraiment borne => resolution inconnue |
969 | return Precision::Parametric(R3d); |
970 | } |
a1eb3afd |
971 | Handle(Geom_ConicalSurface) S (Handle(Geom_ConicalSurface)::DownCast (mySurface)); |
7fd59977 |
972 | Handle(Geom_Curve) C = S->VIso(myVLast); |
c5f3a425 |
973 | const Standard_Real Rayon1 = Handle(Geom_Circle)::DownCast (C)->Radius(); |
7fd59977 |
974 | C = S->VIso(myVFirst); |
c5f3a425 |
975 | const Standard_Real Rayon2 = Handle(Geom_Circle)::DownCast (C)->Radius(); |
7fd59977 |
976 | const Standard_Real R = (Rayon1 > Rayon2)? Rayon1 : Rayon2; |
977 | return (R>Precision::Confusion()? (R3d / R) : 0.); |
978 | } |
979 | case GeomAbs_Plane: |
980 | { |
981 | return R3d; |
982 | } |
983 | case GeomAbs_BezierSurface: |
984 | { |
985 | Standard_Real Ures,Vres; |
c5f3a425 |
986 | Handle(Geom_BezierSurface)::DownCast (mySurface)->Resolution(R3d,Ures,Vres); |
7fd59977 |
987 | return Ures; |
988 | } |
989 | case GeomAbs_BSplineSurface: |
990 | { |
991 | Standard_Real Ures,Vres; |
c5f3a425 |
992 | Handle(Geom_BSplineSurface)::DownCast (mySurface)->Resolution(R3d,Ures,Vres); |
7fd59977 |
993 | return Ures; |
994 | } |
995 | case GeomAbs_OffsetSurface: |
996 | { |
c5f3a425 |
997 | Handle(Geom_Surface) base = Handle(Geom_OffsetSurface)::DownCast (mySurface)->BasisSurface(); |
7fd59977 |
998 | GeomAdaptor_Surface gabase(base,myUFirst,myULast,myVFirst,myVLast); |
999 | return gabase.UResolution(R3d); |
1000 | } |
1001 | default: return Precision::Parametric(R3d); |
1002 | } |
1003 | |
1004 | if ( Res <= 1.) |
1005 | return 2.*ASin(Res); |
1006 | |
c6541a0c |
1007 | return 2.*M_PI; |
7fd59977 |
1008 | } |
1009 | |
1010 | //======================================================================= |
1011 | //function : VResolution |
1012 | //purpose : |
1013 | //======================================================================= |
1014 | |
1015 | Standard_Real GeomAdaptor_Surface::VResolution(const Standard_Real R3d) const |
1016 | { |
1017 | Standard_Real Res = 0.; |
1018 | |
1019 | switch (mySurfaceType) |
1020 | { |
1021 | case GeomAbs_SurfaceOfRevolution: |
1022 | { |
1023 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
1024 | (Handle(Geom_SurfaceOfRevolution)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
1025 | return myBasisCurve.Resolution(R3d); |
1026 | } |
1027 | case GeomAbs_Torus: |
1028 | { |
a1eb3afd |
1029 | Handle(Geom_ToroidalSurface) S (Handle(Geom_ToroidalSurface)::DownCast (mySurface)); |
7fd59977 |
1030 | const Standard_Real R = S->MinorRadius(); |
1031 | if(R>Precision::Confusion()) |
1032 | Res = R3d/(2.*R); |
1033 | break; |
1034 | } |
1035 | case GeomAbs_Sphere: |
1036 | { |
a1eb3afd |
1037 | Handle(Geom_SphericalSurface) S (Handle(Geom_SphericalSurface)::DownCast (mySurface)); |
7fd59977 |
1038 | const Standard_Real R = S->Radius(); |
1039 | if(R>Precision::Confusion()) |
1040 | Res = R3d/(2.*R); |
1041 | break; |
1042 | } |
1043 | case GeomAbs_SurfaceOfExtrusion: |
1044 | case GeomAbs_Cylinder: |
1045 | case GeomAbs_Cone: |
1046 | case GeomAbs_Plane: |
1047 | { |
1048 | return R3d; |
1049 | } |
1050 | case GeomAbs_BezierSurface: |
1051 | { |
1052 | Standard_Real Ures,Vres; |
c5f3a425 |
1053 | Handle(Geom_BezierSurface)::DownCast (mySurface)->Resolution(R3d,Ures,Vres); |
7fd59977 |
1054 | return Vres; |
1055 | } |
1056 | case GeomAbs_BSplineSurface: |
1057 | { |
1058 | Standard_Real Ures,Vres; |
c5f3a425 |
1059 | Handle(Geom_BSplineSurface)::DownCast (mySurface)->Resolution(R3d,Ures,Vres); |
7fd59977 |
1060 | return Vres; |
1061 | } |
1062 | case GeomAbs_OffsetSurface: |
1063 | { |
c5f3a425 |
1064 | Handle(Geom_Surface) base = Handle(Geom_OffsetSurface)::DownCast (mySurface)->BasisSurface(); |
7fd59977 |
1065 | GeomAdaptor_Surface gabase(base,myUFirst,myULast,myVFirst,myVLast); |
1066 | return gabase.VResolution(R3d); |
1067 | } |
1068 | default: return Precision::Parametric(R3d); |
1069 | } |
1070 | |
1071 | if ( Res <= 1.) |
1072 | return 2.*ASin(Res); |
1073 | |
c6541a0c |
1074 | return 2.*M_PI; |
7fd59977 |
1075 | } |
1076 | |
1077 | //======================================================================= |
1078 | //function : Plane |
1079 | //purpose : |
1080 | //======================================================================= |
1081 | |
1082 | gp_Pln GeomAdaptor_Surface::Plane() const |
1083 | { |
1084 | if (mySurfaceType != GeomAbs_Plane) |
1085 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::Plane"); |
c5f3a425 |
1086 | return Handle(Geom_Plane)::DownCast (mySurface)->Pln(); |
7fd59977 |
1087 | } |
1088 | |
1089 | //======================================================================= |
1090 | //function : Cylinder |
1091 | //purpose : |
1092 | //======================================================================= |
1093 | |
1094 | gp_Cylinder GeomAdaptor_Surface::Cylinder() const |
1095 | { |
1096 | if (mySurfaceType != GeomAbs_Cylinder) |
1097 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::Cylinder"); |
c5f3a425 |
1098 | return Handle(Geom_CylindricalSurface)::DownCast (mySurface)->Cylinder(); |
7fd59977 |
1099 | } |
1100 | |
1101 | //======================================================================= |
1102 | //function : Cone |
1103 | //purpose : |
1104 | //======================================================================= |
1105 | |
1106 | gp_Cone GeomAdaptor_Surface::Cone() const |
1107 | { |
1108 | if (mySurfaceType != GeomAbs_Cone) |
1109 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::Cone"); |
c5f3a425 |
1110 | return Handle(Geom_ConicalSurface)::DownCast (mySurface)->Cone(); |
7fd59977 |
1111 | } |
1112 | |
1113 | //======================================================================= |
1114 | //function : Sphere |
1115 | //purpose : |
1116 | //======================================================================= |
1117 | |
1118 | gp_Sphere GeomAdaptor_Surface::Sphere() const |
1119 | { |
1120 | if (mySurfaceType != GeomAbs_Sphere) |
1121 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::Sphere"); |
c5f3a425 |
1122 | return Handle(Geom_SphericalSurface)::DownCast (mySurface)->Sphere(); |
7fd59977 |
1123 | } |
1124 | |
1125 | //======================================================================= |
1126 | //function : Torus |
1127 | //purpose : |
1128 | //======================================================================= |
1129 | |
1130 | gp_Torus GeomAdaptor_Surface::Torus() const |
1131 | { |
1132 | if (mySurfaceType != GeomAbs_Torus) |
1133 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::Torus"); |
c5f3a425 |
1134 | return Handle(Geom_ToroidalSurface)::DownCast (mySurface)->Torus(); |
7fd59977 |
1135 | } |
1136 | |
1137 | //======================================================================= |
1138 | //function : UDegree |
1139 | //purpose : |
1140 | //======================================================================= |
1141 | |
1142 | Standard_Integer GeomAdaptor_Surface::UDegree() const |
1143 | { |
1144 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1145 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->UDegree(); |
7fd59977 |
1146 | if ( mySurfaceType == GeomAbs_BezierSurface) |
c5f3a425 |
1147 | return Handle(Geom_BezierSurface)::DownCast (mySurface)->UDegree(); |
7fd59977 |
1148 | if ( mySurfaceType == GeomAbs_SurfaceOfExtrusion) |
1149 | { |
1150 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
1151 | (Handle(Geom_SurfaceOfLinearExtrusion)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
1152 | return myBasisCurve.Degree(); |
1153 | } |
1154 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::UDegree"); |
1155 | return 0; |
1156 | } |
1157 | |
1158 | //======================================================================= |
1159 | //function : NbUPoles |
1160 | //purpose : |
1161 | //======================================================================= |
1162 | |
1163 | Standard_Integer GeomAdaptor_Surface::NbUPoles() const |
1164 | { |
1165 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1166 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->NbUPoles(); |
7fd59977 |
1167 | if ( mySurfaceType == GeomAbs_BezierSurface) |
c5f3a425 |
1168 | return Handle(Geom_BezierSurface)::DownCast (mySurface)->NbUPoles(); |
7fd59977 |
1169 | if ( mySurfaceType == GeomAbs_SurfaceOfExtrusion) |
1170 | { |
1171 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
1172 | (Handle(Geom_SurfaceOfLinearExtrusion)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
1173 | return myBasisCurve.NbPoles(); |
1174 | } |
1175 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::NbUPoles"); |
1176 | return 0; |
1177 | } |
1178 | |
1179 | //======================================================================= |
1180 | //function : VDegree |
1181 | //purpose : |
1182 | //======================================================================= |
1183 | |
1184 | Standard_Integer GeomAdaptor_Surface::VDegree() const |
1185 | { |
1186 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1187 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->VDegree(); |
7fd59977 |
1188 | if ( mySurfaceType == GeomAbs_BezierSurface) |
c5f3a425 |
1189 | return Handle(Geom_BezierSurface)::DownCast (mySurface)->VDegree(); |
7fd59977 |
1190 | if ( mySurfaceType == GeomAbs_SurfaceOfRevolution) |
1191 | { |
1192 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
1193 | (Handle(Geom_SurfaceOfRevolution)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
1194 | return myBasisCurve.Degree(); |
1195 | } |
1196 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::VDegree"); |
1197 | return 0; |
1198 | } |
1199 | |
1200 | //======================================================================= |
1201 | //function : NbVPoles |
1202 | //purpose : |
1203 | //======================================================================= |
1204 | |
1205 | Standard_Integer GeomAdaptor_Surface::NbVPoles() const |
1206 | { |
1207 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1208 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->NbVPoles(); |
7fd59977 |
1209 | if ( mySurfaceType == GeomAbs_BezierSurface) |
c5f3a425 |
1210 | return Handle(Geom_BezierSurface)::DownCast (mySurface)->NbVPoles(); |
7fd59977 |
1211 | if ( mySurfaceType == GeomAbs_SurfaceOfRevolution) |
1212 | { |
1213 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
1214 | (Handle(Geom_SurfaceOfRevolution)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
1215 | return myBasisCurve.NbPoles(); |
1216 | } |
1217 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::NbVPoles"); |
1218 | return 0; |
1219 | } |
1220 | |
1221 | //======================================================================= |
1222 | //function : NbUKnots |
1223 | //purpose : |
1224 | //======================================================================= |
1225 | |
1226 | Standard_Integer GeomAdaptor_Surface::NbUKnots() const |
1227 | { |
1228 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1229 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->NbUKnots(); |
7fd59977 |
1230 | if ( mySurfaceType == GeomAbs_SurfaceOfExtrusion) |
1231 | { |
1232 | GeomAdaptor_Curve myBasisCurve |
c5f3a425 |
1233 | (Handle(Geom_SurfaceOfLinearExtrusion)::DownCast (mySurface)->BasisCurve(),myUFirst,myULast); |
7fd59977 |
1234 | return myBasisCurve.NbKnots(); |
1235 | } |
1236 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::NbUKnots"); |
1237 | return 0; |
1238 | } |
1239 | |
1240 | //======================================================================= |
1241 | //function : NbVKnots |
1242 | //purpose : |
1243 | //======================================================================= |
1244 | |
1245 | Standard_Integer GeomAdaptor_Surface::NbVKnots() const |
1246 | { |
1247 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1248 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->NbVKnots(); |
7fd59977 |
1249 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::NbVKnots"); |
1250 | return 0; |
1251 | } |
1252 | //======================================================================= |
1253 | //function : IsURational |
1254 | //purpose : |
1255 | //======================================================================= |
1256 | |
1257 | Standard_Boolean GeomAdaptor_Surface::IsURational() const |
1258 | { |
1259 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1260 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->IsURational(); |
7fd59977 |
1261 | if (mySurfaceType == GeomAbs_BezierSurface) |
c5f3a425 |
1262 | return Handle(Geom_BezierSurface)::DownCast (mySurface)->IsURational(); |
7fd59977 |
1263 | return Standard_False; |
1264 | } |
1265 | |
1266 | //======================================================================= |
1267 | //function : IsVRational |
1268 | //purpose : |
1269 | //======================================================================= |
1270 | |
1271 | Standard_Boolean GeomAdaptor_Surface::IsVRational() const |
1272 | { |
1273 | if (mySurfaceType == GeomAbs_BSplineSurface) |
c5f3a425 |
1274 | return Handle(Geom_BSplineSurface)::DownCast (mySurface)->IsVRational(); |
7fd59977 |
1275 | if (mySurfaceType == GeomAbs_BezierSurface) |
c5f3a425 |
1276 | return Handle(Geom_BezierSurface)::DownCast (mySurface)->IsVRational(); |
7fd59977 |
1277 | return Standard_False; |
1278 | } |
1279 | |
1280 | //======================================================================= |
1281 | //function : Bezier |
1282 | //purpose : |
1283 | //======================================================================= |
1284 | |
1285 | Handle(Geom_BezierSurface) GeomAdaptor_Surface::Bezier() const |
1286 | { |
1287 | if (mySurfaceType != GeomAbs_BezierSurface) |
1288 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::Bezier"); |
c5f3a425 |
1289 | return Handle(Geom_BezierSurface)::DownCast (mySurface); |
7fd59977 |
1290 | } |
1291 | |
1292 | //======================================================================= |
1293 | //function : BSpline |
1294 | //purpose : |
1295 | //======================================================================= |
1296 | |
1297 | Handle(Geom_BSplineSurface) GeomAdaptor_Surface::BSpline() const |
1298 | { |
1299 | if (mySurfaceType != GeomAbs_BSplineSurface) |
1300 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::BSpline"); |
c5f3a425 |
1301 | return Handle(Geom_BSplineSurface)::DownCast (mySurface); |
7fd59977 |
1302 | } |
1303 | |
1304 | //======================================================================= |
1305 | //function : AxeOfRevolution |
1306 | //purpose : |
1307 | //======================================================================= |
1308 | |
1309 | gp_Ax1 GeomAdaptor_Surface::AxeOfRevolution() const |
1310 | { |
1311 | if (mySurfaceType != GeomAbs_SurfaceOfRevolution) |
1312 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::AxeOfRevolution"); |
c5f3a425 |
1313 | return Handle(Geom_SurfaceOfRevolution)::DownCast (mySurface)->Axis(); |
7fd59977 |
1314 | } |
1315 | |
1316 | //======================================================================= |
1317 | //function : Direction |
1318 | //purpose : |
1319 | //======================================================================= |
1320 | |
1321 | gp_Dir GeomAdaptor_Surface::Direction() const |
1322 | { |
1323 | if (mySurfaceType != GeomAbs_SurfaceOfExtrusion) |
1324 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::Direction"); |
c5f3a425 |
1325 | return Handle(Geom_SurfaceOfLinearExtrusion)::DownCast (mySurface)->Direction(); |
7fd59977 |
1326 | } |
1327 | |
1328 | //======================================================================= |
1329 | //function : BasisCurve |
1330 | //purpose : |
1331 | //======================================================================= |
1332 | |
1333 | Handle(Adaptor3d_HCurve) GeomAdaptor_Surface::BasisCurve() const |
1334 | { |
1335 | Handle(Geom_Curve) C; |
1336 | if (mySurfaceType == GeomAbs_SurfaceOfExtrusion) |
c5f3a425 |
1337 | C = Handle(Geom_SurfaceOfLinearExtrusion)::DownCast (mySurface)->BasisCurve(); |
7fd59977 |
1338 | else if (mySurfaceType == GeomAbs_SurfaceOfRevolution) |
c5f3a425 |
1339 | C = Handle(Geom_SurfaceOfRevolution)::DownCast (mySurface)->BasisCurve(); |
7fd59977 |
1340 | else |
1341 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::BasisCurve"); |
1342 | return Handle(GeomAdaptor_HCurve)(new GeomAdaptor_HCurve(C)); |
1343 | } |
1344 | |
1345 | //======================================================================= |
1346 | //function : BasisSurface |
1347 | //purpose : |
1348 | //======================================================================= |
1349 | |
1350 | Handle(Adaptor3d_HSurface) GeomAdaptor_Surface::BasisSurface() const |
1351 | { |
1352 | if (mySurfaceType != GeomAbs_OffsetSurface) |
1353 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::BasisSurface"); |
1354 | return new GeomAdaptor_HSurface |
c5f3a425 |
1355 | (Handle(Geom_OffsetSurface)::DownCast (mySurface)->BasisSurface(), |
7fd59977 |
1356 | myUFirst,myULast,myVFirst,myVLast); |
1357 | } |
1358 | |
1359 | //======================================================================= |
1360 | //function : OffsetValue |
1361 | //purpose : |
1362 | //======================================================================= |
1363 | |
1364 | Standard_Real GeomAdaptor_Surface::OffsetValue() const |
1365 | { |
1366 | if (mySurfaceType != GeomAbs_OffsetSurface) |
1367 | Standard_NoSuchObject::Raise("GeomAdaptor_Surface::BasisSurface"); |
c5f3a425 |
1368 | return Handle(Geom_OffsetSurface)::DownCast (mySurface)->Offset(); |
7fd59977 |
1369 | } |
1370 | |
1371 | //======================================================================= |
1372 | //function : IfUVBound <private> |
1373 | //purpose : locates U,V parameters if U,V =First, Last, |
1374 | // processes the finding span and returns the |
1375 | // parameters for LocalDi |
1376 | //======================================================================= |
1377 | |
1378 | Standard_Boolean GeomAdaptor_Surface::IfUVBound(const Standard_Real U, |
1379 | const Standard_Real V, |
1380 | Standard_Integer& IOutDeb, |
1381 | Standard_Integer& IOutFin, |
1382 | Standard_Integer& IOutVDeb, |
1383 | Standard_Integer& IOutVFin, |
1384 | const Standard_Integer USide, |
1385 | const Standard_Integer VSide) const |
1386 | { |
1387 | Standard_Integer Ideb,Ifin; |
6b84c3f7 |
1388 | Handle(Geom_BSplineSurface) myBspl = Handle(Geom_BSplineSurface)::DownCast(mySurface); |
1389 | Standard_Integer anUFKIndx = myBspl->FirstUKnotIndex(), |
470ebb43 |
1390 | anULKIndx = myBspl->LastUKnotIndex(), |
1391 | aVFKIndx = myBspl->FirstVKnotIndex(), aVLKIndx = myBspl->LastVKnotIndex(); |
7fd59977 |
1392 | myBspl->LocateU(U, PosTol, Ideb, Ifin, Standard_False); |
1393 | Standard_Boolean Local = (Ideb == Ifin); |
470ebb43 |
1394 | Span(USide,Ideb,Ifin,Ideb,Ifin,anUFKIndx,anULKIndx); |
7fd59977 |
1395 | Standard_Integer IVdeb,IVfin; |
1396 | myBspl->LocateV(V, PosTol, IVdeb, IVfin, Standard_False); |
1397 | if(IVdeb == IVfin) Local = Standard_True; |
470ebb43 |
1398 | Span(VSide,IVdeb,IVfin,IVdeb,IVfin,aVFKIndx,aVLKIndx); |
7fd59977 |
1399 | |
1400 | IOutDeb=Ideb; IOutFin=Ifin; |
1401 | IOutVDeb=IVdeb; IOutVFin=IVfin; |
1402 | |
1403 | return Local; |
1404 | } |
7fd59977 |
1405 | //======================================================================= |
1406 | //function : Span <private> |
1407 | //purpose : locates U,V parameters if U=UFirst or U=ULast, |
1408 | // processes the finding span and returns the |
1409 | // parameters for LocalDi |
1410 | //======================================================================= |
1411 | |
1412 | void GeomAdaptor_Surface::Span(const Standard_Integer Side, |
1413 | const Standard_Integer Ideb, |
1414 | const Standard_Integer Ifin, |
1415 | Standard_Integer& OutIdeb, |
1416 | Standard_Integer& OutIfin, |
470ebb43 |
1417 | const Standard_Integer theFKIndx, |
1418 | const Standard_Integer theLKIndx) const |
7fd59977 |
1419 | { |
1420 | if(Ideb!=Ifin)//not a knot |
1421 | { |
470ebb43 |
1422 | if(Ideb<theFKIndx) { OutIdeb=theFKIndx; OutIfin=theFKIndx+1; } |
1423 | else if(Ifin>theLKIndx) { OutIdeb=theLKIndx-1; OutIfin=theLKIndx; } |
1424 | else if(Ideb>=(theLKIndx-1)) { OutIdeb=theLKIndx-1; OutIfin=theLKIndx; } |
1425 | else if(Ifin<=theFKIndx+1) { OutIdeb=theFKIndx; OutIfin=theFKIndx+1; } |
7fd59977 |
1426 | else if(Ideb>Ifin) { OutIdeb=Ifin-1; OutIfin=Ifin; } |
1427 | else { OutIdeb=Ideb; OutIfin=Ifin; } |
1428 | } |
1429 | else |
1430 | { |
470ebb43 |
1431 | if(Ideb<=theFKIndx){ OutIdeb=theFKIndx; OutIfin=theFKIndx+1;}//first knot |
1432 | else if(Ifin>=theLKIndx) { OutIdeb=theLKIndx-1;OutIfin=theLKIndx;}//last knot |
7fd59977 |
1433 | else |
1434 | { |
1435 | if(Side==-1){OutIdeb=Ideb-1; OutIfin=Ifin;} |
1436 | else {OutIdeb=Ideb; OutIfin=Ifin+1;} |
1437 | } |
1438 | } |
1439 | } |