]> SALOME platform Git repositories - modules/geom.git/blob - src/GEOMAlgo/GEOMAlgo_FinderShapeOn2.cxx
Salome HOME
GDD: Check that the PipeTShape is a block if it is required by user
[modules/geom.git] / src / GEOMAlgo / GEOMAlgo_FinderShapeOn2.cxx
1 //  Copyright (C) 2007-2010  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 //  Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 //  This library is free software; you can redistribute it and/or
7 //  modify it under the terms of the GNU Lesser General Public
8 //  License as published by the Free Software Foundation; either
9 //  version 2.1 of the License.
10 //
11 //  This library is distributed in the hope that it will be useful,
12 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
13 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 //  Lesser General Public License for more details.
15 //
16 //  You should have received a copy of the GNU Lesser General Public
17 //  License along with this library; if not, write to the Free Software
18 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 //  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 // File:        GEOMAlgo_FinderShapeOn1.cxx
24 // Created:     Fri Mar  4 10:31:06 2005
25 // Author:      Peter KURNEV
26 //              <pkv@irinox>
27 //
28 #include <GEOMAlgo_FinderShapeOn2.ixx>
29 #include <math.h>
30
31
32 #include <Precision.hxx>
33 #include <TColStd_Array1OfInteger.hxx>
34 #include <TColStd_MapOfInteger.hxx>
35
36 #include <gp_Trsf.hxx>
37 #include <gp_Cylinder.hxx>
38 #include <gp_Pnt.hxx>
39
40 #include <TColgp_Array1OfPnt.hxx>
41
42 #include <Poly_Array1OfTriangle.hxx>
43 #include <Poly_Triangle.hxx>
44 #include <Poly_PolygonOnTriangulation.hxx>
45 #include <Poly_Triangulation.hxx>
46 #include <Poly_Polygon3D.hxx>
47
48 #include <Geom_Curve.hxx>
49 #include <Geom_Surface.hxx>
50 #include <GeomAdaptor_Surface.hxx>
51 #include <GeomAbs_SurfaceType.hxx>
52 #include <GeomAdaptor_Curve.hxx>
53 #include <GeomAbs_CurveType.hxx>
54
55 #include <TopAbs_State.hxx>
56
57 #include <TopLoc_Location.hxx>
58 #include <TopoDS.hxx>
59 #include <TopoDS_Shape.hxx>
60 #include <TopoDS_Vertex.hxx>
61 #include <TopoDS_Face.hxx>
62 #include <TopoDS_Edge.hxx>
63
64 #include <TopExp.hxx>
65 #include <TopExp_Explorer.hxx>
66
67 #include <TopTools_IndexedMapOfShape.hxx>
68
69 #include <BRep_Tool.hxx>
70 #include <BRepLib_MakeEdge.hxx>
71
72 #include <GEOMAlgo_ListIteratorOfListOfPnt.hxx>
73
74 #include <GEOMAlgo_SurfaceTools.hxx>
75 #include <GEOMAlgo_StateCollector.hxx>
76 #include <GEOMAlgo_FinderShapeOn.hxx>
77
78 #include <GEOMAlgo_PassKey.hxx>
79 #include <GEOMAlgo_DataMapOfPassKeyInteger.hxx>
80 #include <GEOMAlgo_DataMapIteratorOfDataMapOfPassKeyInteger.hxx>
81
82 //=======================================================================
83 //function : GEOMAlgo_FinderShapeOn1
84 //purpose  : 
85 //=======================================================================
86   GEOMAlgo_FinderShapeOn2::GEOMAlgo_FinderShapeOn2()
87 :
88   GEOMAlgo_ShapeAlgo()
89 {
90   myTolerance=0.0001;
91   myShapeType=TopAbs_VERTEX;
92   myState=GEOMAlgo_ST_UNKNOWN;
93   myNbPntsMin=3;
94   myNbPntsMax=0;
95 }
96 //=======================================================================
97 //function : ~
98 //purpose  : 
99 //=======================================================================
100   GEOMAlgo_FinderShapeOn2::~GEOMAlgo_FinderShapeOn2()
101 {
102 }
103 //=======================================================================
104 //function : SetClsf
105 //purpose  : 
106 //=======================================================================
107   void GEOMAlgo_FinderShapeOn2::SetClsf(const Handle(GEOMAlgo_Clsf)& aClsf)
108 {
109   myClsf=aClsf;
110 }
111 //=======================================================================
112 //function : Clsf
113 //purpose  : 
114 //=======================================================================
115   const Handle(GEOMAlgo_Clsf)& GEOMAlgo_FinderShapeOn2::Clsf() const
116 {
117   return myClsf;
118 }
119 //=======================================================================
120 //function : SetShapeType
121 //purpose  : 
122 //=======================================================================
123   void GEOMAlgo_FinderShapeOn2::SetShapeType(const TopAbs_ShapeEnum aType)
124 {
125   myShapeType=aType;
126 }
127 //=======================================================================
128 //function : ShapeType
129 //purpose  : 
130 //=======================================================================
131   TopAbs_ShapeEnum GEOMAlgo_FinderShapeOn2::ShapeType()const
132 {
133   return myShapeType;
134 }
135 //=======================================================================
136 //function : SetState
137 //purpose  : 
138 //=======================================================================
139   void GEOMAlgo_FinderShapeOn2::SetState(const GEOMAlgo_State aState)
140 {
141   myState=aState;
142 }
143 //=======================================================================
144 //function : State
145 //purpose  : 
146 //=======================================================================
147   GEOMAlgo_State GEOMAlgo_FinderShapeOn2::State() const
148 {
149   return myState;
150 }
151 //=======================================================================
152 //function : SetNbPntsMin
153 //purpose  : 
154 //=======================================================================
155   void GEOMAlgo_FinderShapeOn2::SetNbPntsMin(const Standard_Integer aNb)
156 {
157   myNbPntsMin=aNb;
158 }
159 //=======================================================================
160 //function : NbPntsMin
161 //purpose  : 
162 //=======================================================================
163   Standard_Integer GEOMAlgo_FinderShapeOn2::NbPntsMin()const
164 {
165   return myNbPntsMin;
166 }
167 //=======================================================================
168 //function : SetNbPntsMax
169 //purpose  : 
170 //=======================================================================
171   void GEOMAlgo_FinderShapeOn2::SetNbPntsMax(const Standard_Integer aNb)
172 {
173   myNbPntsMax=aNb;
174 }
175 //=======================================================================
176 //function : NbPntsMax
177 //purpose  : 
178 //=======================================================================
179   Standard_Integer GEOMAlgo_FinderShapeOn2::NbPntsMax()const
180 {
181   return myNbPntsMax;
182 }
183 //=======================================================================
184 // function: MSS
185 // purpose: 
186 //=======================================================================
187   const GEOMAlgo_IndexedDataMapOfShapeState& GEOMAlgo_FinderShapeOn2::MSS() const
188 {
189   return myMSS;
190 }
191 //=======================================================================
192 // function: Shapes
193 // purpose: 
194 //=======================================================================
195   const TopTools_ListOfShape& GEOMAlgo_FinderShapeOn2::Shapes() const
196 {
197   Standard_Integer i, aNb;
198   TopTools_ListOfShape* pL;
199   //
200   pL=(TopTools_ListOfShape*) &myLS;
201   pL->Clear();
202   //
203   aNb=myMSS.Extent();
204   for (i=1; i<=aNb; ++i) {
205     const TopoDS_Shape& aS=myMSS.FindKey(i);
206     if (aS.ShapeType()==myShapeType) {
207       pL->Append(aS);
208     }
209   }
210   return myLS;
211 }
212 //=======================================================================
213 //function : Perform
214 //purpose  : 
215 //=======================================================================
216   void GEOMAlgo_FinderShapeOn2::Perform()
217 {
218   myErrorStatus=0;
219   myWarningStatus=0;
220   myLS.Clear();
221   myMSS.Clear();
222   //
223   CheckData();
224   if(myErrorStatus) {
225     return;
226   }
227   //
228   myClsf->SetTolerance(myTolerance);
229   //
230   // 1
231   ProcessVertices();
232   if(myErrorStatus) {
233     return;
234   }
235   if (myShapeType==TopAbs_VERTEX) {
236     return;
237   }
238   //
239   // 2
240   ProcessEdges();
241   if(myErrorStatus) {
242     return;
243   }
244   if (myShapeType==TopAbs_EDGE) {
245     return;
246   }
247   //
248   // 3
249   ProcessFaces();
250   if(myErrorStatus) {
251     return;
252   }
253   if (myShapeType==TopAbs_FACE) {
254     return;
255   }
256   //
257   // 4
258   ProcessSolids(); 
259   //
260 }
261 //=======================================================================
262 //function : CheckData
263 //purpose  : 
264 //=======================================================================
265   void GEOMAlgo_FinderShapeOn2::CheckData()
266 {
267   Standard_Integer iErr;
268   //
269   myErrorStatus=0;
270   //
271   if(myClsf.IsNull()) {
272     myErrorStatus=10; // myClsf=NULL
273     return;
274   }
275   //
276   myClsf->CheckData();
277   iErr=myClsf->ErrorStatus();
278   if (iErr) {
279     myErrorStatus=41; // invalid data for classifier
280     return;
281   }
282   //
283   if (myShape.IsNull()) {
284     myErrorStatus=11; // myShape=NULL
285     return;
286   }
287   //
288   if (!(myShapeType==TopAbs_VERTEX ||
289         myShapeType==TopAbs_EDGE ||
290         myShapeType==TopAbs_FACE ||
291         myShapeType==TopAbs_SOLID)) {
292     myErrorStatus=12; // unallowed subshape type
293     return;
294   }
295   //
296   if (myState==GEOMAlgo_ST_UNKNOWN || 
297       myState==GEOMAlgo_ST_INOUT) {
298     myErrorStatus=13; // unallowed state type
299     return;
300   }
301 }
302 //=======================================================================
303 //function : ProcessVertices
304 //purpose  : 
305 //=======================================================================
306   void GEOMAlgo_FinderShapeOn2::ProcessVertices()
307 {
308   myErrorStatus=0;
309   //
310   Standard_Boolean bIsConformState;
311   Standard_Integer i, aNb, iErr;
312   gp_Pnt aP;
313   TopTools_IndexedMapOfShape aM;
314   TopAbs_State aSt;
315   //
316   TopExp::MapShapes(myShape, TopAbs_VERTEX, aM);
317   aNb=aM.Extent();
318   for (i=1; i<=aNb; ++i) {
319     const TopoDS_Vertex& aV=TopoDS::Vertex(aM(i));
320     aP=BRep_Tool::Pnt(aV);
321     //
322     myClsf->SetPnt(aP);
323     myClsf->Perform();
324     iErr=myClsf->ErrorStatus();
325     if (iErr) {
326       myErrorStatus=40; // point can not be classified
327       return;
328     }
329     //
330     aSt=myClsf->State();
331     bIsConformState=GEOMAlgo_SurfaceTools::IsConformState(aSt, myState);
332     //
333     if (myShapeType==TopAbs_VERTEX){
334       if (bIsConformState) {
335         myMSS.Add(aV, aSt);
336       }
337     }
338     else if (bIsConformState || aSt==TopAbs_ON) {
339       myMSS.Add(aV, aSt);
340     }
341   }
342 }
343 //=======================================================================
344 //function : ProcessEdges
345 //purpose  : 
346 //=======================================================================
347   void GEOMAlgo_FinderShapeOn2::ProcessEdges()
348 {
349   myErrorStatus=0;
350   //
351   Standard_Boolean bIsConformState, bIsToBreak;
352   Standard_Integer i, aNb, iCnt, iErr;
353   TopAbs_State aSt;
354   TopTools_IndexedMapOfShape aM;
355   TopExp_Explorer aExp;
356   GEOMAlgo_ListIteratorOfListOfPnt aIt;
357   //
358   TopExp::MapShapes(myShape, TopAbs_EDGE, aM);
359   aNb=aM.Extent();
360   for (i=1; i<=aNb; ++i) {
361     GEOMAlgo_ListOfPnt aLP;
362     GEOMAlgo_StateCollector aSC;
363     //
364     const TopoDS_Edge& aE=TopoDS::Edge(aM(i));
365     //
366     aExp.Init(aE, TopAbs_VERTEX);
367     for (; aExp.More(); aExp.Next()) {
368       const TopoDS_Shape& aV=aExp.Current();
369       //
370       bIsConformState=myMSS.Contains(aV);
371       if (!bIsConformState) {
372         break;// vertex has non-conformed state
373       }
374       else {
375         aSt=myMSS.FindFromKey(aV);
376         aSC.AppendState(aSt);
377       }
378     }
379     //
380     if (!bIsConformState) {
381       continue; // vertex has non-conformed state,skip edge
382     }
383     //
384     if (BRep_Tool::Degenerated(aE)) {
385       myMSS.Add(aE, aSt); 
386       continue;
387     }
388     //
389     if (myState==GEOMAlgo_ST_ON) {
390       Standard_Boolean bCanBeON;
391       Standard_Real aT1, aT2;
392       Handle(Geom_Curve) aC; 
393       //
394       aC=BRep_Tool::Curve(aE, aT1, aT2);
395       bCanBeON=myClsf->CanBeON(aC);
396       if(!bCanBeON) {
397         continue;
398       }
399     }
400     //
401     InnerPoints(aE, aLP);
402     if (myErrorStatus) {
403       return;  
404     }
405     //
406     bIsConformState=Standard_True;
407     aIt.Initialize(aLP);
408     for (iCnt=0; aIt.More(); aIt.Next(), ++iCnt) {
409       if (myNbPntsMax) {
410         if (iCnt > myNbPntsMax) {
411           break;
412         }
413       }
414       //
415       const gp_Pnt& aP=aIt.Value();
416       //
417       myClsf->SetPnt(aP);
418       myClsf->Perform();
419       iErr=myClsf->ErrorStatus();
420       if (iErr) {
421         myErrorStatus=40; // point can not be classified
422         return;
423       }
424       //
425       aSt=myClsf->State();
426       //
427       bIsToBreak=aSC.AppendState(aSt);
428       if (bIsToBreak) {
429         break;
430       }
431     }
432     //
433     aSt=aSC.State();    
434     //
435     bIsConformState=GEOMAlgo_SurfaceTools::IsConformState(aSt, myState);
436     if (myShapeType==TopAbs_EDGE) {
437       if (bIsConformState) {
438         myMSS.Add(aE, aSt);
439       }
440     }
441     else if (bIsConformState || aSt==TopAbs_ON) {
442       myMSS.Add(aE, aSt);
443     }
444   } //  for (i=1; i<=aNb; ++i) next edge 
445 }
446 //=======================================================================
447 //function : ProcessFaces
448 //purpose  : 
449 //=======================================================================
450   void GEOMAlgo_FinderShapeOn2::ProcessFaces()
451 {
452   myErrorStatus=0;
453   //
454   Standard_Boolean bIsConformState, bIsToBreak, bCanBeON;
455   Standard_Integer i, aNbF, iCnt, iErr;
456   TopAbs_State aSt;
457   TopTools_IndexedMapOfShape aM;
458   TopExp_Explorer aExp;
459   GEOMAlgo_ListIteratorOfListOfPnt aIt;
460   //
461   TopExp::MapShapes(myShape, TopAbs_FACE, aM);
462   aNbF=aM.Extent();
463   for (i=1; i<=aNbF; ++i) {
464     GEOMAlgo_StateCollector aSC;
465     GEOMAlgo_ListOfPnt aLP;
466     //
467     const TopoDS_Face& aF=TopoDS::Face(aM(i));
468     //
469     if (myState==GEOMAlgo_ST_ON) {
470       Handle(Geom_Surface) aS;
471       //
472       aS=BRep_Tool::Surface(aF);
473       bCanBeON=myClsf->CanBeON(aS);
474       if(!bCanBeON) {
475         continue;
476       }
477     }
478     //
479     aExp.Init(aF, TopAbs_EDGE);
480     for (; aExp.More(); aExp.Next()) {
481       const TopoDS_Shape& aE=aExp.Current();
482       bIsConformState=myMSS.Contains(aE);
483       if (!bIsConformState) {
484         break;// edge has non-conformed state
485       }
486       else {
487         aSt=myMSS.FindFromKey(aE);
488         aSC.AppendState(aSt);
489       }
490     }
491     //
492     if (!bIsConformState) {
493       continue; // edge has non-conformed state,skip face
494     }
495     //
496     InnerPoints(aF, aLP);
497     if (myErrorStatus) {
498       return;
499     }
500     //
501     bIsConformState=Standard_True;
502     aIt.Initialize(aLP);
503     for (iCnt=0; aIt.More(); aIt.Next(), ++iCnt) {
504       if (myNbPntsMax) {
505         if (iCnt > myNbPntsMax) {
506           break;
507         }
508       }
509       //
510       const gp_Pnt& aP=aIt.Value();
511       //
512       myClsf->SetPnt(aP);
513       myClsf->Perform();
514       iErr=myClsf->ErrorStatus();
515       if (iErr) {
516         myErrorStatus=40; // point can not be classified
517         return;
518       }
519       //
520       aSt=myClsf->State();
521       //
522       bIsToBreak=aSC.AppendState(aSt);
523       if (bIsToBreak) {
524         break;
525       }
526     }
527     //
528     aSt=aSC.State();    
529     //
530     bIsConformState=GEOMAlgo_SurfaceTools::IsConformState(aSt, myState);
531     if (myShapeType==TopAbs_FACE) {
532       if (bIsConformState) {
533         myMSS.Add(aF, aSt);
534       }
535     }
536     else if (bIsConformState || aSt==TopAbs_ON) {
537       myMSS.Add(aF, aSt);
538     }
539   }//  for (i=1; i<=aNb; ++i) next face 
540 }
541 //=======================================================================
542 //function : ProcessSolids
543 //purpose  : 
544 //=======================================================================
545   void GEOMAlgo_FinderShapeOn2::ProcessSolids()
546 {
547   myErrorStatus=0;
548   //
549   Standard_Boolean bIsConformState;
550   Standard_Integer i, aNbS, j, aNbF;
551   TopTools_IndexedMapOfShape aM, aMF;
552   TopAbs_State aSt;
553   //
554   TopExp::MapShapes(myShape, TopAbs_SOLID, aM);
555   aNbS=aM.Extent();
556   for (i=1; i<=aNbS; ++i) {
557     GEOMAlgo_StateCollector aSC;
558     //
559     const TopoDS_Shape& aSd=aM(i);
560     aMF.Clear();
561     TopExp::MapShapes(aSd, TopAbs_FACE, aMF);
562     aNbF=aMF.Extent();
563     for (j=1; j<=aNbF; ++j) {
564       const TopoDS_Shape& aF=aMF(j);
565       bIsConformState=myMSS.Contains(aF);
566       if (!bIsConformState) {
567         break;// face has non-conformed state
568       }
569       else {
570         aSt=myMSS.FindFromKey(aF);
571         aSC.AppendState(aSt);
572       }
573     }
574     //
575     if (!bIsConformState) {
576       continue; // face has non-conformed state,skip solid
577     }
578     //
579     aSt=aSC.State();    
580     //
581     bIsConformState=GEOMAlgo_SurfaceTools::IsConformState(aSt, myState);
582     if (bIsConformState) {
583       myMSS.Add(aSd, aSt);
584     }
585   }
586 }
587 //
588 //=======================================================================
589 //function : InnerPoints
590 //purpose  : 
591 //=======================================================================
592   void GEOMAlgo_FinderShapeOn2::InnerPoints(const TopoDS_Face& aF,
593                                           GEOMAlgo_ListOfPnt& aLP)
594 {
595   myErrorStatus=0;
596   //
597   Standard_Integer j, j1, j2, k, n[4], aNbLinks, aNx, aNb, iCnt;//, aNbMax, *pIds;
598   TopLoc_Location aLoc;
599   Handle(Poly_Triangulation) aTRF;
600   TColStd_MapOfInteger aMBN;
601   GEOMAlgo_DataMapOfPassKeyInteger aMPKI;
602   GEOMAlgo_DataMapIteratorOfDataMapOfPassKeyInteger aIt;
603   gp_Pnt aP, aP1, aP2;
604   //
605   aLP.Clear();
606   //
607   aTRF=BRep_Tool::Triangulation(aF, aLoc);
608   if (aTRF.IsNull()) {
609     if (!GEOMAlgo_FinderShapeOn::BuildTriangulation(aF)) {
610       myWarningStatus=20; // no triangulation found
611       return;
612     }
613     aTRF=BRep_Tool::Triangulation(aF, aLoc);
614   }
615   //
616   const gp_Trsf& aTrsf=aLoc.Transformation();
617   const Poly_Array1OfTriangle& aTrs=aTRF->Triangles();
618   const TColgp_Array1OfPnt& aNodes=aTRF->Nodes();
619   //
620   // map link/nbtriangles
621   j1=aTrs.Lower();
622   j2=aTrs.Upper();
623   for (j=j1; j<=j2; ++j) {
624     const Poly_Triangle& aTr=aTrs(j);
625     aTr.Get(n[0], n[1], n[2]);
626     n[3]=n[0];
627     for (k=0; k<3; ++k) {
628       GEOMAlgo_PassKey aPK;
629       //
630       aPK.SetIds(n[k], n[k+1]);
631       if (aMPKI.IsBound(aPK)) {
632         Standard_Integer& iCntX=aMPKI.ChangeFind(aPK);
633         ++iCntX;
634       }
635       else {
636         aMPKI.Bind(aPK, 1);
637       }
638     }
639   }
640   //
641   // boundary nodes aMBN
642   aNbLinks=aMPKI.Extent();
643   aIt.Initialize(aMPKI);
644   for (; aIt.More(); aIt.Next()) {
645     iCnt=aIt.Value();
646     if (iCnt==1) {
647       const GEOMAlgo_PassKey& aPK=aIt.Key();
648       //qf
649       /*
650       aNbMax=aPK.NbMax();
651       pIds=(Standard_Integer*)aPK.Key();
652       for (k=1; k<3; ++k) {
653         aNx=*(pIds+aNbMax-k);
654         aMBN.Add(aNx);
655       }
656       */
657       aNx=(Standard_Integer)aPK.Id(1);
658       aMBN.Add(aNx);
659       aNx=(Standard_Integer)aPK.Id(2);
660       aMBN.Add(aNx);
661       //qt
662     }
663   }
664   //
665   // inner nodes=all_nodes - boundary_nodes 
666   j1=aNodes.Lower();
667   j2=aNodes.Upper();
668   for (j=j1; j<=j2; ++j) {
669     if (!aMBN.Contains(j)) {
670       aP=aNodes(j).Transformed(aTrsf);
671       aLP.Append(aP);
672     }
673   }
674   //
675   aNb=aLP.Extent();
676   //
677   if (!aNb && myNbPntsMin) { 
678     // try to fill it yourself
679     Standard_Boolean bIsDone;
680     Standard_Integer aN1, aN2;
681     Handle(Geom_Surface) aS;
682     GeomAdaptor_Surface aGAS;
683     GeomAbs_SurfaceType aType;
684     //
685     aS=BRep_Tool::Surface(aF);
686     aGAS.Load(aS);
687     aType=aGAS.GetType();
688     if (aType==GeomAbs_Plane || aType==GeomAbs_Cylinder) {
689       // inner links
690       aNbLinks=aMPKI.Extent();
691       aIt.Initialize(aMPKI);
692       for (; aIt.More(); aIt.Next()) {
693         iCnt=aIt.Value();
694         if (iCnt>1) {
695           // take the first having occured inner link
696           // and discretize it
697           const GEOMAlgo_PassKey& aPK=aIt.Key();
698           //qf
699           /*
700           aNbMax=aPK.NbMax();
701           pIds=(Standard_Integer*)aPK.Key();
702           aN1=*(pIds+aNbMax-1);
703           aN2=*(pIds+aNbMax-2);
704           */
705           //
706           aN1=(Standard_Integer)aPK.Id(1);
707           aN2=(Standard_Integer)aPK.Id(2);
708           //qt
709           aP1=aNodes(aN1).Transformed(aTrsf);
710           aP2=aNodes(aN2).Transformed(aTrsf);
711           //
712           if (aType==GeomAbs_Cylinder) {
713             Standard_Real aTolSM;
714             gp_Cylinder aCyl;
715             //
716             aTolSM=1.523e-6;//~1.-cos(0.1 deg)
717             aCyl=aGAS.Cylinder();
718             if (!GEOMAlgo_SurfaceTools::IsCoaxial(aP1, aP2, aCyl, aTolSM)) {
719               continue;
720             }
721           }
722           //
723           BRepLib_MakeEdge aBME(aP1, aP2);
724           bIsDone=aBME.IsDone();
725           if (!bIsDone) {
726             myErrorStatus=30; //can not obtain the line fron the link
727             return;
728           }
729           //
730           const TopoDS_Shape& aSx=aBME.Shape();
731           const TopoDS_Edge& aE=TopoDS::Edge(aSx);
732           //
733           InnerPoints(aE, myNbPntsMin, aLP);
734           break;
735         }// if (iCnt>1)
736       }// for (; aIt.More(); aIt.Next())
737     }// if (aType==GeomAbs_Plane || aType==GeomAbs_Cylinder) 
738   }// if (!aNb && myNbPntsMin) { 
739 }
740 //=======================================================================
741 //function : InnerPoints
742 //purpose  : 
743 //=======================================================================
744   void GEOMAlgo_FinderShapeOn2::InnerPoints(const TopoDS_Edge& aE,
745                                           GEOMAlgo_ListOfPnt& aLP)
746 {
747   myErrorStatus=0;
748   //
749   Standard_Integer j, aNbNodes, aIndex, aNb;
750   Handle(Poly_PolygonOnTriangulation) aPTE;
751   Handle(Poly_Triangulation) aTRE;
752   TopLoc_Location aLoc;
753   gp_Pnt aP;
754   //
755   aLP.Clear();
756   BRep_Tool::PolygonOnTriangulation(aE, aPTE, aTRE, aLoc);
757   if (aTRE.IsNull() || aPTE.IsNull()) {
758     Handle(Poly_Polygon3D) aPE = BRep_Tool::Polygon3D(aE, aLoc);
759     if (aPE.IsNull()) {
760       if (!GEOMAlgo_FinderShapeOn::BuildTriangulation(aE)) {
761         myErrorStatus=20; // no triangulation found
762         return;
763       }
764       aPE = BRep_Tool::Polygon3D(aE, aLoc);
765     }
766     const gp_Trsf& aTrsf=aLoc.Transformation();
767     const TColgp_Array1OfPnt& aNodes=aPE->Nodes();
768     //
769     aNbNodes=aPE->NbNodes();
770     Standard_Integer low = aNodes.Lower(), up = aNodes.Upper();
771     for (j=low+1; j<up; ++j) {
772       aP=aNodes(j).Transformed(aTrsf);
773       aLP.Append(aP);
774     }
775   }
776   else {
777     const gp_Trsf& aTrsf=aLoc.Transformation();
778     const TColgp_Array1OfPnt& aNodes=aTRE->Nodes();
779     //
780     aNbNodes=aPTE->NbNodes();
781     const TColStd_Array1OfInteger& aInds=aPTE->Nodes();
782     for (j=2; j<aNbNodes; ++j) {
783       aIndex=aInds(j);
784       aP=aNodes(aIndex).Transformed(aTrsf);
785       aLP.Append(aP);
786     }
787   }
788   //
789   aNb=aLP.Extent();
790   if (!aNb && myNbPntsMin) { 
791     // try to fill it yourself
792     InnerPoints(aE, myNbPntsMin, aLP);
793     aNb=aLP.Extent();
794   }
795 }
796 //=======================================================================
797 //function : InnerPoints
798 //purpose  : 
799 //=======================================================================
800   void GEOMAlgo_FinderShapeOn2::InnerPoints(const TopoDS_Edge& aE,
801                                           const Standard_Integer aNbPntsMin,
802                                           GEOMAlgo_ListOfPnt& aLP)
803 {
804   // try to fill it yourself
805   Standard_Boolean bInf1, bInf2;
806   Standard_Integer j, aNbT;
807   Standard_Real dT, aT, aT1, aT2;
808   gp_Pnt aP;
809   Handle(Geom_Curve) aC3D;
810   //
811   aC3D=BRep_Tool::Curve(aE, aT1, aT2);
812   if (aC3D.IsNull()) {
813     return;
814   }
815   //
816   bInf1=Precision::IsNegativeInfinite(aT1);
817   bInf2=Precision::IsPositiveInfinite(aT2);
818   if (bInf1 || bInf2) {
819     return;
820   }
821   //
822   aNbT=myNbPntsMin+1;
823   dT=(aT2-aT1)/aNbT;
824   for (j=1; j<=aNbPntsMin; ++j) {
825     aT=aT1+j*dT;
826     aC3D->D0(aT, aP);
827     aLP.Append(aP);
828   }
829 }
830
831 //
832 // myErrorStatus :
833 //
834 // 10 -myClsf=NULL
835 // 11 -myShape=NULL
836 // 12 -unallowed type of subshapes 
837 // 13 -unallowed state  
838 // 15 -unallowed surface type
839 // 20- no triangulation found
840 // 30- can not obtain the line from the link
841 // 40- point can not be classified
842 // 41- invalid data for classifier