Ensemble_old.hpp

Go to the documentation of this file.

#pragma once
 
#include "Headers/RandomMatrices.hpp"
#include "Headers/HilbertSpace/HilbertSpace.hpp"
 
class RandomMatrixEnsemble {
    private:
        HilbertSpace const& m_locHSpace;
        GUEgenerator<>        m_GUE;
        mutable unsigned long count;
 
    public:
        RandomMatrixEnsemble(HilbertSpace const& locHSpace, int seed = 0, double stddev = 1.0)
            : m_locHSpace{locHSpace}, m_GUE(seed, stddev), count{0} {}
 
        void discardSamples(int num) {
            if(num <= 0) return;
            count += num;
            m_GUE.discard(num);
        };
        void getNewSample(void) {
            ++count;
            return;
        };
 
        template<typename EigenMatrix, typename Scalar_t>
        void constructHamiltonian(EigenMatrix& mat, SubHilbertSpace<Scalar_t> const& HSpace) {
            m_GUE(mat);
        };
};