gslpp_vector_complex.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef GSLPP_VECTOR_COMPLEX_H
#define GSLPP_VECTOR_COMPLEX_H
#include <iostream>
#ifndef __GSL_VECTOR_COMPLEX_DOUBLE_H__
#include <gsl/gsl_vector_complex_double.h>
#endif
#ifndef GSLPP_COMPLEX_H
#include "gslpp_complex.h"
#endif
#ifndef GSLPP_VECTOR_BASE_H
#include "gslpp_vector_base.h"
#endif
 
namespace gslpp
{
  template <>
  class vector<complex>
  {
      gsl_vector_complex *_vector;
 
    public:
      vector(const size_t& size, const complex& z);
      vector(const size_t& size, const double& a);
      vector(const vector<complex>& v);
      vector(const vector<double>& v);
      vector(const gsl_vector_complex& v);
      vector(const gsl_vector_complex* v);
      ~vector();
      const complex operator()(const size_t& i) const;
//      complex& operator()(const size_t& i);
      vector<complex>& operator=(const vector<complex>& v);
      vector<complex>& operator=(double a);
      void reset();
      void assign(const size_t& i, const complex& z);
      void assign(const size_t& i, const double& a);
      size_t size() const;
      double mod() const;
      vector<complex> conjugate() const;
      vector<double> real() const;
      vector<double> imag() const;
      gsl_vector_complex* as_gsl_type_ptr() const;
      gsl_vector_complex& as_gsl_type();
      const gsl_vector_complex& as_gsl_type() const;
      vector<complex> operator-() const;
      vector<complex> operator+(const vector<complex>& v) const;
      vector<complex> operator-(const vector<complex>& v) const;
      complex operator*(const vector<complex>& v) const;
//       vector<complex> operator^(const vector<complex>& v);
      vector<complex>& operator+=(const vector<complex>& v);
      vector<complex>& operator-=(const vector<complex>& v);
      vector<complex> operator+(const complex& z);
      vector<complex> operator-(const complex& z);
      vector<complex> operator*(const complex& z) const;
      vector<complex> operator/(const complex& z) const;
      vector<complex>& operator+=(const complex& z);
      vector<complex>& operator-=(const complex& z);
      vector<complex>& operator*=(const complex& z);
      vector<complex>& operator/=(const complex& z);
      vector<complex> operator+(const double& a);
      vector<complex> operator-(const double& a);
      vector<complex> operator*(const double& a) const;
      vector<complex> operator/(const double& a) const;
      vector<complex>& operator+=(const double& a);
      vector<complex>& operator-=(const double& a);
      vector<complex>& operator*=(const double& a);
      vector<complex>& operator/=(const double& a);
      bool operator==(const vector<complex>& a) const;
      bool operator!=(const vector<complex>& a) const
      {
          return(!(*this == a));
      }
 
      friend std::ostream& operator<<(std::ostream& output, const vector<complex>& v);
 
      friend vector<complex> operator+(const complex& z, vector<complex> v);
 
      friend vector<complex> operator-(const complex& z, vector<complex> v);
 
      friend vector<complex> operator*(const complex& z, vector<complex> v);
 
      friend vector<complex> operator+(const double& a, vector<complex> v);
 
      friend vector<complex> operator-(const double& a, vector<complex> v);
 
      friend vector<complex> operator*(const double& a, vector<complex> v);
  };
}
#endif