SMesh* TerrainManager::generateCubicMesh(SimpleVolume<MaterialDensityPair44>& volData)
{
      SMesh* sMesh = new SMesh();
      SurfaceMesh<PositionMaterial> vMesh;

//get the vertices and indices
     PolyVox::CubicSurfaceExtractor<SimpleVolume, MaterialDensityPair44 > surfaceEx(&volData, volData.getEnclosingRegion(), &vMesh,true);
     surfaceEx.execute();

     const std::vector<uint32_t>& vecIndices = vMesh.getIndices();
          const std::vector<PositionMaterial>& vecVertices = vMesh.getVertices();

//put the indices into an array

     u16 *arrIndices = new u16[vecIndices.size()];
     for(int i =0; i < vMesh.getNoOfIndices();i++)
     {
      arrIndices[i] =(u16)vecIndices[i];
     }

//put the vertices into an array of 3dvectors

     irr::core::vector3df *arrVertices = new irr::core::vector3df[vecVertices.size() ];
     for(int i =0; i < vecVertices.size();i++)
     {
       
        arrVertices[i] = irr::core::vector3df((f32)vecVertices[i].getPosition().getX(),
        (f32)vecVertices[i].getPosition().getY(),
        (f32)vecVertices[i].getPosition().getZ());
      
     }

    //add the vertices and indices to a buffer and add them to the mesh
     SMeshBuffer *buf = new SMeshBuffer();
     buf->append(arrVertices,vecVertices.size(),arrIndices,vecIndices.size());
     sMesh->addMeshBuffer(buf);
     return sMesh;
}

void TerrainManager::createSphereInVolume(SimpleVolume<MaterialDensityPair44>& volData, float fRadius)
{
      //This vector hold the position of the center of the volume
        Vector3DFloat v3dVolCenter(volData.getWidth() / 2, volData.getHeight() / 2, volData.getDepth() / 2);

        //This three-level for loop iterates over every voxel in the volume
        for (int z = 0; z < volData.getWidth(); z++)
        {
                for (int y = 0; y < volData.getHeight(); y++)
                {
                        for (int x = 0; x < volData.getDepth(); x++)
                        {
                                //Store our current position as a vector...
                                Vector3DFloat v3dCurrentPos(x,y,z);
                                //And compute how far the current position is from the center of the volume
                                float fDistToCenter = (v3dCurrentPos - v3dVolCenter).length();

                                //If the current voxel is less than 'radius' units from the center then we make it solid.
                                if(fDistToCenter <= fRadius)
                                {
                                        //Our new density value
                                        uint8_t uDensity = MaterialDensityPair44::getMaxDensity();

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

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

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