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


}