Struct lattice_qcd_rs::simulation::monte_carlo::overrelaxation::OverrelaxationSweepRotation

pub struct OverrelaxationSweepRotation;

Overrelaxation algorithm using rotation method.

Alone it can’t advance the simulation as it preserved the hamiltonian. You need to use other method with this one. You can look at super::HybridMethodVec and super::HybridMethodCouple.

The algorithm is described https://arxiv.org/abs/hep-lat/0503041 in section 2.1 up to step 2 using \hat X_{NN}.

Example

use lattice_qcd_rs::simulation::{HeatBathSweep, LatticeState, LatticeStateDefault, OverrelaxationSweepRotation};
use rand::SeedableRng;

let rng = rand::rngs::StdRng::seed_from_u64(0); // change with your seed
let mut heat_bath = HeatBathSweep::new(rng);
let mut overrelax = OverrelaxationSweepRotation::default();

let mut state = LatticeStateDefault::<3>::new_cold(1_f64, 8_f64, 4)?; // 1_f64 : size, 8_f64: beta, 4 number of points.
for _ in 0..2 {
    state = state.monte_carlo_step(&mut heat_bath)?;
    state = state.monte_carlo_step(&mut overrelax)?;
    // operation to track the progress or the evolution
}
// operation at the end of the simulation

Implementations

impl OverrelaxationSweepRotation

pub const fn new() -> Self

Create a new Self with an given RNG

Trait Implementations

impl Clone for OverrelaxationSweepRotation

fn clone(&self) -> OverrelaxationSweepRotation

Returns a copy of the value.

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

Performs copy-assignment from source.

impl Debug for OverrelaxationSweepRotation

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

Formats the value using the given formatter.

impl Default for OverrelaxationSweepRotation

fn default() -> Self

Returns the “default value” for a type.

impl<'de> Deserialize<'de> for OverrelaxationSweepRotation

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

Deserialize this value from the given Serde deserializer.

impl Display for OverrelaxationSweepRotation

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

Formats the value using the given formatter.

impl Hash for OverrelaxationSweepRotation

fn hash<__H: Hasher>(&self, state: &mut __H)

Feeds this value into the given Hasher.

fn hash_slice<H>(data: &[Self], state: &mut H)where H: Hasher, Self: Sized,

Feeds a slice of this type into the given Hasher.

impl<const D: usize> MonteCarlo<LatticeStateDefault<D>, D> for OverrelaxationSweepRotation

type Error = Never

Error returned while getting the next element.

fn next_element( &mut self, state: LatticeStateDefault<D> ) -> Result<LatticeStateDefault<D>, Self::Error>

Do one Monte Carlo simulation step.

impl Ord for OverrelaxationSweepRotation

fn cmp(&self, other: &OverrelaxationSweepRotation) -> Ordering

This method returns an Ordering between self and other.

fn max(self, other: Self) -> Selfwhere Self: Sized,

Compares and returns the maximum of two values.

fn min(self, other: Self) -> Selfwhere Self: Sized,

Compares and returns the minimum of two values.

fn clamp(self, min: Self, max: Self) -> Selfwhere Self: Sized + PartialOrd<Self>,

Restrict a value to a certain interval.

impl PartialEq<OverrelaxationSweepRotation> for OverrelaxationSweepRotation

fn eq(&self, other: &OverrelaxationSweepRotation) -> 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 PartialOrd<OverrelaxationSweepRotation> for OverrelaxationSweepRotation

fn partial_cmp(&self, other: &OverrelaxationSweepRotation) -> Option<Ordering>

This method returns an ordering between self and other values if one exists.

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

This method tests less than (for self and other) and is used by the < operator.

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

This method tests less than or equal to (for self and other) and is used by the <= operator.

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

This method tests greater than (for self and other) and is used by the > operator.

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

This method tests greater than or equal to (for self and other) and is used by the >= operator.

impl Serialize for OverrelaxationSweepRotation

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 OverrelaxationSweepRotation

impl Eq for OverrelaxationSweepRotation

impl StructuralEq for OverrelaxationSweepRotation

impl StructuralPartialEq for OverrelaxationSweepRotation