wafrift_evolution/search/
tabu.rs

1use crate::evolution::crossover::mutation::mutate_with_log;
2use crate::evolution::{Chromosome, GenePool, population::random_chromosome};
3use crate::lineage::Lineage;
4use crate::search::{EvalCandidate, SearchAlgorithm, comparable_fitness, fitness_cmp};
5use crate::types::{Budget, EvolutionError, OracleVerdict, SearchStats};
6use rand::rngs::StdRng;
7use serde::{Deserialize, Serialize};
8use std::collections::{HashSet, VecDeque};
9
10/// Tabu search with aspiration criteria.
11#[derive(Debug, Clone, Serialize, Deserialize)]
12pub struct TabuSearch {
13    current: Chromosome,
14    best: Chromosome,
15    gene_pool: GenePool,
16    generation: u32,
17    eval_counter: u64,
18    tabu_list: VecDeque<u64>,
19    tabu_tenure: usize,
20    tabu_set: HashSet<u64>,
21}
22
23impl TabuSearch {
24    #[must_use]
25    pub fn new(tabu_tenure: usize) -> Self {
26        Self {
27            current: Chromosome::new(vec![]),
28            best: Chromosome::new(vec![]),
29            gene_pool: GenePool::default_wafrift(),
30            generation: 0,
31            eval_counter: 0,
32            tabu_list: VecDeque::new(),
33            tabu_tenure,
34            tabu_set: HashSet::new(),
35        }
36    }
37
38    fn neighbor(&self, rng: &mut StdRng) -> Chromosome {
39        let mut child = self.current.clone();
40        let log = mutate_with_log(&mut child, &self.gene_pool, 0.25, rng);
41        child.lineage = Lineage::mutation(&self.current, log, self.generation);
42        child
43    }
44
45    fn add_tabu(&mut self, hash: u64) {
46        if self.tabu_set.insert(hash) {
47            self.tabu_list.push_back(hash);
48        }
49        while self.tabu_list.len() > self.tabu_tenure {
50            if let Some(old) = self.tabu_list.pop_front() {
51                self.tabu_set.remove(&old);
52            }
53        }
54    }
55}
56
57impl Default for TabuSearch {
58    fn default() -> Self {
59        Self::new(20)
60    }
61}
62
63impl SearchAlgorithm for TabuSearch {
64    fn name(&self) -> &'static str {
65        "tabu_search"
66    }
67
68    fn initialize(&mut self, population: Vec<Chromosome>, gene_pool: &GenePool, rng: &mut StdRng) {
69        self.gene_pool = gene_pool.clone();
70        if let Some(best) = population
71            .iter()
72            .max_by(|a, b| fitness_cmp(a.fitness, b.fitness))
73        {
74            self.current = best.clone();
75            self.best = best.clone();
76            self.add_tabu(best.hash());
77        } else {
78            self.current = random_chromosome(gene_pool, rng);
79            self.best = self.current.clone();
80            self.add_tabu(self.current.hash());
81        }
82    }
83
84    fn request_evaluations(&mut self, n: usize, rng: &mut StdRng) -> Vec<EvalCandidate> {
85        let mut out = Vec::with_capacity(n);
86        let mut attempts = 0;
87        while out.len() < n && attempts < n * 10 {
88            attempts += 1;
89            self.eval_counter += 1;
90            let candidate = self.neighbor(rng);
91            let hash = candidate.hash();
92            // Tabu check. The aspiration criterion ("allow a tabu move
93            // if it beats the current best") is intentionally NOT
94            // applied here: candidate has fitness == 0.0 because it
95            // hasn't been evaluated yet, so the comparison would
96            // always be false and the algorithm would deadlock when
97            // every neighbour is tabu. Aspiration belongs in
98            // submit_evaluations, where fitness is real. Removing it
99            // here matches the SimulatedAnnealing/HillClimbing flow.
100            let is_tabu = self.tabu_set.contains(&hash);
101            if !is_tabu {
102                out.push(EvalCandidate {
103                    id: self.eval_counter,
104                    chromosome: candidate,
105                });
106            }
107        }
108        out
109    }
110
111    fn submit_evaluations(&mut self, results: Vec<(u64, OracleVerdict)>) {
112        for (_id, verdict) in results {
113            let mut candidate = self.current.clone();
114            candidate.record_verdict(&verdict);
115            self.add_tabu(candidate.hash());
116            if comparable_fitness(candidate.fitness) >= comparable_fitness(self.current.fitness) {
117                self.current = candidate;
118                if comparable_fitness(self.current.fitness) > comparable_fitness(self.best.fitness) {
119                    self.best = self.current.clone();
120                }
121            }
122        }
123        self.generation += 1;
124    }
125
126    fn should_terminate(&self, stats: &SearchStats, budget: &Budget) -> bool {
127        stats.evaluations >= budget.max_requests
128            || stats.generation >= budget.max_generations
129            || stats.stagnation_counter >= budget.stagnation_limit
130    }
131
132    fn best(&self) -> Option<&Chromosome> {
133        Some(&self.best)
134    }
135
136    fn checkpoint(&self) -> Result<Vec<u8>, EvolutionError> {
137        serde_json::to_vec(self).map_err(|e| EvolutionError::SerializationFailed(e.to_string()))
138    }
139
140    fn restore(&mut self, bytes: &[u8]) -> Result<(), EvolutionError> {
141        *self = serde_json::from_slice(bytes)
142            .map_err(|e| EvolutionError::DeserializationFailed(e.to_string()))?;
143        Ok(())
144    }
145}
146
147#[cfg(test)]
148mod tests {
149    use super::*;
150    use rand::SeedableRng;
151
152    #[test]
153    fn non_finite_fitness_does_not_block_future_updates() {
154        let mut alg = TabuSearch::new(10);
155        let pool = GenePool::default_wafrift();
156        let mut rng = StdRng::seed_from_u64(11);
157        alg.initialize(vec![Chromosome::new(vec![])], &pool, &mut rng);
158
159        alg.submit_evaluations(vec![(
160            1,
161            OracleVerdict {
162                passed: false,
163                status_delta: 0,
164                body_delta: 0,
165                latency_ms: 0,
166                confidence: f64::NAN,
167                triggered_rules: 1,
168            },
169        )]);
170        let best_after_nan = comparable_fitness(alg.best().expect("best must exist").fitness);
171
172        alg.submit_evaluations(vec![(2, OracleVerdict::from_bool(true))]);
173        let best_after_valid = comparable_fitness(alg.best().expect("best must exist").fitness);
174        assert!(best_after_valid > best_after_nan);
175    }
176}
wafrift_evolution/search/tabu.rs

wafrift_evolution/search/
tabu.rs
