rdbench 0.1.0

Reaction-diffusion simulation based on Gray-Scott model.
Documentation 
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103

use clap::Parser;
use clap::ValueEnum;

pub fn parse() -> Args {
    Args::parse()
}

#[derive(Parser, Debug)]
#[command(
    version,
    about = "Benchmark for MPI-IO using Gray-Scott reaction-diffusion system"
)]
pub struct Args {
    // Tile partitioning settings
    #[arg(long, default_value = "0")]
    pub nr_tiles_x: usize,
    #[arg(long, default_value = "0")]
    pub nr_tiles_y: usize,
    #[arg(long, default_value = "128")]
    pub sz_tile_x: usize,
    #[arg(long, default_value = "128")]
    pub sz_tile_y: usize,

    // Simulation settings
    #[arg(long, default_value = "20000")]
    pub steps: usize,
    #[arg(long, default_value = "200")]
    pub interval: usize,
    #[arg(long)]
    pub init_output: bool,
    #[arg(long)]
    pub verbose: bool,

    // I/O settings
    #[arg(long, default_value = "out/o_")]
    pub output: String,
    #[arg(long, value_enum, default_value = "canonical")]
    pub file_layout: FileLayout,
    #[arg(long, value_enum, default_value = "v")]
    pub io_field: IoField,
    #[arg(long)]
    pub collective: bool,
    #[arg(long)]
    pub nosync: bool,
    #[arg(long)]
    pub validate: bool,
    #[arg(long)]
    pub prettify: bool,

    // Gray-Scott model parameters
    #[arg(long, default_value = "0.04")]
    pub param_f: f64,
    #[arg(long, default_value = "0.06075")]
    pub param_k: f64,
    #[arg(long, default_value = "0.2")]
    pub param_dt: f64,
    #[arg(long, default_value = "0.1")]
    pub param_du: f64,
    #[arg(long, default_value = "0.05")]
    pub param_dv: f64,
}

impl Args {
    #[allow(dead_code)]
    pub fn nr_files(&self) -> usize {
        if self.interval == 0 {
            return 0;
        }
        let output_count = match self.io_field {
            IoField::U | IoField::V => 1,
            IoField::Both => 2,
            IoField::None => 0,
        };

        (self.steps / self.interval + if self.init_output { 1 } else { 0 }) * output_count
    }
}

#[derive(Parser, Debug, Clone, Copy, PartialEq, Eq, ValueEnum)]
pub enum FileLayout {
    Canonical,
    Log,
}

#[derive(Parser, Debug, Clone, Copy, PartialEq, Eq, ValueEnum)]
pub enum IoField {
    U,
    V,
    Both,
    None,
}

impl IoField {
    #[allow(dead_code)]
    pub fn should_io_u(&self) -> bool {
        matches!(self, IoField::U | IoField::Both)
    }

    #[allow(dead_code)]
    pub fn should_io_v(&self) -> bool {
        matches!(self, IoField::V | IoField::Both)
    }
}