From 9c9f270183b15ec0f3af45a313ef557302bc9186 Mon Sep 17 00:00:00 2001
From: Mike Taverne <mtaverne@bluemarlin.engits.com>
Date: Fri, 24 Apr 2009 15:45:22 +0200
Subject: [PATCH] -Added Set2Vector + Vector2Set functions to egvtkobject
 -Create DeleteSetOfPoints() function in class Operation

---
 egvtkobject.h        |  19 ++++++
 guismoothsurface.cpp |  54 +++++++++++++++-
 operation.cpp        | 173 +++++++++++++++++++++++++++++++++++++++++++++++++++
 operation.h          |   3 +-
 4 files changed, 247 insertions(+), 2 deletions(-)

diff --git a/egvtkobject.h b/egvtkobject.h
index 8b5cdcd..7346777 100644
--- a/egvtkobject.h
+++ b/egvtkobject.h
@@ -696,6 +696,25 @@ ostream &operator<<(ostream &out, QVector < pair<T1,T2> > & vector)
   out<<"]";
   return(out);
 }
+
+template <class T>
+QVector <T> Set2Vector(QSet <T> a_set, bool a_sort)
+{
+  QVector <T> l_vector(a_set.size());
+  qCopy(a_set.begin(),a_set.end(),l_vector.begin());
+  if(a_sort) qSort(l_vector.begin(),l_vector.end());
+  return(l_vector);
+}
+
+template <class T>
+QSet <T> Vector2Set(QVector <T> a_vector, bool a_sort)
+{
+  QSet <T> l_set(a_vector.size());
+  qCopy(a_vector.begin(),a_vector.end(),l_set.begin());
+  if(a_sort) qSort(l_set.begin(),l_set.end());
+  return(l_set);
+}
+
 ///////////////////////////////////////////
 // ///////////////////////////////////////////
 // /* Here is how we we get QTextStreams that look like iostreams */
diff --git a/guismoothsurface.cpp b/guismoothsurface.cpp
index 28d84ae..1975549 100644
--- a/guismoothsurface.cpp
+++ b/guismoothsurface.cpp
@@ -210,7 +210,7 @@ void GuiSmoothSurface::before()
   ui.SmoothMethod->addItem("Method 10: Super smoothing :)");
   ui.SmoothMethod->addItem("Method 11: Update current mesh density + node types");
   ui.SmoothMethod->addItem("Method 12: Delete all possible points :)");
-  ui.SmoothMethod->addItem("Method 13");
+  ui.SmoothMethod->addItem("Method 13: Delete selected points");
   ui.SmoothMethod->addItem("Method 14");
   ui.SmoothMethod->addItem("Method 15");
   ui.SmoothMethod->addItem("Method 16");
@@ -1286,6 +1286,58 @@ void GuiSmoothSurface::operate()
     
   }
   //////////////////////////////////////////////////////////////////////////////////////////////
+  else if(ui.SmoothMethod->currentIndex()==13)// Delete selected points
+  {
+    QSet<int> bcs;
+    getSelectedItems(ui.listWidget, bcs);
+    CreateSpecialMapping toto;
+    
+    SetConvergence(ui.doubleSpinBox_Convergence->value());
+    SetFeatureEdgeSmoothing(ui.checkBox_FeatureEdgeSmoothing->checkState());
+    SetFeatureAngle(ui.doubleSpinBox_FeatureAngle->value());
+    SetEdgeAngle(ui.doubleSpinBox_EdgeAngle->value());
+    SetBoundarySmoothing(ui.checkBox_BoundarySmoothing->checkState());
+    
+    int N_newpoints;
+    int N_newcells;
+    
+    vtkIdType N_points=grid->GetNumberOfPoints();
+    vtkIdType N_cells=grid->GetNumberOfCells();
+    
+    QVector <VertexMeshDensity> VMDvector=GetSet();
+//     cout<<"VMDvector="<<VMDvector<<endl;
+    for(unsigned int i=0;i<VMDvector.size();i++)
+    {
+      cout<<"VMDvector["<<i<<"].nodeset="<<VMDvector[i].nodeset<<endl;
+      int N_newpoints=0;
+      int N_newcells=0;
+      DeleteSetOfPoints(grid, VMDvector[i].nodeset, N_newpoints, N_newcells);
+    }
+    
+/*    vector <vtkIdType> nodeId_vector(3);
+    
+    for(vtkIdType i)*/
+    
+    
+/*    bool Global_DelResult=true;
+    while(Global_DelResult)
+    {
+      Global_DelResult=false;
+      vtkIdType DeadNode=0;
+      while(DeadNode<grid->GetNumberOfPoints())
+      {
+        bool Local_DelResult=true;
+        while(Local_DelResult)
+        {
+          Local_DelResult=DeletePoint_2(grid,DeadNode,N_newpoints,N_newcells);
+          if(Local_DelResult) Global_DelResult=true;
+        }
+        DeadNode++;
+      }
+    }*/
+    
+  }
+  //////////////////////////////////////////////////////////////////////////////////////////////
   else
   {
     cout<<"UNKNOWN METHOD"<<endl;
diff --git a/operation.cpp b/operation.cpp
index 51e283a..41816c9 100644
--- a/operation.cpp
+++ b/operation.cpp
@@ -36,6 +36,7 @@
 #include "vtkFloatArray.h"
 #include "vtkMath.h"
 #include <vtkCharArray.h>
+#include "egvtkobject.h"
 
 #include "geometrytools.h"
 using namespace GeometryTools;
@@ -1994,6 +1995,178 @@ bool Operation::DeletePoint_2(vtkUnstructuredGrid *src, vtkIdType DeadNode, int&
 }
 //End of DeletePoint_2
 
+bool Operation::DeleteSetOfPoints(vtkUnstructuredGrid *src, QSet <vtkIdType> DeadNodes, int& N_newpoints, int& N_newcells)
+{
+  QVector <vtkIdType> DeadNode_vector=Set2Vector(DeadNodes,false);
+  
+  getAllSurfaceCells(cells,src);
+//   getNodesFromCells(cells, nodes, src);
+  setGrid(src);
+  setCells(cells);
+  UpdateNodeType();
+  
+  //src grid info
+  int N_points=src->GetNumberOfPoints();
+  int N_cells=src->GetNumberOfCells();
+  
+  QSet <vtkIdType> DeadCells;
+  QSet <vtkIdType> MutatedCells;
+  QSet <vtkIdType> MutilatedCells;
+  QVector <vtkIdType> SnapPoint(DeadNode_vector.size());
+  
+  for(int i=0;i<DeadNode_vector.size();i++)
+  {
+//   N_newpoints=-1;
+//   N_newcells=0;
+    
+    if(DeadNode_vector[i]<0 || DeadNode_vector[i]>=N_points)
+    {
+      cout<<"Warning: Point out of range: DeadNode_vector[i]="<<DeadNode_vector[i]<<" N_points="<<N_points<<endl;
+      return(false);
+    }
+    
+    if(DebugLevel>10) {
+      cout<<"BEFORE FINDSNAPPOINT"<<endl;
+      cout<<"N_points="<<N_points<<endl;
+      cout<<"N_cells="<<N_cells<<endl;
+      cout<<"N_newpoints="<<N_newpoints<<endl;
+      cout<<"N_newcells="<<N_newcells<<endl;
+    }
+    SnapPoint[i]=FindSnapPoint(src,DeadNode_vector[i],DeadCells,MutatedCells,MutilatedCells, N_newpoints, N_newcells);
+    
+    if(DebugLevel>0) cout<<"===>DeadNode_vector[i]="<<DeadNode_vector[i]<<" moving to SNAPPOINT="<<SnapPoint[i]<<" DebugLevel="<<DebugLevel<<endl;
+    if(SnapPoint[i]<0) {cout<<"Sorry no possible SnapPoint found."<<endl; return(false);}
+    
+  }
+  //allocate
+  if(DebugLevel>10) {
+    cout<<"BEFORE ALLOCATION"<<endl;
+    cout<<"N_points="<<N_points<<endl;
+    cout<<"N_cells="<<N_cells<<endl;
+    cout<<"N_newpoints="<<N_newpoints<<endl;
+    cout<<"N_newcells="<<N_newcells<<endl;
+  }
+  N_points=src->GetNumberOfPoints();
+  N_cells=src->GetNumberOfCells();
+  
+  if(DebugLevel>47) {
+    cout<<"N_points="<<N_points<<endl;
+    cout<<"N_cells="<<N_cells<<endl;
+    cout<<"N_newpoints="<<N_newpoints<<endl;
+    cout<<"N_newcells="<<N_newcells<<endl;
+  }
+  EG_VTKSP(vtkUnstructuredGrid,dst);
+  allocateGrid(dst,N_cells+N_newcells,N_points+N_newpoints);
+  
+  //vector used to redefine the new point IDs
+  QVector <vtkIdType> OffSet(N_points);
+  
+  //copy undead points
+  vtkIdType dst_id_node=0;
+  for (vtkIdType src_id_node = 0; src_id_node < N_points; src_id_node++) {//loop through src points
+    if(!DeadNode_vector.contains(src_id_node))//if the node isn't dead, copy it
+    {
+      vec3_t x;
+      src->GetPoints()->GetPoint(src_id_node, x.data());
+      dst->GetPoints()->SetPoint(dst_id_node, x.data());
+      copyNodeData(src, src_id_node, dst, dst_id_node);
+      OffSet[src_id_node]=src_id_node-dst_id_node;
+      dst_id_node++;
+    }
+    else
+    {
+      if(DebugLevel>0) cout<<"src_id_node="<<src_id_node<<" dst_id_node="<<dst_id_node<<endl;
+    }
+  };
+  if(DebugLevel>10) {
+    cout<<"DeadCells="<<DeadCells<<endl;
+    cout<<"MutatedCells="<<MutatedCells<<endl;
+    cout<<"MutilatedCells="<<MutilatedCells<<endl;
+  }
+  //Copy undead cells
+  for (vtkIdType id_cell = 0; id_cell < src->GetNumberOfCells(); ++id_cell) {//loop through src cells
+    if(!DeadCells.contains(id_cell))//if the cell isn't dead
+    {
+      vtkIdType src_N_pts, *src_pts;
+      vtkIdType dst_N_pts, *dst_pts;
+      src->GetCellPoints(id_cell, src_N_pts, src_pts);
+      
+      vtkIdType type_cell = src->GetCellType(id_cell);
+      if(DebugLevel>10) cout<<"-->id_cell="<<id_cell<<endl;
+      if(DebugLevel>10) for(int i=0;i<src_N_pts;i++) cout<<"src_pts["<<i<<"]="<<src_pts[i]<<endl;
+//       src->GetCellPoints(id_cell, dst_N_pts, dst_pts);
+      dst_N_pts=src_N_pts;
+      dst_pts=new vtkIdType[dst_N_pts];
+      if(MutatedCells.contains(id_cell))//mutated cell
+      {
+        if(DebugLevel>10) cout<<"processing mutated cell "<<id_cell<<endl;
+        for(int i=0;i<src_N_pts;i++)
+        {
+          int DeadIndex = DeadNode_vector.indexOf(src_pts[i]);
+          if(DeadIndex!=-1) {
+            if(DebugLevel>10) {
+              cout<<"SnapPoint="<<SnapPoint[DeadIndex]<<endl;
+              cout<<"OffSet[SnapPoint]="<<OffSet[SnapPoint[DeadIndex]]<<endl;
+              cout<<"src_pts["<<i<<"]="<<src_pts[i]<<endl;
+            }
+            dst_pts[i]=SnapPoint[DeadIndex]-OffSet[SnapPoint[DeadIndex]];
+          }
+          else dst_pts[i]=src_pts[i]-OffSet[src_pts[i]];
+        }
+        if(DebugLevel>10) cout<<"--->dst_pts:"<<endl;
+        if(DebugLevel>10) for(int i=0;i<dst_N_pts;i++) cout<<"dst_pts["<<i<<"]="<<dst_pts[i]<<endl;
+        
+      }
+      else if(MutilatedCells.contains(id_cell))//mutilated cell (ex: square becoming triangle)
+      {
+        cout<<"FATAL ERROR: Quads not supported yet."<<endl;EG_BUG;
+        
+        if(DebugLevel>10) cout<<"processing mutilated cell "<<id_cell<<endl;
+        
+        if(type_cell==VTK_QUAD) {
+          type_cell=VTK_TRIANGLE;
+          dst_N_pts-=1;
+        }
+        else {cout<<"FATAL ERROR: Unknown mutilated cell detected! It is not a quad! Potential xenomorph infestation!"<<endl;EG_BUG;}
+        //merge points
+        int j=0;
+        for(int i=0;i<src_N_pts;i++)
+        {
+/*          if(src_pts[i]==SnapPoint) { dst_pts[j]=SnapPoint-OffSet[SnapPoint];j++; }//SnapPoint
+          else if(src_pts[i]!=DeadNode_vector[i]) { dst_pts[j]=src_pts[i]-OffSet[src_pts[i]];j++; }//pre-snap/dead + post-snap/dead*/
+          //do nothing in case of DeadNode_vector[i]
+        }
+      }
+      else//normal cell
+      {
+        if(DebugLevel>10) cout<<"processing normal cell "<<id_cell<<endl;
+        for(int i=0;i<src_N_pts;i++)
+        {
+          dst_pts[i]=src_pts[i]-OffSet[src_pts[i]];
+        }
+      }
+      //copy the cell
+      vtkIdType id_new_cell = dst->InsertNextCell(type_cell, dst_N_pts, dst_pts);
+      copyCellData(src, id_cell, dst, id_new_cell);
+      if(DebugLevel>10) {
+        cout<<"===Copying cell "<<id_cell<<" to "<<id_new_cell<<"==="<<endl;
+        cout<<"src_pts:"<<endl;
+        for(int i=0;i<src_N_pts;i++) cout<<"src_pts["<<i<<"]="<<src_pts[i]<<endl;
+        cout<<"dst_pts:"<<endl;
+        for(int i=0;i<dst_N_pts;i++) cout<<"dst_pts["<<i<<"]="<<dst_pts[i]<<endl;
+        cout<<"OffSet="<<OffSet<<endl;
+        cout<<"===Copying cell end==="<<endl;
+      }
+      delete dst_pts;
+    }
+  };
+  
+//   cout_grid(cout,dst,true,true,true,true);
+  makeCopy(dst, src);
+  return(true);
+}
+//End of DeleteSetOfPoints
+
 void Operation::TxtSave(QString a_filename)
 {
   cout << a_filename.toAscii().data() << endl;
diff --git a/operation.h b/operation.h
index 60fd385..125671c 100644
--- a/operation.h
+++ b/operation.h
@@ -178,7 +178,8 @@ public: // methods
   
   vtkIdType FindSnapPoint(vtkUnstructuredGrid *src, vtkIdType DeadNode,QSet <vtkIdType> & DeadCells,QSet <vtkIdType> & MutatedCells,QSet <vtkIdType> & MutilatedCells, int& N_newpoints, int& N_newcells);
   bool DeletePoint_2(vtkUnstructuredGrid *src, vtkIdType DeadNode, int& N_newpoints, int& N_newcells);
- 
+  bool DeleteSetOfPoints(vtkUnstructuredGrid *src, QSet <vtkIdType> DeadNodes, int& N_newpoints, int& N_newcells);
+    
   void TxtSave(QString a_filename);
   void DualSave(QString a_filename);
     
-- 
2.11.4.GIT