//////////////////////////////////////////////////////////////////////////////////
////////////////////    Skull Test With Specified Directions    ////////////////
//////////////////////////////////////////////////////////////////////////////////


// Include head files
#include "vtkPolyDataReader.h"
#include <string>
#include <iostream>

#include "vtkRadialThicknessCalculate.h"


int main( int argc, char *argv[] )
{
	// class object initialization
	vtkRadialThicknessCalculate radialThicknessCal;

	// read VTK data 1, which is skull inner surface
	vtkPolyDataReader* reader1 = vtkPolyDataReader::New();
	std::string fileName = "skull_inner.vtk";
	std::cout << "Reading " << fileName.c_str() << std::endl; 
	reader1->SetFileName(fileName.c_str());
	reader1->Update();
	vtkPolyData* inner = reader1->GetOutput();

	// read VTK data 2, which is skull outer surface
	vtkPolyDataReader* reader2 = vtkPolyDataReader::New();
	fileName = "skull_outer.vtk";
	std::cout << "Reading " << fileName.c_str() << std::endl; 
	reader2->SetFileName(fileName.c_str());
	reader2->Update();
	vtkPolyData* outer = reader2->GetOutput();


	printf("\n\n Skull Test With Specified Directions \n\n");
	radialThicknessCal.SetMasterMesh(inner);
	radialThicknessCal.SetSupplementaryMesh(outer);
	
	// set file name containing directions
	radialThicknessCal.SetDirectionsFileName("directions");
	
	// set thickness file name as the output
	radialThicknessCal.SetThicknessFileName("skull_directions");

	radialThicknessCal.StartThicknessCalculate();
	printf("Done!\n\n");


	return 0;
}