diffusionx 0.12.0

A multi-threaded crate for random number generation and stochastic process simulation, with optional GPU acceleration.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42

use diffusionx::simulation::{continuous::Bm, prelude::*};

fn main() {
    println!("===== Example of Brownian Motion Simulation =====");

    // Create Brownian motion object
    let bm = Bm::default();

    // Time settings
    let t_max = 10.0; // Maximum time
    let dt = 0.01; // Time step

    // Generate Brownian motion trajectory
    let (times, positions) = bm.simulate(t_max, dt).unwrap();

    // Print some data points
    println!("Time\tPosition");
    for i in (0..times.len()).step_by(100) {
        println!("{:.2}\t{:.6}", times[i], positions[i]);
    }

    // Visualize trajectory
    #[cfg(feature = "visualize")]
    {
        // Create trajectory and visualize
        let traj = bm.duration(t_max).unwrap();
        let config = PlotConfigBuilder::default()
            .time_step(dt)
            .output_path("tmp/bm.svg")
            .caption("Brownian Motion Trajectory")
            .x_label("t")
            .y_label("X(t)")
            .legend("bm")
            .size((800, 600))
            .backend(PlotterBackend::SVG)
            .build()
            .unwrap();

        traj.plot(&config).unwrap();
        println!("Trajectory image saved to tmp/bm.svg");
    }
}