sy_lib.hpp

Go to the documentation of this file.

/**********************************************************************           
    * sy_lib.hpp -- a library of useful processes in the SY MoC       *
    *                                                                 *
    * Authors:  Hosein Attarzadeh (shan2@kth.se)                      *
    *                                                                 *
    * Purpose: Enriching the SY library.                              *
    *                                                                 *
    * Usage:                                                          *
    *                                                                 *
    * License: BSD3                                                   *
    *******************************************************************/

#ifndef SYLIB_H
#define SYLIB_H

#include "sy_moc.hpp"

namespace ForSyDe
{
    
namespace SY
{


class gaussian : public source<double>
{
public:
    gaussian(sc_module_name name_,      
             const double& gaussVar,    
             const double& gaussMean    
            ) : source(name_, [=](abst_ext<double>& out1, const abst_ext<double>& inp)
                               {
                                   double rnd1,rnd2,G,Q,Q1,Q2;
                                   do
                                   {
                                       rnd1 = ((double)my_rand()) / ((double)2147483647) ;
                                       rnd2 = ((double)my_rand()) / ((double)2147483647) ;

                                       Q1 = 2.0 * rnd1 - 1.0 ;
                                       Q2 = 2.0 * rnd2 - 1.0 ;

                                       Q = Q1 * Q1 + Q2 * Q2 ;
                                   } while (Q > 1.0) ;

                                   G = gaussMean+sqrt(gaussVar)*(sqrt(-2.0*log(Q)/Q)*Q1);
                                   out1 = abst_ext<double>(G);
                               }, abst_ext<double>(0)) {}
private:
    // state variable:  
    bool shiftreg[64];  // boolean array for the LFSR random number generator
    void initialize()
    {
        long int seed = 11206341;
        for(int i=63; i>=0; i--) {  // the LFSR shiftregister is initialized with 0
          if(seed>=pow(2.,(double)i)) {
        shiftreg[i]=true;
        seed-=(long int)pow(2.,(double)i);
          }
          else shiftreg[i]=false;
        } 
    }  

  int my_rand() { 
    bool zw = (((shiftreg[59]==shiftreg[60])==shiftreg[62])==shiftreg[63]); // computing feedback
    for(int i=63; i>0; i--) {
      shiftreg[i] = shiftreg[i-1];  // shifting
    }
    shiftreg[0] = zw;           // writing the feedback bit
    double val = 0.;
    for(int i=0; i<31; i++) {
      if(shiftreg[2*i]) val+=pow(2.,(double)i); // extracting random number
    }
    return((int)floor(val)); 
  }
};


template <template <class> class OIf>
inline gaussian* make_gaussian(const std::string& pName,
    const double& gaussVar,    
    const double& gaussMean,   
    OIf<double>& outS
    )
{
    auto p = new gaussian(pName.c_str(), gaussVar, gaussMean);
    
    (*p).oport1(outS);
    
    return p;
}

}
}
#endif