void createPerlinTerrainChunk(Volume<MaterialDensityPair44>& volData, PolyVox::Region reg)
{
    Perlin perlin(2,2,1,234);
    for(int x = reg.getLowerCorner().getX(); x < reg.getUpperCorner().getX(); x++) {
        for(int y = reg.getLowerCorner().getY(); y < reg.getUpperCorner().getY(); y++) {
            float perlinVal = perlin.Get(x / 254.0f , y / 254.0f);
            perlinVal += 1.0f;
            perlinVal *= 0.5f;
            perlinVal *= 255.0f;

            for(int z = reg.getLowerCorner().getZ(); z < reg.getUpperCorner().getZ(); z++) {
                MaterialDensityPair44 voxel;
                if(z < perlinVal)
                {
                    voxel.setMaterial(245);
                    voxel.setDensity(MaterialDensityPair44::getMaxDensity());
                }
                else
                {
                    voxel.setMaterial(0);
                    voxel.setDensity(MaterialDensityPair44::getMinDensity());
                }

                volData.setVoxelAt(x, y, z, voxel);
            }
        }
    }
}

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;
    volData.setBlockCacheSize(4096);
    volData.m_fLoad = boost::bind(&createPerlinTerrainChunk, boost::ref(volData), _1);
    //createSphereInVolume(volData, 30);
    //createPerlinTerrain(volData, Vector3DInt16(255,255,255));
    //createPerlinVolumeSlow(volData);
    std::cout << "Memory usage: " << volData.calculateSizeInBytes() << std::endl;
    //volData.setBlockCacheSize(8);
    std::cout << "Memory usage: " << volData.calculateSizeInBytes() << std::endl;
    std::cout << "Compression ratio: " << volData.calculateCompressionRatio() << std::endl;

    //Extract the surface
    SurfaceMesh<PositionMaterialNormal> mesh;
    SurfaceExtractor<MaterialDensityPair44> surfaceExtractor(&volData, PolyVox::Region(Vector3DInt16(1,1,1), Vector3DInt16(1022,253,253)), &mesh);
    surfaceExtractor.execute();

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

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