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<2,float> domain({0.0,0.0},{1.0,1.0});
 
    // Here we define the boundary conditions of our problem
    size_t bc[2]={PERIODIC,PERIODIC};
 
    // extended boundary around the domain, and the processor domain
    Ghost<2,float> g(0.01);
    
 
 
    vector_dist<2,float, aggregate<float,float[3],float[3][3]> > vd(4096,domain,bc,g);
 
    // the scalar is the element at position 0 in the aggregate
    const int scalar = 0;
 
    // the vector is the element at position 1 in the aggregate
    const int vector = 1;
 
    // the tensor is the element at position 2 in the aggregate
    const int tensor = 2;
 
 
 
    auto it = vd.getDomainIterator();
 
    while (it.isNext())
    {
        // Particle p
        auto p = it.get();
 
        // we define x, assign a random position between 0.0 and 1.0
        vd.getPos(p)[0] = (float)rand() / RAND_MAX;
 
        // we define y, assign a random position between 0.0 and 1.0
        vd.getPos(p)[1] = (float)rand() / RAND_MAX;
 
        // the scalar property
        vd.template getProp<scalar>(p) = vd.getPos(p)[0];
 
        // The vector
        vd.template getProp<vector>(p)[0] = 1.0;
        vd.template getProp<vector>(p)[1] = 2.0;
        vd.template getProp<vector>(p)[2] = 3.0;
 
        // The tensor
        vd.template getProp<tensor>(p)[0][0] = 1.0;
        vd.template getProp<tensor>(p)[0][1] = 2.0;
        vd.template getProp<tensor>(p)[0][2] = 3.0;
        vd.template getProp<tensor>(p)[1][0] = 4.0;
        vd.template getProp<tensor>(p)[1][1] = 5.0;
        vd.template getProp<tensor>(p)[1][2] = 6.0;
        vd.template getProp<tensor>(p)[2][0] = 7.0;
        vd.template getProp<tensor>(p)[2][1] = 8.0;
        vd.template getProp<tensor>(p)[2][2] = 9.0;
 
        // next particle
        ++it;
    }
 
 
 
    vd.map();
 
 
 
    // write in vtk to see the result before
    vd.write("before_ghost_get");
 
    // Here we synchronize the ghost get only the scalar and the tensor property
    vd.ghost_get<scalar,tensor>();
 
    // write in vtk to see the result after
    vd.write("after_ghost_get");
 
 
 
    vd.ghost_get<vector>(KEEP_PROPERTIES);
    
    // write in vtk to see the result before
    vd.write("after_ghost_get_kp");
 
 
 
    it = vd.getGhostIterator();
 
    while (it.isNext())
    {
        // Particle p
        auto p = it.get();
 
        // we shift down he particles
        vd.getPos(p)[0] -= 1.0;
 
        // we shift
        vd.getPos(p)[1] -= 1.0;
 
        // The vector
        vd.template getProp<vector>(p)[0] = -1.0;
        vd.template getProp<vector>(p)[1] = -2.0;
        vd.template getProp<vector>(p)[2] = -3.0;
 
        // next particle
        ++it;
    }
 
    vd.ghost_get<vector>(KEEP_PROPERTIES | NO_POSITION);
 
    // write in vtk to see the result before
    vd.write("after_ghost_get_kp_no_pos");
 
    vd.ghost_get<vector>(KEEP_PROPERTIES);
 
    vd.write("after_ghost_get_kp_no_pos_restore");
 
 
 
    vd.ghost_put<add_,vector>();
 
    // write in vtk to see the result before
    vd.write("after_ghost_put");
 
 
 
    openfpm_finalize();
 
 
}