vector_dist_operators_unit_tests.cpp

/*
 * vector_dist_operators_unit_tests.hpp
 *
 *  Created on: Jun 11, 2016
 *      Author: i-bird
 */

#ifndef OPENFPM_NUMERICS_SRC_OPERATORS_VECTOR_VECTOR_DIST_OPERATORS_UNIT_TESTS_HPP_
#define OPENFPM_NUMERICS_SRC_OPERATORS_VECTOR_VECTOR_DIST_OPERATORS_UNIT_TESTS_HPP_

#define BOOST_TEST_DYN_LINK
#include <boost/test/unit_test.hpp>

#include "tests/vector_dist_operators_tests_util.hpp"



BOOST_AUTO_TEST_SUITE( vector_dist_operators_test )

BOOST_AUTO_TEST_CASE( vector_dist_operators_test )
{
    if (create_vcluster().getProcessingUnits() > 3)
        return;

    Box<3,float> box({0.0,0.0,0.0},{1.0,1.0,1.0});

    // Boundary conditions
    size_t bc[3]={PERIODIC,PERIODIC,PERIODIC};

    // ghost
    Ghost<3,float> ghost(0.05);

    // vector type
    typedef vector_dist<3,float,aggregate<float,float,float,VectorS<3,float>,VectorS<3,float>,VectorS<3,float>>> vtype;

    vector_dist<3,float,aggregate<float,float,float,VectorS<3,float>,VectorS<3,float>,VectorS<3,float>>> vd(100,box,bc,ghost);

    check_all_expressions<comp_host>::check(vd);
}


void reset(vector_dist<3,float,aggregate<float,float,float,float,float,float,float,float>> & v1)
{
    auto it = v1.getDomainIterator();

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

        v1.template getProp<0>(p) = -1.0;
        v1.template getProp<1>(p) = -1.0;
        v1.template getProp<2>(p) = -1.0;
        v1.template getProp<3>(p) = -1.0;
        v1.template getProp<4>(p) = -1.0;
        v1.template getProp<5>(p) = -1.0;
        v1.template getProp<6>(p) = -1.0;
        v1.template getProp<7>(p) = -1.0;

        ++it;
    }
}

template<unsigned int i> void loop_check(vector_dist<3,float,aggregate<float,float,float,float,float,float,float,float>> & v1, float check)
{
    auto it = v1.getDomainIterator();
    bool ret = true;

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

        ret &= v1.template getProp<i>(p) == check;

        ++it;
    }

    BOOST_REQUIRE_EQUAL(ret,true);
}

void check(vector_dist<3,float,aggregate<float,float,float,float,float,float,float,float>> & v1, size_t i)
{
    float check = -1.0;

    if (0 < i) check = 0.0;
    else check = -1.0;
    loop_check<0>(v1,check);

    if (1 < i) check = 1.0;
    else check = -1.0;
    loop_check<1>(v1,check);

    if (2 < i) check = 2.0;
    else check = -1.0;
    loop_check<2>(v1,check);

    if (3 < i) check = 3.0;
    else check = -1.0;
    loop_check<3>(v1,check);

    if (4 < i) check = 4.0;
    else check = -1.0;
    loop_check<4>(v1,check);

    if (5 < i) check = 5.0;
    else check = -1.0;
    loop_check<5>(v1,check);

    if (6 < i) check = 6.0;
    else check = -1.0;
    loop_check<6>(v1,check);

    if (7 < i) check = 7.0;
    else check = -1.0;
    loop_check<7>(v1,check);
}

BOOST_AUTO_TEST_CASE( vector_dist_operators_assign_test )
{
    if (create_vcluster().getProcessingUnits() > 3)
        return;

    Box<3,float> box({0.0,0.0,0.0},{1.0,1.0,1.0});

    // Boundary conditions
    size_t bc[3]={PERIODIC,PERIODIC,PERIODIC};

    // ghost
    Ghost<3,float> ghost(0.05);

    vector_dist<3,float,aggregate<float,float,float,float,float,float,float,float>> vd(100,box,bc,ghost);

    auto v1 = getV<0>(vd);
    auto v2 = getV<1>(vd);
    auto v3 = getV<2>(vd);
    auto v4 = getV<3>(vd);
    auto v5 = getV<4>(vd);
    auto v6 = getV<5>(vd);
    auto v7 = getV<6>(vd);
    auto v8 = getV<7>(vd);

    auto v_pos = getV<POS_PROP>(vd);

    reset(vd);

    assign(v1,0.0,
           v2,1.0);

    check(vd,2);

    reset(vd);

    assign(v1,0.0,
           v2,1.0,
           v3,2.0);

    check(vd,3);

    reset(vd);

    assign(v1,0.0,
           v2,1.0,
           v3,2.0,
           v4,3.0);

    check(vd,4);

    reset(vd);

    assign(v1,0.0,
           v2,1.0,
           v3,2.0,
           v4,3.0,
           v5,4.0);

    check(vd,5);

    reset(vd);

    assign(v1,0.0,
           v2,1.0,
           v3,2.0,
           v4,3.0,
           v5,4.0,
           v6,5.0);

    check(vd,6);

    reset(vd);

    assign(v1,0.0,
           v2,1.0,
           v3,2.0,
           v4,3.0,
           v5,4.0,
           v6,5.0,
           v7,6.0);

    check(vd,7);

    reset(vd);

    assign(v1,0.0,
           v2,1.0,
           v3,2.0,
           v4,3.0,
           v5,4.0,
           v6,5.0,
           v7,6.0,
           v8,7.0);

    check(vd,8);

    // Self assign;
    assign(v_pos,v_pos,
           v_pos,v_pos);
}

BOOST_AUTO_TEST_SUITE_END()


#endif /* OPENFPM_NUMERICS_SRC_OPERATORS_VECTOR_VECTOR_DIST_OPERATORS_UNIT_TESTS_HPP_ */