template<class IV >

void vtkfig::DataSet::SetSimplexGridBoundaryCells ( const IV & cells )

inline

Set boundary cells for simplex grids.

Template Parameters

IV	Vector class counting from zero with member functions size() and operator[]. std::vector will work.

Parameters

cr	Cell region numbers
name	Name of function

Definition at line 393 of file vtkfigDataSet.h.

   {
     if (this->boundary_data==NULL)
       this->boundary_data=vtkSmartPointer<vtkUnstructuredGrid>::New();
 
     auto budata=vtkUnstructuredGrid::SafeDownCast(this->boundary_data);
     budata->Reset();
 
     auto bgridpoints = vtkSmartPointer<vtkPoints>::New();
     budata->SetPoints(bgridpoints);
     
     auto udata=vtkUnstructuredGrid::SafeDownCast(this->data);
     auto pdata=udata->GetPoints();
     
     int np=pdata->GetNumberOfPoints();
     std::vector<int>pmask(np);
     for (size_t i=0;i<pmask.size(); i++) 
     {
       pmask[i]=-1;
     }
 
 
 
     int ip=0;
     for (int icell=0;icell<cells.size();icell+=this->spacedim)
     {
       for (int id=0;id<this->spacedim;id++)
       {
         if (pmask[cells[icell+id]]==-1)
         {
           pmask[cells[icell+id]]=ip++;
           double point[3];
           pdata->GetPoint(cells[icell+id],point);
           bgridpoints->InsertNextPoint(point[0],point[1],point[2]);
         }
 
       }
       if (this->spacedim==2)
       {
         vtkIdType        c[2]={pmask[cells[icell+0]],pmask[cells[icell+1]]};
         budata->InsertNextCell(VTK_LINE,2,c);
       }
       else
       {
         vtkIdType        c[3]={pmask[cells[icell+0]],pmask[cells[icell+1]],pmask[cells[icell+2]]};
         budata->InsertNextCell(VTK_TRIANGLE,3,c);
       }
 
     }
 
   }