c07838e613ef5ccab83f01eb8734fb116a3b06dc
[occt.git] / src / FilletSurf / FilletSurf_InternalBuilder.cxx
1 // Copyright (c) 1997-1999 Matra Datavision
2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
3 //
4 // This file is part of Open CASCADE Technology software library.
5 //
6 // This library is free software; you can redistribute it and/or modify it under
7 // the terms of the GNU Lesser General Public License version 2.1 as published
8 // by the Free Software Foundation, with special exception defined in the file
9 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
10 // distribution for complete text of the license and disclaimer of any warranty.
11 //
12 // Alternatively, this file may be used under the terms of Open CASCADE
13 // commercial license or contractual agreement.
14
15
16 #include <Adaptor3d_TopolTool.hxx>
17 #include <BRep_Tool.hxx>
18 #include <BRepAdaptor_HCurve.hxx>
19 #include <BRepAdaptor_HCurve2d.hxx>
20 #include <BRepAdaptor_HSurface.hxx>
21 #include <BRepBlend_ConstRad.hxx>
22 #include <BRepBlend_ConstRadInv.hxx>
23 #include <BRepBlend_Line.hxx>
24 #include <BRepExtrema_DistShapeShape.hxx>
25 #include <ChFi3d_Builder_0.hxx>
26 #include <ChFi3d_FilBuilder.hxx>
27 #include <ChFiDS_CircSection.hxx>
28 #include <ChFiDS_CommonPoint.hxx>
29 #include <ChFiDS_FaceInterference.hxx>
30 #include <ChFiDS_FilSpine.hxx>
31 #include <ChFiDS_HData.hxx>
32 #include <ChFiDS_HElSpine.hxx>
33 #include <ChFiDS_ListIteratorOfListOfStripe.hxx>
34 #include <ChFiDS_ListOfStripe.hxx>
35 #include <ChFiDS_SecHArray1.hxx>
36 #include <ChFiDS_SequenceOfSurfData.hxx>
37 #include <ChFiDS_Spine.hxx>
38 #include <ChFiDS_Stripe.hxx>
39 #include <ChFiDS_SurfData.hxx>
40 #include <ElSLib.hxx>
41 #include <FilletSurf_ErrorTypeStatus.hxx>
42 #include <FilletSurf_InternalBuilder.hxx>
43 #include <FilletSurf_StatusDone.hxx>
44 #include <FilletSurf_StatusType.hxx>
45 #include <Geom2d_Curve.hxx>
46 #include <Geom_Circle.hxx>
47 #include <Geom_Curve.hxx>
48 #include <Geom_Plane.hxx>
49 #include <Geom_Surface.hxx>
50 #include <Geom_TrimmedCurve.hxx>
51 #include <GeomAbs_Shape.hxx>
52 #include <GeomAdaptor_HCurve.hxx>
53 #include <GeomAdaptor_HSurface.hxx>
54 #include <GeomAdaptor_Surface.hxx>
55 #include <gp_Circ.hxx>
56 #include <gp_Pln.hxx>
57 #include <IntCurveSurface_HInter.hxx>
58 #include <IntCurveSurface_IntersectionPoint.hxx>
59 #include <Precision.hxx>
60 #include <Standard_ConstructionError.hxx>
61 #include <TopoDS.hxx>
62 #include <TopoDS_Face.hxx>
63 #include <TopoDS_Shape.hxx>
64 #include <TopOpeBRepDS_HDataStructure.hxx>
65 #include <TopTools_ListIteratorOfListOfShape.hxx>
66
67 static Standard_Boolean isinlist(const TopoDS_Shape&         E,
68                                  const TopTools_ListOfShape& L){
69   TopTools_ListIteratorOfListOfShape It;
70   for (It.Initialize(L); It.More(); It.Next()){
71     if(E.IsSame(It.Value())) return 1;
72   }
73   return 0;
74 }
75
76 static Standard_Boolean IntPlanEdge(Handle(BRepAdaptor_HCurve)& Ed,
77                                     const gp_Pln&               P,
78                                     Standard_Real&              w,
79                                     const Standard_Real         tol3d)
80 {
81   Standard_Boolean done = 0;
82   Standard_Real f = Ed->FirstParameter();
83   Standard_Real l = Ed->LastParameter();
84   gp_Pnt Or = P.Location();
85   Handle(Geom_Plane) Pln = new Geom_Plane(P);
86   Handle(GeomAdaptor_HSurface) 
87     Plan = new GeomAdaptor_HSurface(GeomAdaptor_Surface(Pln));
88
89   IntCurveSurface_HInter Intersection;
90   Standard_Integer nbp ,iip;
91   IntCurveSurface_IntersectionPoint IP;
92   Standard_Real dist = RealLast();
93   
94   Intersection.Perform(Ed,Plan);
95   
96   if(Intersection.IsDone()){
97     nbp = Intersection.NbPoints();
98     for (iip = 1; iip <= nbp; iip++){
99       IP = Intersection.Point(iip);
100       gp_Pnt pint = IP.Pnt();
101       Standard_Real d = pint.Distance(Or);
102       if(d<dist){
103         done = 1;
104         w = IP.W();
105         dist = d;
106       }
107     }
108   }
109   gp_Pnt pdeb = Ed->Value(f);
110   gp_Pnt pfin = Ed->Value(l);
111   Standard_Real u,v;
112   //check if the extremities are not solution
113   ElSLib::Parameters(P,pdeb,u,v);
114   gp_Pnt projdeb = ElSLib::Value(u,v,P);
115   Standard_Real dprojdeb = pdeb.Distance(projdeb);
116   if(dprojdeb<tol3d){
117     Standard_Real d = pdeb.Distance(Or);
118     if(d<dist){
119       done = 1;
120       w = f;
121       dist = d;
122     }
123   }
124   ElSLib::Parameters(P,pfin,u,v);
125   gp_Pnt projfin = ElSLib::Value(u,v,P);
126   Standard_Real dprojfin = pfin.Distance(projfin);
127   if(dprojfin<tol3d){
128     Standard_Real d = pfin.Distance(Or);
129     if(d<dist){
130       done = 1;
131       w = l;
132       dist = d;
133     }
134   }
135   return done;
136 }
137
138 static Standard_Boolean ComputeEdgeParameter(const Handle(ChFiDS_Spine)&    Spine, 
139                                              const Standard_Integer         ind,
140                                              const Standard_Real            pelsp,
141                                              Standard_Real&                 ped,
142                                              const Standard_Real            tol3d)
143 {
144   Handle(ChFiDS_HElSpine) Guide = Spine->ElSpine(ind);
145   gp_Pnt P; gp_Vec V;
146   Guide->D1(pelsp,P,V);
147   gp_Pln pln(P,V);
148   Handle(BRepAdaptor_HCurve) ed = new BRepAdaptor_HCurve();
149   ed->ChangeCurve() = Spine->CurrentElementarySpine(ind);
150   return IntPlanEdge(ed,pln,ped,tol3d);
151 }
152                                              
153 //=======================================================================
154 //function : FilletSurf_InternalBuilder
155 //purpose  : 
156 //=======================================================================
157
158 FilletSurf_InternalBuilder::FilletSurf_InternalBuilder
159 (const TopoDS_Shape&      S, 
160  const ChFi3d_FilletShape FShape, 
161  const Standard_Real      Ta,
162  const Standard_Real      Tapp3d,
163  const Standard_Real      Tapp2d):ChFi3d_FilBuilder(S,FShape,Ta)
164 {
165   SetParams(Ta,Tapp3d,Tapp2d,
166             Tapp3d,Tapp2d,1.e-3);
167   SetContinuity(GeomAbs_C2,Ta);
168 }
169
170 //=======================================================================
171 //function : Add
172 //purpose  : creation of spine on a set of edges
173 // 
174 //  0 : no problem  
175 //  1 : empty list 
176 //  2 : non g1 edges
177 //  3 : non G1 adjacent faces
178 //  4 : edge is not on the shape
179 //  5 : edge is not alive 
180 //=======================================================================
181
182 Standard_Integer  FilletSurf_InternalBuilder::Add(const TopTools_ListOfShape& E, 
183                                                   const Standard_Real R)
184 {
185   if(E.IsEmpty()) return 1;
186   TopTools_ListIteratorOfListOfShape It;
187   for(It.Initialize(E); It.More(); It.Next()){
188     TopoDS_Edge cured = TopoDS::Edge(It.Value());
189     if(cured.IsNull()) return 4;
190     if(!myEFMap.Contains(cured)) return 4;
191     //check if the edge is a fracture edge
192     TopoDS_Face ff1,ff2;  
193     for(It.Initialize(myEFMap(cured));It.More();It.Next()){  
194       if (ff1.IsNull()) {
195         ff1 = TopoDS::Face(It.Value());
196       }
197       else {
198         ff2 = TopoDS::Face(It.Value());
199         if(!ff2.IsSame(ff1)){
200           break;
201         }
202       }
203     }  
204     if(ff1.IsNull() || ff2.IsNull()) return 5;
205     if(ff1.IsSame(ff2)) return 5;
206     if(BRep_Tool::Continuity(cured,ff1,ff2) != GeomAbs_C0) return 5;
207   }  
208   TopoDS_Edge ed = TopoDS::Edge(E.First());
209   ed.Orientation(TopAbs_FORWARD);
210   ChFi3d_FilBuilder::Add(R,ed);
211   Handle(ChFiDS_Stripe) st = myListStripe.First();
212   Handle(ChFiDS_Spine)& sp = st->ChangeSpine();
213   Standard_Boolean periodic = sp->IsPeriodic();
214   
215   //It is checked if edges of list E are in the contour,
216   //the edges that arenot in the list are removed from the contour,
217   //it is checked that the remainder is monoblock.
218  
219   for(It.Initialize(E); It.More(); It.Next()){
220     TopoDS_Edge cured = TopoDS::Edge(It.Value());
221     if(!Contains(cured)) return 2;
222   }  
223         
224   Handle(ChFiDS_FilSpine) newsp = new ChFiDS_FilSpine();
225   Standard_Boolean debut = 0;
226   Standard_Integer premierquinyestpas = 0;
227   Standard_Integer yatrou = 0;
228   for(Standard_Integer i = 1; i <= sp->NbEdges(); i++){
229     TopoDS_Edge cured = sp->Edges(i);
230     if(isinlist(cured,E)){
231       debut = 1;
232       if(premierquinyestpas) {
233         yatrou = 1;
234         break;
235       }
236       newsp->SetEdges(cured);
237     }
238     else if(debut && !premierquinyestpas) premierquinyestpas = i;
239   }
240   if(!periodic && yatrou) return 2;
241   if(periodic && yatrou){
242     Standard_Boolean vraitrou = 0, aLocalDebut = 0;
243     for(Standard_Integer i = sp->NbEdges(); i > yatrou; i--){
244       TopoDS_Edge cured = sp->Edges(i);
245       if(isinlist(cured,E)){
246         if(vraitrou) return 2;
247         newsp->PutInFirst(cured);
248       }
249       else if(aLocalDebut) vraitrou = 1;
250       aLocalDebut = 1;
251     }
252   }
253
254   if(newsp->NbEdges() != sp->NbEdges()) {
255     newsp->Load();
256     newsp->SetRadius(R);
257     sp = newsp;
258   }
259
260   //ElSpine is immediately constructed
261   Handle(ChFiDS_HElSpine) hels =  new ChFiDS_HElSpine();
262   gp_Vec TFirst,TLast;
263   gp_Pnt PFirst,PLast;
264   sp->D1(sp->FirstParameter(),PFirst,TFirst);
265   sp->D1(sp->LastParameter(),PLast,TLast);
266   hels->ChangeCurve().FirstParameter(sp->FirstParameter());
267   hels->ChangeCurve().SetFirstPointAndTgt(PFirst,TFirst);
268   hels->ChangeCurve().LastParameter(sp->LastParameter());
269   hels->ChangeCurve().SetLastPointAndTgt(PLast,TLast);
270   ChFi3d_PerformElSpine(hels,sp,myConti,tolesp);
271   sp->AppendElSpine(hels);
272   sp->SplitDone(Standard_True);
273   return 0;
274 }
275
276 //=======================================================================
277 //function : Perform
278 //purpose  : 
279 //=======================================================================
280
281 void FilletSurf_InternalBuilder::Perform()
282 {
283   //PerformSetOfSurfOnElSpine is enough.
284   
285   Handle(ChFiDS_Stripe) Stripe = myListStripe.First();
286   Handle(ChFiDS_HData)&  HData  = Stripe->ChangeSetOfSurfData();
287   HData =  new ChFiDS_HData();
288   Handle(ChFiDS_Spine)& Spine = Stripe->ChangeSpine();
289   TopAbs_Orientation RefOr1,RefOr2;
290   Standard_Integer RefChoix;
291   StripeOrientations(Spine,RefOr1,RefOr2,RefChoix);
292   Stripe->OrientationOnFace1(RefOr1);
293   Stripe->OrientationOnFace2(RefOr2);
294   Stripe->Choix(RefChoix);
295   PerformSetOfKGen(Stripe,0);
296 }
297
298 //=======================================================================
299 //function : PerformSurf
300 //purpose  : 
301 //=======================================================================
302
303 Standard_Boolean
304  FilletSurf_InternalBuilder::PerformSurf(ChFiDS_SequenceOfSurfData& SeqData, 
305                                          const Handle(ChFiDS_HElSpine)& Guide, 
306                                          const Handle(ChFiDS_Spine)& Spine, 
307                                          const Standard_Integer Choix, 
308                                          const Handle(BRepAdaptor_HSurface)& S1, 
309                                          const Handle(Adaptor3d_TopolTool)& I1, 
310                                          const Handle(BRepAdaptor_HSurface)& S2, 
311                                          const Handle(Adaptor3d_TopolTool)& I2, 
312                                          const Standard_Real MaxStep, 
313                                          const Standard_Real Fleche, 
314                                          const Standard_Real TolGuide, 
315                                          Standard_Real& First, 
316                                          Standard_Real& Last, 
317                                          const Standard_Boolean Inside, 
318                                          const Standard_Boolean Appro, 
319                                          const Standard_Boolean Forward, 
320                                          const Standard_Boolean RecOnS1, 
321                                          const Standard_Boolean RecOnS2, 
322                                          const math_Vector& Soldep, 
323                                          Standard_Boolean& Intf, 
324                                          Standard_Boolean& Intl)
325 {
326   Handle(ChFiDS_SurfData) Data = SeqData(1);
327   Handle(ChFiDS_FilSpine) fsp = Handle(ChFiDS_FilSpine)::DownCast(Spine);
328   if(fsp.IsNull()) Standard_ConstructionError::Raise
329     ("PerformSurf : this is not the spine of a fillet");
330   Handle(BRepBlend_Line) lin;
331   TopAbs_Orientation Or = S1->ChangeSurface().Face().Orientation();
332   if(!fsp->IsConstant()) Standard_ConstructionError::Raise
333     ("PerformSurf : no variable radiuses");
334   // Standard_Boolean maybesingular; //pour scinder les Surfdata singulieres 
335   
336   BRepBlend_ConstRad Func(S1,S2,Guide);
337   BRepBlend_ConstRadInv FInv(S1,S2,Guide);
338   Func.Set(fsp->Radius(),Choix);
339   FInv.Set(fsp->Radius(),Choix);
340   switch (GetFilletShape()) {
341   case ChFi3d_Rational:
342     Func.Set(BlendFunc_Rational);
343     break;
344   case ChFi3d_QuasiAngular:
345     Func.Set(BlendFunc_QuasiAngular);
346     break;
347   case ChFi3d_Polynomial:
348     Func.Set(BlendFunc_Polynomial);
349   }
350   Standard_Real PFirst = First;
351   done = SimulData(Data,Guide,lin,S1,I1,
352                    S2,I2,Func,FInv,PFirst,MaxStep,Fleche,
353                    TolGuide,First,Last,Inside,Appro,Forward,Soldep,
354                    20,RecOnS1,RecOnS2);
355   if(!done) return Standard_False;
356   if(lin->StartPointOnFirst().NbPointOnRst() !=0){
357     ChFi3d_FilCommonPoint(lin->StartPointOnFirst(),lin->TransitionOnS1(),
358                           Standard_True, Data->ChangeVertexFirstOnS1(), tolesp);
359   }
360   if(lin->EndPointOnFirst().NbPointOnRst() !=0){
361     ChFi3d_FilCommonPoint(lin->EndPointOnFirst(),lin->TransitionOnS1(),
362                           Standard_False,Data->ChangeVertexLastOnS1(), tolesp);
363   }
364   if(lin->StartPointOnSecond().NbPointOnRst() !=0){
365     ChFi3d_FilCommonPoint(lin->StartPointOnSecond(),lin->TransitionOnS2(),
366                           Standard_True, Data->ChangeVertexFirstOnS2(), tolesp);
367   }
368   if(lin->EndPointOnSecond().NbPointOnRst() !=0){
369     ChFi3d_FilCommonPoint(lin->EndPointOnSecond(),lin->TransitionOnS2(),
370                           Standard_False, Data->ChangeVertexLastOnS2(), tolesp);
371   }
372   done = CompleteData(Data,Func,lin,S1,S2,Or,0,0,0,0);
373   if(!done)  Standard_Failure::Raise("PerformSurf : Failed approximation!");
374 //  maybesingular = (Func.GetMinimalDistance()<=100*tolapp3d);
375   Standard_Boolean ok = 0;
376   if(!Forward){
377     Intf = 0;
378     const ChFiDS_CommonPoint& cpf1 = Data->VertexFirstOnS1();
379     if(cpf1.IsOnArc()){
380       TopoDS_Face F1 = S1->ChangeSurface().Face();
381       TopoDS_Face bid;
382       ok = Intf = !SearchFace(Spine,cpf1,F1,bid);
383     }
384     const ChFiDS_CommonPoint& cpf2 = Data->VertexFirstOnS2();
385     if(cpf2.IsOnArc() && !ok){
386       TopoDS_Face F2 = S2->ChangeSurface().Face();
387       TopoDS_Face bid;
388       Intf = !SearchFace(Spine,cpf2,F2,bid);
389     }
390   }
391   Intl = 0;
392   ok = 0;
393   const ChFiDS_CommonPoint& cpl1 = Data->VertexLastOnS1();
394   if(cpl1.IsOnArc()){
395     TopoDS_Face F1 = S1->ChangeSurface().Face();
396     TopoDS_Face bid;
397     ok = Intl = !SearchFace(Spine,cpl1,F1,bid);
398   }
399   const ChFiDS_CommonPoint& cpl2 = Data->VertexLastOnS2();
400   if(cpl2.IsOnArc() && !ok){
401     TopoDS_Face F2 = S2->ChangeSurface().Face();
402     TopoDS_Face bid;
403     Intl = !SearchFace(Spine,cpl2,F2,bid);
404   }
405   Data->FirstSpineParam(First);
406   Data->LastSpineParam (Last);
407
408 //  if (maybesingular) SplitSurf(SeqData, lin);
409 //  Necessite de trimer resultats : A faire
410   return Standard_True;
411 }
412
413 void FilletSurf_InternalBuilder::PerformSurf(ChFiDS_SequenceOfSurfData&, const Handle(ChFiDS_HElSpine)&, const Handle(ChFiDS_Spine)&, const int, const Handle(BRepAdaptor_HSurface)&, const Handle(Adaptor3d_TopolTool)&, const Handle(BRepAdaptor_HCurve2d)&, const Handle(BRepAdaptor_HSurface)&, const Handle(BRepAdaptor_HCurve2d)&, unsigned int&, const Handle(BRepAdaptor_HSurface)&, const Handle(Adaptor3d_TopolTool)&, const TopAbs_Orientation, const double, const double, const double, double&, double&, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const math_Vector&)
414 {
415  Standard_DomainError::Raise("BlendFunc_CSConstRad::Section : Not implemented");
416 }
417
418 void FilletSurf_InternalBuilder::PerformSurf(ChFiDS_SequenceOfSurfData&, const Handle(ChFiDS_HElSpine)&, const Handle(ChFiDS_Spine)&, const int, const Handle(BRepAdaptor_HSurface)&, const Handle(Adaptor3d_TopolTool)&, const Handle(BRepAdaptor_HCurve2d)&, const Handle(BRepAdaptor_HSurface)&, const Handle(BRepAdaptor_HCurve2d)&, unsigned int&, const TopAbs_Orientation, const Handle(BRepAdaptor_HSurface)&, const Handle(Adaptor3d_TopolTool)&, const Handle(BRepAdaptor_HCurve2d)&, const Handle(BRepAdaptor_HSurface)&, const Handle(BRepAdaptor_HCurve2d)&, unsigned int&, const TopAbs_Orientation, const double, const double, const double, double&, double&, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const math_Vector&)
419 {
420  Standard_DomainError::Raise("BlendFunc_CSConstRad::Section : Not implemented");
421 }
422
423 void FilletSurf_InternalBuilder::PerformSurf(ChFiDS_SequenceOfSurfData&, const Handle(ChFiDS_HElSpine)&, const Handle(ChFiDS_Spine)&, const int, const Handle(BRepAdaptor_HSurface)&, const Handle(Adaptor3d_TopolTool)&, const TopAbs_Orientation, const Handle(BRepAdaptor_HSurface)&, const Handle(Adaptor3d_TopolTool)&, const Handle(BRepAdaptor_HCurve2d)&, const Handle(BRepAdaptor_HSurface)&, const Handle(BRepAdaptor_HCurve2d)&, unsigned int&, const double, const double, const double, double&, double&, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const unsigned int, const math_Vector&)
424 {
425  Standard_DomainError::Raise("BlendFunc_CSConstRad::Section : Not implemented");
426 }
427
428 Standard_Boolean  FilletSurf_InternalBuilder::Done() const
429 {
430  return done;
431 }
432 //=======================================================================
433 //function : NbSurface
434 //purpose  :  gives the number of NUBS surfaces  of the Fillet
435 //=======================================================================
436
437 Standard_Integer FilletSurf_InternalBuilder::NbSurface() const 
438 {
439   return myListStripe.First()->SetOfSurfData()->Length();
440 }
441
442 //=======================================================================
443 //function : SurfaceFillet
444 //purpose  : gives the NUBS surface of index Index
445 //=======================================================================
446
447 const Handle(Geom_Surface)& FilletSurf_InternalBuilder::SurfaceFillet(const Standard_Integer Index) const 
448 {
449   Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->Surf();
450
451   return myDS->Surface(isurf).Surface();
452 }
453
454  //=======================================================================
455 //function : TolApp3d
456 //purpose  :  gives the 3d tolerance reached during approximation 
457 //=======================================================================
458
459 Standard_Real  FilletSurf_InternalBuilder::TolApp3d(const Standard_Integer Index) const 
460 {
461   Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->Surf();
462
463   return myDS->Surface(isurf).Tolerance();
464 }
465
466 //=======================================================================
467 //function : SupportFace1 
468 //purpose  : gives the first support  face relative to SurfaceFillet(Index)
469 //=======================================================================
470 const TopoDS_Face& FilletSurf_InternalBuilder::SupportFace1(const Standard_Integer Index) const
471 {
472   Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->IndexOfS1();
473
474   return TopoDS::Face(myDS->Shape(isurf));
475 }
476 //=======================================================================
477 //function : SupportFace2
478 //purpose  : gives the second support face relative to SurfaceFillet(Index)
479 //=======================================================================
480 const TopoDS_Face& FilletSurf_InternalBuilder::SupportFace2(const Standard_Integer Index) const 
481 {
482   Standard_Integer isurf = myListStripe.First()->SetOfSurfData()->Value(Index)->IndexOfS2();
483
484   return TopoDS::Face(myDS->Shape(isurf));
485 }
486 //===============================================================================
487 //function : CurveOnSup1 
488 //purpose  :  gives  the 3d curve  of SurfaceFillet(Index)  on SupportFace1(Index)
489 //===============================================================================
490
491 const Handle(Geom_Curve)& FilletSurf_InternalBuilder::CurveOnFace1(const Standard_Integer Index) const 
492 {
493   Standard_Integer icurv = myListStripe.First()->SetOfSurfData()->Value(Index)->
494     InterferenceOnS1().LineIndex();
495   return myDS->Curve(icurv).Curve();
496 }
497
498 //=======================================================================
499 //function : CurveOnSup2
500 //purpose  : gives the 3d  curve of  SurfaceFillet(Index) on SupportFace2(Index
501 //=======================================================================
502
503 const Handle(Geom_Curve)& FilletSurf_InternalBuilder::CurveOnFace2(const Standard_Integer Index) const 
504 {
505   Standard_Integer icurv=myListStripe.First()->SetOfSurfData()->Value(Index)->
506     InterferenceOnS2().LineIndex();
507   return myDS->Curve(icurv).Curve();
508 }
509 //=======================================================================
510 //function : PCurveOnFace1
511 //purpose  : gives the  PCurve associated to CurvOnSup1(Index)  on the support face
512 //=======================================================================
513 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurveOnFace1(const Standard_Integer Index) const 
514 {
515   return myListStripe.First()->SetOfSurfData()->Value(Index)->
516     InterferenceOnS1().PCurveOnFace();
517
518 }
519
520 //=======================================================================
521 //function : PCurveOnFillet1
522 //purpose  : gives the PCurve associated to CurveOnFace1(Index) on the Fillet
523 //=======================================================================
524
525 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurve1OnFillet(const Standard_Integer Index) const 
526 {return myListStripe.First()->SetOfSurfData()->Value(Index)->
527     InterferenceOnS1().PCurveOnSurf();
528  
529 }
530
531 //=======================================================================
532 //function : PCurveOnFace2
533 //purpose  : gives the  PCurve associated to CurvOnSup2(Index)  on the support face
534 //=======================================================================
535 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurveOnFace2(const Standard_Integer Index) const 
536 {
537   return myListStripe.First()->SetOfSurfData()->Value(Index)->
538     InterferenceOnS2().PCurveOnFace();
539 }
540
541 //=======================================================================
542 //function : PCurveOnFillet2
543 //purpose  : gives the PCurve associated to CurveOnFace2(Index) on the Fillet
544 //=======================================================================
545 const Handle(Geom2d_Curve)& FilletSurf_InternalBuilder::PCurve2OnFillet(const Standard_Integer Index) const 
546 {return myListStripe.First()->SetOfSurfData()->Value(Index)->
547     InterferenceOnS2().PCurveOnSurf();
548   
549 }
550
551 //=======================================================================
552 //function : FirstParameter 
553 //purpose  : gives the parameter of the fillet  on the first edge
554 //=======================================================================
555
556 Standard_Real FilletSurf_InternalBuilder::FirstParameter() const
557 {
558   const Handle(ChFiDS_Stripe)& st = myListStripe.First();
559   const Handle(ChFiDS_Spine)& sp = st->Spine();
560   const Handle(ChFiDS_SurfData)& sd = st->SetOfSurfData()->Value(1);
561   Standard_Real p = sd->FirstSpineParam();
562   Standard_Integer ind = 1;
563   if(sp->IsPeriodic()) ind = sp->Index(p);
564   Standard_Real ep;
565   if(ComputeEdgeParameter(sp,ind,p,ep,tolesp)) return ep;
566   return 0.0;
567 }
568 //=======================================================================
569 //function : LastParameter
570 //purpose  :  gives the parameter of the fillet  on the last edge
571 //=======================================================================
572 Standard_Real FilletSurf_InternalBuilder::LastParameter() const
573 {
574   const Handle(ChFiDS_Stripe)& st = myListStripe.First();
575   const Handle(ChFiDS_Spine)& sp = st->Spine();
576   const Handle(ChFiDS_SurfData)& sd = st->SetOfSurfData()->Value(NbSurface());
577   Standard_Real p = sd->LastSpineParam();
578   Standard_Integer ind = sp->NbEdges();
579   if(sp->IsPeriodic()) ind = sp->Index(p);
580   Standard_Real ep;
581   if(ComputeEdgeParameter(sp,ind,p,ep,tolesp)) return ep;
582   return 0.0;
583 }
584
585 //=======================================================================
586 //function : StatusStartSection
587 //purpose  :  returns: 
588 //       twoExtremityonEdge: each extremity of  start section of the Fillet is
589 //                        on the edge of  the corresponding support face.  
590 //       OneExtremityOnEdge:  only one  of  the extremities of  start  section  of the  Fillet 
591 //                           is on the  edge of the corresponding support face.  
592 //       NoExtremityOnEdge:  any extremity of  start  section  of the fillet is  on  
593 //                           the edge  of   the  corresponding support face.
594 //=======================================================================
595
596
597 FilletSurf_StatusType  FilletSurf_InternalBuilder::StartSectionStatus() const 
598
599  
600    Standard_Boolean isonedge1 = myListStripe.First()->SetOfSurfData()->Value(1)->
601                                 VertexFirstOnS1().IsOnArc();
602    Standard_Boolean isonedge2=  myListStripe.First()->SetOfSurfData()->Value(1)->
603                                 VertexFirstOnS2().IsOnArc();;
604
605   if (isonedge1 && isonedge2) 
606   {return FilletSurf_TwoExtremityOnEdge;}
607   else if ((!isonedge1)&& (!isonedge2)) 
608   {return FilletSurf_NoExtremityOnEdge ;}
609   else 
610     {return FilletSurf_OneExtremityOnEdge;}
611 }
612 //=======================================================================
613 //function : StatusEndSection()
614 //purpose  :  returns: 
615 //            twoExtremityonEdge: each extremity of  end section of the Fillet is
616 //                                on the edge of  the corresponding support face.  
617 //            OneExtremityOnEdge: only one  of  the extremities of  end section  of the  Fillet 
618 //                                is on the  edge of the corresponding support face.  
619 //            NoExtremityOnEdge:  any extremity of   end  section  of the fillet is  on  
620 //                                the edge  of   the  corresponding support face.
621 //=======================================================================
622 FilletSurf_StatusType  FilletSurf_InternalBuilder::EndSectionStatus() const 
623
624    Standard_Boolean isonedge1 = myListStripe.First()->SetOfSurfData()->Value(NbSurface())->
625                                 VertexLastOnS1().IsOnArc();
626    Standard_Boolean isonedge2=  myListStripe.First()->SetOfSurfData()->Value(NbSurface())->
627                                 VertexLastOnS2().IsOnArc();
628
629   if (isonedge1 && isonedge2) 
630   { return FilletSurf_TwoExtremityOnEdge;}
631   else if ((!isonedge1)&& (!isonedge2)) 
632    { return FilletSurf_NoExtremityOnEdge;}
633   else 
634     { return FilletSurf_OneExtremityOnEdge;}
635 }
636
637
638 //=======================================================================
639 //function : Simulate
640 //purpose  : computes only the sections used in the computation of the fillet
641 //=======================================================================
642 void FilletSurf_InternalBuilder::Simulate()
643 {
644  //ChFi3d_FilBuilder::Simulate(1);
645  Handle(ChFiDS_Stripe) Stripe = myListStripe.First();
646   Handle(ChFiDS_HData)&  HData  = Stripe->ChangeSetOfSurfData();
647   HData =  new ChFiDS_HData();
648   Handle(ChFiDS_Spine)& Spine = Stripe->ChangeSpine();
649   TopAbs_Orientation RefOr1,RefOr2;
650   Standard_Integer RefChoix;
651   StripeOrientations(Spine,RefOr1,RefOr2,RefChoix);
652   Stripe->OrientationOnFace1(RefOr1);
653   Stripe->OrientationOnFace2(RefOr2);
654   Stripe->Choix(RefChoix);
655   PerformSetOfKGen(Stripe,1);
656 }
657
658
659 //=======================================================================
660 //function : NbSection
661 //purpose  :  gives the number of sections relative to SurfaceFillet(IndexSurf)
662 //=======================================================================
663 Standard_Integer FilletSurf_InternalBuilder::NbSection(const Standard_Integer IndexSurf) const 
664 {
665 return Sect(1,IndexSurf)->Length();
666 }
667
668 //=======================================================================
669 //function : Section
670 //purpose  : gives the   arc of circle corresponding    to section number 
671 //           IndexSec  of  SurfaceFillet(IndexSurf)  (The   basis curve  of the 
672 //           trimmed curve is a Geom_Circle) 
673 //=======================================================================
674 void FilletSurf_InternalBuilder::Section(const Standard_Integer IndexSurf, const Standard_Integer IndexSec, Handle(Geom_TrimmedCurve)& Circ) const 
675 {
676  gp_Circ c;
677  Standard_Real deb,fin; 
678  Sect(1,IndexSurf)->Value(IndexSec).Get(c,deb,fin);
679  Handle(Geom_Circle) Gc = new Geom_Circle(c);
680  Circ= new Geom_TrimmedCurve(Gc,deb,fin);
681 }
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696