compilation on Linux

[modules/hydro.git] / src / HYDROData / HYDROData_BSplineOperation.cxx

diff --git a/src/HYDROData/HYDROData_BSplineOperation.cxx b/src/HYDROData/HYDROData_BSplineOperation.cxx
index 293dbd442b6c6032e802982bbf1b17e46f87ffca..f3d451f1654dc1beb5207e6869a0a1e60bb1ede4 100644 (file)
--- a/src/HYDROData/HYDROData_BSplineOperation.cxx
+++ b/src/HYDROData/HYDROData_BSplineOperation.cxx
@@ -1,3 +1,24 @@
+// Copyright (C) 2007-2015  CEA/DEN, EDF R&D, OPEN CASCADE
+//
+// Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
+// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
 
 #include <HYDROData_BSplineOperation.h>
 
@@ -11,15 +32,30 @@
 
 HYDROData_BSplineOperation::HYDROData_BSplineOperation(
   const NCollection_Sequence<gp_XYZ>& thePoints,
-  const bool                          theIsClosed )
+  const bool                          theIsClosed,
+  const double                        theTolerance )
 {
-  // fill array for algorithm by the received coordinates
+  // skip equal points
   int aNbPoints = thePoints.Size();
+  NCollection_Sequence<gp_XYZ> aPoints;
+  if ( aNbPoints > 0 ) {
+    gp_XYZ aPrevPoint = thePoints.Value( 1 );
+    aPoints.Append( aPrevPoint );
+    for( int i = 2 ; i <= aNbPoints; ++i )
+    {
+      gp_XYZ aPoint( thePoints.Value( i ) );
+      if ( !aPoint.IsEqual( aPrevPoint, theTolerance ) )
+        aPoints.Append( aPoint );
+      aPrevPoint = aPoint;
+    }
+  }
 
+  // fill array for algorithm by the received coordinates
+  aNbPoints = aPoints.Size();
   Handle(TColgp_HArray1OfPnt) aHCurvePoints = new TColgp_HArray1OfPnt( 1, aNbPoints );
   for ( int i = 1; i <= aNbPoints; i++ )
   {
-    gp_Pnt aPnt( thePoints.Value( i ) );
+    gp_Pnt aPnt( aPoints.Value( i ) );
     aHCurvePoints->SetValue( i, aPnt );
   }