1 // Copyright (C) 2007-2014 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
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, or (at your option) any later version.
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.
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
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 #include "GEOMImpl_PolylineDriver.hxx"
25 #include "GEOMImpl_ICurveParametric.hxx"
26 #include "GEOMImpl_ICurvesOperations.hxx"
27 #include "GEOMImpl_IPolyline.hxx"
28 #include "GEOMImpl_IPolyline2D.hxx"
29 #include "GEOMImpl_Types.hxx"
30 #include "GEOM_Function.hxx"
31 #include <GEOMUtils.hxx>
32 #include <Sketcher_Utils.hxx>
34 #include <BRepBuilderAPI_MakePolygon.hxx>
35 #include <BRepBuilderAPI_MakeVertex.hxx>
36 #include <BRep_Builder.hxx>
37 #include <BRep_Tool.hxx>
38 #include <Precision.hxx>
39 #include <TColgp_Array1OfPnt.hxx>
43 #include <TopoDS_Shape.hxx>
44 #include <TopoDS_Vertex.hxx>
45 #include <TopoDS_Wire.hxx>
49 //=======================================================================
52 //=======================================================================
53 const Standard_GUID& GEOMImpl_PolylineDriver::GetID()
55 static Standard_GUID aPolylineDriver("FF1BBB31-5D14-4df2-980B-3A668264EA16");
56 return aPolylineDriver;
60 //=======================================================================
61 //function : GEOMImpl_PolylineDriver
63 //=======================================================================
64 GEOMImpl_PolylineDriver::GEOMImpl_PolylineDriver()
68 //=======================================================================
69 //function : MakePolyline2D
71 //=======================================================================
72 Standard_Integer GEOMImpl_PolylineDriver::MakePolyline2D
73 (TFunction_Logbook& log) const
75 if (Label().IsNull()) {
79 Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
80 GEOMImpl_IPolyline2D aCI(aFunction);
81 Standard_Integer aType = aFunction->GetType();
85 Handle(TColStd_HArray1OfExtendedString) aNames = aCI.GetNames();
86 Handle(TColStd_HArray1OfByte) aTypes = aCI.GetTypes();
87 Handle(TColStd_HArray1OfByte) aClosedFlags = aCI.GetClosedFlags();
88 std::list <std::list <double> > aCoords;
91 aCI.GetCoords(aCoords);
93 // Check the data validity
94 if (aNames.IsNull()) {
98 Standard_Integer aNbSections = aNames->Length();
100 if (aTypes.IsNull() || aNbSections != aTypes->Length()) {
104 if (aClosedFlags.IsNull() || aNbSections != aClosedFlags->Length()) {
108 if (aNbSections != aCoords.size()) {
112 if (aType == POLYLINE2D_PLN_COORDS) {
113 Handle(TColStd_HArray1OfReal) aPlaneCoords = aCI.GetWorkingPlaneDbls();
115 if (aPlaneCoords.IsNull()) {
119 if (aPlaneCoords->Length() != 9) {
123 Standard_Integer i = aPlaneCoords->Lower();
124 gp_Pnt aOrigin(aPlaneCoords->Value(i), aPlaneCoords->Value(i + 1),
125 aPlaneCoords->Value(i + 2));
126 gp_Dir aDirZ(aPlaneCoords->Value(i + 3), aPlaneCoords->Value(i + 4),
127 aPlaneCoords->Value(i + 5));
128 gp_Dir aDirX(aPlaneCoords->Value(i + 6), aPlaneCoords->Value(i + 7),
129 aPlaneCoords->Value(i + 8));
130 aWPlane = gp_Ax3(aOrigin, aDirZ, aDirX);
131 } else if (aType == POLYLINE2D_PLN_OBJECT) {
132 Handle(GEOM_Function) aRefFace = aCI.GetWorkingPlane();
133 TopoDS_Shape aShape = aRefFace->GetValue();
135 aWPlane = GEOMUtils::GetPosition(aShape);
140 // Construct a shape.
141 Standard_Integer iN = aNames->Lower();
142 Standard_Integer iT = aTypes->Lower();
143 Standard_Integer iC = aClosedFlags->Lower();
144 std::list <std::list <double> >::const_iterator anIter = aCoords.begin();
145 BRep_Builder aBuilder;
146 Standard_Boolean isEmpty = Standard_True;
148 if (aNbSections > 1) {
149 aBuilder.MakeCompound(TopoDS::Compound(aShape));
152 for (; anIter != aCoords.end(); ++anIter, ++iN, ++iT, ++iC) {
153 Standard_Integer aType = aTypes->Value(iT);
154 TopoDS_Shape aSection;
156 if (aType == GEOMImpl_ICurvesOperations::Polyline) {
157 aSection = Sketcher_Utils::MakePolyline
158 (*anIter, aClosedFlags->Value(iC), aWPlane);
159 } else if (aType == GEOMImpl_ICurvesOperations::Interpolation) {
160 aSection = Sketcher_Utils::MakeInterpolation
161 (*anIter, aClosedFlags->Value(iC), aWPlane);
164 if (aNbSections > 1) {
165 // There are multiple sections.
166 if (aSection.IsNull() == Standard_False) {
167 aBuilder.Add(aShape, aSection);
168 isEmpty = Standard_False;
171 // There is only one section.
172 isEmpty = aSection.IsNull();
181 aFunction->SetValue(aShape);
182 log.SetTouched(Label());
187 //=======================================================================
190 //=======================================================================
191 Standard_Integer GEOMImpl_PolylineDriver::Execute(TFunction_Logbook& log) const
193 if (Label().IsNull()) return 0;
194 Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
195 Standard_Integer aType = aFunction->GetType();
197 if (aType == POLYLINE2D_PLN_COORDS || aType == POLYLINE2D_PLN_OBJECT) {
198 return MakePolyline2D(log);
201 GEOMImpl_IPolyline aCI (aFunction);
205 if (aType == POLYLINE_POINTS) {
207 bool useCoords = aCI.GetConstructorType() == COORD_CONSTRUCTOR;
208 TColgp_Array1OfPnt points(1, (useCoords ? aCI.GetLength() : 1) );
210 Handle(TColStd_HArray1OfReal) aCoordsArray = aCI.GetCoordinates();
211 int anArrayLength = aCoordsArray->Length();
212 for (int i = 0, j = 1; i <= (anArrayLength-3); i += 3) {
213 gp_Pnt aPnt = gp_Pnt(aCoordsArray->Value(i+1), aCoordsArray->Value(i+2), aCoordsArray->Value(i+3));
214 points.SetValue(j,aPnt);
219 int aLen = aCI.GetLength();
221 BRepBuilderAPI_MakePolygon aMakePoly;
222 for (; ind <= aLen; ind++)
225 aMakePoly.Add(BRepBuilderAPI_MakeVertex(points.Value(ind)));
227 Handle(GEOM_Function) aRefPoint = aCI.GetPoint(ind);
228 TopoDS_Shape aShapePnt = aRefPoint->GetValue();
229 if (aShapePnt.ShapeType() != TopAbs_VERTEX) {
230 Standard_TypeMismatch::Raise
231 ("Polyline creation aborted : arguments are not a vertexes");
234 if (aShapePnt.ShapeType() == TopAbs_VERTEX) {
235 aMakePoly.Add(TopoDS::Vertex(aShapePnt));
236 //if (!aMakePoly.Added()) return 0;
240 // Compare first and last point coordinates and close polyline if it's the same.
244 aV1 = BRepBuilderAPI_MakeVertex(points.Value(1));
246 Handle(GEOM_Function) aFPoint = aCI.GetPoint(1);
247 TopoDS_Shape aFirstPnt = aFPoint->GetValue();
248 aV1 = TopoDS::Vertex(aFirstPnt);
253 aV2 = BRepBuilderAPI_MakeVertex(points.Value(aLen));
255 Handle(GEOM_Function) aLPoint = aCI.GetPoint(aLen);
256 TopoDS_Shape aLastPnt = aLPoint->GetValue();
257 aV2 = TopoDS::Vertex(aLastPnt);
260 if ( (!aV1.IsNull() && !aV2.IsNull() && aV1.IsSame(aV2)) ||
265 if (aMakePoly.IsDone()) {
266 aShape = aMakePoly.Wire();
272 if (aShape.IsNull()) return 0;
274 aFunction->SetValue(aShape);
276 log.SetTouched(Label());
281 //================================================================================
283 * \brief Returns a name of creation operation and names and values of creation parameters
285 //================================================================================
287 bool GEOMImpl_PolylineDriver::
288 GetCreationInformation(std::string& theOperationName,
289 std::vector<GEOM_Param>& theParams)
291 if (Label().IsNull()) return 0;
292 Handle(GEOM_Function) function = GEOM_Function::GetFunction(Label());
294 GEOMImpl_IPolyline aCI( function );
295 GEOMImpl_ICurveParametric aIP( function );
296 Standard_Integer aType = function->GetType();
299 case POLYLINE_POINTS:
300 theOperationName = "CURVE";
301 AddParam( theParams, "Type", "Polyline");
304 AddParam( theParams, "X(t) equation", aIP.GetExprX() );
305 AddParam( theParams, "Y(t) equation", aIP.GetExprY() );
306 AddParam( theParams, "Z(t) equation", aIP.GetExprZ() );
307 AddParam( theParams, "Min t", aIP.GetParamMin() );
308 AddParam( theParams, "Max t", aIP.GetParamMax() );
309 if ( aIP.GetParamNbStep() )
310 AddParam( theParams, "Number of steps", aIP.GetParamNbStep() );
312 AddParam( theParams, "t step", aIP.GetParamStep() );
316 GEOM_Param& pntParam = AddParam( theParams, "Points");
317 if ( aCI.GetConstructorType() == COORD_CONSTRUCTOR )
319 Handle(TColStd_HArray1OfReal) coords = aCI.GetCoordinates();
320 if ( coords->Length() > 3 )
321 pntParam << ( coords->Length() ) / 3 << " points: ";
322 for ( int i = coords->Lower(), nb = coords->Upper(); i <= nb; )
323 pntParam << "( " << coords->Value( i++ )
324 << ", " << coords->Value( i++ )
325 << ", " << coords->Value( i++ ) << " ) ";
329 if ( aCI.GetLength() > 1 )
330 pntParam << aCI.GetLength() << " points: ";
331 for ( int i = 1, nb = aCI.GetLength(); i <= nb; ++i )
332 pntParam << aCI.GetPoint( i ) << " ";
334 AddParam( theParams, "Is closed", aCI.GetIsClosed() );
337 case POLYLINE2D_PLN_COORDS:
338 case POLYLINE2D_PLN_OBJECT:
340 theOperationName = "SKETCH";
342 GEOMImpl_IPolyline2D aP2d(function);
343 Handle(TColStd_HArray1OfExtendedString) aNames = aP2d.GetNames();
345 if (aNames.IsNull() == Standard_False) {
346 if (aNames->Length() == 1) {
347 // This is the single curve. Make its full dump.
348 AddParam(theParams, "Name", aNames->Value(aNames->Lower()));
350 Handle(TColStd_HArray1OfByte) aTypes = aP2d.GetTypes();
352 if (aTypes.IsNull() == Standard_False && aTypes->Length() == 1) {
353 Standard_Integer aType = aTypes->Value(aTypes->Lower());
355 if (aType == GEOMImpl_ICurvesOperations::Polyline) {
356 AddParam(theParams, "Type") << "Polyline";
357 } else if (aType == GEOMImpl_ICurvesOperations::Interpolation) {
358 AddParam(theParams, "Type") << "Interpolation";
362 Handle(TColStd_HArray1OfByte) aCloseds = aP2d.GetClosedFlags();
364 if (aCloseds.IsNull() == Standard_False && aCloseds->Length() == 1) {
366 aCloseds->Value(aCloseds->Lower()) ? "Yes" : "No";
368 AddParam(theParams, "Is closed", aYesNo);
371 std::list <std::list <double> > aCoords;
373 aP2d.GetCoords(aCoords);
375 if (aCoords.size() == 1) {
376 AddParam(theParams, "Number of points", aCoords.front().size()/2);
379 // There are more than 1 curve.
380 Standard_Integer aNbCurves = aNames->Length();
382 std::list <std::list <double> > aCoords;
384 AddParam(theParams, "Number of curves", aNbCurves);
385 aP2d.GetCoords(aCoords);
387 Standard_Integer aNbCoords = aCoords.size();
388 std::list <std::list <double> >::const_iterator
389 anIt = aCoords.begin();
391 for (i = 0; i < aNbCurves; i++) {
392 TCollection_AsciiString aName("Curve ");
393 TCollection_ExtendedString
394 aValue(aNames->Value(aNames->Lower() + i));
396 aName.AssignCat(i + 1);
398 if (anIt != aCoords.end()) {
399 aValue.AssignCat(" (");
400 aValue.AssignCat(Standard_Integer(anIt->size()));
401 aValue.AssignCat(" points)");
405 AddParam(theParams, aName.ToCString(), aValue);
410 if (aType == POLYLINE2D_PLN_COORDS) {
411 Handle(TColStd_HArray1OfReal) aPln = aP2d.GetWorkingPlaneDbls();
413 if (aPln.IsNull() == Standard_False && aPln->Length() == 9) {
414 Standard_Integer i = aPln->Lower();
416 AddParam( theParams, "Origin")
417 << aPln->Value(i) << " "
418 << aPln->Value(i + 1) << " "
419 << aPln->Value(i + 2);
420 AddParam( theParams, "OZ")
421 << aPln->Value(i + 3) << " "
422 << aPln->Value(i + 4) << " "
423 << aPln->Value(i + 5);
424 AddParam( theParams, "OX")
425 << aPln->Value(i + 6) << " "
426 << aPln->Value(i + 7) << " "
427 << aPln->Value(i + 8);
430 AddParam(theParams, "Working plane", aP2d.GetWorkingPlane(), "XOY");
441 IMPLEMENT_STANDARD_HANDLE (GEOMImpl_PolylineDriver,GEOM_BaseDriver);
442 IMPLEMENT_STANDARD_RTTIEXT (GEOMImpl_PolylineDriver,GEOM_BaseDriver);