main.cpp

#include "Grid/grid_dist_id.hpp"
#include "data_type/aggregate.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();


}