DrawParticles_unit_tests.hpp

/*
 * DrawParticles_unit_tests.hpp
 *
 *  Created on: Jan 3, 2017
 *      Author: i-bird
 */
 
#ifndef OPENFPM_NUMERICS_SRC_DRAW_DRAWPARTICLES_UNIT_TESTS_HPP_
#define OPENFPM_NUMERICS_SRC_DRAW_DRAWPARTICLES_UNIT_TESTS_HPP_
 
#include "PointIterator.hpp"
#include "PointIteratorSkin.hpp"
#include "DrawParticles.hpp"
 
BOOST_AUTO_TEST_SUITE( draw_particles )
 
BOOST_AUTO_TEST_CASE(point_iterator)
{
 
    size_t sz[] = {23,27,20};
 
    Box<3,double> domain({-1.2,0.5,-0.6},{1.0,3.1,1.3});
    Box<3,double> sub_domain({-0.15,0.75,0.15},{1.05,1.15,1.05});
    Box<3,double> sub_real({-0.1,0.80,0.2},{1.0,1.1,1.0});
 
    size_t sz_sub[] = {12,4,9};
 
    // Boundary conditions
    size_t bc[3]={NON_PERIODIC,NON_PERIODIC,NON_PERIODIC};
 
    // ghost, big enough to contain the interaction radius
    Ghost<3,double> ghost(0.01);
 
    vector_dist<3,double,aggregate<double>> vd(0,domain,bc,ghost);
 
    auto p = DrawParticles::DrawBox(vd,sz,domain,sub_domain);
 
    size_t cnt = 0;
 
    bool good = true;
    while (p.isNext())
    {
        vd.add();
 
        vd.getLastPos()[0] = p.get().get(0);
        vd.getLastPos()[1] = p.get().get(1);
        vd.getLastPos()[2] = p.get().get(2);
 
        good &= sub_domain.isInside(p.get());
 
        cnt++;
 
        ++p;
    }
 
 
    Vcluster<> & v_cl = create_vcluster();
 
    v_cl.sum(cnt);
    v_cl.execute();
 
    BOOST_REQUIRE_EQUAL(cnt,sz_sub[0]*sz_sub[1]*sz_sub[2]);
    BOOST_REQUIRE_EQUAL(good,true);
 
    Box<3,double> marg = p.getBoxMargins();
 
    BOOST_REQUIRE_CLOSE(marg.getLow(0), sub_real.getLow(0) ,0.0001);
    BOOST_REQUIRE_CLOSE(marg.getLow(1), sub_real.getLow(1) ,0.0001);
    BOOST_REQUIRE_CLOSE(marg.getLow(2), sub_real.getLow(2) ,0.0001);
 
    BOOST_REQUIRE_CLOSE(marg.getHigh(0), sub_real.getHigh(0) ,0.0001);
    BOOST_REQUIRE_CLOSE(marg.getHigh(1), sub_real.getHigh(1) ,0.0001);
    BOOST_REQUIRE_CLOSE(marg.getHigh(2), sub_real.getHigh(2) ,0.0001);
}
 
BOOST_AUTO_TEST_CASE(point_iterator_skin)
{
    size_t sz[] = {23,27,20};
 
    Box<3,double> domain({-1.2,0.5,-0.6},{1.0,3.1,1.3});
    Box<3,double> sub_domainA({-0.15,0.75,0.15},{1.05,1.15,1.05});
    Box<3,double> sub_domainB({-0.25,0.65,0.05},{0.95,1.05,1.05});
 
    // Boundary conditions
    size_t bc[3]={NON_PERIODIC,NON_PERIODIC,NON_PERIODIC};
 
    // ghost, big enough to contain the interaction radius
    Ghost<3,double> ghost(0.01);
 
    vector_dist<3,double,aggregate<double>> vd(0,domain,bc,ghost);
 
    auto p = DrawParticles::DrawSkin(vd,sz,domain,sub_domainA, sub_domainB);
 
    size_t cnt = 0;
 
    bool good = true;
    while (p.isNext())
    {
        vd.add();
 
        vd.getLastPos()[0] = p.get().get(0);
        vd.getLastPos()[1] = p.get().get(1);
        vd.getLastPos()[2] = p.get().get(2);
 
        if (!((p.get().get(0) > -0.21 && p.get().get(0) < -0.19) || (p.get().get(1) > 0.69 && p.get().get(1) < 0.71) || (p.get().get(2) > 0.09 && p.get().get(2) < 0.11)) )
            good &= false;
 
        if (sub_domainA.isInside(p.get()))
            good &= false;
 
        cnt++;
 
        ++p;
    }
 
    Vcluster<> & v_cl = create_vcluster();
 
    v_cl.sum(cnt);
    v_cl.execute();
 
    BOOST_REQUIRE_EQUAL(good,true);
}
 
BOOST_AUTO_TEST_SUITE_END()
 
#endif /* OPENFPM_NUMERICS_SRC_DRAW_DRAWPARTICLES_UNIT_TESTS_HPP_ */