misc.h

Go to the documentation of this file.

 
#ifndef MISC_H
#define MISC_H
 
#include <complex>
#include <string>
#include <limits>
#include <itpp/itexports.h>
 
namespace std
{

template <class T>
std::ostream& operator<<(std::ostream &os, const std::complex<T> &x)
{
  os << x.real();
  ios::fmtflags saved_format = os.setf(ios::showpos);
  os << x.imag();
  os.setf(saved_format, ios::showpos);
  return os << 'i';
}

template <class T>
std::istream& operator>>(std::istream &is, std::complex<T> &x)
{
  T re, im;
  char c;
  is >> c;
  if (c == '(') {
    is >> re >> c;
    if (c == ',') {
      is >> im >> c;
      if (c == ')') {
        x = complex<T>(re, im);
      }
      else {
        is.setstate(ios_base::failbit);
      }
    }
    else if (c == ')') {
      x = complex<T>(re, T(0));
    }
    else {
      is.setstate(ios_base::failbit);
    }
  }
  else {
    is.putback(c);
    is >> re;
    if (!is.eof() && ((c = static_cast<char>(is.peek())) == '+' || c == '-')) {
      is >> im >> c;
      if (c == 'i') {
        x = complex<T>(re, im);
      }
      else {
        is.setstate(ios_base::failbit);
      }
    }
    else {
      x = complex<T>(re, T(0));
    }
  }
  return is;
}
 
} // namespace std
 

namespace itpp
{

const double pi = 3.14159265358979323846;

const double m_2pi = 2 * pi;

const double eps = std::numeric_limits<double>::epsilon();


inline bool is_int(double x)
{
  double dummy;
  return (modf(x, &dummy) == 0.0);
}

inline bool is_even(int x) { return ((x&1) == 0); }

ITPP_EXPORT std::string itpp_version();

ITPP_EXPORT bool is_bigendian();

inline bool check_big_endianness() { return is_bigendian(); }

 
} //namespace itpp
 
 
#endif // #ifndef MISC_H