model.h

Go to the documentation of this file.

#pragma once
 
#include <vector>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <iostream>
#include <random>
 
#ifdef _OPENMP
  #include <omp.h>
#endif
 
#include "grid/grid.h"
 
#include <cuda_runtime.h>
#include <curand.h>
#include <curand_kernel.h>
 
#define CUDA_CALL(x) do { if((x)!=cudaSuccess) { \
    printf("Error at %s:%d\n",__FILE__,__LINE__);\
    return;}} while(0)
 
#define CURAND_CALL(x) do { if((x)!=CURAND_STATUS_SUCCESS) { \
    printf("Error at %s:%d\n",__FILE__,__LINE__);\
    return;}} while(0)
 
 
namespace HB
{
  class Model
  {
   public:
    Model() = default;
    /****************************************************************************/
    Model(float J, float3 B) : mJ(J), mB(B) {};
 
    virtual ~Model() = default;
    
    /****************************************************************************/
    virtual float calcEnergy(const Grid<short int> &grid) = 0;
 
    /****************************************************************************/
    virtual float calcMagnetization(const Grid<short int> &grid) = 0;
    
    /****************************************************************************/
    float3 getBField() const {return mB;}
 
    /****************************************************************************/
    float getCouplingConstant() const {return mJ;}
    
    /****************************************************************************/
    void setBField(const float3 B) {mB = B;}
 
    /****************************************************************************/
    void setCouplingConstant(const float J) {mJ = J;}
 
    /****************************************************************************/
    virtual float calcEnergy(const Grid<short int> &grid, const dim3 index) = 0;
 
    /****************************************************************************/
    virtual void flip(const dim3 index, Grid<short int> &grid) = 0;
 
    virtual void simulate() const = 0;
    
    float3 mB; 
    float mJ; 
};
}