Domain_icells_cart_unit_test.cu

#define BOOST_TEST_DYN_LINK
#include <boost/test/unit_test.hpp>
#include "Space/Shape/Box.hpp"
#include "Decomposition/Domain_icells_cart.hpp"
#include "VCluster/VCluster.hpp"
#include "VTKWriter/VTKWriter.hpp"
#include "Grid/iterators/grid_skin_iterator.hpp"

BOOST_AUTO_TEST_SUITE( domain_icells_cart )

BOOST_AUTO_TEST_CASE( domain_icells_use )
{
#if 0

    domain_icell_calculator<3,float,memory_traits_inte,CudaMemory> dcc;

    openfpm::vector_gpu<SpaceBox<3,float>> domain_proc;

    Box<3,float> box1({0.0,0.0,0.0},{0.5,0.5,0.5});
    Box<3,float> box2({0.5,0.25,0.25},{0.75,0.5,0.5});

    domain_proc.add(box1);
    domain_proc.add(box2);

    openfpm::vector_gpu<Box<3,float>> ie_ghost;

    // first sub ghost area

    Box<3,float> g1({0.0,0.0,0.0},{0.1,0.5,0.5});
    Box<3,float> g2({0.0,0.0,0.4},{0.5,0.5,0.5});
    Box<3,float> g3({0.0,0.0,0.0},{0.5,0.1,0.5});
    Box<3,float> g4({0.0,0.4,0.0},{0.5,0.5,0.5});
    Box<3,float> g5({0.0,0.0,0.0},{0.5,0.5,0.1});
    Box<3,float> g6_1({0.4,0.0,0.0},{0.5,0.25,0.5});
    Box<3,float> g6_2({0.4,0.0,0.0},{0.5,0.5,0.25});

    // Second sub ghost area

    Box<3,float> g7({0.5,0.25,0.25},{0.6,0.5,0.5});
    Box<3,float> g8({0.65,0.25,0.25},{0.75,0.5,0.5});
    Box<3,float> g9({0.5,0.25,0.25},{0.75,0.35,0.25});
    Box<3,float> g10({0.5,0.4,0.25},{0.75,0.5,0.5});
    Box<3,float> g11({0.5,0.25,0.25},{0.75,0.5,0.35});
    Box<3,float> g12({0.5,0.25,0.4},{0.75,0.5,0.5});

    ie_ghost.add(g1);
    ie_ghost.add(g2);
    ie_ghost.add(g3);
    ie_ghost.add(g4);
    ie_ghost.add(g5);
    ie_ghost.add(g6_1);
    ie_ghost.add(g6_2);

    ie_ghost.add(g7);
    ie_ghost.add(g8);
    ie_ghost.add(g9);
    ie_ghost.add(g10);
    ie_ghost.add(g11);
    ie_ghost.add(g12);

    Box<3,float> pbox({0.0,0.0,0.0},{0.75,0.5,0.5});

    Ghost<3,float> g(0.1);

    auto & v_cl = create_vcluster<CudaMemory>();

    dcc.CalculateInternalCells(v_cl,ie_ghost,domain_proc,pbox,0.1,g);

    // Check

    auto & ic = dcc.getIcells();
    auto & dc = dcc.getDcells();

    ic.template deviceToHost<0>();
    dc.template deviceToHost<0>();

    openfpm::vector<unsigned int> icheck;
    openfpm::vector<unsigned int> dcheck;

    grid_sm<3,void> gr = dcc.getGrid();

    Box<3,size_t> b1({2,2,2},{6,6,6});

    size_t bc[3] = {NON_PERIODIC,NON_PERIODIC,NON_PERIODIC};

    grid_skin_iterator_bc<3> it_skin(gr,b1,b1,bc);

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

        if (p.get(0) == 6 && p.get(1) == 5 && p.get(2) == 5)
        {}
        else
        {icheck.add(gr.LinId(p));}

        ++it_skin;
    }

    Box<3,size_t> b3({7,4,4},{9,6,6});
    grid_key_dx_iterator_sub<3> it_sub(gr,b3.getKP1(),b3.getKP2());

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

        icheck.add(gr.LinId(p));

        ++it_sub;
    }

    // dcheck has only the internal part of the first domain + 6,5,5

    Box<3,size_t> b1i({3,3,3},{5,5,5});
    grid_key_dx_iterator_sub<3> it_sub2(gr,b1i.getKP1(),b1i.getKP2());

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

        dcheck.add(gr.LinId(p));

        ++it_sub2;
    }

    grid_key_dx<3> kk({6,5,5});
    dcheck.add(gr.LinId(kk));

    icheck.sort();
    dcheck.sort();

    for (size_t i = 0 ; i < icheck.size() ; i++)
    {BOOST_REQUIRE_EQUAL(icheck.template get<0>(i),ic.template get<0>(i));}

    for (size_t i = 0 ; i < dcheck.size() ; i++)
    {BOOST_REQUIRE_EQUAL(dcheck.template get<0>(i),dc.template get<0>(i));}

    #endif
}


BOOST_AUTO_TEST_SUITE_END()