ecrs 0.1.0-beta.4

Evolutionary computation tools & algorithms
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 itertools::Itertools;
use serde::Serialize;

use crate::aco::probe::Probe;
use crate::aco::FMatrix;
use crate::aco::Solution;

#[derive(Serialize)]
#[doc(hidden)]
struct BestSolutionRecord {
    from: usize,
    to: usize,
    iter: usize,
}

#[derive(Serialize)]
#[doc(hidden)]
struct FMatrixRecord {
    from: usize,
    to: usize,
    iter: usize,
    val: f64,
}

/// Probe implementation for writing algorithm output into a csv file.
///
/// Data is flushed only after algorithm ends.
pub struct CsvProbe {
    iteration: usize,
    best_sols: Vec<BestSolutionRecord>,
    pher: Vec<FMatrixRecord>,
    best_sol: Solution,
}

impl CsvProbe {
    pub fn new() -> CsvProbe {
        CsvProbe {
            iteration: 0,
            best_sols: vec![],
            pher: vec![],
            best_sol: Solution::default(),
        }
    }

    #[doc(hidden)]
    fn flush(&mut self) {
        let mut wrt = csv::WriterBuilder::new().from_path("best.csv").unwrap();

        for record in self.best_sols.iter() {
            wrt.serialize(record).expect("Could not serialize record");
        }
        wrt.flush().expect("Could not flush data");

        let mut wrt = csv::WriterBuilder::new().from_path("pheromone.csv").unwrap();

        for record in self.pher.iter() {
            wrt.serialize(record).expect("Could not serialize record");
        }
        wrt.flush().expect("Could not flush data");
    }
}

impl Probe<FMatrix> for CsvProbe {
    fn on_pheromone_update(&mut self, _old_pheromone: &FMatrix, new_pheromone: &FMatrix) {
        for (i, row) in new_pheromone.row_iter().enumerate() {
            for (j, val) in row.iter().enumerate() {
                self.pher.push(FMatrixRecord {
                    from: i,
                    to: j,
                    iter: self.iteration,
                    val: *val,
                })
            }
        }
    }

    fn on_current_best(&mut self, best: &Solution) {
        self.best_sol = best.clone();
    }

    fn on_iteration_start(&mut self) {
        self.iteration += 1;
    }

    fn on_iteration_end(&mut self) {
        for (i, j) in self
            .best_sol
            .path
            .iter()
            .circular_tuple_windows::<(&usize, &usize)>()
        {
            self.best_sols.push(BestSolutionRecord {
                from: *i,
                to: *j,
                iter: self.iteration,
            })
        }
    }

    fn on_end(&mut self) {
        self.flush();
    }
}