text: /*=========================================================================
   *
   *  Copyright Insight Software Consortium
   *
   *  Licensed under the Apache License, Version 2.0 (the "License");
   *  you may not use this file except in compliance with the License.
   *  You may obtain a copy of the License at
   *
   *         http://www.apache.org/licenses/LICENSE-2.0.txt
   *
   *  Unless required by applicable law or agreed to in writing, software
   *  distributed under the License is distributed on an "AS IS" BASIS,
   *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
   *  See the License for the specific language governing permissions and
   *  limitations under the License.
   *
   *=========================================================================*/
  #include "itkQuadEdgeMesh.h"
  #include "itkVTKPolyDataReader.h"
  
  #include "itkQuadEdgeMeshExtendedTraits.h"
  #include "itkDiscreteGaussianCurvatureQuadEdgeMeshFilter.h"
  #include "itkQuadEdgeMeshScalarDataVTKPolyDataWriter.h"
  
  int itkDiscreteGaussianCurvatureQuadEdgeMeshFilterTest( int argc, char* argv[] )
  {
    if( argc < 2 )
      {
      std::cout <<"*** GaussianCurvature ***" <<std::endl;
      std::cout <<"This example requires at least one argument:" <<std::endl;
      std::cout <<" 1- FileName" <<std::endl;
      return EXIT_FAILURE;
      }
  
    const unsigned int Dimension = 3;
    typedef double CoordType;
  
    typedef itk::QuadEdgeMeshExtendedTraits <
      CoordType,
      Dimension,
      2,
      CoordType,
      CoordType,
      CoordType,
      bool,
      bool > Traits;
  
    typedef itk::QuadEdgeMesh< CoordType, Dimension, Traits > MeshType;
    typedef itk::DiscreteGaussianCurvatureQuadEdgeMeshFilter<
      MeshType, MeshType > CurvatureFilterType;
  
    typedef itk::VTKPolyDataReader< MeshType > ReaderType;
  
    ReaderType::Pointer reader = ReaderType::New( );
    reader->SetFileName( argv[1] );
    try
      {
      reader->Update( );
      }
    catch( itk::ExceptionObject & excp )
      {
      std::cerr << "Exception thrown while reading the input file " << std::endl;
      std::cerr << excp << std::endl;
      return EXIT_FAILURE;
      }
  
    MeshType::Pointer mesh = reader->GetOutput();
  
    CurvatureFilterType::Pointer gaussian_curvature = CurvatureFilterType::New();
    gaussian_curvature->SetInput( mesh );
    gaussian_curvature->Update();
  
    MeshType::Pointer output = gaussian_curvature->GetOutput();
  
    typedef itk::QuadEdgeMeshScalarDataVTKPolyDataWriter< MeshType > WriterType;
    WriterType::Pointer writer = WriterType::New();
    writer->SetInput( output );
    writer->SetFileName( "gaussian_curvature.vtk" );
    writer->Update();
  
    // ** PRINT **
    std::cout << gaussian_curvature;
    return EXIT_SUCCESS;
  }
  