main.cpp

#include "Grid/grid_dist_id.hpp"
 
int main(int argc, char* argv[])
{
 
    // Initialize the library
    openfpm_init(&argc,&argv);
 
    // 3D physical domain
    Box<3,float> domain({0.0,0.0,0.0},{1.0,1.0,1.0});
    
    // Grid size on eaxh dimension
    size_t sz[3] = {100,100,100};
 
    // Ghost part
    Ghost<3,float> g(0.1);
    
 
 
    grid_dist_id<3, float, aggregate<float[3]>> g_dist(sz,domain,g);
    
 
 
    // Get the iterator (No ghost)
    auto dom = g_dist.getDomainIterator();
    
    // Counter
    size_t count = 0;
    
    // Iterate over all the grid points
    while (dom.isNext())
    {
 
 
        // local grid key from iterator
        auto key = dom.get();
        
 
 
        auto key_g = g_dist.getGKey(key);
 
 
 
        g_dist.template get<0>(key)[0] = key_g.get(0);
        g_dist.template get<0>(key)[1] = key_g.get(1);
        g_dist.template get<0>(key)[2] = key_g.get(2);
        
        // Count the points
        count++;
 
 
 
        // next point
        ++dom;
    }
 
 
 
    g_dist.template ghost_get<0>();
    
 
 
    // Get the VCluster object
    Vcluster<> & vcl = create_vcluster();
 
    // queue an operation of sum for the counter count
    vcl.sum(count);
 
    // execute the operation
    vcl.execute();
    
    // only master output
    if (vcl.getProcessUnitID() == 0)
      std::cout << "Number of points: " << count << "\n";
 
 
 
    g_dist.write("output");
    
 
 
    g_dist.getDecomposition().write("out_dec");
    
 
 
    openfpm_finalize();
 
 
}