utilities.h

Go to the documentation of this file.

/***************************************************************************/
#ifndef AURESERVOIR_UTILITIES_H__
#define AURESERVOIR_UTILITIES_H__

// external includes
#include <flens/flens.h>
#include <string>
#include <ctime>
#include <iostream>
#include <iomanip>
#include <sstream>
#include <cstdlib>
// #include <complex>

#include "auexcept.h"
#include "denormal.h"

namespace aureservoir
{

using std::string;
using std::complex;
using flens::_;

template<typename T = float>
struct SPMatrix
{
  typedef flens::SparseGeMatrix<flens::CRS<T> > Type;
};

template<typename T = float>
struct DEMatrix
{
  typedef flens::GeMatrix<flens::FullStorage<T, flens::ColMajor> > Type;
};

template<typename T = float>
struct DEVector
{
  typedef flens::DenseVector<flens::Array<T> > Type;
};

template<typename T = float>
struct CDEVector
{
  typedef flens::DenseVector<flens::Array< complex<T> > > Type;
};


template <typename T=float>
class Rand
{
 public:

  static void initSeed()
  { srand(time(0)); }

  static T uniform(float min=-1, float max=1)
  {
   T tmp = std::rand() / (T(RAND_MAX)+1); // between [0|1)
   return tmp*(max-min) + min;
  }

  static void uniform(typename DEVector<T>::Type &vec, float min=-1, float max=1)
  {
    for(int i=0; i<vec.length(); ++i)
      vec.data()[i] = std::rand() / (T(RAND_MAX)+1); // between [0|1)

    vec *= (max-min);
    vec += min;
  }
};

// template<>
// static dcplx Rand<dcplx>::uniform(float min, float max)
// {
//   double real = std::rand() / (double(RAND_MAX)+1);
//   double imag = std::rand() / (double(RAND_MAX)+1);
//   dcplx tmp(real,imag);
//   return tmp*(dcplx(max)-dcplx(min)) + dcplx(min);
// }


inline double stringToDouble(const string& s)
  throw(AUExcept)
{
  std::istringstream is(s);
  double val;
  char c;

  // throw error if conversion is wrong or if we still have values
  if( !(is >> val) || is.get(c) )
    throw AUExcept("stringToDouble: could not convert parameter: " + s);

  return val;
}

inline int stringToInt(const string& s)
  throw(AUExcept)
{
  std::istringstream is(s);
  int val;
  char c;

  // throw error if conversion is wrong or if we still have values
  if( !(is >> val) || is.get(c) )
    throw AUExcept("stringToInt: could not convert parameter: " + s);

  return val;
}

} // end of namespace aureservoir

#endif // AURESERVOIR_UTILITIES_H__

---