field_class.hpp

Go to the documentation of this file.

 
#ifndef FIELD_CLASS_INCLUDED
#define FIELD_CLASS_INCLUDED
 
#include <iostream>
#include <math.h>
#include <random>
#include "Eigen/LU"
 
using namespace Eigen;
using std::cout;
using std::endl;
using std::cerr;
 
const std::complex<double> IM(0, 1.0);
 
 
enum class Distribution {
  Zero_conf     ,
  Gauss_unit    ,
  Gauss_L       ,
  Gauss_MT_unit ,
  Gauss_MT_L    
};
 
 
 
class Field {
protected:
  int Li;
  int num_f;
  MatrixXcd field;
 
  void setgauss(const int num_col,
                const double mean = 0.0, const double dev = 1.0);
  void setgaussl(const int num_col,
                 const double mean = 0.0);
  void setgaussmt(const int num_col,
                  const double mean = 0.0, const double dev = 1.0);
  void setgaussmtl(const int num_col,
                   const double mean = 0.0);
  void setconf(const int num_col, const Distribution n,
               const double mean = 0.0, const double dev = 1.0);
  void setconf(const VectorXcd &v, const int num_field = 1);
  void setconf(const VectorXd &v,  const int num_field);
  void setconf(const MatrixXcd &m);
  void setconf(const MatrixXcd &m, const int num_field);
 
  MatrixXcd vector2matrix(const int n) const;
  MatrixXcd field2matrix(const int n) const;
  MatrixXcd dfield2matrix(const int n) const;
 
public:
  explicit Field() {};
  explicit Field(const int n1, const int n2) {Li = n1; num_f = n2;};
  Field(const int n1, const int n2, const Distribution n3,
        const double mean = 0.0, const double dev = 1.0) {
    Li = n1; num_f = n2; setconf(num_f, n3, mean, dev);
  };
  Field(const int n1, const int n2, const int n3, const Distribution n4,
        const double mean = 0.0, const double dev = 1.0) {
    Li = n1; num_f = n2; setconf(n3, n4, mean, dev);
  };
  Field(const VectorXcd &v, const int num_field = 1) {setconf(v, num_field);};
  Field(const VectorXd &v, const int num_field) {setconf(v, num_field);};
  Field(const MatrixXcd &m) {setconf(m);};
  Field(const MatrixXcd &m, const int num_field) {setconf(m, num_field);};
  virtual ~Field() {};
 
  static Field Zero(const int n1, const int n2) {
    return Field {n1, n2, Distribution::Zero_conf}; };
  static Field Zero(const int n1, const int n2, const int n3) {
    return Field {n1, n2, n3, Distribution::Zero_conf}; };
  static Field Gauss(const int n1, const int n2,
                     const double mean = 0.0, const double dev = 1.0) {
    return Field {n1, n2, Distribution::Gauss_unit, mean, dev}; };
  static Field Gauss(const int n1, const int n2, const int n3,
                     const double mean = 0.0, const double dev = 1.0) {
    return Field {n1, n2, n3, Distribution::Gauss_unit, mean, dev}; };
  static Field GaussL(const int n1, const int n2,
                      const double mean = 0.0) {
    return Field {n1, n2, Distribution::Gauss_L, mean}; };
  static Field GaussL(const int n1, const int n2, const int n3,
                      const double mean = 0.0) {
    return Field {n1, n2, n3, Distribution::Gauss_L, mean}; };
  static Field GaussMT(const int n1, const int n2,
                       const double mean = 0.0, const double dev = 1.0) {
    return Field {n1, n2, Distribution::Gauss_MT_unit, mean, dev}; };
  static Field GaussMT(const int n1, const int n2, const int n3,
                       const double mean = 0.0, const double dev = 1.0) {
    return Field {n1, n2, n3, Distribution::Gauss_MT_unit, mean, dev}; };
  static Field GaussMTL(const int n1, const int n2,
                        const double mean = 0.0) {
    return Field {n1, n2, Distribution::Gauss_MT_L, mean}; };
  static Field GaussMTL(const int n1, const int n2, const int n3,
                        const double mean = 0.0) {
    return Field {n1, n2, n3, Distribution::Gauss_MT_L, mean}; };
  static Field Vector(const VectorXcd &v, const int num_field = 1) {
    return Field {v, num_field}; };
  static Field Matrix(const MatrixXcd &m) {return Field {m};};
  static Field Matrix(const MatrixXcd &m, const int num_field) {
    return Field {m, num_field}; };
 
  MatrixXcd conf() const {return field;};
  VectorXcd conf(const int n) const {return field.col(n);};
  VectorXd conf_real(const int n) const {return field.col(n).real();};
  VectorXd conf_imag(const int n) const {return field.col(n).imag();};
 
  virtual bool operator==(const Field &f) {return (Li==f.Li) && (num_f==f.num_f);};
  virtual bool operator!=(const Field &f) {return !(*this==f);};
 
  virtual Field &operator*=(const double n) { field *= n; return *this; };
  virtual Field &operator/=(const double n) { field *= 1.0/n; return *this;};
 
  Field conv(const int n1, const int n2) const;
  Field conjconv(const int n) const;
  Field conv_pw(const int pw, const int n = 0) const;
 
  Field combine_with(const Field &f);
 
  virtual void show() const;
};
 
 
#endif