int main(int argc, char *argv[])
{
   //Create and show the Qt OpenGL window
   QApplication app(argc, argv);
   OpenGLWidget openGLWidget(0);
   openGLWidget.show();

   //Create an empty volume and then place a sphere in it
   Volume<MaterialDensityPair44> volData(PolyVox::Region(Vector3DInt32(0,0,0), Vector3DInt32(63, 63, 63)));
   //createSphereInVolume(volData, 30);

   //Our new density value
   uint8_t uDensity = MaterialDensityPair44::getMaxDensity();

   //Get the old voxel
   MaterialDensityPair44 voxel;

   //Modify the density
   voxel.setDensity(uDensity);
   voxel.setMaterial(1);

   //Wrte the voxel value into the volume   
   volData.setVoxelAt(32,32,32, voxel);
   volData.setVoxelAt(32,32,33, voxel);
   volData.setVoxelAt(32,33,32, voxel);
   volData.setVoxelAt(32,33,33, voxel);
   volData.setVoxelAt(33,32,32, voxel);
   volData.setVoxelAt(33,32,33, voxel);
   volData.setVoxelAt(33,33,32, voxel);
   volData.setVoxelAt(33,33,33, voxel);

   //Extract the surface
   /*SurfaceMesh<PositionMaterialNormal> mesh;
   CubicSurfaceExtractorWithNormals<MaterialDensityPair44> surfaceExtractor(&volData, PolyVox::Region(Vector3DInt32(32,32,32), Vector3DInt32(33,33,33)), &mesh);
   surfaceExtractor.execute();*/

   SurfaceMesh<PositionMaterial> mesh;
   CubicSurfaceExtractor<MaterialDensityPair44> surfaceExtractor(&volData, PolyVox::Region(Vector3DInt32(32,32,32), Vector3DInt32(33,33,33)), &mesh);
   //CubicSurfaceExtractor<MaterialDensityPair44> surfaceExtractor(&volData, volData.getEnclosingRegion(), &mesh);
   surfaceExtractor.execute();

   std::cout << "Vertices = " << mesh.getNoOfVertices() << ", Indices = " << mesh.getNoOfIndices() << ", Triangles = " << mesh.getNoOfIndices()/3 << std::endl;

   //Pass the surface to the OpenGL window
   //openGLWidget.setSurfaceMeshToRender(mesh);

   //Run the message pump.
   return app.exec();
}

#include "MaterialDensityPair.h"
#include "CubicSurfaceExtractor.h"
#include "CubicSurfaceExtractorWithNormals.h"
#include "SurfaceMesh.h"
#include "Volume.h"

//#define USE_NORMALS

using namespace PolyVox;

int main(int argc, char *argv[])
{
   typedef MaterialDensityPair<uint8_t, 4, 4> Material44;

   // Create a volume
   Volume<Material44> volData(PolyVox::Region(Vector3DInt32(0,0,0), Vector3DInt32(63, 63, 63)));
   
   // Set up a voxel material
   uint8_t density = Material44::getMaxDensity();
   Material44 voxel;
   voxel.setDensity(density);
   voxel.setMaterial(1);
   
   // Set some voxels just inside the edge of the region to be extracted. Don't set any voxels on the other side of the region border
   volData.setVoxelAt(0, 0, 31, voxel);
   volData.setVoxelAt(1, 0, 31, voxel);
   volData.setVoxelAt(2, 0, 31, voxel);
   volData.setVoxelAt(3, 0, 31, voxel);

   // Extract the surface mesh from Region0 (contains voxels)
#ifdef USE_NORMALS
   SurfaceMesh<PositionMaterialNormal> meshRegion0;
   CubicSurfaceExtractorWithNormals<Material44> surfaceExtractorRegion0(&volData, PolyVox::Region(Vector3DInt32(0, 0, 0), Vector3DInt32(63, 63, 31)), &meshRegion0);
#else
   SurfaceMesh<PositionMaterial> meshRegion0;
   CubicSurfaceExtractor<Material44> surfaceExtractorRegion0(&volData, PolyVox::Region(Vector3DInt32(0, 0, 0), Vector3DInt32(63, 63, 31)), &meshRegion0);
#endif
   surfaceExtractorRegion0.execute();

   // Extract the surface mesh from Region1 (does not contain voxels)
#ifdef USE_NORMALS
   SurfaceMesh<PositionMaterialNormal> meshRegion1;
   CubicSurfaceExtractorWithNormals<Material44> surfaceExtractorRegion1(&volData, PolyVox::Region(Vector3DInt32(0, 0, 32), Vector3DInt32(63, 63, 63)), &meshRegion1);
#else
   SurfaceMesh<PositionMaterial> meshRegion1;
   CubicSurfaceExtractor<Material44> surfaceExtractorRegion1(&volData, PolyVox::Region(Vector3DInt32(0, 0, 32), Vector3DInt32(63, 63, 63)), &meshRegion1);
#endif
   surfaceExtractorRegion1.execute();

   // Check vertices on meshRegion0 - should be 20 vertices
   std::cout << "meshRegion0 Vertices = " << meshRegion0.getNoOfVertices() << ", Indices = " << meshRegion0.getNoOfIndices() << ", Triangles = " << meshRegion0.getNoOfIndices()/3 << std::endl;

   // Check vertices on meshRegion1 - should be 0 vertices
   std::cout << "meshRegion1 Vertices = " << meshRegion1.getNoOfVertices() << ", Indices = " << meshRegion1.getNoOfIndices() << ", Triangles = " << meshRegion1.getNoOfIndices()/3 << std::endl;

   return 0;
}