// Alternatively, this file may be used under the terms of Open CASCADE
// commercial license or contractual agreement.
-#include <stdio.h>
-
-#include <BRepFill_PipeShell.ixx>
#include <BRep_Builder.hxx>
#include <BRep_Tool.hxx>
-#include <TopExp.hxx>
-#include <TopTools_SequenceOfShape.hxx>
-#include <TopoDS.hxx>
-#include <TopoDS_Shell.hxx>
-#include <TopoDS_Solid.hxx>
-#include <TopoDS_Iterator.hxx>
-#include <TopLoc_Location.hxx>
-
-#include <BRepLib_MakeEdge.hxx>
-#include <BRepLib_MakeFace.hxx>
#include <BRepAdaptor_HCompCurve.hxx>
+#include <BRepBuilderAPI_Copy.hxx>
+#include <BRepBuilderAPI_Transform.hxx>
#include <BRepClass3d_SolidClassifier.hxx>
-
#include <BRepFill.hxx>
-#include <BRepFill_Sweep.hxx>
-#include <BRepFill_SectionPlacement.hxx>
-#include <BRepFill_Edge3DLaw.hxx>
#include <BRepFill_ACRLaw.hxx>
-#include <BRepFill_EdgeOnSurfLaw.hxx>
-#include <BRepFill_ShapeLaw.hxx>
#include <BRepFill_CompatibleWires.hxx>
-#include <BRepFill_NSections.hxx>
#include <BRepFill_DataMapOfShapeHArray2OfShape.hxx>
-#include <TColStd_HArray1OfReal.hxx>
-
-#include <GeomFill_TrihedronLaw.hxx>
+#include <BRepFill_Edge3DLaw.hxx>
+#include <BRepFill_EdgeOnSurfLaw.hxx>
+#include <BRepFill_LocationLaw.hxx>
+#include <BRepFill_NSections.hxx>
+#include <BRepFill_PipeShell.hxx>
+#include <BRepFill_Section.hxx>
+#include <BRepFill_SectionLaw.hxx>
+#include <BRepFill_SectionPlacement.hxx>
+#include <BRepFill_ShapeLaw.hxx>
+#include <BRepFill_Sweep.hxx>
+#include <BRepGProp.hxx>
+#include <BRepLib_MakeEdge.hxx>
+#include <BRepLib_MakeFace.hxx>
+#include <GeomAdaptor_HCurve.hxx>
+#include <GeomAdaptor_HSurface.hxx>
+#include <GeomFill_ConstantBiNormal.hxx>
#include <GeomFill_CorrectedFrenet.hxx>
-#include <GeomFill_Frenet.hxx>
+#include <GeomFill_CurveAndTrihedron.hxx>
#include <GeomFill_DiscreteTrihedron.hxx>
#include <GeomFill_Fixed.hxx>
-#include <GeomFill_ConstantBiNormal.hxx>
-#include <GeomFill_SectionLaw.hxx>
-#include <GeomFill_CurveAndTrihedron.hxx>
+#include <GeomFill_Frenet.hxx>
#include <GeomFill_GuideTrihedronAC.hxx>
#include <GeomFill_GuideTrihedronPlan.hxx>
#include <GeomFill_LocationGuide.hxx>
-
-//Specification Guide
-#include <GeomAdaptor_HCurve.hxx>
-
-#include <gp_Trsf.hxx>
+#include <GeomFill_SectionLaw.hxx>
+#include <GeomFill_TrihedronLaw.hxx>
+#include <gp_Ax2.hxx>
#include <gp_Dir.hxx>
+#include <gp_Trsf.hxx>
#include <gp_Vec.hxx>
-#include <Precision.hxx>
-
-#include <Standard_NotImplemented.hxx>
-#include <Standard_ConstructionError.hxx>
-#include <StdFail_NotDone.hxx>
-
-#include <BRepBuilderAPI_Copy.hxx>
-#include <BRepBuilderAPI_Transform.hxx>
-
#include <GProp_GProps.hxx>
-#include <BRepGProp.hxx>
-#include <GeomAdaptor_HSurface.hxx>
#include <IntCurveSurface_HInter.hxx>
#include <IntCurveSurface_IntersectionPoint.hxx>
-#include <TColgp_HArray1OfPnt2d.hxx>
+#include <Law_Function.hxx>
#include <Law_Interpol.hxx>
+#include <Precision.hxx>
+#include <Standard_ConstructionError.hxx>
+#include <Standard_DomainError.hxx>
+#include <Standard_NotImplemented.hxx>
+#include <Standard_Type.hxx>
+#include <StdFail_NotDone.hxx>
+#include <TColgp_HArray1OfPnt2d.hxx>
+#include <TColStd_HArray1OfReal.hxx>
+#include <TColStd_SequenceOfInteger.hxx>
+#include <TopExp.hxx>
+#include <TopLoc_Location.hxx>
+#include <TopoDS.hxx>
+#include <TopoDS_Iterator.hxx>
+#include <TopoDS_Shape.hxx>
+#include <TopoDS_Shell.hxx>
+#include <TopoDS_Solid.hxx>
+#include <TopoDS_Vertex.hxx>
+#include <TopoDS_Wire.hxx>
+#include <TopTools_SequenceOfShape.hxx>
+#include <stdio.h>
+IMPLEMENT_STANDARD_RTTIEXT(BRepFill_PipeShell,MMgt_TShared)
+
+//Specification Guide
#ifdef DRAW
#include <Draw.hxx>
#include <DrawTrSurf.hxx>
#include <TopTools_ListIteratorOfListOfShape.hxx>
#include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
+#include <TopTools_DataMapOfIntegerShape.hxx>
#include <TopoDS_Compound.hxx>
-static Standard_Boolean UpdateMap(const TopoDS_Shape& theKey,
- const TopoDS_Shape& theValue,
- TopTools_DataMapOfShapeListOfShape& theMap);
static Standard_Boolean BuildBoundaries(const BRepFill_Sweep& theSweep,
const Handle(BRepFill_SectionLaw)& theSection,
BRepFill_CompatibleWires CW(SSh);
CW.SetPercent(0.1);
CW.Perform();
- if (!CW.IsDone()) StdFail_NotDone::Raise("Uncompatible wires");
+ if (!CW.IsDone()) throw StdFail_NotDone("Uncompatible wires");
GeomFill_SequenceOfTrsf EmptyTrsfs;
Handle(BRepFill_NSections) SL = new (BRepFill_NSections) (CW.Shape(),EmptyTrsfs,SR,0.,1.);
Standard_Real US = p1/(p1+p2);
SetTolerance();
myMaxDegree = 11;
- myMaxSegments = 30;
+ myMaxSegments = 100;
// Attention to closed non-declared wire !
if (!mySpine.Closed()) {
BRepFill_CompatibleWires CW(Seq);
CW.SetPercent(0.1);
CW.Perform();
- if (!CW.IsDone()) StdFail_NotDone::Raise("Uncompatible wires");
+ if (!CW.IsDone()) throw StdFail_NotDone("Uncompatible wires");
TheGuide = TopoDS::Wire(CW.Shape().Value(2));
}
else if (GuideClose) {
const Standard_Boolean WithContact,
const Standard_Boolean WithCorrection)
{
- Delete(Profile); // No duplication
+ DeleteProfile(Profile); // No duplication
if (myIsAutomaticLaw)
{
mySeq.Clear();
//function : Delete
//purpose : Delete a section
//=======================================================================
- void BRepFill_PipeShell::Delete(const TopoDS_Shape& Profile)
+ void BRepFill_PipeShell::DeleteProfile(const TopoDS_Shape& Profile)
{
Standard_Boolean isVertex = (Profile.ShapeType() == TopAbs_VERTEX);
TopTools_MapOfShape Dummy;
BRepFill_DataMapOfShapeHArray2OfShape Dummy2;
BRepFill_DataMapOfShapeHArray2OfShape Dummy3;
- MkSw.Build(Dummy, Dummy2, Dummy3, myTransition, theContinuity);
+ MkSw.Build(Dummy, Dummy2, Dummy3, myTransition, theContinuity,
+ GeomFill_Location, myMaxDegree, myMaxSegments);
myStatus = myLocation->GetStatus();
Ok = (MkSw.IsDone() && (myStatus == GeomFill_PipeOk));
if (Ok) {
myShape = MkSw.Shape();
+ myErrorOnSurf = MkSw.ErrorOnSurface();
TopoDS_Shape aBottomWire = myFirst;
TopoDS_Shape aTopWire = myLast;
Standard_Boolean BRepFill_PipeShell::MakeSolid()
{
if (myShape.IsNull())
- StdFail_NotDone::Raise("PipeShell is not built");
+ throw StdFail_NotDone("PipeShell is not built");
Standard_Boolean B = myShape.Closed();
BRep_Builder BS;
return myShape;
}
+//=======================================================================
+//function : ErrorOnSurface
+//purpose :
+//=======================================================================
+
+Standard_Real BRepFill_PipeShell::ErrorOnSurface() const
+{
+ return myErrorOnSurf;
+}
+
//=======================================================================
//function : FirstShape
//purpose : Return the start section
void BRepFill_PipeShell::Generated(const TopoDS_Shape& theShape,
TopTools_ListOfShape& theList)
{
- // Standard_NotImplemented::Raise("Generated:Pas Fait");
+ // throw Standard_NotImplemented("Generated:Pas Fait");
+
theList.Clear();
if(myGenMap.IsBound(theShape)) {
theList = myGenMap.Find(theShape);
- }
+ }
}
//=======================================================================
//=======================================================================
void BRepFill_PipeShell::Prepare()
{
+ WSeq.Clear();
+ myEdgeNewEdges.Clear();
+
TopoDS_Wire theSect;
- if (!IsReady()) StdFail_NotDone::Raise("PipeShell");
+ if (!IsReady()) throw StdFail_NotDone("PipeShell");
if (!myLocation.IsNull() && !mySection.IsNull()) return; // It is ready
//Check set of section for right configuration of punctual sections
Standard_Integer i;
- TopoDS_Iterator iter;;
+ TopoDS_Iterator iter;
for (i = 2; i <= mySeq.Length()-1; i++)
{
Standard_Boolean wdeg = Standard_True;
wdeg = wdeg && (BRep_Tool::Degenerated(anEdge));
}
if (wdeg)
- Standard_Failure::Raise("Wrong usage of punctual sections");
+ throw Standard_Failure("Wrong usage of punctual sections");
}
if (mySeq.Length() <= 2)
{
wdeg = wdeg && (BRep_Tool::Degenerated(anEdge));
}
if (wdeg)
- Standard_Failure::Raise("Wrong usage of punctual sections");
+ throw Standard_Failure("Wrong usage of punctual sections");
}
// Construction of the law of location
}
default :
{ // Not planned!
- Standard_ConstructionError::Raise("PipeShell");
+ throw Standard_ConstructionError("PipeShell");
}
}
}
else
mySection = new BRepFill_ShapeLaw(aLocalShape);
// mySection = new (BRepFill_ShapeLaw) (TopoDS::Wire(theSect));
- }
- else
+
+ WSeq.Append(theSect);
+ //Simple case of single section
+ myIndOfSec.Append(1);
+ TopoDS_Iterator itw(theSect);
+ for (; itw.More(); itw.Next())
{
- TColStd_SequenceOfReal Param;
- TopTools_SequenceOfShape WSeq;
- GeomFill_SequenceOfTrsf Transformations;
- //WSeq.Clear();
- //Param.Clear();
- Standard_Integer NbL = myLocation->NbLaw();
- gp_Trsf aTrsf;
- Standard_Real V1, V2, param;
- myLocation->CurvilinearBounds(NbL, V1, V2);
- V1 = 0.;
- Standard_Integer ideb = 0, ifin = 0;
-// for (Standard_Integer iseq=1;iseq<=mySeq.Length();iseq++) {
- Standard_Integer iseq;
- for (iseq=1;iseq<=mySeq.Length();iseq++) {
- Place(mySeq(iseq), theSect, aTrsf, param);
- Param.Append(param);
- WSeq.Append(theSect);
-// WSeq.Append(TopoDS::Wire(theSect));
- Transformations.Append(aTrsf);
- if (param==V1) ideb = iseq;
- if (param==V2) ifin = iseq;
+ const TopoDS_Shape& anEdge = itw.Value();
+ TopTools_ListOfShape Elist;
+ Elist.Append(anEdge);
+ myEdgeNewEdges.Bind(anEdge, Elist);
+ }
+ ///////////////////////////////
+ }
+ else
+ {
+ TColStd_SequenceOfReal Param;
+ TColStd_SequenceOfInteger IndSec;
+ GeomFill_SequenceOfTrsf Transformations;
+ Standard_Integer NbL = myLocation->NbLaw();
+ gp_Trsf aTrsf;
+ Standard_Real V1, V2, param;
+ myLocation->CurvilinearBounds(NbL, V1, V2);
+ V1 = 0.;
+ Standard_Integer ideb = 0, ifin = 0;
+ Standard_Integer iseq;
+ for (iseq = 1; iseq <= mySeq.Length(); iseq++) {
+ IndSec.Append(iseq);
+ Place(mySeq(iseq), theSect, aTrsf, param);
+ Param.Append(param);
+ WSeq.Append(theSect);
+ Transformations.Append(aTrsf);
+ if (param==V1) ideb = iseq;
+ if (param==V2) ifin = iseq;
+ }
+
+
+ // looping sections ?
+ if (myLocation->IsClosed()) {
+ if (ideb>0) {
+ // place the initial section at the final position
+ Param.Append(V2);
+ WSeq.Append(WSeq(ideb));
}
-
-
- // looping sections ?
- if (myLocation->IsClosed()) {
- if (ideb>0) {
- // place the initial section at the final position
- Param.Append(V2);
- WSeq.Append(WSeq(ideb));
+ else if (ifin>0) {
+ // place the final section at the initial position
+ Param.Append(V1);
+ WSeq.Append(WSeq(ifin));
+ }
+ else {
+ // it is necessary to find a medium section to impose by V1 and by V2
+ Standard_Real pmin = RealLast(), pmax = RealFirst();
+ TopoDS_Wire Wmin, Wmax;
+ for (iseq = 1; iseq <= WSeq.Length(); iseq++) {
+ if (Param.Value(iseq)<pmin) {
+ pmin = Param.Value(iseq);
+ Wmin = TopoDS::Wire(WSeq.Value(iseq));
+ }
+ if (Param.Value(iseq)>pmax) {
+ pmax = Param.Value(iseq);
+ Wmax = TopoDS::Wire(WSeq.Value(iseq));
+ }
}
- else if (ifin>0) {
- // place the final section at the initial position
+ // medium section between Wmin and Wmax
+ TopoDS_Wire Wres;
+ Standard_Real dmin = Abs(pmin-V1);
+ Standard_Real dmax = Abs(pmax-V2);
+ if (ComputeSection(Wmin,Wmax,dmin,dmax,Wres)) {
+ // impose section Wres at the beginning and the end
Param.Append(V1);
- WSeq.Append(WSeq(ifin));
+ WSeq.Append(Wres);
+ IndSec.Append(WSeq.Length());
+ Param.Append(V2);
+ WSeq.Append(Wres);
+ IndSec.Append(WSeq.Length());
}
- else {
- // it is necessary to find a medium section to impose by V1 and by V2
- Standard_Real pmin = Param.Value(1), pmax = Param.Value(1);
- TopoDS_Wire Wmin = TopoDS::Wire(WSeq.Value(1)), Wmax;
- for (iseq=2;iseq<=WSeq.Length();iseq++) {
- if (Param.Value(iseq)<pmin) {
- pmin = Param.Value(iseq);
- Wmin = TopoDS::Wire(WSeq.Value(iseq));
- }
- if (Param.Value(iseq)>pmax) {
- pmax = Param.Value(iseq);
- Wmax = TopoDS::Wire(WSeq.Value(iseq));
- }
- }
- // medium section between Wmin and Wmax
- TopoDS_Wire Wres;
- Standard_Real dmin = Abs(pmin-V1);
- Standard_Real dmax = Abs(pmax-V2);
- if (ComputeSection(Wmin,Wmax,dmin,dmax,Wres)) {
- // impose section Wres at the beginning and the end
- Param.Append(V1);
- WSeq.Append(Wres);
- Param.Append(V2);
- WSeq.Append(Wres);
-
+ }
+ }
+
+ // parse sections by increasing parameter
+ Standard_Boolean play_again = Standard_True;
+ while (play_again) {
+ play_again = Standard_False;
+ for (iseq=1;iseq<=WSeq.Length();iseq++) {
+ for (Standard_Integer jseq=iseq+1;jseq<=WSeq.Length();jseq++) {
+ if (Param.Value(iseq) > Param.Value(jseq)) {
+ Param.Exchange(iseq,jseq);
+ WSeq.Exchange(iseq,jseq);
+ IndSec.Exchange(iseq,jseq);
+ play_again = Standard_True;
}
-
}
}
-
- // parse sections by increasing parameter
- Standard_Boolean play_again = Standard_True;
- while (play_again) {
- play_again = Standard_False;
- for (iseq=1;iseq<=WSeq.Length();iseq++) {
- for (Standard_Integer jseq=iseq+1;jseq<=WSeq.Length();jseq++) {
- if (Param.Value(iseq)>Param.Value(jseq)) {
- Param.Exchange(iseq,jseq);
- WSeq.Exchange(iseq,jseq);
- play_again = Standard_True;
- }
- }
- }
- }
-
+ }
+ //Fill the array of real indices of sections
+ for (Standard_Integer ii = 1; ii <= mySeq.Length(); ii++)
+ for (Standard_Integer jj = 1; jj <= IndSec.Length(); jj++)
+ if (IndSec(jj) == ii)
+ {
+ myIndOfSec.Append(jj);
+ break;
+ }
+
#ifdef DRAW
- if ( Affich) {
- char* name = new char[100];
- Standard_Integer NBSECT = 0;
- for (Standard_Integer i=1;i<=WSeq.Length();i++) {
- NBSECT++;
- sprintf(name,"WSeq_%d",NBSECT);
- DBRep::Set(name,TopoDS::Wire(WSeq.Value(i)));
+ if ( Affich) {
+ char* name = new char[100];
+ Standard_Integer NBSECT = 0;
+ for (Standard_Integer i=1;i<=WSeq.Length();i++) {
+ NBSECT++;
+ sprintf(name,"WSeq_%d",NBSECT);
+ DBRep::Set(name,TopoDS::Wire(WSeq.Value(i)));
+ }
}
- }
#endif
-
-
-
- // Calculate work sections
- TopTools_SequenceOfShape WorkingSections;
- WorkingSections.Clear();
- TopTools_DataMapOfShapeListOfShape WorkingMap;
- WorkingMap.Clear();
- BRepFill_CompatibleWires Georges(WSeq);
- Georges.SetPercent(0.1);
- Georges.Perform(Standard_False);
- if (Georges.IsDone()) {
- WorkingSections = Georges.Shape();
- WorkingMap = Georges.Generated();
- }
- else {
- Standard_ConstructionError::Raise("PipeShell : uncompatible wires");
+
+
+
+ // Calculate work sections
+ TopTools_SequenceOfShape WorkingSections;
+ WorkingSections.Clear();
+ TopTools_DataMapOfShapeListOfShape WorkingMap;
+ BRepFill_CompatibleWires Georges(WSeq);
+ Georges.SetPercent(0.1);
+ Georges.Perform(Standard_False);
+ if (Georges.IsDone()) {
+ WorkingSections = Georges.Shape();
+ WorkingMap = Georges.Generated();
+ //For each sub-edge of each section
+ //we save its splits
+ for (Standard_Integer ii = 1; ii <= WSeq.Length(); ii++)
+ {
+ TopExp_Explorer Explo(WSeq(ii), TopAbs_EDGE);
+ for (; Explo.More(); Explo.Next())
+ {
+ const TopoDS_Edge& anEdge = TopoDS::Edge(Explo.Current());
+ TopTools_ListOfShape aNewEdges = Georges.GeneratedShapes(anEdge);
+ myEdgeNewEdges.Bind(anEdge, aNewEdges);
+ }
}
- mySection = new (BRepFill_NSections) (WorkingSections,Transformations,Param,V1,V2);
-
- }// else
+ }
+ else {
+ throw Standard_ConstructionError("PipeShell : uncompatible wires");
+ }
+ mySection = new (BRepFill_NSections) (WorkingSections,Transformations,Param,V1,V2);
+
+ }// else
// modify the law of location if contact
if ( (myTrihedron == GeomFill_IsGuidePlanWithContact)
}
myStatus = myLocation->GetStatus();
+ if (!mySection->IsDone())
+ myStatus = GeomFill_PipeNotOk;
}
//=======================================================================
//=======================================================================
//function : BuildHistory
-//purpose : Builds history for edges of spine,
-// for built bottom shape of sweep,
-// for boundary vertices of bottom shape of sweep,
-// for boundary profiles
+//purpose : Builds history for edges and vertices
+// of sections
//=======================================================================
void BRepFill_PipeShell::BuildHistory(const BRepFill_Sweep& theSweep)
{
- Handle(TopTools_HArray2OfShape) aFaces = theSweep.SubShape();
- Handle(TopTools_HArray2OfShape) aVEdges = theSweep.Sections();
- Handle(TopTools_HArray2OfShape) aUEdges = theSweep.InterFaces();
- Standard_Integer i = 0, j = 0;
- Standard_Boolean bPrevModified = Standard_False;
-
- for(i = 1; i <= mySection->NbLaw(); i++) {
- if((!aVEdges->Value(i, 1).IsNull()) && (aVEdges->Value(i, 1).ShapeType() == TopAbs_FACE)) {
- bPrevModified = Standard_True;
- break;
- }
- }
-
- for(j = myLocation->NbLaw(); j >= 1; j--) {
- TopTools_ListOfShape aListOfFace;
-
- if(bPrevModified) {
- for(i = 1; i <= mySection->NbLaw(); i++) {
- Standard_Integer lessindex = j + 1;
- lessindex = (lessindex > myLocation->NbLaw()) ? 1 : lessindex;
-
- if((!aVEdges->Value(i, lessindex).IsNull()) && (aVEdges->Value(i, lessindex).ShapeType() == TopAbs_FACE)) {
- aListOfFace.Append(aVEdges->Value(i, lessindex));
- const TopoDS_Shape& aBottomEdge = aVEdges->Value(i, 1);
-
- if((!aBottomEdge.IsNull()) && (aBottomEdge.ShapeType() == TopAbs_EDGE)) {
- UpdateMap(aBottomEdge, aVEdges->Value(i, lessindex), myGenMap);
- }
- }
- }
+ //Filling of <myGenMap>
+ const Handle(TopTools_HArray2OfShape)& anUEdges = theSweep.InterFaces();
+ BRep_Builder BB;
+
+ TopTools_DataMapOfIntegerShape IndWireMap;
+
+ Standard_Integer indw, inde;
+ TopoDS_Iterator itw;
+ for (indw = 1; indw <= mySeq.Length(); indw++)
+ {
+ const TopoDS_Wire& aSection = mySeq(indw).Wire();
+ Standard_Boolean IsPunctual = mySeq(indw).IsPunctual();
+ if (IsPunctual)
+ {
+ //for punctual sections (first or last)
+ //we take all the wires generated along the path
+ TopExp_Explorer Explo(aSection, TopAbs_VERTEX);
+ const TopoDS_Shape& VerSection = Explo.Current();
+ TopTools_ListOfShape Elist;
+ for (Standard_Integer i = 1; i <= anUEdges->UpperRow(); i++)
+ for (Standard_Integer j = 1; j <= anUEdges->UpperCol(); j++)
+ Elist.Append(anUEdges->Value(i,j));
+ myGenMap.Bind(VerSection, Elist);
+ continue;
}
- bPrevModified = Standard_False;
-
- for(i = 1; i <= mySection->NbLaw(); i++) {
- if((!aVEdges->Value(i, j).IsNull()) && (aVEdges->Value(i, j).ShapeType() == TopAbs_FACE)) {
- aListOfFace.Append(aVEdges->Value(i, j));
- bPrevModified = Standard_True;
-
- const TopoDS_Shape& aBottomEdge = aVEdges->Value(i, 1);
-
- if((!aBottomEdge.IsNull()) && (aBottomEdge.ShapeType() == TopAbs_EDGE)) {
- UpdateMap(aBottomEdge, aVEdges->Value(i, j), myGenMap);
- }
- }
-
- if(aFaces->Value(i, j).ShapeType() == TopAbs_FACE) {
- aListOfFace.Append(aFaces->Value(i, j));
- const TopoDS_Shape& aBottomEdge = aVEdges->Value(i, 1);
-
- if((!aBottomEdge.IsNull()) && (aBottomEdge.ShapeType() == TopAbs_EDGE)) {
- UpdateMap(aBottomEdge, aFaces->Value(i, j), myGenMap);
- }
+ //Take the real index of section on the path
+ Standard_Integer IndOfW = myIndOfSec(indw);
+ const TopoDS_Wire& theWire = TopoDS::Wire(WSeq(IndOfW));
+ BRepTools_WireExplorer wexp_sec(aSection);
+ for (inde = 1; wexp_sec.More(); wexp_sec.Next())
+ {
+ const TopoDS_Edge& anEdge = TopoDS::Edge(wexp_sec.Current());
+ if (BRep_Tool::Degenerated(anEdge))
+ continue;
+
+ TopoDS_Shell aShell;
+ BB.MakeShell(aShell);
+ TopoDS_Vertex aVertex [2];
+ TopExp::Vertices(anEdge, aVertex[0], aVertex[1]);
+ Standard_Integer SignOfAnEdge =
+ (anEdge.Orientation() == TopAbs_FORWARD)? 1 : -1;
+
+ //For each non-degenerated inde-th edge of <aSection>
+ //we find inde-th edge in <theWire>
+ TopoDS_Edge theEdge;
+ BRepTools_WireExplorer wexp(theWire);
+ for (Standard_Integer i = 1; wexp.More(); wexp.Next())
+ {
+ theEdge = TopoDS::Edge(wexp.Current());
+ if (BRep_Tool::Degenerated(anEdge))
+ continue;
+ if (i == inde)
+ break;
+ i++;
}
- }
- if(!myGenMap.IsBound(myLocation->Edge(j)))
- myGenMap.Bind(myLocation->Edge(j), aListOfFace);
- else
- myGenMap.ChangeFind(myLocation->Edge(j)).Append(aListOfFace);
-
- // build history for free booundaries.begin
- if(!mySection->IsUClosed()) {
- TopoDS_Compound aFaceComp;
- BRep_Builder aB;
- aB.MakeCompound(aFaceComp);
- TopTools_ListIteratorOfListOfShape anIt(aListOfFace);
-
- for(; anIt.More(); anIt.Next()) {
- aB.Add(aFaceComp, anIt.Value());
+ //Take the list of splits for <theEdge>
+ const TopTools_ListOfShape& NewEdges = myEdgeNewEdges(theEdge);
+ Standard_Integer SignOfANewEdge = 0, SignOfIndex = 0;
+ TopTools_ListIteratorOfListOfShape iter(NewEdges);
+ for (; iter.More(); iter.Next())
+ {
+ const TopoDS_Edge& aNewEdge = TopoDS::Edge(iter.Value());
+ SignOfANewEdge = (aNewEdge.Orientation() == TopAbs_FORWARD)? 1 : -1;
+ Standard_Integer anIndE = mySection->IndexOfEdge(aNewEdge);
+ SignOfIndex = (anIndE > 0)? 1 : -1;
+ anIndE = Abs(anIndE);
+ //For an edge generated shape is a "tape" -
+ //a shell usually containing this edge and
+ //passing from beginning of path to its end
+ TopoDS_Shape aTape = theSweep.Tape(anIndE);
+ TopoDS_Iterator itsh(aTape);
+ for (; itsh.More(); itsh.Next())
+ BB.Add(aShell, itsh.Value());
}
- TopTools_IndexedDataMapOfShapeListOfShape aMapEF;
- TopExp::MapShapesAndAncestors(aFaceComp, TopAbs_EDGE, TopAbs_FACE, aMapEF);
- Standard_Integer eit = 0;
-
- for(eit = aUEdges->LowerRow(); eit <= aUEdges->UpperRow(); eit++) {
- const TopoDS_Shape& aShape = aUEdges->Value(eit, j);
-
- if(aMapEF.Contains(aShape)) {
- const TopTools_ListOfShape& aList = aMapEF.FindFromKey(aShape);
- if(aList.Extent() < 2) {
- UpdateMap(myLocation->Edge(j), aShape, myGenMap);
-
- TopoDS_Shape aGenVertex;
- TopTools_IndexedDataMapOfShapeListOfShape aMapVE;
-
- for(i = 1; i <= mySection->NbLaw(); i++) {
- const TopoDS_Shape& aBottomEdge = aVEdges->Value(i, aVEdges->LowerCol());
-
- if((!aBottomEdge.IsNull()) && (aBottomEdge.ShapeType() == TopAbs_EDGE)) {
- TopExp::MapShapesAndAncestors(aBottomEdge, TopAbs_VERTEX, TopAbs_EDGE, aMapVE);
- }
- }
- const TopoDS_Shape& aFreeEdge = aUEdges->Value(eit, aUEdges->LowerCol());
- TopExp::MapShapesAndAncestors(aFreeEdge, TopAbs_VERTEX, TopAbs_EDGE, aMapVE);
- TopExp_Explorer anExpV(aFreeEdge, TopAbs_VERTEX);
-
- for(; anExpV.More(); anExpV.Next()) {
- if(aMapVE.Contains(anExpV.Current())) {
- const TopTools_ListOfShape& aListOfV = aMapVE.FindFromKey(anExpV.Current());
-
- if(aListOfV.Extent() >= 2) {
- aGenVertex = anExpV.Current();
- }
- }
- }
-
- if(!aGenVertex.IsNull()) {
- UpdateMap(aGenVertex, aShape, myGenMap);
- }
- }
- }
+ //Processing of vertices of <anEdge>
+ //We should choose right index in <anUEdges>
+ //for each vertex of edge
+ Standard_Integer ToReverse = SignOfAnEdge * SignOfANewEdge * SignOfIndex;
+ Standard_Integer UIndex [2];
+ UIndex[0] = Abs(mySection->IndexOfEdge(NewEdges.First()));
+ UIndex[1] = Abs(mySection->IndexOfEdge(NewEdges.Last())) + ToReverse;
+ if (ToReverse == -1)
+ {
+ UIndex[0]++;
+ UIndex[1]++;
}
- // end for(eit = aUEdges->LowerRow...
- }
- // build history for free booundaries.end
- }
-
- // build history for boundary section wires. begin
-
- if(!mySeq.IsEmpty()) {
- Standard_Integer iseq;
- TopoDS_Wire aSect;
- gp_Trsf aTrsf;
- Standard_Real param = 0., aparmin = RealLast(), aparmax = -RealLast();
- Standard_Integer ideb = 1, ifin = mySeq.Length();
-
- for (iseq = 1;iseq <= mySeq.Length(); iseq++) {
- Place(mySeq(iseq), aSect, aTrsf, param);
-
- if(param < aparmin) {
- ideb = iseq;
- aparmin = param;
+ if (mySection->IsUClosed())
+ {
+ if (UIndex[0] > mySection->NbLaw())
+ UIndex[0] = 1;
+ if (UIndex[1] > mySection->NbLaw())
+ UIndex[1] = 1;
}
+ //if (SignOfAnEdge * SignOfANewEdge == -1)
+ if (SignOfAnEdge == -1 ||
+ SignOfANewEdge == -1)
+ { Standard_Integer Tmp = UIndex[0]; UIndex[0] = UIndex[1]; UIndex[1] = Tmp; }
+
+ TopTools_IndexedDataMapOfShapeListOfShape VEmap;
+ TopExp::MapShapesAndAncestors(aShell, TopAbs_VERTEX, TopAbs_EDGE, VEmap);
+ for (Standard_Integer kk = 0; kk < 2; kk++)
+ {
+ if (myGenMap.IsBound(aVertex[kk]))
+ continue;
+ if (IndWireMap.IsBound(UIndex[kk]))
+ {
+ TopTools_ListOfShape Wlist;
+ Wlist.Append(IndWireMap(UIndex[kk]));
+ myGenMap.Bind(aVertex[kk], Wlist);
+ continue;
+ }
+
+ //Collect u-edges
+ TopTools_SequenceOfShape SeqEdges;
+ Standard_Integer jj;
+ for (jj = 1; jj <= anUEdges->UpperCol(); jj++)
+ SeqEdges.Append(anUEdges->Value(UIndex[kk], jj));
+
+ //Assemble the wire ("rail" along the path)
+ //checking for possible holes
+ //(they appear with option "Round Corner")
+ //and filling them
+ //Missed edges are taken from <aShell>
+ TopoDS_Wire aWire;
+ BB.MakeWire(aWire);
+ const TopoDS_Edge& FirstEdge = TopoDS::Edge(SeqEdges(1));
+ if (FirstEdge.IsNull())
+ continue;
+ BB.Add(aWire, FirstEdge);
+ TopoDS_Vertex FirstVertex, CurVertex;
+ TopExp::Vertices(FirstEdge, FirstVertex, CurVertex);
+ TopoDS_Edge CurEdge;
+ for (jj = 2; jj <= SeqEdges.Length(); jj++)
+ {
+ CurEdge = TopoDS::Edge(SeqEdges(jj));
+ TopoDS_Vertex Vfirst, Vlast;
+ TopExp::Vertices(CurEdge, Vfirst, Vlast);
+ if (CurVertex.IsSame(Vfirst))
+ CurVertex = Vlast;
+ else //a hole
+ {
+ const TopTools_ListOfShape& Elist = VEmap.FindFromKey(Vfirst);
+ TopTools_ListIteratorOfListOfShape itl(Elist);
+ for (; itl.More(); itl.Next())
+ {
+ const TopoDS_Edge& Candidate = TopoDS::Edge(itl.Value());
+ if (Candidate.IsSame(CurEdge))
+ continue;
+ TopoDS_Vertex V1, V2;
+ TopExp::Vertices(Candidate, V1, V2);
+ if (V1.IsSame(CurVertex) || V2.IsSame(CurVertex))
+ {
+ BB.Add(aWire, Candidate);
+ break;
+ }
+ }
+ }
+ CurVertex = Vlast;
+ BB.Add(aWire, CurEdge);
+ } //for (jj = 2; jj <= SeqEdges.Length(); jj++)
+ //case of closed wire
+ if (mySection->IsVClosed() &&
+ !CurVertex.IsSame(FirstVertex))
+ {
+ const TopTools_ListOfShape& Elist = VEmap.FindFromKey(CurVertex);
+ TopTools_ListIteratorOfListOfShape itl(Elist);
+ for (; itl.More(); itl.Next())
+ {
+ const TopoDS_Edge& Candidate = TopoDS::Edge(itl.Value());
+ if (Candidate.IsSame(CurEdge))
+ continue;
+ TopoDS_Vertex V1, V2;
+ TopExp::Vertices(Candidate, V1, V2);
+ if (V1.IsSame(FirstVertex) || V2.IsSame(FirstVertex))
+ {
+ BB.Add(aWire, Candidate);
+ break;
+ }
+ }
+ }
+ TopTools_ListOfShape Wlist;
+ Wlist.Append(aWire);
+ myGenMap.Bind(aVertex[kk], Wlist);
+ //Save already built wire with its index
+ IndWireMap.Bind(UIndex[kk], aWire);
+ } //for (Standard_Integer kk = 0; kk < 2; kk++)
+ ////////////////////////////////////
+
+ TopTools_ListOfShape ListShell;
+ ListShell.Append(aShell);
+ myGenMap.Bind(anEdge, ListShell);
+ ////////////////////////
- if(param > aparmax) {
- ifin = iseq;
- aparmax = param;
- }
+ inde++;
}
-
- UpdateMap(mySeq(ideb).Wire(), myFirst, myGenMap);
- UpdateMap(mySeq(ifin).Wire(), myLast, myGenMap);
}
- // build history for boundary section wires. end
}
-// ---------------------------------------------------------------------------------
-// static function: UpdateMap
-// purpose:
-// ---------------------------------------------------------------------------------
-Standard_Boolean UpdateMap(const TopoDS_Shape& theKey,
- const TopoDS_Shape& theValue,
- TopTools_DataMapOfShapeListOfShape& theMap) {
-
- if(!theMap.IsBound(theKey)) {
- TopTools_ListOfShape thelist;
- theMap.Bind(theKey, thelist);
- }
- TopTools_ListOfShape& aList = theMap.ChangeFind(theKey);
- TopTools_ListIteratorOfListOfShape anIt(aList);
- Standard_Boolean found = Standard_False;
-
- for(; anIt.More(); anIt.Next()) {
- if(theValue.IsSame(anIt.Value())) {
- found = Standard_True;
- break;
- }
- }
-
- if(!found)
- aList.Append(theValue);
- return !found;
-}
// ---------------------------------------------------------------------------------
// static function: BuildBoundaries
projects / occt.git / blobdiff