KiThe 0.3.4

A numerical suite for chemical kinetics and thermodynamics, combustion, heat and mass transfer,chemical engeneering. Work in progress. Advices and contributions will be appreciated
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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123

use crate::gui::experimental_kinetics_gui::model::{Colours, PlotModel};
use crate::gui::ivp_common::{IvpTaskState, ivp_result_to_experiment};

#[derive(Debug, Clone)]
pub struct DirectProblemDialogState {
    pub open: bool,
    pub task: IvpTaskState,
    pub output_experiment_id: String,
    pub last_error: Option<String>,
    pub last_success: Option<String>,
}

impl Default for DirectProblemDialogState {
    fn default() -> Self {
        Self::new()
    }
}

impl DirectProblemDialogState {
    pub fn new() -> Self {
        Self {
            open: false,
            task: IvpTaskState::new(),
            output_experiment_id: String::from("ivp_result"),
            last_error: None,
            last_success: None,
        }
    }

    pub fn open(&mut self) {
        self.open = true;
        self.last_error = None;
        self.last_success = None;
    }

    fn resolved_experiment_id(&self, model: &PlotModel) -> String {
        let base = self.output_experiment_id.trim();
        let candidate = if base.is_empty() {
            "ivp_result".to_string()
        } else {
            base.to_string()
        };
        if !model.series.ids().iter().any(|id| id == &candidate) {
            return candidate;
        }

        let mut suffix = 1usize;
        loop {
            let next = format!("{}_{}", candidate, suffix);
            if !model.series.ids().iter().any(|id| id == &next) {
                return next;
            }
            suffix += 1;
        }
    }

    pub fn show(&mut self, ctx: &egui::Context, model: &mut PlotModel) {
        if !self.open {
            return;
        }

        let mut open = self.open;
        egui::Window::new("Direct Problem")
            .open(&mut open)
            .default_size([1100.0, 760.0])
            .show(ctx, |ui| {
                self.task.show_model_table(ui);
                ui.separator();
                self.task.show_solver_combo(ui);
                ui.separator();
                self.task.show_parameter_inputs(ui);
                ui.separator();
                self.task.show_problem_inputs(ui);
                ui.separator();

                ui.horizontal(|ui| {
                    ui.label("output experiment:");
                    ui.text_edit_singleline(&mut self.output_experiment_id);
                });

                self.task.show_summary(ui);
                ui.add_space(10.0);

                if ui.button("Solve IVP").clicked() {
                    match self.task.run_kinetic_model() {
                        Ok(result) => {
                            let experiment_id = self.resolved_experiment_id(model);
                            match ivp_result_to_experiment(&result, experiment_id.clone()) {
                                Ok(experiment) => match model
                                    .push_experiment_and_create_plot(experiment, Colours::Red)
                                {
                                    Ok(()) => {
                                        self.last_error = None;
                                        self.last_success =
                                            Some(format!("Created experiment '{}'", experiment_id));
                                    }
                                    Err(err) => {
                                        self.last_error = Some(format!("{:?}", err));
                                    }
                                },
                                Err(err) => {
                                    self.last_error = Some(format!("{:?}", err));
                                }
                            }
                        }
                        Err(err) => {
                            self.last_error = Some(err);
                        }
                    }
                }

                if let Some(msg) = &self.last_success {
                    ui.separator();
                    ui.label(msg);
                }
                if let Some(err) = &self.last_error {
                    ui.separator();
                    ui.label(format!("Error: {}", err));
                }
            });
        self.open = open;
    }
}