parse-monitors 0.4.0

Giant Magellan Telescope Computation Fluid Dynamics Parser
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

//! Pressure maps
//!
//! Plot the pressure maps on M1 and M2 segments

use parse_monitors::{
    cfd::{self, BaselineTrait},
    pressure::Pressure,
    CFD_YEAR,
};
use rayon::prelude::*;
use std::{path::Path, time::Instant};

trait Config {
    fn configure(cfd_case: cfd::CfdCase<CFD_YEAR>) -> anyhow::Result<(String, Vec<String>)>;
}
impl Config for geotrans::M1 {
    fn configure(cfd_case: cfd::CfdCase<CFD_YEAR>) -> anyhow::Result<(String, Vec<String>)> {
        Ok((
            "M1p.csv.z".to_string(),
            cfd::CfdDataFile::<CFD_YEAR>::M1Pressure
                .glob(cfd_case)?
                .into_iter()
                .map(|p| p.to_str().unwrap().to_string())
                .collect(),
        ))
    }
}
impl Config for geotrans::M2 {
    fn configure(cfd_case: cfd::CfdCase<CFD_YEAR>) -> anyhow::Result<(String, Vec<String>)> {
        Ok((
            "M2p.csv.z".to_string(),
            cfd::CfdDataFile::<CFD_YEAR>::M2Pressure
                .glob(cfd_case)?
                .into_iter()
                .map(|p| p.to_str().unwrap().to_string())
                .collect(),
        ))
    }
}

fn main() -> anyhow::Result<()> {
    type M12 = geotrans::M2;

    cfd::Baseline::<CFD_YEAR>::default()
        .into_iter()
        .collect::<Vec<cfd::CfdCase<CFD_YEAR>>>()
        .into_par_iter()
        .skip(38)
        .take(2)
        .for_each(|cfd_case| {
            let now = Instant::now();
            let case_path = cfd::Baseline::<CFD_YEAR>::path()
                .expect("undefined path to CFD repository")
                .join(cfd_case.to_string());
            let (_geometry, files) = M12::configure(cfd_case).unwrap();

            let _ = files.last().map(|file| {
                let path = Path::new(file);
                let csv_pressure = Pressure::<M12>::decompress(path.to_path_buf()).unwrap();
                //let csv_geometry =
                //    Pressure::<M12>::decompress(path.with_file_name(geometry)).unwrap();
                let mut pressures = Pressure::<M12>::load(csv_pressure).unwrap();
                pressures.pressure_map(case_path);
            });

            println!(
                "{:<32}{}: {:>8}s",
                cfd_case,
                files.last().unwrap(),
                now.elapsed().as_secs()
            );
        });
    Ok(())
}