TestLinearity.hpp

Go to the documentation of this file.

#ifndef MTTOOLBOX_TEST_LINEARITY_HPP
#define MTTOOLBOX_TEST_LINEARITY_HPP

#include <iostream>
#include <iomanip>
#include <MTToolBox/AbstractGenerator.hpp>
#include <MTToolBox/EquidistributionCalculatable.hpp>
#include <MTToolBox/util.hpp>

namespace MTToolBox {
    template<typename U, typename V = U>
    class TestLinearity {
    public:
        bool operator()(const EquidistributionCalculatable<U, V>& generator) {
            EquidistributionCalculatable<U, V> *g1 = generator.clone();
            EquidistributionCalculatable<U, V> *g2 = generator.clone();
            g1->seed(convert<U>(1234U));
            g2->seed(convert<U>(4321U));
            bool result = test1(*g1) && test2(*g1, *g2);
            delete g1;
            delete g2;
            return result;
        }

    private:
        bool test1(EquidistributionCalculatable<U, V>& g1) {
            using namespace std;
            EquidistributionCalculatable<U, V> *g2 = g1.clone();
            bool result = true;
            g2->add(g1);
            for (int i = 0; i < 100; i++) {
                if (!isZero(g2->generate())) {
                    result = false;
                    break;
                }
            }
            delete g2;
#if defined(DEBUG)
            if (result) {
                cout << "test1 passed" << endl;
            } else {
                cout << "test1 failed" << endl;
            }
#endif
            return result;
        }

        bool test2(EquidistributionCalculatable<U, V>& g1,
                   EquidistributionCalculatable<U, V>& g2) {
            using namespace std;
            EquidistributionCalculatable<U, V> *g3 = g2.clone();
            g3->add(g1);
            bool result = true;
            for (int i = 0; i < 100; i++) {
                U res1 = g1.generate();
                U res2 = g2.generate();
                U res3 = g3->generate();
                if (!((res1 ^ res2) == res3)) {
                    result = false;
#if defined(DEBUG)
                    cout << "i,res1,res2,res3 = " << dec << i << ","
                         << hex << res1 << ","
                         << res2 <<"," << res3 << ","
                         << (res1 ^ res2) << endl;
#endif
                    break;
                }
            }
            delete g3;
#if defined(DEBUG)
            if (result) {
                cout << "test2 passed" << endl;
            } else {
                cout << "test2 failed" << endl;
            }
#endif
            return result;
        }

    };
}
#endif // MTTOOLBOX_TEST_LINEARITY_HPP