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