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();
 
 
}