unique-method modified

extensions/VectorOperations.h

@@ -633,14 +633,19 @@ namespace vector_operations {
   inline void unique(std::vector<WorldVector<double> > &vec, double tol, std::vector<unsigned> &ind)
   {
     compareTol<WorldVector<double> > comp(tol);
-    unsigned newVec=0;
-    for(unsigned i=0; i<vec.size();++i) {
-      bool inNew=false;
-      for(unsigned j=0;j<newVec;++j) {
-	inNew=inNew||comp(vec[i],vec[j]);
-	if(inNew) break;
+    size_t newVec = 0;
+    for (size_t i = 0; i < vec.size(); ++i) {
+      bool inNew = false;
+      for (size_t j = 0; j < newVec; ++j) {
+	inNew = inNew || comp(vec[i], vec[j]);
+	if(inNew) 
+	  break;
+      }
+      if (!inNew) {
+	vec[newVec] = vec[i]; 
+	newVec++; 
+	ind.push_back(i);
       }
-      if(!inNew) {vec[newVec]=vec[i]; newVec++; ind.push_back(i);}
     }
     vec.erase(vec.begin()+newVec,vec.end());
   }