piecewise_linear_random_variable.h

Go to the documentation of this file.

//
// Authors:
// Guillermo Casas gcasas@cimne-upc.edu
//
 
#ifndef PIECEWISE_LINEAR_RANDOM_VARIABLE_H
#define PIECEWISE_LINEAR_RANDOM_VARIABLE_H
 
 
// System includes
#include <string>
#include <iostream>
 
// Project includes
#include "random_variable.h"
#include <random>
 
namespace Kratos {
 
class KRATOS_API(DEM_APPLICATION) PiecewiseLinearRandomVariable: public RandomVariable {
 
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(PiecewiseLinearRandomVariable);
 
    PiecewiseLinearRandomVariable();
    PiecewiseLinearRandomVariable(const Parameters rParameters);
    PiecewiseLinearRandomVariable(const Parameters rParameters, const int seed);
 
    double Sample() override;
    double ProbabilityDensity(const double x);
    double GetMean() override;
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "PiecewiseLinearRandomVariable" ;
        return buffer.str();
 
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "PiecewiseLinearRandomVariable";
    }
 
 
    void PrintData(std::ostream& rOStream) const override
    {
 
    }
 
protected:
    void Check() override;
 
 
private:
 
    void CalculateTrapezoidProbabilitiesAndNormalize();
    std::size_t SampleTrapezoidChoice();
    double SampleWithinTrapezoid(const double H, const double B1, const double B2);
    double SampleWithinStandardTrapezoid(const double b);
    double SamplePositiveSlopingStandardTriangle();
    double SampleNegativeSlopingStandardTriangle();
 
    PiecewiseLinearRandomVariable & operator=(PiecewiseLinearRandomVariable const& rOther);
 
    double mRelativeClosenessTolerance = 0.0;
    std::vector<double> mPDFValues;
    std::vector<double> mPDFBreakpoints;
    std::mt19937 mRandomNumberGenerator;
    std::discrete_distribution<int> mTrapezoidsDiscreteDistribution;
 
}; // Class PiecewiseLinearRandomVariable
 
} // namespace Kratos
 
#endif // PIECEWISE_LINEAR_RANDOM_VARIABLE_H defined