StandardNormal.h

Go to the documentation of this file.

#ifndef SIMPT_DEFAULT_MOVE_GEN_STANDARD_NORMAL_H_INCLUDED
#define SIMPT_DEFAULT_MOVE_GEN_STANDARD_NORMAL_H_INCLUDED
/*
 * Copyright 2011-2016 Universiteit Antwerpen
 *
 * Licensed under the EUPL, Version 1.1 or  as soon they will be approved by
 * the European Commission - subsequent versions of the EUPL (the "Licence");
 * You may not use this work except in compliance with the Licence.
 * You may obtain a copy of the Licence at: http://ec.europa.eu/idabc/eupl5
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the Licence is distributed on an "AS IS" basis,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the Licence for the specific language governing
 * permissions and limitations under the Licence.
 */
#include "math/constants.h"
#include "math/RandomEngine.h"
#include "sim/CoreData.h"

#include <trng/normal_dist.hpp>
#include <trng/uniform01_dist.hpp>
#include <array>
#include <cmath>

namespace SimPT_Default {
namespace MoveGenerator {

using namespace SimPT_Sim;
using namespace SimPT_Sim::Util;

class StandardNormal
{
public:
        StandardNormal(const CoreData& cd)
        {
                assert( cd.Check() && "CoreData not ok.");
                Initialize(cd);
        }

        void Initialize(const CoreData& cd)
        {
                m_normal_generator = cd.m_random_engine->GetGenerator(trng::normal_dist<double>(0.0, 1.0));
                m_uniform_generator = cd.m_random_engine->GetGenerator(trng::uniform01_dist<double>());
        }

        std::array<double, 3> operator()()
        {
                // Angles range [0.0, pi], amplitudes [- very large, very large]
                const double angle  = m_uniform_generator() * pi();
                const double ampli  = m_normal_generator();
                return std::array<double, 3>{{ ampli * std::sin(angle), ampli * std::cos(angle), 0.0 }};
        }

private:
        std::function<double()>   m_normal_generator;   
        std::function<double()>   m_uniform_generator;  
};

} // namespace
} // namespace

#endif // end_of_include_guard