projects / occt.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
0024059: Eliminate compiler warning C4701 in MSVC++ with warning level 4
[occt.git] / src / BRepFill / BRepFill_MultiLine.cxx
CommitLineData
b311480e 1// Created on: 1994-11-14
2// Created by: Bruno DUMORTIER
3// Copyright (c) 1994-1999 Matra Datavision
4// Copyright (c) 1999-2012 OPEN CASCADE SAS
5//
6// The content of this file is subject to the Open CASCADE Technology Public
7// License Version 6.5 (the "License"). You may not use the content of this file
8// except in compliance with the License. Please obtain a copy of the License
9// at http://www.opencascade.org and read it completely before using this file.
10//
11// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
12// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
13//
14// The Original Code and all software distributed under the License is
15// distributed on an "AS IS" basis, without warranty of any kind, and the
16// Initial Developer hereby disclaims all such warranties, including without
17// limitation, any warranties of merchantability, fitness for a particular
18// purpose or non-infringement. Please see the License for the specific terms
19// and conditions governing the rights and limitations under the License.
20
7fd59977 21
22
23#include <BRepFill_MultiLine.ixx>
24
25#include <BRepIntCurveSurface_Inter.hxx>
26#include <gp.hxx>
27#include <gp_Lin.hxx>
28#include <gp_Ax3.hxx>
29#include <gp_Lin2d.hxx>
30#include <gp_Circ2d.hxx>
31#include <Precision.hxx>
32#include <BRep_Tool.hxx>
33#include <Geom_Line.hxx>
34#include <Geom_Plane.hxx>
35#include <Geom_TrimmedCurve.hxx>
36#include <Geom2d_TrimmedCurve.hxx>
37#include <Geom_Surface.hxx>
38#include <Geom_RectangularTrimmedSurface.hxx>
39#include <Geom2d_Line.hxx>
40#include <GeomProjLib.hxx>
41#include <Geom2dAPI_ProjectPointOnCurve.hxx>
42#include <Geom2dInt_GInter.hxx>
43#include <IntRes2d_IntersectionPoint.hxx>
44#include <Standard_ConstructionError.hxx>
45#include <IntRes2d_IntersectionSegment.hxx>
46#include <TopExp_Explorer.hxx>
47#include <TopoDS.hxx>
48#include <ElCLib.hxx>
49
50
51#include <GeomAdaptor_Surface.hxx>
52#include <GeomAbs_SurfaceType.hxx>
53
54#ifdef DRAW
55#include <DrawTrSurf.hxx>
56#endif
57
58#ifdef DEB
59static Standard_Boolean AffichCurve = Standard_False;
60static Standard_Integer NbProj = 1;
61#endif
62//POP pour NT
63#include <stdio.h>
64
65//=======================================================================
66//function : isIsoU
67//purpose :
68//=======================================================================
69
70static Standard_Boolean isIsoU(const TopoDS_Face& Face,
71 const TopoDS_Edge& Edge )
72{
73 Handle(Geom2d_Curve) C;
74 Handle(Geom2d_Line) Li;
75 Standard_Real f,l;
76
77 C = BRep_Tool::CurveOnSurface(Edge,Face, f, l);
78 if ( C.IsNull()) {
79 Standard_ConstructionError::Raise ("BRepFill_MultiLine : Edge without PCurve");
80 }
81
82 gp_Dir2d D = C->DN(f,1);
83
84 if (Abs(D.Dot(gp::DX2d())) < Abs(D.Dot(gp::DY2d())))
85 return Standard_True;
86 else
87 return Standard_False;
88}
89
90
91
92
93//=======================================================================
94//function : BRepFill_MultiLine
95//purpose :
96//=======================================================================
97
98BRepFill_MultiLine::BRepFill_MultiLine()
99{
100}
101
102
103//=======================================================================
104//function : BRepFill_MultiLine
105//purpose :
106//=======================================================================
107
108BRepFill_MultiLine::BRepFill_MultiLine(const TopoDS_Face& Face1,
109 const TopoDS_Face& Face2,
110 const TopoDS_Edge& Edge1,
111 const TopoDS_Edge& Edge2,
112 const Standard_Boolean Inv1,
113 const Standard_Boolean Inv2,
114 const Handle(Geom2d_Curve)& Bissec) :
115myFace1(Face1 ),
116myFace2(Face2 ),
117myBis (Bissec),
118myKPart(0)
119{
120 // eval if myedges are IsoU or not
121 myIsoU1 = isIsoU(Face1, Edge1);
122 myIsoU2 = isIsoU(Face2, Edge2);
123
124 // eval myU1, myV1, myU2, myV2;
125 Handle(Geom_Plane) RefPlane;
126 Handle(Geom_Plane) BasisPlane = new Geom_Plane(0.,0.,1.,0.);
127 TopLoc_Location L;
128
129 TopExp_Explorer Exp;
1d47d8d0 130 Standard_Real Umin = 0.,Vmin = 0.,Umax = 0.,Vmax = 0.,U,V;
7fd59977 131 gp_Pnt2d P1,P2;
132 gp_Vec DZ;
133 gp_Pnt P;
134
135 // Result on Face1
136 Standard_Boolean First = Standard_True;
137 for (Exp.Init(myFace1,TopAbs_EDGE);Exp.More(); Exp.Next()) {
138 TopoDS_Edge CurEdge = TopoDS::Edge(Exp.Current());
139 BRep_Tool::UVPoints(CurEdge,myFace1,P1,P2);
140 if ( First) {
141 First = Standard_False;
142 Umin = Min(P1.X(),P2.X());
143 Umax = Max(P1.X(),P2.X());
144
145 Vmin = Min(P1.Y(),P2.Y());
146 Vmax = Max(P1.Y(),P2.Y());
147 }
148 else {
149 U = Min(P1.X(),P2.X());
150 Umin = Min(Umin,U);
151 U = Max(P1.X(),P2.X());
152 Umax = Max(Umax,U);
153
154 V = Min(P1.Y(),P2.Y());
155 Vmin = Min(Vmin,V);
156 V = Max(P1.Y(),P2.Y());
157 Vmax = Max(Vmax,V);
158 }
159 }
160
161
0d969553 162 // return isos in their domain of restriction.
7fd59977 163 Handle(Geom_Curve) UU1, UU2, VV1, VV2;
164 Handle(Geom_Surface) S;
165 S = BRep_Tool::Surface(myFace1,L);
166 if (!L.IsIdentity())
167 S = Handle(Geom_Surface)::DownCast(S->Transformed(L.Transformation()));
168
169 if ( myIsoU1) {
170 if (!BRep_Tool::Degenerated(Edge1) && !Inv1) {
171 UU1 = S->UIso(Umin);
172 GeomAdaptor_Curve Dummy(UU1);
173 if (Dummy.GetType() == GeomAbs_Circle &&
174 Dummy.Circle().Radius() < Precision::Confusion()) {
175 UU1 = S->UIso(Umax);
176 }
177 }
178 else {
179 UU1 = S->UIso(Umax);
180 GeomAdaptor_Curve Dummy(UU1);
181 if (Dummy.GetType() == GeomAbs_Circle &&
182 Dummy.Circle().Radius() < Precision::Confusion()) {
183 UU1 = S->UIso(Umin);
184 }
185 }
186 VV1 = S->VIso(Vmin);
187 }
188 else {
189 if (!BRep_Tool::Degenerated(Edge1) && !Inv1) {
190 UU1 = S->VIso(Vmin);
191 GeomAdaptor_Curve Dummy(UU1);
192 if (Dummy.GetType() == GeomAbs_Circle &&
193 Dummy.Circle().Radius() < Precision::Confusion()) {
194 UU1 = S->VIso(Vmax);
195 }
196 }
197 else {
198 UU1 = S->VIso(Vmax);
199 GeomAdaptor_Curve Dummy(UU1);
200 if (Dummy.GetType() == GeomAbs_Circle &&
201 Dummy.Circle().Radius() < Precision::Confusion()) {
202 UU1 = S->VIso(Vmin);
203 }
204 }
205 VV1 = S->UIso(Umin);
206 }
207
208 if ( myIsoU1) {
209 Standard_Real dummyUmin = Umin, dummyUmax = Umax;
210 Umin = Vmin;
211 Umax = Vmax;
212 Vmin = dummyUmin;
213 Vmax = dummyUmax;
214 }
215
0d969553 216 // try duplication
7fd59977 217 GeomAdaptor_Surface GAS1(S);
218 GeomAbs_SurfaceType Type1 = GAS1.GetType();
219
220 if ( UU1->IsPeriodic()) {
221 ElCLib::AdjustPeriodic(UU1->FirstParameter(),
222 UU1->LastParameter(),
223 Precision::PConfusion(),
224 Umin, Umax);
225 }
226 if ( VV1->IsPeriodic()) {
227 ElCLib::AdjustPeriodic(VV1->FirstParameter(),
228 VV1->LastParameter(),
229 Precision::PConfusion(),
230 Vmin, Vmax);
231 }
232 if (GAS1.GetType() == GeomAbs_Sphere) {
233 if (myIsoU1)
c6541a0c 234 ElCLib::AdjustPeriodic(-M_PI/2.,M_PI/2.,
7fd59977 235 Precision::PConfusion(),
236 Umin, Umax);
237 else
c6541a0c 238 ElCLib::AdjustPeriodic(-M_PI/2.,M_PI/2.,
7fd59977 239 Precision::PConfusion(),
240 Vmin, Vmax);
241 }
0d969553 242 // end try duplication
7fd59977 243
244 myU1 = Geom2dAdaptor_Curve(GeomProjLib::Curve2d(UU1, BasisPlane),
245 Umin, Umax);
246
247 UU1->D1(Umin, P, DZ);
248 RefPlane = new Geom_Plane(gp_Ax3(P,DZ,gp::DZ()));
249
250 myV1 = Geom2dAdaptor_Curve(GeomProjLib::Curve2d(VV1, RefPlane),
251 Vmin, Vmax);
252
253
254 First = Standard_True;
255 for (Exp.Init(myFace2,TopAbs_EDGE);Exp.More(); Exp.Next()) {
256 TopoDS_Edge CurEdge = TopoDS::Edge(Exp.Current());
257 BRep_Tool::UVPoints(CurEdge,myFace2,P1,P2);
258 if ( First) {
259 First = Standard_False;
260 Umin = Min(P1.X(),P2.X());
261 Umax = Max(P1.X(),P2.X());
262
263 Vmin = Min(P1.Y(),P2.Y());
264 Vmax = Max(P1.Y(),P2.Y());
265 }
266 else {
267 U = Min(P1.X(),P2.X());
268 Umin = Min(Umin,U);
269 U = Max(P1.X(),P2.X());
270 Umax = Max(Umax,U);
271
272 V = Min(P1.Y(),P2.Y());
273 Vmin = Min(Vmin,V);
274 V = Max(P1.Y(),P2.Y());
275 Vmax = Max(Vmax,V);
276 }
277 }
278
0d969553 279 // return isos in their domain of restriction.
7fd59977 280 S = BRep_Tool::Surface(myFace2,L);
281
282 if (!L.IsIdentity())
283 S = Handle(Geom_Surface)::DownCast(S->Transformed(L.Transformation()));
284
285 if ( myIsoU2) {
286 if (!BRep_Tool::Degenerated(Edge2) && !Inv2) {
287 UU2 = S->UIso(Umin);
288 GeomAdaptor_Curve Dummy(UU2);
289 if (Dummy.GetType() == GeomAbs_Circle &&
290 Dummy.Circle().Radius() < Precision::Confusion()) {
291 UU2 = S->UIso(Umax);
292 }
293 }
294 else {
295 UU2 = S->UIso(Umax);
296 GeomAdaptor_Curve Dummy(UU2);
297 if (Dummy.GetType() == GeomAbs_Circle &&
298 Dummy.Circle().Radius() < Precision::Confusion()) {
299 UU2 = S->UIso(Umin);
300 }
301 }
302 VV2 = S->VIso(Vmin);
303 }
304 else {
305 if (!BRep_Tool::Degenerated(Edge2) && !Inv2) {
306 UU2 = S->VIso(Vmin);
307 GeomAdaptor_Curve Dummy(UU2);
308 if (Dummy.GetType() == GeomAbs_Circle &&
309 Dummy.Circle().Radius() < Precision::Confusion()) {
310 UU2 = S->VIso(Vmax);
311 }
312 }
313 else {
314 UU2 = S->VIso(Vmax);
315 GeomAdaptor_Curve Dummy(UU2);
316 if (Dummy.GetType() == GeomAbs_Circle &&
317 Dummy.Circle().Radius() < Precision::Confusion()) {
318 UU2 = S->VIso(Vmin);
319 }
320 }
321 VV2 = S->UIso(Umin);
322 }
323
324 if ( myIsoU2) {
325 Standard_Real dummyUmin = Umin, dummyUmax = Umax;
326 Umin = Vmin;
327 Umax = Vmax;
328 Vmin = dummyUmin;
329 Vmax = dummyUmax;
330 }
331
0d969553 332 // try duplication
7fd59977 333 GeomAdaptor_Surface GAS2(S);
334 GeomAbs_SurfaceType Type2 = GAS2.GetType();
335
336 if ( UU2->IsPeriodic()) {
337 ElCLib::AdjustPeriodic(UU2->FirstParameter(),
338 UU2->LastParameter(),
339 Precision::PConfusion(),
340 Umin, Umax);
341 }
342 if ( VV2->IsPeriodic()) {
343 ElCLib::AdjustPeriodic(VV2->FirstParameter(),
344 VV2->LastParameter(),
345 Precision::PConfusion(),
346 Vmin, Vmax);
347 }
348 if (GAS2.GetType() == GeomAbs_Sphere) {
349 if (myIsoU2)
c6541a0c 350 ElCLib::AdjustPeriodic(-M_PI/2.,M_PI/2.,
7fd59977 351 Precision::PConfusion(),
352 Umin, Umax);
353 else
c6541a0c 354 ElCLib::AdjustPeriodic(-M_PI/2.,M_PI/2.,
7fd59977 355 Precision::PConfusion(),
356 Vmin, Vmax);
357 }
0d969553 358 // end try duplication
7fd59977 359
360 myU2 = Geom2dAdaptor_Curve(GeomProjLib::Curve2d(UU2, BasisPlane),
361 Umin, Umax);
362
363 UU2->D1(Umin, P, DZ);
364 RefPlane = new Geom_Plane(gp_Ax3(P,DZ,gp::DZ()));
365 myV2 = Geom2dAdaptor_Curve(GeomProjLib::Curve2d(VV2, RefPlane),
366 Vmin, Vmax);
367
368 // eval if in a particular case.
0d969553
Y		 369 // Particular case if :
370 // 1) - Straight Bissectrice
371 // - Bissectrice orthogonal to the base element.
372 // ==> Iso on 2 faces.
373 // 2) - Straight Bissectrice
374 // - 2 surfaces are planes.
7fd59977 375 myCont = GeomAbs_C0;
376
377 if ( myBis.GetType() == GeomAbs_Line) {
378 Standard_Real DeltaU = myBis.LastParameter() - myBis.FirstParameter();
379 gp_Pnt2d P1 = ValueOnF1(myBis.FirstParameter() + 0.1*DeltaU);
380 gp_Pnt2d P2 = ValueOnF1(myBis.FirstParameter() + 0.9*DeltaU);
381 if ( myIsoU1) {
382 if ( Abs(P1.Y() - P2.Y()) < Precision::Confusion())
383 myKPart = 1;
384 }
385 else {
386 if ( Abs(P1.X() - P2.X()) < Precision::Confusion())
387 myKPart = 1;
388 }
389
390 if ( myKPart == 1)
391 myCont = GeomAbs_G1;
392
393 if ( (Type1 == GeomAbs_Plane) && (Type2 == GeomAbs_Plane)) {
394 myKPart = 2;
395 }
396 }
397}
398
399
400//=======================================================================
401//function : IsParticularCase
402//purpose :
403//=======================================================================
404
405Standard_Boolean BRepFill_MultiLine::IsParticularCase() const
406{
407 return ( myKPart != 0);
408}
409
410
411//=======================================================================
412//function : Curves
413//purpose :
414//=======================================================================
415
416void BRepFill_MultiLine::Curves(Handle(Geom_Curve)& Curve,
417 Handle(Geom2d_Curve)& PCurve1,
418 Handle(Geom2d_Curve)& PCurve2) const
419{
420 if ( myKPart == 1) {
421 gp_Pnt2d P1,P2,PMil;
422 Standard_Real f,l;
423
424 P1 = ValueOnF1(myBis.FirstParameter());
425 P2 = ValueOnF1(myBis.LastParameter());
426
0d969553
Y		 427 // find value of the with medium point
428 // the ends can be degenerated points.
7fd59977 429
430 PMil = ValueOnF1(0.5*(myBis.FirstParameter() + myBis.LastParameter()));
431
432 TopLoc_Location L;
433 Handle(Geom_Surface) S = BRep_Tool::Surface(myFace1,L);
434 if (!L.IsIdentity())
435 S = Handle(Geom_Surface)::DownCast(S->Transformed(L.Transformation()));
436
437 Standard_Boolean Sens;
438 if ( !myIsoU1) {
439 Curve = S->UIso(PMil.X());
440 Sens = P1.Y() < P2.Y();
441 if ( Sens)
442 Curve = new Geom_TrimmedCurve(Curve, P1.Y(), P2.Y(), Sens);
443 else
444 Curve = new Geom_TrimmedCurve(Curve, P2.Y(), P1.Y(), Sens);
445
446 f = Curve->FirstParameter();
447 l = Curve->LastParameter();
448 if ( Sens)
449 PCurve1 = new Geom2d_Line(gp_Pnt2d(PMil.X(),P1.Y() - f), gp::DY2d());
450 else
451 PCurve1 = new Geom2d_Line(gp_Pnt2d(PMil.X(),P1.Y() + f),-gp::DY2d());
452 PCurve1 = new Geom2d_TrimmedCurve( PCurve1 ,f ,l);
453 }
454 else {
455 Curve = S->VIso(PMil.Y());
456 Sens = P1.X() < P2.X();
457 if (Sens)
458 Curve = new Geom_TrimmedCurve(Curve, P1.X(), P2.X(), Sens);
459 else
460 Curve = new Geom_TrimmedCurve(Curve, P2.X(), P1.X(), Sens);
461
462 f = Curve->FirstParameter();
463 l = Curve->LastParameter();
464 if ( Sens)
465 PCurve1 = new Geom2d_Line(gp_Pnt2d(P1.X() - f,PMil.Y()), gp::DX2d());
466 else
467 PCurve1 = new Geom2d_Line(gp_Pnt2d(P1.X() + f,PMil.Y()), -gp::DX2d());
468 PCurve1 = new Geom2d_TrimmedCurve( PCurve1 ,f ,l);
469 }
470
471 P1 = ValueOnF2(myBis.FirstParameter());
472 P2 = ValueOnF2(myBis.LastParameter());
473 PMil = ValueOnF2(0.5*(myBis.FirstParameter() + myBis.LastParameter()));
474
475 if (!myIsoU2) {
476 Sens = P1.Y() < P2.Y();
477 if ( Sens)
478 PCurve2 = new Geom2d_Line(gp_Pnt2d(PMil.X(),(P1.Y() - f)), gp::DY2d());
479 else
480 PCurve2 = new Geom2d_Line(gp_Pnt2d(PMil.X(),(P1.Y() + f)), -gp::DY2d());
481 }
482 else {
483 Sens = P1.X() < P2.X();
484 if ( Sens)
485 PCurve2 = new Geom2d_Line(gp_Pnt2d(P1.X() - f ,PMil.Y()), gp::DX2d());
486 else
487 PCurve2 = new Geom2d_Line(gp_Pnt2d(P1.X() + f ,PMil.Y()), -gp::DX2d());
488 }
489 PCurve2 = new Geom2d_TrimmedCurve( PCurve2 ,f ,l);
490 }
491 else if ( myKPart == 2) {
492 TopLoc_Location L;
493
494 Handle(Geom_Surface) S = BRep_Tool::Surface(myFace1,L);
495 if (!L.IsIdentity())
496 S = Handle(Geom_Surface)::DownCast(S->Transformed(L.Transformation()));
497
498 if (S->IsKind(STANDARD_TYPE(Geom_RectangularTrimmedSurface)))
499 S = Handle(Geom_RectangularTrimmedSurface)::DownCast(S)->BasisSurface();
500
501 Handle(Geom_Plane) Plane = Handle(Geom_Plane)::DownCast(S);
502 // eval the 3d curve corresponding to the bissectrice.
503 gp_Pnt2d P = myBis.Line().Location();
504 gp_Dir2d D = myBis.Line().Direction();
505 Handle(Geom_Line) Line = new Geom_Line(gp_Pnt(P.X(),P.Y(),0.),
506 gp_Dir(D.X(),D.Y(),0.) );
507 Handle(Geom_TrimmedCurve) TLine =
508 new Geom_TrimmedCurve(Line, myBis.FirstParameter(),
509 myBis.LastParameter());
510 Curve = GeomProjLib::ProjectOnPlane(TLine, Plane,
511 gp::DZ(), Standard_False);
512
513#ifdef DRAW
514 if ( AffichCurve) {
515//POP pour NT
516 // char name[100];
517 char* name = new char[100];
518 sprintf(name,"C2_%d",NbProj);
519 DrawTrSurf::Set(name,TLine);
520 sprintf(name,"C3_%d",NbProj);
521 DrawTrSurf::Set(name,Curve);
522 sprintf(name,"SS_%d",NbProj);
523 DrawTrSurf::Set(name,Plane);
524 NbProj++;
525 }
526#endif
527
528 // eval PCurve1
529 PCurve1 = GeomProjLib::Curve2d(Curve,Plane);
530
531 // eval PCurve2
532 S = BRep_Tool::Surface(myFace2,L);
533 if (!L.IsIdentity())
534 S = Handle(Geom_Surface)::DownCast(S->Transformed(L.Transformation()));
535 if (S->IsKind(STANDARD_TYPE(Geom_RectangularTrimmedSurface)))
536 S = Handle(Geom_RectangularTrimmedSurface)::DownCast(S)->BasisSurface();
537 Plane = Handle(Geom_Plane)::DownCast(S);
538 PCurve2 = GeomProjLib::Curve2d(Curve,Plane);
539 }
540}
541
542
543//=======================================================================
544//function : FirstParameter
545//purpose :
546//=======================================================================
547
548Standard_Real BRepFill_MultiLine::FirstParameter() const
549{
550 return myBis.FirstParameter();
551}
552
553
554//=======================================================================
555//function : LastParameter
556//purpose :
557//=======================================================================
558
559Standard_Real BRepFill_MultiLine::LastParameter() const
560{
561 return myBis.LastParameter();
562}
563
564
565//=======================================================================
566//function : Value
567//purpose :
568//=======================================================================
569
570gp_Pnt BRepFill_MultiLine::Value(const Standard_Real U) const
571{
572 Handle(Geom_Surface) S;
573 TopLoc_Location L;
574
575 S = BRep_Tool::Surface(myFace1,L);
576
577 gp_Pnt2d P2d = ValueOnF1(U);
578
579 gp_Pnt P3d = S->Value(P2d.X(), P2d.Y());
580 P3d.Transform(L.Transformation());
581
582 return P3d;
583}
584
585
586//=======================================================================
587//function : ValueOnFace
588//purpose :
589//=======================================================================
590
591static gp_Pnt2d ValueOnFace(const Standard_Real U,
592 const Geom2dAdaptor_Curve& TheBis,
593 const Geom2dAdaptor_Curve& TheU,
594 const Geom2dAdaptor_Curve& TheV,
595 const Standard_Boolean IsIsoU)
596{
597 gp_Pnt2d P = TheBis.Value(U);
598
599 Geom2dAPI_ProjectPointOnCurve Ext(P,TheU.Curve(),
600 TheU.FirstParameter(),
601 TheU.LastParameter());
602#ifdef DRAW
603 if (AffichCurve) {
604 char* TheUname = "TheU";
605 char* PP1name = "PP1";
606 DrawTrSurf::Set(TheUname,TheU.Curve());
607 DrawTrSurf::Set(PP1name,P);
608// DrawTrSurf::Set("TheU",TheU.Curve());
609// DrawTrSurf::Set("PP1",P);
610 }
611#endif
612
613 Standard_Real UU =0., Dist = Precision::Infinite(), D1, D2;
614
615 if ( !Ext.NbPoints() == 0 ) {
616 UU = Ext.LowerDistanceParameter();
617 Dist = Ext.LowerDistance();
618 }
0d969553 619 // Control with `ends`
7fd59977 620 D1 = P.Distance(TheU.Value(TheU.FirstParameter()));
621 D2 = P.Distance(TheU.Value(TheU.LastParameter()));
622
623 if (D1 < Dist || D2 < Dist) {
624 if ( Abs( D1 - D2) < Precision::Confusion()) {
625 if ( TheU.GetType() == GeomAbs_Circle) {
626 gp_Vec2d Axis = TheU.Circle().XAxis().Direction();
627 gp_Vec2d D12d = TheBis.DN(TheBis.FirstParameter(),1);
628 Standard_Real Ang = Axis.Angle(D12d);
629 if ( !TheU.Circle().IsDirect()) Ang = -Ang;
630 UU = ElCLib::InPeriod( Ang, TheU.FirstParameter(),
c6541a0c 631 TheU.FirstParameter() + 2*M_PI);
7fd59977 632 Dist = TheU.Circle().Radius();
633 }
634 else {
635#ifdef DEB
636 cout << "MultiLine : D1 = D2 and the Curve is not a circle" << endl;
637 cout << " ---> ValueOnFace failed at parameter U = " << U << endl;
638#endif
639 Standard_ConstructionError::Raise("BRepFill_MultiLine: ValueOnFace");
640 }
641 }
642 else if ( D1 < D2) {
643 Dist = D1;
644 UU = TheU.FirstParameter();
645 }
646 else {
647 Dist = D2;
648 UU = TheU.LastParameter();
649 }
650 }
651
652 Standard_Real Tol = Precision::Confusion();
653 Standard_Real VV;
654
655 gp_Pnt2d PF = TheV.Value(TheV.FirstParameter());
656 gp_Pnt2d PL = TheV.Value(TheV.LastParameter());
657
658 if (Abs(Dist - Abs(PF.Y())) < Tol) {
659 VV = TheV.FirstParameter();
660 }
661 else if (Abs(Dist - Abs(PL.Y())) < Tol) {
662 VV = TheV.LastParameter();
663 }
664 else {
0d969553 665 // test if the curve is at the side `negative Y`.
7fd59977 666 if ( Min( PF.Y(),PL.Y()) < -Tol) Dist = -Dist;
667
668 Handle(Geom2d_Line) Line
669 = new Geom2d_Line(gp_Pnt2d(0., Dist), gp::DX2d());
670
671#ifdef DRAW
672 if (AffichCurve) {
673 static Standard_CString aTheV = "TheV" ;
674 DrawTrSurf::Set(aTheV,TheV.Curve());
675 static Standard_CString aLINF1 = "LINF1" ;
676 DrawTrSurf::Set(aLINF1,Line);
677 }
678#endif
679
680 Geom2dAdaptor_Curve Cu1 = TheV;
681 Geom2dAdaptor_Curve Cu2( Line);
682
683 Standard_Real TolConf = 0.;
684
685 Geom2dInt_GInter Intersector(Cu1,Cu2,TolConf,Tol);
686
687 if ( !Intersector.IsDone()) {
688#ifdef DEB
689 cout << "Intersector not done" << endl;
690 cout << " ---> ValueonFace failed at parameter U = " << U << endl;
691#endif
692 return gp_Pnt2d(0.,0.);
693 }
694 else {
695 if ( Intersector.NbPoints() > 0) {
696 VV = Intersector.Point(1).ParamOnFirst();
697 }
698 else if ( Intersector.NbSegments() > 0) {
699 IntRes2d_IntersectionSegment Seg = Intersector.Segment(1);
700 Standard_Real VS1 = Seg.FirstPoint().ParamOnFirst();
701 Standard_Real VS2 = Seg.LastPoint().ParamOnFirst();
702 gp_Pnt2d PS1 = TheV.Value(VS1);
703 gp_Pnt2d PS2 = TheV.Value(VS2);
704 Standard_Real Alp = (Dist - PS1.Y())/(PS2.Y() - PS1.Y());
705 VV = Alp*(VS2 - VS1) + VS1;
706 }
707 else {
708#ifdef DEB
709 cout << "Intersector done, but no points found" << endl;
710 cout << " ---> ValueonFace failed at parameter U = " << U << endl;
711#endif
712 if (Abs(Dist - PL.Y()) < Abs(Dist - PF.Y()))
713 VV = TheV.LastParameter();
714 else
715 VV = TheV.FirstParameter();
716 }
717 }
718 }
719
720 if ( IsIsoU)
721 return gp_Pnt2d(VV,UU);
722 else
723 return gp_Pnt2d(UU,VV);
724}
725
726//=======================================================================
727//function : ValueOnF1
728//purpose :
729//=======================================================================
730
731gp_Pnt2d BRepFill_MultiLine::ValueOnF1(const Standard_Real U) const
732{
733 return ValueOnFace(U,myBis,myU1,myV1,myIsoU1);
734}
735
736
737//=======================================================================
738//function : ValueOnF2
739//purpose :
740//=======================================================================
741
742gp_Pnt2d BRepFill_MultiLine::ValueOnF2(const Standard_Real U) const
743{
744 return ValueOnFace(U,myBis,myU2,myV2,myIsoU2);
745}
746
747//=======================================================================
748//function : Value3dOnF1OnF2
749//purpose :
750//=======================================================================
751
752void BRepFill_MultiLine::Value3dOnF1OnF2(const Standard_Real U,
753 gp_Pnt& P3d,
754 gp_Pnt2d& PF1,
755 gp_Pnt2d& PF2)
756const
757{
758 PF1 = ValueOnFace(U,myBis,myU1,myV1,myIsoU1);
759 PF2 = ValueOnFace(U,myBis,myU2,myV2,myIsoU2);
760
761 Handle(Geom_Surface) S;
762 TopLoc_Location L;
763
764 S = BRep_Tool::Surface(myFace1,L);
765 P3d = S->Value(PF1.X(), PF1.Y());
766 P3d.Transform(L.Transformation());
767}
768
769//=======================================================================
770//function : Continuity
771//purpose :
772//=======================================================================
773
774GeomAbs_Shape BRepFill_MultiLine::Continuity() const
775{
776 return myCont;
777}
Open CASCADE Technology repository