Kernels_unit_tests.hpp

/*
 * Kernels_unit_tests.hpp
 *
 *  Created on: Feb 17, 2016
 *      Author: i-bird
 */

#ifndef OPENFPM_NUMERICS_SRC_PSE_KERNELS_UNIT_TESTS_HPP_
#define OPENFPM_NUMERICS_SRC_PSE_KERNELS_UNIT_TESTS_HPP_

#include "PSE/Kernels_test_util.hpp"
#ifdef HAVE_LIBQUADMATH
#include <boost/multiprecision/float128.hpp>
#endif

BOOST_AUTO_TEST_SUITE( pse_kernels_unit_tests )

BOOST_AUTO_TEST_CASE( pse_ker )
{
    Vcluster & v_cl = create_vcluster();

    // This test is not made to run in parallel
    if (v_cl.getProcessingUnits() > 1)
        return;

    openfpm::vector<openfpm::vector<double>> y;
    openfpm::vector<openfpm::vector<double>> y_res;

    // Load the result of the test

#ifdef HAVE_LIBQUADMATH
    y_res.load("test/PSE_convergence");
#else
    y_res.load("test/PSE_convergence_osx");
#endif

    // Every time increase the number of particles by 2
    for (size_t i = 250 ; i <= 2097152000 ; i*=2)
    {
        y.add();

        PSEError err;


#ifdef HAVE_LIBQUADMATH

        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,2>>(i,2,err);
        y.last().add(err.linf_error);

        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,2>>(i,4,err);
        y.last().add(err.linf_error);
#endif

        PSE_test<double,Lap_PSE<1,double,2>>(i,2,err);
        y.last().add(err.linf_error);

        PSE_test<double,Lap_PSE<1,double,2>>(i,4,err);
        y.last().add(err.linf_error);

        PSE_test<float,Lap_PSE<1,float,2>>(i,2,err);
        y.last().add(err.linf_error);

        PSE_test<float,Lap_PSE<1,float,2>>(i,4,err);
        y.last().add(err.linf_error);


#ifdef HAVE_LIBQUADMATH

        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,4>>(i,2,err);
        y.last().add(err.linf_error);

        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,4>>(i,4,err);
        y.last().add(err.linf_error);



        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,6>>(i,2,err);
        y.last().add(err.linf_error);

        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,6>>(i,4,err);
        y.last().add(err.linf_error);




        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,8>>(i,8,err);
        y.last().add(err.linf_error);

        PSE_test<boost::multiprecision::float128,Lap_PSE<1,boost::multiprecision::float128,8>>(i,16,err);
        y.last().add(err.linf_error);

#endif
    }

    // Check the result
    for (size_t i = 0 ; i < y.size(); i++)
    {
        for (size_t j = 0 ; j < y.get(i).size(); j++)
        {
            double c1 = y.get(i).get(j);
            double c2 = y_res.get(i).get(j);

#ifdef HAVE_LIBQUADMATH

            // In divergent mode the system is too sensitive
            // to compiler/hardware differences disable them
            if (j != 4 && j != 5)
            {BOOST_REQUIRE_CLOSE(c1,c2,3.0);}

#else

            // In divergent mode the system is too sensitive
            // to compiler/hardware differences disable them
            if (j != 2 && j != 3)
            {BOOST_REQUIRE_CLOSE(c1,c2,3.0);}

#endif
        }
    }
}

BOOST_AUTO_TEST_SUITE_END()

#endif /* OPENFPM_NUMERICS_SRC_PSE_KERNELS_UNIT_TESTS_HPP_ */