CommonInline.cxx

// Copyright (C) 2001-2012 Vivien Mallet
//
// This file is part of the linear-algebra library Seldon,
// http://seldon.sourceforge.net/.
//
// Seldon is free software; you can redistribute it and/or modify it under the
// terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 2.1 of the License, or (at your option)
// any later version.
//
// Seldon is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
// FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
// more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with Seldon. If not, see http://www.gnu.org/licenses/.
 
 
#ifndef SELDON_FILE_COMMON_INLINE_CXX
 
 
namespace std
{
  template<class T>
  inline T conjugate(const T& x)
  {
    return x;
  }
 
  template<class T>
  inline complex<T> conjugate(const complex<T>& x)
  {
    return conj(x);
  }
 
  template<class T>
  inline T realpart(const T& x)
  {
    return x;
  }
 
  template<class T>
  inline T realpart(const complex<T>& x)
  {
    return real(x);
  }
}
 
namespace Seldon
{
 
  template <class T>
  inline void PrintArray(T* v, int lgth)
  {
    for (int k = 0; k < lgth - 1; k++)
      std::cout << v[k] << " | ";
    std::cout << v[lgth - 1] << std::endl;
  }
 
 
  inline Str::Str()
  {
  }
 
 
 
  inline Str::Str(const Str& s)
  {
    string s_input = s;
    output_ << s_input;
  }
 
 
  inline Str::operator std::string() const
  {
    return output_.str();
  }
 
 
 
  template <class T>
  inline Str& Str::operator << (const T& input)
  {
    output_ << input;
    return *this;
  }
 
 
 
  template <class T>
  inline Str operator + (const Str& s, const T& input)
  {
    string s_input = s;
    Str output;
    output << s_input << input;
    return output;
  }
 
 
  inline ostream& operator << (ostream& out, Str& in)
  {
    string output = in;
    out << output;
    return out;
  }
 
 
  inline ostream& operator << (ostream& out, Str in)
  {
    string output = in;
    out << output;
    return out;
  }
 
 
 
  template<typename T>
  inline std::string to_str(const T& input)
  {
    std::ostringstream output;
    output.precision(cout.precision());
    output << input;
    return output.str();
  }
 
 
 
  template <class T>
  inline void to_num(std::string s, T& num)
  {
    std::istringstream str(s);
    str >> num;
  }
 
 
 
  template <class T>
  inline T to_num(std::string s)
  {
    T num;
    std::istringstream str(s);
    str >> num;
    return num;
  }
 
 
 
  template <class T>
  inline void SetComplexZero(T& number)
  {
    number = T(0);
  }
 
 
 
  template <class T>
  inline void SetComplexZero(std::complex<T>& number)
  {
    number = std::complex<T>(T(0), T(0));
  }
 
 
 
  template <class T>
  inline void SetComplexOne(T& number)
  {
    number = T(1);
  }
 
 
 
  template <class T>
  inline void SetComplexOne(std::complex<T>& number)
  {
    number = complex<T>(T(1), T(0));
  }
 
 
 
  template <class T>
  inline void SetComplexReal(int n, T& number)
  {
    number = T(n);
  }
 
 
 
  template <class T>
  inline void SetComplexReal(int n, std::complex<T>& number)
  {
    number = std::complex<T>(n, 0);
  }
 
 
 
  template <class T>
  inline void SetComplexReal(bool n, std::complex<T>& number)
  {
    number = std::complex<T>(n, 0);
  }
 
 
  inline void SetComplexReal(int x, std::complex<int>& number)
  {
    number = std::complex<int>(x, 0);
  }
 
 
 
  template <class T>
  inline void SetComplexReal(const T& x, std::complex<T>& number)
  {
    number = std::complex<T>(x, 0);
  }
 
 
 
  template <class T>
  inline void SetComplexReal(const std::complex<T>& x, T& number)
  {
    number = std::realpart(x);
  }
 
 
 
  template <class T0, class T1>
  inline void SetComplexReal(const T0& x, T1& number)
  {
    number = x;
  }
 
 
  template <class T>
  inline bool IsComplexNumber(const T& number)
  {
    return false;
  }
 
 
  template <class T>
  inline bool IsComplexNumber(const std::complex<T>& number)
  {
    return true;
  }
 
 
  template<class T>
  inline T ComplexAbs(const T& val)
  {
    return abs(val);
  }
 
 
  template<class T>
  inline T ComplexAbs(const std::complex<T>& val)
  {
#if defined(SELDON_WITH_BLAS) && !defined(SELDON_WITH_LAPACK)
    // we choose Blas convention
    return abs(real(val)) + abs(imag(val));
#else
    // otherwise usual modulus
    return abs(val);
#endif
  }
 
 
  template<class T>
  inline T absSquare(const std::complex<T>& z)
  {
    // more optimal than real(z * conj(z))
    return real(z)*real(z) + imag(z)*imag(z);
  }
 
 
  template<class T>
  inline T absSquare(const T& z)
  {
    return z*z;
  }
 
}  // namespace Seldon.
 
#define SELDON_FILE_COMMON_INLINE_CXX
#endif