#include "itkSWaveletTest.h"

#include <algorithm>
#include <fstream>


int main( int argc, char * argv [] )
{

  if (argc != 2) {
    cerr<<"format: itkSWaveletTest levelOfSubdivision"<<endl;
    exit(-1);
  }


  SphereMeshSourceType::Pointer  mySphereMeshSource = SphereMeshSourceType::New();

  PointType center; 
  center.Fill( 0.0 );

  VectorType scale;
  scale.Fill( 1.0 );
  
  int n = atoi(argv[1]);
  if (n>=7 && n<0) 
    {
      std::cerr<<"Finest resolution can't be less than 0, and best to be no larger than 6..."<<std::endl;
      exit(-1);
    }

  mySphereMeshSource->SetCenter( center );
  mySphereMeshSource->SetResolution( n );
  mySphereMeshSource->SetScale( scale );

  mySphereMeshSource->Modified();

  try
    {
      mySphereMeshSource->Update();
    }
  catch( itk::ExceptionObject & excp )
    {
      std::cerr << "Error during Update() " << std::endl;
      std::cerr << excp << std::endl;
    }

  // testing task 1:
  std::cerr<<"Testing task 1:"<<std::endl; 
  std::cerr<<"    Given different resolution in the command line, generate and visulize subdivided mesh at different resolutions."<<std::endl;

  MeshType::Pointer myMesh = mySphereMeshSource->GetOutput();

  displayITKmesh(myMesh);

  return EXIT_SUCCESS;
}