ExtensiveMoments

class ExtensiveMoments : public feasst::Analyze

Accumulate extensive moments for derivatives from fluctuations. This is currently implemented for a grand canonical ensemble:

\(N_i^j N_k^m U^p\)

where \(i,k\) are particle types, and the powers \(j,m,p\) are collected from 0 to a maximum order cutoff.

Arguments

  • max_order: maximum order cutoff for the powers (default: 3).

  • Stepper arguments.

Public Functions

std::string header(const MonteCarlo &mc) const

Return the header for writing.

void initialize(MonteCarlo *mc)

Initialize and precompute before trials.

void update(const MonteCarlo &mc)

Perform update action.

std::string write(const MonteCarlo &mc)

Write the moments by serialization of a 5d accumulator. If checkpoints do not work, this file can be read using feasst_deserialize

const Accumulator &moments(const int p, const int m, const int k, const int j, const int i) const

Return the extensive moments.