grid_dist_id_util_tests.hpp

/*
 * grid_dist_id_util_tests.hpp
 *
 *  Created on: Oct 15, 2017
 *      Author: i-bird
 */
 
#ifndef SRC_GRID_TESTS_GRID_DIST_ID_UTIL_TESTS_HPP_
#define SRC_GRID_TESTS_GRID_DIST_ID_UTIL_TESTS_HPP_
 
 
static void print_test(std::string test, size_t sz)
{
    if (create_vcluster().getProcessUnitID() == 0)
        std::cout << test << " " << sz << "\n";
}
 
 
static void Test2D_core(grid_dist_id<2, float, aggregate<double>> & g_dist, const size_t (& sz)[2], size_t k)
{
    // check the consistency of the decomposition
    bool val = g_dist.getDecomposition().check_consistency();
    BOOST_REQUIRE_EQUAL(val,true);
 
    // Grid sm
    grid_sm<2,void> info(sz);
 
    // get the domain iterator
    size_t count = 0;
 
    auto dom = g_dist.getDomainIterator();
 
    while (dom.isNext())
    {
        auto key = dom.get();
        auto key_g = g_dist.getGKey(key);
 
        g_dist.template get<0>(key) = info.LinId(key_g);
 
        // Count the point
        count++;
 
        ++dom;
    }
 
 
    // Get the virtual cluster machine
    Vcluster<> & vcl = g_dist.getVC();
 
    // reduce
    vcl.sum(count);
    vcl.execute();
 
    // Check
    BOOST_REQUIRE_EQUAL(count,(size_t)k*k);
 
    auto dom2 = g_dist.getDomainIterator();
 
    grid_key_dx<2> start = dom2.getStart();
    grid_key_dx<2> stop = dom2.getStop();
 
    BOOST_REQUIRE_EQUAL((long int)stop.get(0),(long int)g_dist.size(0)-1);
    BOOST_REQUIRE_EQUAL((long int)stop.get(1),(long int)g_dist.size(1)-1);
 
    BOOST_REQUIRE_EQUAL(start.get(0),0);
    BOOST_REQUIRE_EQUAL(start.get(1),0);
 
    bool match = true;
 
    // check that the grid store the correct information
    while (dom2.isNext())
    {
        auto key = dom2.get();
        auto key_g = g_dist.getGKey(key);
 
        match &= (g_dist.template get<0>(key) == info.LinId(key_g))?true:false;
 
        ++dom2;
    }
 
    BOOST_REQUIRE_EQUAL(match,true);
 
    g_dist.template ghost_get<0>();
 
    // check that the communication is correctly completed
 
    auto domg = g_dist.getDomainGhostIterator();
 
    // check that the grid with the ghost part store the correct information
    while (domg.isNext())
    {
        auto key = domg.get();
        auto key_g = g_dist.getGKey(key);
 
        // In this case the boundary condition are non periodic
        if (g_dist.isInside(key_g))
        {
            match &= (g_dist.template get<0>(key) == info.LinId(key_g))?true:false;
        }
 
        ++domg;
    }
 
    BOOST_REQUIRE_EQUAL(match,true);
}
 
#endif /* SRC_GRID_TESTS_GRID_DIST_ID_UTIL_TESTS_HPP_ */