gslpp_complex.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

#include <iostream>
#ifndef GSLPP_COMPLEX
#define GSLPP_COMPLEX
#ifndef __GSL_COMPLEX_H__
# include <gsl/gsl_complex.h>
#endif
#ifndef __GSL_COMPLEX_MATH_H__
# include <gsl/gsl_complex_math.h>
#endif

namespace gslpp
{
  class complex
  {
      gsl_complex _complex;
    public:
      static const complex& i();
      complex();
      
      complex(const double& real, const double& imag, bool polar = false);
      complex(const complex& z);
      complex(const double& a);
      complex(const gsl_complex* z);
      complex(const gsl_complex& z);
      virtual ~complex();
      bool is_real() const;
      bool is_imag() const;
      void assign(const double& real, const double& imag, bool polar);
      const double& real() const;
      const double& imag() const;
      double& real();
      double& imag();
      double arg() const;
      double abs() const;
      double abs2() const;
      double log_of_abs() const;
      complex& operator=(const complex& z);
      complex& operator=(const double& x);
      bool operator==(const complex& z1) const;
      bool operator!=(const complex& z1) const;
      complex operator-() const;
      complex operator+(const complex& z1) const;
      complex operator-(const complex& z1) const;
      complex operator*(const complex& z1) const;
      complex operator/(const complex& z1) const;
      complex& operator+=(const complex& z1);
      complex& operator-=(const complex& z1);
      complex& operator*=(const complex& z1);
      complex& operator/=(const complex& z1);
      complex operator+(const double& a) const;
      complex operator-(const double& a) const;
      complex operator*(const double& a) const;
      complex operator/(const double& a) const;
      complex& operator+=(const double& a);
      complex& operator-=(const double& a);
      complex& operator*=(const double& a);
      complex& operator/=(const double& a);
      complex conjugate() const;
      complex inverse() const;
      gsl_complex* as_gsl_type_ptr() const;
      gsl_complex& as_gsl_type();
      const gsl_complex& as_gsl_type() const;
      operator gsl_complex& ();
      operator const gsl_complex& () const;

      friend std::ostream& operator<<(std::ostream& output, const complex& z);
      friend complex operator+(const double& x1, const complex& z2);

      friend complex operator-(const double& x1, const complex& z2);

      friend complex operator*(const double& x1, const complex& z2);

      friend complex operator/(const double& x1,const complex& z2);
      friend complex exp(const complex& z);
      friend complex log(const complex& z);
      friend complex log10(const complex& z);
      friend complex log(const complex& z,
                          const complex& b);
      friend complex dilog(const complex& z);
      friend complex sqrt(const complex& z);
      friend complex pow(const complex& z1,
                          const complex& z2);
      friend complex pow(const complex& z,
                          const double  x);
      friend complex sin(const complex& z);
      friend complex cos(const complex& z);
      friend complex tan(const complex& z);
      friend complex sec(const complex& z);
      friend complex csc(const complex& z);
      friend complex cot(const complex& z);
      friend complex arcsin(const complex& z);
      friend complex arccos(const complex& z);
      friend complex arctan(const complex& z);
      friend complex arcsec(const complex& z);
      friend complex arccsc(const complex& z);
      friend complex arccot(const complex& z);
      friend complex sinh(const complex& z);
      friend complex cosh(const complex& z);
      friend complex tanh(const complex& z);
      friend complex sech(const complex& z);
      friend complex csch(const complex& z);
      friend complex coth(const complex& z);
      friend complex arcsinh(const complex& z);
      friend complex arccosh(const complex& z);
      friend complex arctanh(const complex& z);
      friend complex arcsech(const complex& z);
      friend complex arccsch(const complex& z);
      friend complex arccoth(const complex& z);
  };
}
#endif