appli/examples/example_SimplexMeshDeleteEdge.cxx

   1 #include <cmath>
   2 #include <cstdlib>
   3 #include <iostream>
   4 #include <string>
   5
   6 #include <cpm/DataStructures/QuadEdgeMesh.h>
   7 #include <cpm/DataStructures/SimplexMesh.h>
   8 #include <cpm/Algorithms/QuadEdge/MeshToDualFilter.h>
   9 #include <cpm/IO/MeshReader.h>
  10
  11 #include <vtkActor.h>
  12 #include <vtkCallbackCommand.h>
  13 #include <vtkProperty.h>
  14 #include <vtkRenderer.h>
  15 #include <vtkRenderWindow.h>
  16 #include <vtkRenderWindowInteractor.h>
  17 #include <vtkSmartPointer.h>
  18 #include <cpm/VTK/MeshMapper.h>
  19
  20 // -------------------------------------------------------------------------
  21 const unsigned int Dimension = 3;
  22 const unsigned int Order = 2;
  23
  24 typedef float TScalar;
  25 typedef cpm::DataStructures::QuadEdgeMesh< TScalar, Dimension > TMesh;
  26 typedef cpm::DataStructures::SimplexMesh< TScalar, Order, Dimension > TSimplex;
  27 typedef cpm::VTK::MeshMapper< TMesh > TMeshMapper;
  28
  29 // -------------------------------------------------------------------------
  30 int main( int argc, char* argv[] )
  31 {
  32   if( argc < 2 )
  33   {
  34     std::cerr
  35       << "Usage: " << argv[ 0 ]
  36       << " input_mesh"
  37       << std::endl;
  38     return( 1 );
  39
  40   } // fi
  41
  42   typedef cpm::IO::MeshReader< TMesh > TReader;
  43   TReader::Pointer reader = TReader::New( );
  44   reader->SetFileName( argv[ 1 ] );
  45
  46   typedef
  47     cpm::Algorithms::QuadEdge::MeshToDualFilter< TMesh, TSimplex >
  48     TDualFilter;
  49   TDualFilter::Pointer dual_filter = TDualFilter::New( );
  50   dual_filter->SetInput( reader->GetOutput( ) );
  51   dual_filter->Update( );
  52   TSimplex::Pointer simplex = dual_filter->GetOutput( );
  53
  54   std::cout
  55     << " --> Read " << simplex->GetNumberOfPoints( )
  56     << " simplex points."
  57     << std::endl;
  58
  59   bool finish = false;
  60   while( !finish )
  61   {
  62     std::cout
  63       << " --> Delete an edge between a and b (-1 to finish) : "
  64       << std::flush;
  65     int a, b;
  66     std::cin >> a;
  67     if( a == -1 )
  68     {
  69       finish = true;
  70       continue;
  71
  72     } // fi
  73     std::cin >> b;
  74     if( simplex->OperatorDeleteEdge( a, b ) )
  75       std::cout << "Deleted" << std::endl;
  76     else
  77       std::cout << "Not present" << std::endl;
  78
  79   } // elihw
  80
  81   // Map mesh
  82   vtkSmartPointer< TMeshMapper > mapper =
  83     vtkSmartPointer< TMeshMapper >::New( );
  84   mapper->SetInputData( reader->GetOutput( ) );
  85
  86   // Map dual mesh
  87   vtkSmartPointer< TMeshMapper > dual_mapper =
  88     vtkSmartPointer< TMeshMapper >::New( );
  89   dual_mapper->SetInputData( dual_filter->GetOutput( ) );
  90
  91   // Create actor
  92   vtkSmartPointer< vtkActor > actor =
  93     vtkSmartPointer< vtkActor >::New( );
  94   actor->SetMapper( mapper );
  95   actor->GetProperty( )->SetColor( 1, 1, 0 );
  96   actor->GetProperty( )->SetOpacity( 0.5 );
  97
  98   // Create dual actor
  99   vtkSmartPointer< vtkActor > dual_actor =
 100     vtkSmartPointer< vtkActor >::New( );
 101   dual_actor->SetMapper( dual_mapper );
 102   dual_actor->GetProperty( )->SetColor( 1, 0, 0 );
 103   dual_actor->GetProperty( )->SetOpacity( 1 );
 104   dual_actor->GetProperty( )->SetLineWidth( 2 );
 105   dual_actor->GetProperty( )->SetRepresentationToWireframe( );
 106
 107   // Configure visualization objects
 108   vtkSmartPointer< vtkRenderer > renderer =
 109     vtkSmartPointer< vtkRenderer >::New( );
 110   renderer->SetBackground( 0.1, 0.3, 0.5 );
 111
 112   vtkSmartPointer< vtkRenderWindow > window =
 113     vtkSmartPointer< vtkRenderWindow >::New( );
 114   window->AddRenderer( renderer );
 115   window->SetSize( 800, 800 );
 116
 117   // Set up the interaction
 118   vtkSmartPointer< vtkRenderWindowInteractor > interactor =
 119     vtkSmartPointer< vtkRenderWindowInteractor >::New( );
 120   window->SetInteractor( interactor );
 121
 122   // Associate actors
 123   //renderer->AddActor( actor );
 124   renderer->AddActor( dual_actor );
 125
 126   // Begin interaction
 127   window->Render( );
 128   interactor->Start( );
 129
 130   return( 0 );
 131 }
 132
 133 // eof - $RCSfile$