Struct ChirikovStandardMap 

pub struct ChirikovStandardMap {
    pub k: f64,
}

The Chirikov–Taylor standard map (area-preserving):

p_{n+1} = p_n + K sin(θ_n)   (mod 2π)
θ_{n+1} = θ_n + p_{n+1}      (mod 2π)

K is the stochasticity parameter; chaos sets in for K ≳ 1.

Fields

k: f64

Stochasticity parameter K.

Implementations

impl ChirikovStandardMap

pub fn new(k: f64) -> Self

Create a Chirikov standard map with parameter K.

pub fn apply(&self, theta: f64, p: f64) -> (f64, f64)

Apply one step. Both theta and p are taken modulo 2π. Returns (theta_new, p_new).

pub fn orbit_2d(&self, theta0: f64, p0: f64, n: usize) -> Vec<(f64, f64)>

Compute an orbit for n steps.

pub fn lyapunov_exponent(&self, theta0: f64, p0: f64, n: usize) -> f64

Estimate the Lyapunov exponent via the method of Benettin et al.

Trait Implementations

impl Clone for ChirikovStandardMap

fn clone(&self) -> ChirikovStandardMap

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for ChirikovStandardMap

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

Formats the value using the given formatter.

impl DiscreteMapIterate for ChirikovStandardMap

fn dim(&self) -> usize

State dimension.

fn step(&self, state: &[f64]) -> Vec<f64>

Advance the state by one step. state must have length dim().

fn orbit(&self, init: &[f64], n: usize) -> Vec<Vec<f64>>

Compute an orbit of length n starting from init.