main.cpp

#include "Vector/vector_dist.hpp"


int main(int argc, char* argv[])
{


    // initialize the library
    openfpm_init(&argc,&argv);

    // Here we define our domain a 2D box with internals from 0 to 1.0 for x and y
    Box<3,float> domain({0.0,0.0,0.0},{22.0,5.0,5.0});

    // Here we define the boundary conditions of our problem
    size_t bc[3]={PERIODIC,PERIODIC,PERIODIC};

    // extended boundary around the domain, and the processor domain
    Ghost<3,float> g(0.01);
    


    vector_dist<3,float, aggregate<float[3],float[3],float[3],float[3],float[3]> > vd(4096,domain,bc,g);



    auto it = vd.getDomainIterator();

    while (it.isNext())
    {
        auto key = it.get();

        // we define x, assign a random position between 0.0 and 1.0
        vd.getPos(key)[0] = 22.0*((float)rand() / RAND_MAX);

        // we define y, assign a random position between 0.0 and 1.0
        vd.getPos(key)[1] = 5.0*((float)rand() / RAND_MAX);

        // we define y, assign a random position between 0.0 and 1.0
        vd.getPos(key)[1] = 5.0*((float)rand() / RAND_MAX);

        // next particle
        ++it;
    }



    vd.map();



    // Get a particle iterator
    it = vd.getDomainIterator();

    // For each particle ...
    while (it.isNext())
    {
        // ... p
        auto p = it.get();

        // we set the properties of the particle p

        vd.template getProp<0>(p)[0] = 1.0;
        vd.template getProp<0>(p)[1] = 1.0;
        vd.template getProp<0>(p)[2] = 1.0;

        vd.template getProp<1>(p)[0] = 2.0;
        vd.template getProp<1>(p)[1] = 2.0;
        vd.template getProp<1>(p)[2] = 2.0;

        vd.template getProp<2>(p)[0] = 3.0;
        vd.template getProp<2>(p)[1] = 3.0;
        vd.template getProp<2>(p)[2] = 3.0;

        vd.template getProp<3>(p)[0] = 4.0;
        vd.template getProp<3>(p)[1] = 4.0;
        vd.template getProp<3>(p)[2] = 4.0;

        vd.template getProp<4>(p)[0] = 5.0;
        vd.template getProp<4>(p)[1] = 5.0;
        vd.template getProp<4>(p)[2] = 5.0;


        // next particle
        ++it;
    }


    
    vd.save("particles_save.hdf5");



    vector_dist<3,float, aggregate<float[3],float[3],float[3],float[3],float[3]> > vd2(vd.getDecomposition(),0);

    // load the particles on another vector
    vd2.load("particles_save.hdf5");



    vector_dist<3,float, aggregate<float> > vd3(vd.getDecomposition(),0);

    auto it2 = vd2.getDomainIterator();

    while (it2.isNext())
    {
        auto p = it2.get();

        float magn_vort = sqrt(vd2.getProp<0>(p)[0]*vd2.getProp<0>(p)[0] +
                               vd2.getProp<0>(p)[1]*vd2.getProp<0>(p)[1] +
                               vd2.getProp<0>(p)[2]*vd2.getProp<0>(p)[2]);

        if (magn_vort < 0.1)
        {
            ++it2;
            continue;
        }

        vd3.add();

        vd3.getLastPos()[0] = vd2.getPos(p)[0];
        vd3.getLastPos()[1] = vd2.getPos(p)[1];
        vd3.getLastPos()[2] = vd2.getPos(p)[2];

        vd3.template getLastProp<0>() = magn_vort;

        ++it2;
    }

    // We write the vtk file out from vd3
    vd3.write("output", VTK_WRITER | FORMAT_BINARY );
    vd3.save("particles_post_process_2");



    openfpm_finalize();


}