gslpp_vector_double.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef GSLPP_VECTOR_DOUBLE_H
#define GSLPP_VECTOR_DOUBLE_H
#include <iostream>
#ifndef __GSL_VECTOR_DOUBLE_H__
# include <gsl/gsl_vector_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<double>
  {
    gsl_vector *_vector;
 
    public:
      vector(const size_t& size);
      vector(const size_t& size, const double& a);
      vector(const vector<double>& v);
      vector(const gsl_vector& v);
      vector(const gsl_vector* v);
      ~vector();
      const double& operator()(const size_t& i) const;
      double& operator()(const size_t& i);
      void reset();
      vector<double>& operator=(const vector<double>& v);
      size_t size() const;
      double mod() const;
      double max() const;
      double min() const;
      gsl_vector* as_gsl_type_ptr() const;
      gsl_vector& as_gsl_type();
      const gsl_vector& as_gsl_type() const;
      vector<double> operator-() const;
      vector<double> operator+(const vector<double>& v) const;
      vector<complex> operator+(const vector<complex>& v) const;
      vector<double> operator-(const vector<double>& v) const;
      vector<complex> operator-(const vector<complex>& v) const;
      double operator*(const vector<double>& v) const;
      complex operator*(const vector<complex>& v) const;
//       vector<double> operator^(const vector<double>& v);
      vector<double>& operator+=(const vector<double>& v);
      vector<double>& operator-=(const vector<double>& v);
      vector<double> operator+(const double& a) const;
      vector<double> operator-(const double& a) const;
      vector<double> operator*(const double& a) const;
      vector<double> operator/(const double& a) const;
      vector<double>& operator+=(const double& a);
      vector<double>& operator-=(const double& a);
      vector<double>& operator*=(const double& a);
      vector<double>& operator/=(const double& a);
      vector<complex> operator+(const complex& z) const;
      vector<complex> operator-(const complex& z) const;
      vector<complex> operator*(const complex& z) const;
      vector<complex> operator/(const complex& z) const;
 
      bool operator==(const vector<double>& a) const;
      bool operator!=(const vector<double>& a) const
      {
          return(!(*this == a));
      }
 
      friend std::ostream& operator<<(std::ostream& output, const vector<double>& v);
      friend vector<double> operator+(const double& a, vector<double> v);
 
      friend vector<double> operator-(const double& a, vector<double>  v);
 
      friend vector<double> operator*(const double& a, vector<double> v);
 
      friend vector<complex> operator+(const complex& z, vector<double> v);
 
      friend vector<complex> operator-(const complex& z, vector<double> v);
 
      friend vector<complex> operator*(const complex& z, vector<double> v);
  };
}
 
#endif