PolyVox::Volume<PolyVox::MaterialDensityPair44> m_VolumeData(64, 64, 64);

       Vector3 v3dVolCenter(m_VolumeData.getWidth() / 2, m_VolumeData.getHeight() / 2, m_VolumeData.getDepth() / 2);

       //This three-level for loop iterates over every voxel in the volume
       for (int z = 0; z < m_VolumeData.getWidth(); z++)
       {
          for (int y = 0; y < m_VolumeData.getHeight(); y++)
          {
             for (int x = 0; x < m_VolumeData.getDepth(); x++)
             {
                //Store our current position as a vector...
                Vector3 v3dCurrentPos(x,y,z);
                //And compute how far the current position is from the center of the volume
                if (IsInRange(v3dCurrentPos, v3dCurrentPos, 5))
                {
                   //Our new density value
                   uint8_t uDensity = PolyVox::MaterialDensityPair44::getMaxDensity();

                   //Get the old voxel
                   PolyVox::MaterialDensityPair44 voxel = m_VolumeData.getVoxelAt(x,y,z);

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

                   //Wrte the voxel value into the volume
                   m_VolumeData.setVoxelAt(x, y, z, voxel);
                }
             }
          }
       }

        PolyVox::SurfaceMesh<PolyVox::PositionMaterialNormal> mesh;
       PolyVox::CubicSurfaceExtractorWithNormals<PolyVox::MaterialDensityPair44> surfaceExtractor(&m_VolumeData, m_VolumeData.getEnclosingRegion(), &mesh);

       surfaceExtractor.execute();

       const vector<uint32_t>& vecIndices = mesh.getIndices();
       const vector<PolyVox::PositionMaterialNormal>& vecVertices = mesh.getVertices();

       m_VertexBuffer.Create(sizeof(PolyVox::PositionMaterialNormal) * mesh.getNoOfVertices());

       m_IndexBuffer.Create(mesh.getNoOfIndices() * sizeof(DWORD), 8UL, D3DFMT_INDEX32);

       void* pVoid;
       m_IndexBuffer.Lock(&pVoid);

       const void* pIndices = static_cast<const void*>(&(vecIndices[0]));      

   

      memcpy((pVoid), pIndices, sizeof(uint32_t) * vecIndices.size());
   

       m_IndexBuffer.Unlock();

       const void* pVertices = static_cast<const void*>(&(vecVertices[0]));   

       m_VertexBuffer.Lock(&pVoid);

       count = 0;

       memcpy((pVoid), pVertices, sizeof(PolyVox::PositionMaterialNormal) * vecVertices.size());



       m_VertexBuffer.Unlock();

       m_NumVertices = vecVertices.size();
       m_NumFaces = vecIndices.size() / 3;

bool IsInRange(CONST D3DXVECTOR3& pos1, CONST D3DXVECTOR3& pos2, float range)
   {
      Vector3 diffVec = pos2 - pos1;

      float distance = diffVec.x * diffVec.x + diffVec.y * diffVec.y + diffVec.z * diffVec.z;

      range = range * range;

      if (distance <= range)
         return true;

      return false;
   }

-      mesh   {m_Region={...} m_iTimeStamp=-1 m_iNoOfLod0Tris=-858993460 ...}   PolyVox::SurfaceMesh<PolyVox::PositionMaterialNormal>
-      m_Region   {m_v3dLowerCorner={...} m_v3dUpperCorner={...} dummy='d' }   PolyVox::Region
-      m_v3dLowerCorner   {m_tElements=0x006ae504 }   PolyVox::Vector<3,short>
-      m_tElements   0x006ae504   short [3]
      [0]   0   short
      [1]   0   short
      [2]   0   short
-      m_v3dUpperCorner   {m_tElements=0x006ae50a }   PolyVox::Vector<3,short>
-      m_tElements   0x006ae50a   short [3]
      [0]   63   short
      [1]   63   short
      [2]   63   short
      dummy   100 'd'   unsigned char
      m_iTimeStamp   -1   int
      m_iNoOfLod0Tris   -858993460   int
+      m_vecTriangleIndices   [0]()   std::vector<unsigned int,std::allocator<unsigned int> >
+      m_vecVertices   [12288]({position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=0.00000000 },{position={...} normal={...} material=??? },...,...)   std::vector<PolyVox::PositionMaterialNormal,std::allocator<PolyVox::PositionMaterialNormal> >
-      m_vecLodRecords   [1]({beginIndex=0 endIndex=0 })   std::vector<PolyVox::LodRecord,std::allocator<PolyVox::LodRecord> >
      [size]   1   unsigned long
      [capacity]   1   unsigned long
+      [0]   {beginIndex=0 endIndex=0 }   PolyVox::LodRecord
-      m_mapUsedMaterials   [0]()   std::set<unsigned char,std::less<unsigned char>,std::allocator<unsigned char> >
      [comp]   less   std::less<unsigned char>
      pIndices   0xcccccccc   const void *
      count   -858993460   int

if (IsInRange(v3dCurrentPos, v3dCurrentPos, 5))