randomT_test.cpp

Go to the documentation of this file.

/** \file randomT_test.cpp
 */
#include "../../catch2/catch.hpp"
 
#include <Eigen/Dense>
 
#define MX_NO_ERROR_REPORTS
 
#include "../../../include/math/randomT.hpp"
 
/** Verify compilation and basic operation of randomT with std::distributions.
 * Basic tests include verification that seeding works.  Note that there is the very slight
 * possibility that normal operation could return two consecutive variates with the same value
 * which would cause some of the below tests to fail.  Will probably never happen . . .
 *
 * \anchor tests_math_randomT_basic
 */
SCENARIO( "Verify compilation and basic operation of randomT with std::distributions", "[math::randomT]" )
{
    GIVEN( "a uniform distribution is desired" )
    {
        WHEN( "two double random numbers from same generator requested" )
        {
            mx::math::randomT<double, std::mt19937_64, std::uniform_real_distribution<double>> uniDist;
 
            double r1 = uniDist;
            double r2 = uniDist;
 
            REQUIRE( r1 != r2 );
        }
 
        WHEN( "two double random numbers requested from different generators with same seed" )
        {
            mx::math::randomT<double, std::mt19937_64, std::uniform_real_distribution<double>> uniDist1( false );
            mx::math::randomT<double, std::mt19937_64, std::uniform_real_distribution<double>> uniDist2( false );
 
            uniDist1.seed( 10 );
            uniDist2.seed( 10 );
            double r1 = uniDist1;
            double r2 = uniDist2;
 
            REQUIRE( r1 == r2 );
        }
 
        WHEN( "two double random numbers requested from different generators with different seed" )
        {
            mx::math::randomT<double, std::mt19937_64, std::uniform_real_distribution<double>> uniDist1( false );
            mx::math::randomT<double, std::mt19937_64, std::uniform_real_distribution<double>> uniDist2( false );
 
            uniDist1.seed( 10 );
            uniDist2.seed( 11 );
            double r1 = uniDist1;
            double r2 = uniDist2;
 
            REQUIRE( r1 != r2 );
        }
    }
 
    GIVEN( "a normal distribution is desired" )
    {
        WHEN( "two double random numbers from same generator requested" )
        {
            mx::math::randomT<double, std::mt19937_64, std::normal_distribution<double>> normDist;
 
            double r1 = normDist;
            double r2 = normDist;
 
            REQUIRE( r1 != r2 );
        }
 
        WHEN( "two double random numbers requested from different generators with same seed" )
        {
            mx::math::randomT<double, std::mt19937_64, std::normal_distribution<double>> normDist1( false );
            mx::math::randomT<double, std::mt19937_64, std::normal_distribution<double>> normDist2( false );
 
            normDist1.seed( 10 );
            normDist2.seed( 10 );
            double r1 = normDist1;
            double r2 = normDist2;
 
            REQUIRE( r1 == r2 );
        }
 
        WHEN( "two double random numbers requested from different generators with different seed" )
        {
            mx::math::randomT<double, std::mt19937_64, std::normal_distribution<double>> normDist1( false );
            mx::math::randomT<double, std::mt19937_64, std::normal_distribution<double>> normDist2( false );
 
            normDist1.seed( 10 );
            normDist2.seed( 11 );
            double r1 = normDist1;
            double r2 = normDist2;
 
            REQUIRE( r1 != r2 );
        }
    }
}
 
/** Verify compilation and basic operation of randomT with the Lapace distribution
 * Basic tests include verification that seeding works.  Note that there is the very slight
 * possibility that normal operation could return two consecutive variates with the same value
 * which would cause some of the below tests to fail.  Will probably never happen . . .
 *
 * \anchor tests_math_randomT_basic_laplace
 */
SCENARIO( "Verify compilation and basic operation of randomT with the Lapace distribution",
          "[math::laplace_distribution]" )
{
    GIVEN( "a laplace distribution is desired" )
    {
        WHEN( "two double random numbers from same generator requested" )
        {
            mx::math::randomT<double, std::mt19937_64, mx::math::laplace_distribution<double>> lapDist;
 
            double r1 = lapDist;
            double r2 = lapDist;
 
            REQUIRE( r1 != r2 );
        }
 
        WHEN( "two double random numbers requested from different generators with same seed" )
        {
            mx::math::randomT<double, std::mt19937_64, mx::math::laplace_distribution<double>> lapDist1( false );
            mx::math::randomT<double, std::mt19937_64, mx::math::laplace_distribution<double>> lapDist2( false );
 
            lapDist1.seed( 10 );
            lapDist2.seed( 10 );
            double r1 = lapDist1;
            double r2 = lapDist2;
 
            REQUIRE( r1 == r2 );
        }
 
        WHEN( "two double random numbers requested from different generators with different seed" )
        {
            mx::math::randomT<double, std::mt19937_64, mx::math::laplace_distribution<double>> lapDist1( false );
            mx::math::randomT<double, std::mt19937_64, mx::math::laplace_distribution<double>> lapDist2( false );
 
            lapDist1.seed( 10 );
            lapDist2.seed( 11 );
            double r1 = lapDist1;
            double r2 = lapDist2;
 
            REQUIRE( r1 != r2 );
        }
    }
}