diagram.h

Go to the documentation of this file.

 
#pragma once
 
#include <list>
#include <random>
#include <chrono>
 
#define EPSILON 1e-10  //theshold for floating point comparison
 
class Diagram_core
{
 
    protected:
 
    double _beta;                 
    int _s0;                      
    double _H;                    
    double _GAMMA;                
    std::list<double> _vertices;  
 
 
    void assert_parameters_validity(double beta, int s0, double H, double GAMMA, std::list<double> vertices) const;
 
 
    public:
 
    Diagram_core(double beta, int s0, double H, double GAMMA, std::list<double> vertices=std::list<double>() );
 
    bool operator==(const Diagram_core & other) const;
 
 
    bool operator!=(const Diagram_core & other) const;
 
    double operator/(const Diagram_core & other) const;
 
    double sum_deltatau() const;
 
 
    double value() const;
 
    size_t order() const; 
 
    double get_beta() const;
 
    int get_s0() const;
 
    double get_H() const;
 
    double get_GAMMA() const;
 
    std::list<double> get_vertices() const;
 
 
    double acceptance_rate_add(double tau1, double tau2, double tau2max, double new_segment_spin) const;
 
    double acceptance_rate_remove(double tau1, double tau2, double tau2max, double segment_toberemoved_spin) const;
 
    double acceptance_rate_flip() const;
 
    bool attempt_add_segment(double RN1, double RN2, double RNacc);
 
    bool attempt_remove_segment(double RN1, double RNacc);
 
    bool attempt_spin_flip(double RNacc);
 
 
};
 
 
class Diagram: public Diagram_core
{
 
    private:
        std::uniform_real_distribution<double> _uniform_dist; 
        std::mt19937 _mt_generator;                           
 
 
    public:
 
    Diagram(double beta, int s0, double H, double GAMMA, 
        std::list<double> _vertices=std::list<double>(), 
        unsigned int seed = std::chrono::system_clock::now().time_since_epoch().count());
 
    using Diagram_core::operator/ ;
    //using Diagram_core::operator== ;
 
 
    bool attempt_add_segment();
 
    bool attempt_remove_segment(); 
 
    bool attempt_spin_flip();
 
    void reset_diagram(double beta, int s0, double H, double GAMMA, 
        std::list<double> vertices=std::list<double>(), 
        unsigned int seed = std::chrono::system_clock::now().time_since_epoch().count());
 
};
 
 
bool lists_are_float_equal(const std::list<double>& list1, const std::list<double>& list2, double epsilon);