Ensemble.hpp

Go to the documentation of this file.

#include "../ShortRange_mBody/Ensemble.cuh"
#include <Eigen/Dense>
 
template<typename Scalar_t = Complex_t<Real_t>>
class RandomMatrixEnsemble {
    private:
        using MatrixType = typename HamiltonianEnsemble::MatrixType;
        mutable ScalableOperator<MatrixType> m_sample;
 
    public:
        RandomMatrixEnsemble(HilbertSpace<int> const& locHSpace) {
            debug_constructor_printf(1);
            debug_print("\tlocHSpace.dim()=" << locHSpace.dim()
                                             << ", m_dLocMat.ptr()=" << m_dLocMat.ptr());
        }
 
        void discardSamples(int num) {
            if(num <= 0) return;
            auto& hEnsemble = HamiltonianEnsemble::get_instance(2);
            hEnsemble.discard(num);
        };
        void getNewSample(void) {
            auto& hEnsemble = HamiltonianEnsemble::get_instance(2);
            m_sample        = hEnsemble.sample();
        };
 
        template<typename EigenMatrix, class TotalSpace, typename Scalar_,
                 typename std::enable_if_t<on_GPU_v<EigenMatrix>>* = nullptr>
        void constructHamiltonian(EigenMatrix&                         mat,
                                  SubSpace<TotalSpace, Scalar_> const& subSpace) const {
            mat = m_sample.construct_globalOp(subSpace.totalSpace().sysSize());
        }
};