main.cpp

#include <stddef.h>
#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())
    {
        auto key = it.get();

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

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

        // next particle
        ++it;
    }



    vd.map();



    //Counter we use it later
    size_t cnt = 0;

    // 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
        
         // the scalar property
        vd.template getProp<scalar>(p) = 1.0;

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

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

        // increment the counter
        cnt++;

        // next particle
        ++it;
    }


    
    auto & v_cl = create_vcluster();
    v_cl.sum(cnt);
    v_cl.execute();
    


    openfpm::vector<std::string> names({"scalar","vector","tensor"});
    vd.setPropNames(names);

    // save vtk format (vtk is always the default)
    vd.write("particles");

    // save in vtk binary format
    vd.write("particles_bin",VTK_WRITER | FORMAT_BINARY);

    // save in vtk format with time
    vd.write("particles_with_time","time=1.234");

    // save in vtk format with time
    vd.write("particles_with_time_bin","time=1.234",VTK_WRITER | FORMAT_BINARY);



    openfpm_finalize();


}