1 // Created on: 1996-03-07
2 // Created by: Jean Yves LEBEY
3 // Copyright (c) 1996-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
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.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
18 #include <BRep_Tool.hxx>
19 #include <BRepClass3d_SolidExplorer.hxx>
20 #include <BRepTopAdaptor_FClass2d.hxx>
22 #include <Standard_NoSuchObject.hxx>
23 #include <TCollection_AsciiString.hxx>
25 #include <TopoDS_Edge.hxx>
26 #include <TopoDS_Face.hxx>
27 #include <TopoDS_Shape.hxx>
28 #include <TopoDS_Vertex.hxx>
29 #include <TopOpeBRepBuild_Builder.hxx>
30 #include <TopOpeBRepBuild_define.hxx>
31 #include <TopOpeBRepBuild_EdgeBuilder.hxx>
32 #include <TopOpeBRepBuild_FaceBuilder.hxx>
33 #include <TopOpeBRepBuild_FuseFace.hxx>
34 #include <TopOpeBRepBuild_GTopo.hxx>
35 #include <TopOpeBRepBuild_HBuilder.hxx>
36 #include <TopOpeBRepBuild_PaveSet.hxx>
37 #include <TopOpeBRepBuild_ShapeSet.hxx>
38 #include <TopOpeBRepBuild_ShellFaceSet.hxx>
39 #include <TopOpeBRepBuild_SolidBuilder.hxx>
40 #include <TopOpeBRepBuild_WireEdgeSet.hxx>
41 #include <TopOpeBRepDS_BuildTool.hxx>
42 #include <TopOpeBRepDS_CurveIterator.hxx>
43 #include <TopOpeBRepDS_EXPORT.hxx>
44 #include <TopOpeBRepDS_HDataStructure.hxx>
45 #include <TopOpeBRepDS_PointIterator.hxx>
46 #include <TopOpeBRepDS_samdom.hxx>
47 #include <TopOpeBRepDS_SurfaceIterator.hxx>
48 #include <TopOpeBRepTool_EXPORT.hxx>
49 #include <TopOpeBRepTool_SC.hxx>
50 #include <TopOpeBRepTool_ShapeExplorer.hxx>
53 #define DEBSHASET(sarg,meth,shaset,str) TCollection_AsciiString sarg((meth));(sarg)=(sarg)+(shaset).DEBNumber()+(str);
54 Standard_EXPORT Standard_Boolean TopOpeBRepBuild_GetcontextNOFUFA();
55 Standard_EXPORT void debffsfs(const Standard_Integer i){std::cout<<"+++ debffsfs "<<i<<std::endl;}
56 Standard_EXPORT void debffflo(const Standard_Integer i){std::cout<<"+++ debffflo "<<i<<std::endl;}
59 static Standard_Boolean STATIC_motheropedef = Standard_False; static TopOpeBRepBuild_GTopo STATIC_Gmotherope;
60 Standard_EXPORT void FUN_setmotherope(const TopOpeBRepBuild_GTopo& G) {STATIC_Gmotherope = G;STATIC_motheropedef = Standard_True;}
61 Standard_EXPORT void FUN_unsetmotherope() { STATIC_motheropedef = Standard_False; }
62 Standard_EXPORT Standard_Boolean FUN_ismotheropedef() { return STATIC_motheropedef; }
63 Standard_EXPORT const TopOpeBRepBuild_GTopo& FUN_motherope() { return STATIC_Gmotherope; }
65 //Standard_IMPORT extern Standard_Boolean GLOBAL_classifysplitedge;
66 Standard_EXPORTEXTERN Standard_Boolean GLOBAL_classifysplitedge;
67 //Standard_IMPORT extern Standard_Boolean GLOBAL_revownsplfacori;
68 Standard_EXPORTEXTERN Standard_Boolean GLOBAL_revownsplfacori;
69 Standard_EXPORT void FUNBUILD_ANCESTORRANKPREPARE
70 (TopOpeBRepBuild_Builder& B,const TopTools_ListOfShape& LF1,const TopTools_ListOfShape& LF2,const TopOpeBRepDS_Config c1,const TopOpeBRepDS_Config c2);
71 Standard_EXPORT void FUNBUILD_ANCESTORRANKGET(TopOpeBRepBuild_Builder& B,const TopoDS_Shape& f,Standard_Boolean& of1,Standard_Boolean& of2);
73 static Standard_Integer FUN_getAncestorFsp(TopOpeBRepBuild_Builder& B,TopOpeBRepTool_ShapeClassifier& SC,const TopTools_ListOfShape& LF,const TopoDS_Shape& fsp,Standard_Boolean& p3ddef,gp_Pnt& p3d);
74 static Standard_Integer FUN_getAncestorFsp(TopOpeBRepBuild_Builder& B,TopOpeBRepTool_ShapeClassifier& SC,const TopTools_ListOfShape& LF1,const TopTools_ListOfShape& LF2,const TopoDS_Shape& fsp);
76 //static void FUN_getAncestorFsp(const Handle(TopOpeBRepDS_HDataStructure)& HDS,TopOpeBRepTool_ShapeClassifier& SC,const TopTools_ListOfShape& LF1,const TopTools_ListOfShape& LF2,const TopTools_ListOfShape& spFOR,
77 // TopTools_DataMapOfShapeInteger* SplitAnc);
80 static Standard_Integer FUN_getAncestorFsp(TopOpeBRepBuild_Builder& B,TopOpeBRepTool_ShapeClassifier& SC,const TopTools_ListOfShape& LF,const TopoDS_Shape& fsp,Standard_Boolean& p3ddef,gp_Pnt& p3d)
82 const TopOpeBRepDS_DataStructure& BDS = B.DataStructure()->DS(); // How to do static <--> const
84 TopTools_ListIteratorOfListOfShape itf(LF);
85 for (; itf.More(); itf.Next()){
86 const TopoDS_Face& f = TopoDS::Face(itf.Value());
87 TopAbs_State st = SC.StateShapeShape(fsp,f,1);
88 if ((st == TopAbs_UNKNOWN) || (st == TopAbs_OUT)) continue;
89 if (st == TopAbs_ON) {
91 Standard_Boolean ok = BRepClass3d_SolidExplorer::FindAPointInTheFace(TopoDS::Face(fsp),p3d);
93 p3ddef = Standard_True;
95 gp_Pnt2d p2d;Standard_Real dd = 0.; Standard_Boolean ok = FUN_tool_projPonF(p3d,f,p2d,dd);
97 Standard_Real tolf = BRep_Tool::Tolerance(f)*1.e1;
98 if (dd > tolf) return 0;
99 Standard_Real TolClass = 1e-8; BRepTopAdaptor_FClass2d FClass2d(f,TolClass);
100 st = FClass2d.Perform(p2d);
102 if (st == TopAbs_IN) {
103 Standard_Integer ianc = BDS.Shape(f);
110 static Standard_Integer FUN_getAncestorFsp(TopOpeBRepBuild_Builder& B,TopOpeBRepTool_ShapeClassifier& SC,
111 const TopTools_ListOfShape& LF1,const TopTools_ListOfShape& LF2,const TopoDS_Shape& fsp)
113 // IMPORTANT : fsp is split IN/OUT of F1,so it has only one ancestor face
114 // LF1 faces sdm with LF2
115 const TopOpeBRepDS_DataStructure& BDS = B.DataStructure()->DS();
116 Standard_Boolean of1,of2; FUNBUILD_ANCESTORRANKGET(B,fsp,of1,of2);
117 Standard_Integer rkfsp = 0;
118 if (of1 && !of2) rkfsp = 1;
119 else if (of2 && !of1) rkfsp = 2;
120 Standard_Boolean unk = (rkfsp == 0);
122 Standard_Integer rkf1 = BDS.AncestorRank(LF1.First());
123 Standard_Integer rkf2 = BDS.AncestorRank(LF2.First());
124 Standard_Boolean p3ddef = Standard_False; gp_Pnt p3d;
126 Standard_Boolean ison1 = (rkf1 == rkfsp);
127 Standard_Boolean ison2 = (rkf2 == rkfsp);
129 Standard_Integer ianc1=0,ianc2=0;
131 ianc1 = FUN_getAncestorFsp(B,SC,LF1,fsp,p3ddef,p3d);
132 if (ison1) return ianc1;
135 ianc2 = FUN_getAncestorFsp(B,SC,LF2,fsp,p3ddef,p3d);
136 if (ison2) return ianc2;
138 if (ianc1 + ianc2 >0) {
139 if (ianc1 == 0) return ianc2;
140 else if (ianc2 == 0) return ianc1;
141 else return 0; // fsp has 2 ancestor faces
146 Standard_EXPORT TopTools_DataMapOfShapeInteger* GLOBAL_SplitAnc = NULL; //xpu260598
148 static void FUN_getAncestorFsp(TopOpeBRepBuild_Builder& B,TopOpeBRepTool_ShapeClassifier& SC,
149 const TopTools_ListOfShape& LF1,const TopTools_ListOfShape& LF2,
150 const TopoDS_Shape& FOR,TopTools_DataMapOfShapeInteger* SplitAnc)
152 if (SplitAnc == NULL) return;
154 Standard_Boolean issplitIN = B.IsSplit(FOR,TopAbs_IN);
155 Standard_Boolean issplitOU = B.IsSplit(FOR,TopAbs_OUT);
156 if (!issplitIN && !issplitOU) return;
158 TopTools_ListOfShape spFOR;
159 if (issplitIN) FDS_copy(B.Splits(FOR,TopAbs_IN),spFOR);
160 if (issplitOU) FDS_copy(B.Splits(FOR,TopAbs_OUT),spFOR);
162 for (TopTools_ListIteratorOfListOfShape itsp(spFOR); itsp.More(); itsp.Next()){
163 const TopoDS_Shape& fsp = itsp.Value();
165 Standard_Boolean isbound = SplitAnc->IsBound(fsp);
166 if (isbound) continue;
168 Standard_Integer ianc = FUN_getAncestorFsp(B,SC,LF1,LF2,fsp);
170 SplitAnc->Bind(fsp,ianc);
175 Standard_EXPORT TopTools_ListOfShape* GLOBAL_lfr1 = NULL;
176 Standard_EXPORT Standard_Boolean GLOBAL_lfrtoprocess = Standard_False;
177 //Standard_IMPORT extern TopTools_ListOfShape* GLOBAL_lfr1;
178 //Standard_IMPORT TopTools_ListOfShape* GLOBAL_lfr1;
179 //Standard_IMPORT extern Standard_Boolean GLOBAL_lfrtoprocess;
180 //Standard_IMPORT Standard_Boolean GLOBAL_lfrtoprocess;
182 //=======================================================================
183 //function : GFillFaceSFS
185 //=======================================================================
186 void TopOpeBRepBuild_Builder::GFillFaceSFS(const TopoDS_Shape& FOR,const TopTools_ListOfShape& LSO2,const TopOpeBRepBuild_GTopo& Gin,TopOpeBRepBuild_ShellFaceSet& SFS)
188 TopAbs_State TB1,TB2; Gin.StatesON(TB1,TB2);
191 Standard_Integer iF; Standard_Boolean tSPS = GtraceSPS(FOR,iF);
192 if(tSPS)std::cout<<std::endl;
194 const TopOpeBRepDS_DataStructure& BDS = myDataStructure->DS();
195 Standard_Boolean tosplit = GToSplit(FOR,TB1);
196 Standard_Boolean tomerge = GToMerge(FOR);
198 // Standard_Integer iFOR = BDS.Shape(FOR);
200 Standard_Integer rkFOR = BDS.AncestorRank(FOR);
204 GdumpSHASTA(FOR,TB1,"--- GFillFaceSFS START ");
205 std::cout<<" tosplit "<<tosplit<<" tomerge "<<tomerge<<std::endl;
210 TopoDS_Shape FF = FOR; FF.Orientation(TopAbs_FORWARD);
211 Standard_Boolean hsd = myDataStructure->HasSameDomain(FOR); //xpu280598
212 GLOBAL_lfrtoprocess = Standard_False;
214 if ( tosplit && tomerge ) {
216 // merge des faces SameDomain
218 // on effectue le merge ssi FOR est la reference de ses faces SameDomain
219 Standard_Integer iref = myDataStructure->SameDomainReference(FOR);
220 const TopoDS_Shape& fref = myDataStructure->Shape(iref);
221 Standard_Boolean isref = FOR.IsSame(fref);
223 Standard_Boolean makemerge = isref;
224 // xpu280199 PRO16958 : f10=ref(f6), GToMerge(f10)=false
225 Standard_Boolean makemergeref = GToMerge(fref);
226 makemerge = makemerge || (!makemergeref && (rkFOR == 1));
232 GdumpSHASTA(FOR,TB1,"[[[[[[[[[[[[[[[[[[[[[[[[[[ GFillFaceSFS makemerge START ");
233 std::cout<<std::endl;
237 Standard_Boolean performfufa = Standard_True;
239 performfufa = !TopOpeBRepBuild_GetcontextNOFUFA();
242 GLOBAL_lfrtoprocess = Standard_True;
243 if (GLOBAL_lfrtoprocess) {
244 if (GLOBAL_lfr1==NULL) GLOBAL_lfr1=(TopTools_ListOfShape*)new TopTools_ListOfShape();
245 GLOBAL_lfr1->Clear();
249 // xpu280598 : Filling up GLOBAL_SplitAnc = {(fsp,ifanc)}
250 // . fsp = spIN/OU(fanc),
251 // . fanc hsdm is the unique ancestor face
252 if ( GLOBAL_SplitAnc == NULL ) GLOBAL_SplitAnc = (TopTools_DataMapOfShapeInteger*) new TopTools_DataMapOfShapeInteger();
253 GLOBAL_SplitAnc->Clear();
255 TopTools_ListOfShape LFSO,LFDO,LFSO1,LFDO1,LFSO2,LFDO2;
256 GFindSamDomSODO(FF,LFSO,LFDO); // -980617
257 // FDSSDM_sordor(FF,LFSO,LFDO);
258 Standard_Integer rankF=GShapeRank(FF),rankX=(rankF)?((rankF==1)?2:1):0;
259 GFindSameRank(LFSO,rankF,LFSO1);
260 GFindSameRank(LFDO,rankF,LFDO1);
261 GFindSameRank(LFSO,rankX,LFSO2);
262 GFindSameRank(LFDO,rankX,LFDO2);
263 TopTools_ListOfShape LF1,LF2;
265 TopOpeBRepBuild_GTopo GM;
267 Standard_Boolean dodo;
268 Standard_Integer l1,l2;
270 // WES : toutes les faces de meme orientation topologique
271 LF1 = LFSO1; // NYI pointeurs
272 LF2 = LFSO2; // NYI pointeurs
273 l1 = LF1.Extent(); l2 = LF2.Extent();
274 dodo = (l1!=0) && (l2!=0);
276 FUN_unsetmotherope(); // +12/07
279 GM.ChangeConfig(TopOpeBRepDS_SAMEORIENTED,TopOpeBRepDS_SAMEORIENTED);
283 GdumpSAMDOM(LF1, (char *) "LF1 (LFSO1) : ");
284 GdumpSAMDOM(LF2, (char *) "LF2 (LFSO2) : ");
287 GLOBAL_classifysplitedge = Standard_True;
288 GFillFacesWESMakeFaces(LF1,LF2,LSO2,GM);
289 GLOBAL_classifysplitedge = Standard_False;
291 GLOBAL_revownsplfacori = Standard_True;
292 FUNBUILD_ANCESTORRANKPREPARE
293 (*this,LF1,LF2,TopOpeBRepDS_SAMEORIENTED,TopOpeBRepDS_SAMEORIENTED);
294 if (hsd) FUN_getAncestorFsp((*this),myShapeClassifier,
295 LF1,LF2,FOR,GLOBAL_SplitAnc); //xpu280598
297 // GLOBAL_lfrtoprocess = t
298 // ==> on ne stocke PAS les faces 'startelement' dans le SFS
299 // mais dans GLOBAL_lfr1,
300 // GLOBAL_lfrtoprocess = f
301 // ==> on stocke normalement dans le SFS et pas dans GLOBAL_lfr1
302 // NYI : + argument a la methode GSplitFaceSFS ?? a voir
304 // ici : GLOBAL_lfrtoprocess = t
305 if (GLOBAL_lfr1==NULL) GLOBAL_lfr1=(TopTools_ListOfShape*)new TopTools_ListOfShape(); //flo150998
306 GLOBAL_lfr1->Clear();
307 GSplitFaceSFS(FOR,LSO2,GM,SFS);
308 GLOBAL_lfrtoprocess = Standard_False;
309 GLOBAL_revownsplfacori = Standard_False;
312 // WES : FOR + faces d'orientation topologique opposee
313 LF1 = LFSO1; // NYI pointeurs
314 LF2 = LFDO2; // NYI pointeurs
315 l1 = LF1.Extent(); l2 = LF2.Extent();
316 dodo = (l1!=0) && (l2!=0);
319 GM = Gin; // OUT,OUT-> OUT,IN + IN OUT
321 NTB = TopAbs_IN; // NTB = (TB == TopAbs_OUT) ? TopAbs_IN : TopAbs_OUT;
322 GM.ChangeValue(TB,TopAbs_ON,Standard_False);
323 GM.ChangeValue(NTB,TopAbs_ON,Standard_True);
324 GM.ChangeConfig(TopOpeBRepDS_SAMEORIENTED,TopOpeBRepDS_DIFFORIENTED);
325 FUN_setmotherope(GM); // +12/07
329 TopAbs_State TB11,TB21; GM.StatesON(TB11,TB21);
330 std::cout<<std::endl;
331 std::cout <<"@@@@"<<std::endl;
332 std::cout <<"@@@@@@@@ partie 1 : ";TopAbs::Print(TB11,std::cout);std::cout<<" ";TopAbs::Print(TB21,std::cout);std::cout<<std::endl;
333 std::cout <<"@@@@"<<std::endl;
334 GdumpSAMDOM(LF1, (char *) "LF1 (LFSO1) : ");
335 GdumpSAMDOM(LF2, (char *) "LF2 (LFDO2) : ");
336 std::cout<<std::endl;
339 GLOBAL_classifysplitedge = Standard_True;
340 GFillFacesWESMakeFaces(LF1,LF2,LSO2,GM);
341 GLOBAL_classifysplitedge = Standard_False;
343 GLOBAL_revownsplfacori = Standard_True;
344 FUNBUILD_ANCESTORRANKPREPARE
345 (*this,LF1,LF2,TopOpeBRepDS_SAMEORIENTED,TopOpeBRepDS_DIFFORIENTED);
346 if (hsd) FUN_getAncestorFsp((*this),myShapeClassifier,
347 LF1,LF2,FOR,GLOBAL_SplitAnc); //xpu280598
350 Standard_Boolean issplitIN = IsSplit(FOR,TopAbs_IN);
352 const TopTools_ListOfShape& spFOR = Splits(FOR,TopAbs_IN);
353 TopTools_ListOfShape spFORcopy; FDS_copy(spFOR,spFORcopy);
355 TopTools_ListIteratorOfListOfShape it(LF1);
356 for (; it.More(); it.Next()) {
357 const TopoDS_Shape& f = it.Value();
358 Standard_Boolean issplit = IsSplit(f,TopAbs_IN);
359 if ( issplit ) ChangeSplit(f,TopAbs_IN).Clear();
362 for (; it.More(); it.Next()) {
363 const TopoDS_Shape& f = it.Value();
364 Standard_Boolean issplit = IsSplit(f,TopAbs_IN);
365 if ( issplit ) ChangeSplit(f,TopAbs_IN).Clear();
367 ChangeSplit(FOR,TopAbs_IN).Append(spFORcopy); // keep split for reference
371 GSplitFaceSFS(FOR,LSO2,GM,SFS);
372 GLOBAL_revownsplfacori = Standard_False;
377 GM = Gin; // OUT,OUT-> OUT,IN + IN OUT
379 // NTB = (TB == TopAbs_OUT) ? TopAbs_IN : TopAbs_OUT;
381 GM.ChangeValue(TopAbs_ON,TB,Standard_False);
382 GM.ChangeValue(TopAbs_ON,NTB,Standard_True);
383 GM.ChangeConfig(TopOpeBRepDS_SAMEORIENTED,TopOpeBRepDS_DIFFORIENTED);
384 FUN_setmotherope(GM); // +12/07
388 TopAbs_State TB12,TB22; GM.StatesON(TB12,TB22);
389 std::cout<<std::endl;
390 std::cout <<"@@@@"<<std::endl;
391 std::cout <<"@@@@@@@@ partie 2 : ";TopAbs::Print(TB12,std::cout);std::cout<<" ";TopAbs::Print(TB22,std::cout);std::cout<<std::endl;
392 std::cout <<"@@@@"<<std::endl;
393 GdumpSAMDOM(LF1, (char *) "LF1 (LFSO1) : ");
394 GdumpSAMDOM(LF2, (char *) "LF2 (LFDO2) : ");
395 std::cout<<std::endl;
398 GLOBAL_classifysplitedge = Standard_True;
399 GFillFacesWESMakeFaces(LF1,LF2,LSO2,GM);
400 GLOBAL_classifysplitedge = Standard_False;
402 GLOBAL_revownsplfacori = Standard_True;
403 FUNBUILD_ANCESTORRANKPREPARE
404 (*this,LF1,LF2,TopOpeBRepDS_SAMEORIENTED,TopOpeBRepDS_DIFFORIENTED);
405 if (hsd) FUN_getAncestorFsp((*this),myShapeClassifier,
406 LF1,LF2,FOR,GLOBAL_SplitAnc); //xpu280598
407 GSplitFaceSFS(FOR,LSO2,GM,SFS);
408 GLOBAL_revownsplfacori = Standard_False;
412 FUN_unsetmotherope(); // +12/07
416 GdumpSHASTA(FOR,TB1,"]]]]]]]]]]]]]]]]]]]]]]]]]] GFillFaceSFS makemerge END ");
417 std::cout<<std::endl;
421 GLOBAL_SplitAnc->Clear(); // xpu280598
424 SFS.ChangeStartShapes().Extent();
427 if(tSPS) debffflo(iF);
429 // const TopTools_ListOfShape& lou = Splits(FF,TopAbs_OUT); Standard_Integer nou = lou.Extent();
430 // const TopTools_ListOfShape& lin = Splits(FF,TopAbs_IN); Standard_Integer nin = lin.Extent();
431 // GCopyList(lou,*GLOBAL_lfr1);
432 // GCopyList(lin,*GLOBAL_lfr1);
434 // Standard_Integer nlfr1 = GLOBAL_lfr1->Extent();
437 // NYI : Builder += methode pour le process fufa
438 TopOpeBRepBuild_FuseFace fufa; TopTools_ListOfShape ldum; Standard_Integer addinternal = 1; // disparition
439 fufa.Init(ldum,*GLOBAL_lfr1,addinternal);
441 Standard_Boolean isdone = fufa.IsDone();
444 // Standard_Boolean ismodified = fufa.IsModified();
446 const TopTools_ListOfShape& lfr2 = fufa.LFuseFace();
448 // const TopTools_ListOfShape& lfr2 = *GLOBAL_lfr1
449 // les faces remplacantes
450 for (TopTools_ListIteratorOfListOfShape itlfr2(lfr2);itlfr2.More();itlfr2.Next()) {
451 const TopoDS_Shape& flfr2 = itlfr2.Value();
455 DEBSHASET(ss,"--- FillFaceSFS apres fufa",SFS," AddStartElement SFS+ face ");
456 GdumpSHA(flfr2,(Standard_Address)ss.ToCString());
457 std::cout<<" ";TopAbs::Print(TB1,std::cout)<<" : 1 face ";
458 TopAbs::Print(flfr2.Orientation(),std::cout); std::cout<<std::endl;
461 SFS.AddStartElement(flfr2);
463 } // performfufa (context)
466 } // tosplit && tomerge
468 else if ( tosplit && !tomerge ) {
469 GSplitFace(FOR,Gin,LSO2);
470 GSplitFaceSFS(FOR,LSO2,Gin,SFS);
472 else if ( !tosplit && !tomerge ) {
473 GSplitFaceSFS(FOR,LSO2,Gin,SFS);
480 GdumpSHASTA(FOR,TB1,"--- GFillFaceSFS END ");std::cout<<std::endl;