IsingModel.hpp

Go to the documentation of this file.

 
#ifndef __ISING_MODEL__
#define __ISING_MODEL__
 
#include "constants.hpp"
#include "data_type.hpp"
#include "typedefs.hpp"
#include "utils.hpp"
 
#include <armadillo>
#include <random>
#include <unordered_map>
 
#define INDEX(I, N) (I + N) % N
 
// Indeces for the neighbor matrix.
#define UP 0
#define LEFT 0
#define DOWN 1
#define RIGHT 1
 
class IsingModel {
private:
    friend class IsingModelTest;
    arma::Mat<int> lattice;
 
    arma::Mat<int> neighbors;
 
    std::unordered_map<int, double> energy_diff;
 
    double T;
 
    int L;
 
    int E;
 
    int M;
 
    void initialize_lattice();
 
    void initialize_neighbors();
 
    void initialize_energy_diff();
 
    void initialize_magnetization();
 
    void initialize_energy();
 
    IsingModel();
 
public:
    IsingModel(int L, double T);
 
    IsingModel(int L, double T, int val);
 
    data_t Metropolis();
 
    int get_E();
 
    int get_M();
};
 
#endif