MicroCanonicalAverage Class Reference

Calculate the microcanonical averages with respect to a given sorted vector 'eigVal'. More...

#include <MicroCanonicalAverage.hpp>

Collaboration diagram for MicroCanonicalAverage:

Public Member Functions
template<typename InVector_t >
__host__ __device__	MicroCanonicalAverage (InVector_t const &eigVal, double const shellWidth)
__host__ void	status (void) const
	Show the current status of the instance.
__host__ __device__ int	dim () const
__host__ __device__ double	shellWidth () const
__host__ __device__ int	dimShell (int index) const
__host__ __device__ Vector_t< int > &	dimShell ()
__host__ __device__ Vector_t< int > const &	dimShell () const
__host__ __device__ Vector_t< double > &	eigVal ()
__host__ __device__ Vector_t< double > const &	eigVal () const
template<class VectorType >
__host__ __device__ void	operator() (VectorType &res, VectorType const &expVal) const

Private Types
template<class T >
using	Vector_t = Eigen::VectorX< T >

Private Attributes
int	m_dim
double	m_shellWidth
Vector_t< double >	m_eigVal
Vector_t< int >	m_idMin
Vector_t< int >	m_idMax
Vector_t< int >	m_dimShell

Detailed Description

Calculate the microcanonical averages with respect to a given sorted vector 'eigVal'.

Template Parameters

int
Container

Member Typedef Documentation

Vector_t

template<class T >

using MicroCanonicalAverage::Vector_t = Eigen::VectorX<T>

private

Constructor & Destructor Documentation

MicroCanonicalAverage()

template<typename InVector_t >

__host__ __device__ MicroCanonicalAverage::MicroCanonicalAverage ( InVector_t const &  eigVal, double const  shellWidth  )

inline

            : m_dim(eigVal.size()), m_shellWidth(shellWidth), m_eigVal{eigVal} {
            m_idMin.resize(m_dim);
            m_idMax.resize(m_dim);
            m_dimShell.resize(m_dim);
 
            // clang-format off
            #pragma omp parallel for
            // clang-format on
            for(auto j = 0; j < m_dim; ++j) {
                int idMin, idMax;
                for(idMin = j; idMin >= 0 && eigVal[j] - eigVal[idMin] <= shellWidth; --idMin) {};
                m_idMin[j] = (++idMin);
                for(idMax = j; idMax < m_dim && eigVal[idMax] - eigVal[j] <= shellWidth; ++idMax) {};
                m_idMax[j]    = (--idMax);
                m_dimShell[j] = idMax - idMin + 1;
            }
        }

Member Function Documentation

dim()

__host__ __device__ int MicroCanonicalAverage::dim ( ) const

inline

{ return m_dim; }

dimShell() [1/3]

__host__ __device__ Vector_t< int > & MicroCanonicalAverage::dimShell ( )

inline

{ return m_dimShell; }

dimShell() [2/3]

__host__ __device__ Vector_t< int > const & MicroCanonicalAverage::dimShell ( ) const

inline

{ return m_dimShell; }

dimShell() [3/3]

__host__ __device__ int MicroCanonicalAverage::dimShell ( int  index ) const

inline

{ return m_dimShell[index]; }

eigVal() [1/2]

__host__ __device__ Vector_t< double > & MicroCanonicalAverage::eigVal ( )

inline

{ return m_eigVal; }

eigVal() [2/2]

__host__ __device__ Vector_t< double > const & MicroCanonicalAverage::eigVal ( ) const

inline

{ return m_eigVal; }

operator()()

template<class VectorType >

__host__ __device__ void MicroCanonicalAverage::operator() ( VectorType &  res, VectorType const &  expVal  ) const

inline

                                                                                             {
            res.resize(m_dim);
            // clang-format off
            #pragma omp parallel for
            // clang-format on
            for(int j = 0; j < m_dim; ++j) {
                res[j] = 0;
                for(auto id = m_idMin[j]; id <= m_idMax[j]; ++id) res[j] += expVal[id];
                res[j] /= double(m_dimShell[j]);
            }
        }

shellWidth()

__host__ __device__ double MicroCanonicalAverage::shellWidth ( ) const

inline

{ return m_shellWidth; }

status()

__host__ void MicroCanonicalAverage::status ( void  ) const

inline

Show the current status of the instance.

• Mean: Mean of each elements
• Stddev: Standard deviation of diagonal elements
• Count: Number of random numbers generated by the instance so far

                                         {
            std::cout << "MicroCanonicalAverage.status()"
                      << "\n"
                      << "\tDim:   " << this->dim() << "\n"
                      << "\tShellWidth: " << this->shellWidth() << std::endl;
        }

Member Data Documentation

m_dim

int MicroCanonicalAverage::m_dim

private

m_dimShell

Vector_t<int> MicroCanonicalAverage::m_dimShell

private

m_eigVal

Vector_t<double> MicroCanonicalAverage::m_eigVal

private

m_idMax

Vector_t<int> MicroCanonicalAverage::m_idMax

private

m_idMin

Vector_t<int> MicroCanonicalAverage::m_idMin

private

m_shellWidth

double MicroCanonicalAverage::m_shellWidth

private

---

The documentation for this class was generated from the following file:

• /Users/shoki/GitHub/Locality/Headers/StatMech/MicroCanonicalAverage.hpp