xy.h

Go to the documentation of this file.

#ifndef XY_H
#define XY_H
 
#include <model.h>
 
 
class XY : public Model
{
 public:
  /******************************************************************************************/
  XY(double beta, double J, double B, unsigned int Nx, unsigned int Ny, unsigned int Nz, bool hc = true) : Model(beta, J, B, Nx, Ny, Nz)
  {
    if(hc) hot_start();
    else  cold_start();
  };
  
 
  //functions
  double compute_energy() final;
  void  hot_start() final;
  void cold_start() final;
 
 protected:
  double delta_energy(unsigned int) final;
  void flip(unsigned int, gsl_rng*) final;
  
  
 private:
  vector<double> mgrid; 
};
 
#endif