IsingModel.hpp

Go to the documentation of this file.

 
#ifndef __ISING_MODEL__
#define __ISING_MODEL__
 
#include "data_type.hpp"
#include "utils.hpp"
 
#include <armadillo>
#include <cstdint>
#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;
 
    double energy_diff[17];
 
    double T;
 
    int L;
 
    int64_t E;
 
    int64_t M;
 
    std::mt19937 engine;
 
    void initialize_engine();
 
    void initialize_lattice();
 
    void initialize_lattice(int val);
 
    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();
};
 
#endif