vector_dist_util_unit_tests.hpp

/*
 * vector_dist_util_unit_tests.hpp
 *
 *  Created on: Feb 14, 2018
 *      Author: i-bird
 */

#ifndef SRC_VECTOR_TESTS_VECTOR_DIST_UTIL_UNIT_TESTS_HPP_
#define SRC_VECTOR_TESTS_VECTOR_DIST_UTIL_UNIT_TESTS_HPP_


template<unsigned int dim,typename vector_dist> inline void count_local_n_local(vector_dist & vd, vector_dist_iterator & it, size_t (& bc)[dim] , Box<dim,float> & box, Box<dim,float> & dom_ext, size_t & l_cnt, size_t & nl_cnt, size_t & n_out)
{
    const CartDecomposition<dim,float> & ct = vd.getDecomposition();

    while (it.isNext())
    {
        auto key = it.get();
        // Check if it is in the domain
        if (box.isInsideNP(vd.getPos(key)) == true)
        {
            // Check if local
            if (ct.isLocalBC(vd.getPos(key),bc) == true)
                l_cnt++;
            else
                nl_cnt++;
        }
        else
        {
            nl_cnt++;
        }

        Point<dim,float> xp = vd.getPos(key);

        // Check that all particles are inside the Domain + Ghost part
        if (dom_ext.isInside(xp) == false)
                n_out++;

        ++it;
    }
}



#endif /* SRC_VECTOR_TESTS_VECTOR_DIST_UTIL_UNIT_TESTS_HPP_ */