vector_dist_operators_unit_tests.cu

#include "config.h"
#define BOOST_TEST_DYN_LINK
#include <boost/test/unit_test.hpp>
 
#include "VCluster/VCluster.hpp"
#include "Vector/vector_dist.hpp"
#include "Operators/Vector/vector_dist_operators.hpp"
#include "Operators/Vector/tests/vector_dist_operators_tests_util.hpp"
 
 
BOOST_AUTO_TEST_SUITE( vector_dist_operators_test_gpu )
 
BOOST_AUTO_TEST_CASE(vector_dist_operators_list_ker_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.2);
 
    vector_dist_gpu<3,float,aggregate<float,float,float,VectorS<3,float>,VectorS<3,float>,VectorS<3,float>>> vd(127*create_vcluster().size(),box,bc,ghost);
 
    typedef decltype(vd.toKernel()) vector_dist_kernel;
 
    vector_dist_ker_list<vector_dist_kernel> & vdkl = vd.private_get_vector_dist_ker_list();
 
    BOOST_REQUIRE_EQUAL(vdkl.n_entry(),0);
 
    {
        auto vA = getV<A,comp_dev>(vd);
 
        BOOST_REQUIRE_EQUAL(vdkl.n_entry(),1);
 
        {
            auto vB = getV<B,comp_dev>(vd);
            auto vC = getV<C,comp_dev>(vd);
 
            auto vVA = getV<VA,comp_dev>(vd);
            auto vVB = getV<VB,comp_dev>(vd);
            auto vVC = getV<VC,comp_dev>(vd);
 
            BOOST_REQUIRE_EQUAL(vdkl.n_entry(),6);
 
            // fill vd with some value
            fill_values<comp_dev>(vd);
 
            vd.map(RUN_ON_DEVICE);
 
            // check that the entry has been updated
            BOOST_REQUIRE_EQUAL(vdkl.check(vd.toKernel()),true);
 
            vd.template ghost_get<0,1,2,3,4,5>(RUN_ON_DEVICE);
 
            // check that the entry has been updated
            BOOST_REQUIRE_EQUAL(vdkl.check(vd.toKernel()),true);
        }
 
        BOOST_REQUIRE_EQUAL(vdkl.n_entry(),1);
    }
 
    BOOST_REQUIRE_EQUAL(vdkl.n_entry(),0);
}
 
BOOST_AUTO_TEST_CASE( vector_dist_operators_gpu_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_gpu<3,float,aggregate<float,float,float,VectorS<3,float>,VectorS<3,float>,VectorS<3,float>>> vtype;
 
    vector_dist_gpu<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);
    check_all_expressions<comp_dev>::check(vd);
}
 
BOOST_AUTO_TEST_SUITE_END()