Struct lattice_qcd_rs::simulation::monte_carlo::metropolis_hastings::MetropolisHastingsDiagnostic

pub struct MetropolisHastingsDiagnostic { /* private fields */ }

Metropolis Hastings algorithm with diagnostics. Very slow, use MetropolisHastingsDeltaDiagnostic instead.

Similar to MetropolisHastingsDiagnostic but with diagnostic information.

Note that this method does not do a sweep but change random link matrix, for a sweep there is super::MetropolisHastingsSweep.

Example

see example of super::McWrapper

Implementations

impl MetropolisHastingsDiagnostic

pub fn new(number_of_update: usize, spread: Real) -> Option<Self>

spread should be between 0 and 1 both not included and number_of_update should be greater than 0. 0.1_f64 is a good choice for this parameter.

number_of_update is the number of times a link matrix is randomly changed. spread is the spread factor for the random matrix change ( used in su3::random_su3_close_to_unity).

pub const fn prob_replace_last(&self) -> Real

Get the last probably of acceptance of the random change.

pub const fn has_replace_last(&self) -> bool

Get if last step has accepted the replacement.

pub const fn number_of_update(&self) -> usize

Get the number of updates per steps.

pub const fn spread(&self) -> Real

Get the spread parameter.

Trait Implementations

impl Clone for MetropolisHastingsDiagnostic

fn clone(&self) -> MetropolisHastingsDiagnostic

Returns a copy of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for MetropolisHastingsDiagnostic

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl Default for MetropolisHastingsDiagnostic

fn default() -> Self

Returns the “default value” for a type.

impl<'de> Deserialize<'de> for MetropolisHastingsDiagnostic

fn deserialize<__D>(__deserializer: __D) -> Result<Self, __D::Error>where __D: Deserializer<'de>,

Deserialize this value from the given Serde deserializer.

impl Display for MetropolisHastingsDiagnostic

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl<State, const D: usize> MonteCarloDefault<State, D> for MetropolisHastingsDiagnosticwhere State: LatticeState<D> + LatticeStateNew<D>,

type Error = <State as LatticeStateNew<D>>::Error

Error returned while getting the next element.

fn potential_next_element<Rng>( &mut self, state: &State, rng: &mut Rng ) -> Result<State, Self::Error>where Rng: Rng + ?Sized,

Generate a radom element from the previous element (like a Markov chain).

fn next_element_default<Rng>( &mut self, state: State, rng: &mut Rng ) -> Result<State, Self::Error>where Rng: Rng + ?Sized,

Get the next element in the chain either the old state or a new one replacing it.

fn probability_of_replacement(old_state: &State, new_state: &State) -> Real

probability of the next element to replace the current one.

impl PartialEq<MetropolisHastingsDiagnostic> for MetropolisHastingsDiagnostic

fn eq(&self, other: &MetropolisHastingsDiagnostic) -> bool

This method tests for self and other values to be equal, and is used by ==.

fn ne(&self, other: &Rhs) -> bool

This method tests for !=. The default implementation is almost always sufficient, and should not be overridden without very good reason.

impl Serialize for MetropolisHastingsDiagnostic

fn serialize<__S>(&self, __serializer: __S) -> Result<__S::Ok, __S::Error>where __S: Serializer,

Serialize this value into the given Serde serializer.

impl Copy for MetropolisHastingsDiagnostic

impl StructuralPartialEq for MetropolisHastingsDiagnostic