Struct lattice_qcd_rs::integrator::symplectic_euler_rayon::SymplecticEulerRayon

pub struct SymplecticEulerRayon;

Basic symplectic Euler integrator using Rayon.

It is slightly faster than super::SymplecticEuler but use slightly more memory.

Example

use lattice_qcd_rs::integrator::{SymplecticEulerRayon, SymplecticIntegrator};
use lattice_qcd_rs::simulation::{
    LatticeStateDefault, LatticeStateEFSyncDefault, LatticeStateWithEField,
};
use rand::SeedableRng;

let mut rng = rand::rngs::StdRng::seed_from_u64(0); // change with your seed
let state1 = LatticeStateEFSyncDefault::new_random_e_state(
    LatticeStateDefault::<3>::new_determinist(1_f64, 2_f64, 4, &mut rng)?,
    &mut rng,
);
let integrator = SymplecticEulerRayon::default();
let state2 = integrator.integrate_symplectic(&state1, 0.000_001_f64)?;
let h = state1.hamiltonian_total();
let h2 = state2.hamiltonian_total();
println!("The error on the Hamiltonian is {}", h - h2);

Implementations

impl SymplecticEulerRayon

pub const fn new() -> Self

Create a new SymplecticEulerRayon

Trait Implementations

impl Clone for SymplecticEulerRayon

fn clone(&self) -> SymplecticEulerRayon

Returns a copy of the value.

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

Performs copy-assignment from source.

impl Debug for SymplecticEulerRayon

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

Formats the value using the given formatter.

impl Default for SymplecticEulerRayon

fn default() -> Self

Identical to SymplecticEulerRayon::new.

impl<'de> Deserialize<'de> for SymplecticEulerRayon

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

Deserialize this value from the given Serde deserializer.

impl Display for SymplecticEulerRayon

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

Formats the value using the given formatter.

impl Hash for SymplecticEulerRayon

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 Ord for SymplecticEulerRayon

fn cmp(&self, other: &SymplecticEulerRayon) -> 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<SymplecticEulerRayon> for SymplecticEulerRayon

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

fn partial_cmp(&self, other: &SymplecticEulerRayon) -> 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 SymplecticEulerRayon

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

Serialize this value into the given Serde serializer.

impl<State, const D: usize> SymplecticIntegrator<State, SimulationStateLeap<State, D>, D> for SymplecticEulerRayonwhere State: SimulationStateSynchronous<D> + LatticeStateWithEField<D> + LatticeStateWithEFieldNew<D>,

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

Type of error returned by the Integrator.

fn integrate_sync_sync( &self, l: &State, delta_t: Real ) -> Result<State, Self::Error>

Integrate a sync state to a sync state by advancing the link matrices and the electrical field simultaneously.

fn integrate_leap_leap( &self, l: &SimulationStateLeap<State, D>, delta_t: Real ) -> Result<SimulationStateLeap<State, D>, Self::Error>

Integrate a leap state to a leap state using leap frog algorithm.

fn integrate_sync_leap( &self, l: &State, delta_t: Real ) -> Result<SimulationStateLeap<State, D>, Self::Error>

Integrate a sync state to a leap state by doing a half step for the conjugate momenta.

fn integrate_leap_sync( &self, l: &SimulationStateLeap<State, D>, delta_t: Real ) -> Result<State, Self::Error>

Integrate a leap state to a sync state by finishing doing a step for the position and finishing the half step for the conjugate momenta.

fn integrate_symplectic( &self, l: &State, delta_t: Real ) -> Result<State, Self::Error>

Integrate a Sync state by going to leap and then back to sync. This is the symplectic method of integration, which should conserve approximately the hamiltonian

impl Copy for SymplecticEulerRayon

impl Eq for SymplecticEulerRayon

impl StructuralEq for SymplecticEulerRayon

impl StructuralPartialEq for SymplecticEulerRayon