field_potential.hpp

Go to the documentation of this file.

 
#ifndef FIELD_POTENTIAL_INCLUDED
#define FIELD_POTENTIAL_INCLUDED
 
#include "field_class.hpp"
 
namespace MatrixComp {
  extern MatrixXcd ResizeMatrix (const MatrixXcd &m,
                                 const int rw, const int cl);
  extern MatrixXcd DeleteZeroLine (const MatrixXcd &m);
  extern double DiffLU (const MatrixXcd &m,
                        const PartialPivLU<MatrixXcd> &lu);
  extern int SignDet (const MatrixXcd &m, const int debug = 0);
};
 
 
 
class Potential :public Field {
// protected:
  void setpotential(const MatrixXcd &m);
  void setpotential(const MatrixXcd &m, const int num_field);
 
  MatrixXcd jacobian() const;
 
  MatrixXcd unifield_tilde_type() const;
  MatrixXcd unifield_jacobian() const;
  MatrixXcd dualfield_jacobian() const;
  MatrixXcd triplefield_jacobian1() const;
  MatrixXcd triplefield_jacobian2() const;
  MatrixXcd multifield_jacobian() const;
 
public:
  explicit Potential() {};
  explicit Potential(const int n1, const int n2,
                     const Distribution n3 = Distribution::Zero_conf)
    : Field(n1, n2) {
    int num_p = int(num_f*(num_f+1)/2); setconf(num_p, n3); };
  explicit Potential(const MatrixXcd &m) {
    setpotential(m);
  };
  explicit Potential(const MatrixXcd &m, const int num_field) {
    setpotential(m, num_field);
  };
  virtual ~Potential() {};
 
  virtual bool operator==(const Potential &p) {return (Li==p.Li) && (num_f == p.num_f);};
  virtual bool operator!=(const Potential &p) {return !(*this==p);};
 
  int sign_det(const int debug = 0) const;
  int sign_det4multifield (const int debug = 0) const;
  virtual void show() const;
};
 
 
 
class PotentialNR: public Potential {
  MatrixXcd field2;
  MatrixXcd field3;
 
public:
  explicit PotentialNR() {};
  explicit PotentialNR(const int n1, const int n2,
                       const Distribution n3 = Distribution::Zero_conf)
    : Potential(n1, n2, n3) {
    int N = (Li+1)*(Li+1);
    field2 = MatrixXcd::Zero(N, num_f);
    field3 = MatrixXcd::Zero(N, num_f);
  };
  explicit PotentialNR(const int num_field,
                       const MatrixXcd &m1,
                       const MatrixXcd &m2, const MatrixXcd &m3)
    : Potential(m1, num_field) {
    int N = m1.rows();
    try {
      if (m2.rows()!=N || m3.rows()!=N
          || m2.cols()!=num_f || m3.cols()!=num_f) throw 1;
      field2 = m2; field3 = m3;
    } catch (const int i) {
      field2 = MatrixComp::ResizeMatrix(m2, N, num_f);
      field3 = MatrixComp::ResizeMatrix(m3, N, num_f);
    }
  };
  ~PotentialNR() {};
 
  bool operator==(const PotentialNR &p) {return (Li==p.Li) && (num_f == p.num_f);};
  bool operator!=(const PotentialNR &p) {return !(*this==p);};
 
  VectorXd nr_pvec() const;
  MatrixXd nr_mat() const;
 
  VectorXcd nrerr_pvec() const;
 
  void show() const;
};
 
#endif