using Arc<Mutex<Vec<Human>>>

src/human.rs

@@ -14,20 +14,20 @@ pub struct Human {
     pub x: i32,
     pub y: i32,
 }
-impl Human {
+// impl Human {
-    // humans with state "Normal"
+//     // humans with state "Normal"
-    pub fn new(pos_x: i32, pos_y: i32) -> Self {
+//     pub fn new(pos_x: i32, pos_y: i32) -> Self {
-        Self {
+//         Self {
-            present_state: State::Normal,
+//             present_state: State::Normal,
-            x: pos_x,
+//             x: pos_x,
-            y: pos_y,
+//             y: pos_y,
-        }
+//         }
-    }
+//     }
-    // pub fn new_empty() -> Self{
+//     // pub fn new_empty() -> Self{
-    //     Self{
+//     //     Self{
-    //         present_state:State::Normal,
+//     //         present_state:State::Normal,
-    //         x : 0,
+//     //         x : 0,
-    //         y : 0,
+//     //         y : 0,
-    //     }
+//     //     }
-    // }
+//     // }
-}
+// }

src/main.rs

@@ -17,18 +17,13 @@ fn main() {
     let term = Term::stdout();
     term.write_line("********** Rusty Propagation (Console) 2022 **********")
         .expect("Oops Looks like we have a problem here...");
-    term.write_line("Press any key to start the propagation")
-        .expect("Oops Looks like we have a problem here...");
 
     let disease = Disease::new(20, 10, 5, String::from("Covid 44"));
     let mut population = Population::new(20, 10, 5, 1000, 1000, disease);
     //population.change_disease(disease);
-    println!("Before Filling");
-    //population.display();
-    population.generate();
     println!("After Filling");
     //population.display();
-    println!("After Propagation");
+    //population.display();
     let mut stats: [i32; 4];
     // = [0,0,0,0];
     let mut counter: u32 = 0;

src/population.rs

@@ -1,4 +1,8 @@
+use std::sync::Arc;
+use std::sync::Mutex;
 use crate::prelude::*;
+use std::thread;
+
 
 #[derive(Debug)]
 pub struct Point {
@@ -10,11 +14,12 @@ pub struct Population {
     pub start_infected_ratio: i32,
     pub start_immune_ratio: i32,
     pub start_dead_ratio: i32,
-    pub humans: Vec<Human>,
+    pub humans: Arc<Mutex<Vec<Human>>>,
     pub width: i32,
     pub height: i32,
     pub age: i32,
     pub plague: Disease,
+    pub size: usize,
 }
 
 pub fn human_idx(x: i32, y: i32, width: i32) -> usize {
@@ -30,12 +35,37 @@ impl Population {
         height: i32,
         plague: Disease,
     ) -> Self {
-        let mut the_humans: Vec<Human> = vec![Human{x: 0, y: 0, present_state: State::Normal}; (width * height) as usize];
+        let mut rng = rand::thread_rng();
+
+        let size: usize = (width * height) as usize;
+
+        let the_humans_arc = Arc::new(Mutex::new(vec![
+            Human {
+                x: 0,
+                y: 0,
+                present_state: State::Normal
+            };
+            size
+        ]));
+        let the_humans = Arc::clone(&the_humans_arc);
         for x in 0..width {
             for y in 0..height {
                 let idx = human_idx(x, y, width);
-                the_humans[idx].x = x;
+                let mut present_state = State::Normal;
-                the_humans[idx].y = y;
+                if (start_infected_ratio > 0)
+                    && (rng.gen_range(0..CORRECTED_PERCENTAGE) <= start_infected_ratio)
+                {
+                    present_state = State::Infected;
+                } else if (start_immune_ratio > 0)
+                    && (rng.gen_range(0..CORRECTED_PERCENTAGE) <= start_immune_ratio)
+                {
+                    present_state = State::Immune;
+                } else if (start_dead_ratio > 0)
+                    && (rng.gen_range(0..CORRECTED_PERCENTAGE) <= start_dead_ratio)
+                {
+                    present_state = State::Dead;
+                }
+                the_humans.lock().unwrap()[idx] = Human{x: x, y: y, present_state: present_state};
             }
         }
         Self {
@@ -46,36 +76,13 @@ impl Population {
             height: height,
             plague: plague,
             age: 0,
-            humans: the_humans,
+            humans: the_humans_arc,
+            size: size,
         }
     }
     // pub fn change_disease(&mut self, plague:Disease){
     //     self.plague = plague;
     // }
-    pub fn generate(&mut self) {
-        //The ratios will not be exact, for example someone who wants 100% infected 100% immune and 100% dead, he will have 100% dead because they are overwriting each others
-        //Maybe consider limiting the total to not exceed 100 in the view
-        //Other thing, there will always be more of the last one because someone who is already infected for example could be then put to immune or dead
-        //One solution to this issue would be to have if else else statements but in this case for example 20% would be lower with the last because its 20% on the remaining population and not of the all
-        //In other words I did it that way but it can be changed just its not the right method to have perfect ratios
-        let mut rng = rand::thread_rng();
-
-        for x in self.humans.iter_mut() {
-            if (self.start_infected_ratio) > 0
-                && (rng.gen_range(0..CORRECTED_PERCENTAGE) <= self.start_infected_ratio as i32)
-            {
-                x.present_state = State::Infected;
-            } else if self.start_immune_ratio > 0
-                && rng.gen_range(0..CORRECTED_PERCENTAGE) <= self.start_immune_ratio as i32
-            {
-                x.present_state = State::Immune;
-            } else if self.start_dead_ratio > 0
-                && rng.gen_range(0..CORRECTED_PERCENTAGE) <= self.start_dead_ratio as i32
-            {
-                x.present_state = State::Dead;
-            }
-        }
-    }
 
     fn is_inside(&self, pos: &Point) -> bool {
         if pos.x >= 0 && pos.x < self.width && pos.y >= 0 && pos.y < self.height {
@@ -85,31 +92,44 @@ impl Population {
         }
     }
 
-    fn push_if_inside(&self, point_list: &mut Vec<Point>, point: Point) {
+    fn is_inside_and_infected(&self, point: Point) -> bool {
+        let the_humans_arc = Arc::clone(&self.humans);
         if self.is_inside(&point) {
-            point_list.push(point);
+            let idx = human_idx(point.x, point.y, self.width);
+            let humans = the_humans_arc.lock().unwrap();
+            if humans[idx].present_state == State::Infected {
+                roll(self.plague.infection_rate)
+            } else {
+                false
+            }
+        } else {
+            false
         }
     }
 
     pub fn propagate(&mut self) -> [i32; 4] {
         let mut people_to_check: Vec<Point> =
-            Vec::with_capacity((self.width * self.height) as usize);
+            Vec::with_capacity(self.size);
+        let mut possible_infected: Vec<Point> =
+            Vec::with_capacity(self.size);
         let mut people_to_infect: Vec<Point> =
-            Vec::with_capacity((self.width * self.height) as usize);
+            Vec::with_capacity(self.size);
         let mut people_to_cure: Vec<Point> =
-            Vec::with_capacity((self.width * self.height) as usize);
+            Vec::with_capacity(self.size);
         let mut people_to_kill: Vec<Point> =
-            Vec::with_capacity((self.width * self.height) as usize);
+            Vec::with_capacity(self.size);
         let mut stats: [i32; 4] = [0, 0, 0, 0];
         // stats[0] Normal stats[1] Infected stats[2] Immune stats[3] Dead
 
-        for h in self.humans.iter() {
+        let humans = Arc::clone(&self.humans);
+        for h in humans.lock().unwrap().iter() {
             match h.present_state {
                 State::Normal => {
+                    possible_infected.push(Point{ x: h.x, y: h.y});
                     stats[0] += 1;
                 }
                 State::Infected => {
-                    people_to_check.push(Point{ x: h.x, y: h.y });
+                    people_to_check.push(Point { x: h.x, y: h.y });
                     stats[1] += 1;
                 }
                 State::Immune => {
@@ -133,120 +153,123 @@ impl Population {
                 if roll(self.plague.death_rate) {
                     //cheks if the man dies
                     people_to_kill.push(Point { x: pos.x, y: pos.y });
-                } else {
-                    let mut possible_infections: Vec<Point> = Vec::with_capacity(8);
-                    // Vec::new();
-                    //possible_infections.push(Point{x:pos.x,y:pos.y});
-
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x - 1,
-                            y: pos.y - 1,
-                        },
-                    ); //Top Left
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x,
-                            y: pos.y - 1,
-                        },
-                    ); //Top
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x + 1,
-                            y: pos.y - 1,
-                        },
-                    ); //Top Right
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x - 1,
-                            y: pos.y,
-                        },
-                    ); //Left
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x + 1,
-                            y: pos.y,
-                        },
-                    ); //Right
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x - 1,
-                            y: pos.y + 1,
-                        },
-                    ); //Bottom Left
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x,
-                            y: pos.y + 1,
-                        },
-                    ); //Bottom
-                    self.push_if_inside(
-                        &mut possible_infections,
-                        Point {
-                            x: pos.x + 1,
-                            y: pos.y + 1,
-                        },
-                    ); //Bottom Right
-                    for poss_infected_pos in possible_infections.iter() {
-                        //possible_infections.iter().map(|poss_infected_pos|{
-                        let inf_idx =
-                            human_idx(poss_infected_pos.x, poss_infected_pos.y, self.width);
-                        if self.humans[inf_idx].present_state == State::Normal {
-                            if roll(self.plague.infection_rate) {
-                                people_to_infect.push(Point {
-                                    x: poss_infected_pos.x,
-                                    y: poss_infected_pos.y,
-                                });
-                            }
-                        }
-                    }
                 }
             }
         }
-
+        for pos in possible_infected.iter() {
-        for infected_position in people_to_infect.iter() {
+            let infected: bool = self.is_inside_and_infected(
-            // println!("To infect: {:?}", infected_position);
+                Point {
-            //people_to_infect.iter().map(|infected_position|{
+                    x: pos.x - 1,
-            let infected_index = human_idx(infected_position.x, infected_position.y, self.width);
+                    y: pos.y - 1,
-            // let _ = infected_position.x;
+                },
-            //DEBUG
+            ) || //Top Left
-            //println!("x: {} y: {} index: {}",infected_position.x,infected_position.y,infected_index);
+            self.is_inside_and_infected(
-            self.humans[infected_index].present_state = State::Infected;
+                Point {
-            //DEBUG
+                    x: pos.x,
-            //println!("Infected someone");
+                    y: pos.y - 1,
+                },
+            ) || //Top
+            self.is_inside_and_infected(
+                Point {
+                    x: pos.x + 1,
+                    y: pos.y - 1,
+                },
+            ) || //Top Right
+            self.is_inside_and_infected(
+                Point {
+                    x: pos.x - 1,
+                    y: pos.y,
+                },
+            ) || //Left
+            self.is_inside_and_infected(
+                Point {
+                    x: pos.x + 1,
+                    y: pos.y,
+                },
+            ) || //Right
+            self.is_inside_and_infected(
+                Point {
+                    x: pos.x - 1,
+                    y: pos.y + 1,
+                },
+            ) || //Bottom Left
+            self.is_inside_and_infected(
+                Point {
+                    x: pos.x,
+                    y: pos.y + 1,
+                },
+            ) || //Bottom
+            self.is_inside_and_infected(
+                Point {
+                    x: pos.x + 1,
+                    y: pos.y + 1,
+                },
+            ); //Bottom Right
+            if infected {
+                people_to_infect.push(Point { x: pos.x, y: pos.y });
+            }
         }
 
-        for cured_position in people_to_cure.iter() {
+//        for infected_position in people_to_infect.iter() {
-            //people_to_cure.iter().map(|cured_position|{
+//            // println!("To infect: {:?}", infected_position);
-            let cured_index = human_idx(cured_position.x, cured_position.y, self.width);
+//            //people_to_infect.iter().map(|infected_position|{
-            if self.humans[cured_index].present_state != State::Infected {
+//            let infected_index = human_idx(infected_position.x, infected_position.y, self.width);
-                println!("not infected");
+//            // let _ = infected_position.x;
-            }
+//            //DEBUG
-            self.humans[cured_index].present_state = State::Immune;
+//            //println!("x: {} y: {} index: {}",infected_position.x,infected_position.y,infected_index);
-            //DEBUG
+//            self.humans[infected_index].present_state = State::Infected;
-            //println!("Cured someone");
+//            //DEBUG
+//            //println!("Infected someone");
+//        }
+
+        let mut threads = vec![];
+        {
+            let humans = Arc::clone(&self.humans);
+            let width = self.width;
+            threads.push(thread::spawn(move || {
+                for infected_position in people_to_infect.iter() {
+                    let infected_index = human_idx(infected_position.x, infected_position.y, width);
+                    humans.lock().unwrap()[infected_index].present_state = State::Infected;
+                }
+            }));
+        }
+
+        {
+            let humans = Arc::clone(&self.humans);
+            let width = self.width;
+            threads.push(thread::spawn(move || {
+                for cured_position in people_to_cure.iter() {
+                    //people_to_cure.iter().map(|cured_position|{
+                    let cured_index = human_idx(cured_position.x, cured_position.y, width);
+                    if humans.lock().unwrap()[cured_index].present_state != State::Infected {
+                        println!("not infected");
+                    } else {
+                        humans.lock().unwrap()[cured_index].present_state = State::Immune;
+                    }
+                    //DEBUG
+                    //println!("Cured someone");
+                }
+            }));
+        }
+
+        for t in threads {
+            t.join().unwrap();
         }
 
         for dead_position in people_to_kill.iter() {
+            let humans = Arc::clone(&self.humans);
             //people_to_kill.iter().map(|dead_position|{
             let dead_index = human_idx(dead_position.x, dead_position.y, self.width);
-            if self.humans[dead_index].present_state == State::Dead {
+            if humans.lock().unwrap()[dead_index].present_state == State::Dead {
                 // println!("Already dead");
             } else {
-                self.humans[dead_index].present_state = State::Dead;
+                humans.lock().unwrap()[dead_index].present_state = State::Dead;
             }
             //DEBUG
         }
         assert_eq!(
             stats[0] + stats[1] + stats[2] + stats[3],
-            self.humans.len() as i32
+            self.size as i32
         );
         stats
     }
@@ -343,7 +366,11 @@ mod tests {
         let mut stats: Stats;
 
         for _ in 0..10 {
-            humans.push(Human::new(0, 0));
+            humans.push(Human {
+                present_state: State::Normal,
+                x: 0,
+                y: 0,
+            });
         }
         stats = humans_stats(&humans);
         assert_eq!(stats.normal, 10);
@@ -369,29 +396,22 @@ mod tests {
         let disease = Disease::new(20, 10, 5, String::from("Covid 44"));
         let (width, height) = (5, 7);
         let population = Population::new(20, 10, 5, 5, 7, disease);
-        assert_eq!(population.humans.len(), 5 * 7);
+        let humans = Arc::clone(&population.humans);
-        for human in population.humans.iter() {
+        assert_eq!(humans.lock().unwrap().len(), 5 * 7);
-            assert!(
+        for h in humans.lock().unwrap().iter() {
-                human.present_state == State::Normal,
-                "all humans should be normal"
-            );
-        }
-        for h in population.humans.iter() {
             let idx = human_idx(h.x, h.y, width);
-            assert_eq!(population.humans[idx].x, h.x, "coordinates should match");
+            assert_eq!(humans.lock().unwrap()[idx].x, h.x, "coordinates should match");
-            assert_eq!(population.humans[idx].y, h.y, "coordinates should match");
+            assert_eq!(humans.lock().unwrap()[idx].y, h.y, "coordinates should match");
         }
-        assert_eq!(population.humans.len(), (width * height) as usize);
+        assert_eq!(humans.lock().unwrap().len(), (width * height) as usize);
-
     }
 
     #[test]
     fn population_gen() {
         let disease = Disease::new(20, 10, 5, String::from("Covid 44"));
         let (width, height) = (5, 7);
-        let mut population = Population::new(20, 10, 5, 5, 7, disease);
+        let population = Population::new(20, 10, 5, 5, 7, disease);
 
-        population.generate();
         let stats: Stats = humans_stats(&population.humans);
         println!("Stats: {:?}", stats);
 
@@ -410,28 +430,24 @@ mod tests {
 
         disease = Disease::new(0, 0, 0, String::from("Test"));
         population = Population::new(0, 0, 0, width, height, disease);
-        population.generate();
         stats = humans_stats(&population.humans);
         println!("should be normal: {:?}", stats);
         assert_eq!(stats.normal, width * height);
 
         disease = Disease::new(0, 0, 0, String::from("Test"));
         population = Population::new(100, 0, 0, width, height, disease);
-        population.generate();
         stats = humans_stats(&population.humans);
         println!("should be infected: {:?}", stats);
         assert_eq!(stats.infected, width * height);
 
         disease = Disease::new(0, 0, 0, String::from("Test"));
         population = Population::new(0, 100, 0, width, height, disease);
-        population.generate();
         stats = humans_stats(&population.humans);
         println!("should be immune: {:?}", stats);
         assert_eq!(stats.immune, width * height);
 
         disease = Disease::new(0, 0, 0, String::from("Test"));
         population = Population::new(0, 0, 100, width, height, disease);
-        population.generate();
         stats = humans_stats(&population.humans);
         println!("should be dead: {:?}", stats);
         assert_eq!(stats.dead, width * height);
@@ -457,16 +473,10 @@ mod tests {
         let mut stats: Stats;
         let mut propagate_stats: [i32; 4];
 
-        // start with normal population
+        // infect every one
         stats = humans_stats(&population.humans);
         println!("stats after init: {:?}", stats);
-        assert_eq!(stats.normal, 100);
+        assert_eq!(stats.infected, 100, "everybody should be infected");
-
-        // infect every one
-        population.generate();
-        stats = humans_stats(&population.humans);
-        println!("stats after generate: {:?}", stats);
-        assert_eq!(stats.infected, 100);
 
         // kill every one
         propagate_stats = population.propagate();
@@ -499,15 +509,51 @@ mod tests {
 
         // start with normal population
         population.humans = vec![
-            Human{present_state: State::Normal, x: 0, y: 0},
+            Human {
-            Human{present_state: State::Normal, x: 1, y: 0},
+                present_state: State::Normal,
-            Human{present_state: State::Normal, x: 2, y: 0},
+                x: 0,
-            Human{present_state: State::Normal, x: 0, y: 1},
+                y: 0,
-            Human{present_state: State::Infected, x: 1, y: 1},
+            },
-            Human{present_state: State::Normal, x: 2, y: 1},
+            Human {
-            Human{present_state: State::Normal, x: 0, y: 2},
+                present_state: State::Normal,
-            Human{present_state: State::Normal, x: 1, y: 2},
+                x: 1,
-            Human{present_state: State::Normal, x: 2, y: 2},
+                y: 0,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 2,
+                y: 0,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 0,
+                y: 1,
+            },
+            Human {
+                present_state: State::Infected,
+                x: 1,
+                y: 1,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 2,
+                y: 1,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 0,
+                y: 2,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 1,
+                y: 2,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 2,
+                y: 2,
+            },
         ];
         stats = humans_stats(&population.humans);
         println!("stats after init: {:?}", stats);
@@ -544,44 +590,113 @@ mod tests {
 
         // start with normal population
         population.humans = vec![
-            Human{present_state: State::Normal, x: 0, y: 0},
+            Human {
-            Human{present_state: State::Normal, x: 1, y: 0},
+                present_state: State::Normal,
-            Human{present_state: State::Normal, x: 2, y: 0},
+                x: 0,
-            Human{present_state: State::Normal, x: 0, y: 1},
+                y: 0,
-            Human{present_state: State::Infected, x: 1, y: 1},
+            },
-            Human{present_state: State::Normal, x: 2, y: 1},
+            Human {
-            Human{present_state: State::Normal, x: 0, y: 2},
+                present_state: State::Normal,
-            Human{present_state: State::Normal, x: 1, y: 2},
+                x: 1,
-            Human{present_state: State::Normal, x: 2, y: 2},
+                y: 0,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 2,
+                y: 0,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 0,
+                y: 1,
+            },
+            Human {
+                present_state: State::Infected,
+                x: 1,
+                y: 1,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 2,
+                y: 1,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 0,
+                y: 2,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 1,
+                y: 2,
+            },
+            Human {
+                present_state: State::Normal,
+                x: 2,
+                y: 2,
+            },
         ];
         stats = humans_stats(&population.humans);
         println!("stats after init: {:?}", stats);
         assert_eq!(stats.normal, 8);
 
-        // kill every one
+        // infect every one
         propagate_stats = population.propagate();
         stats = humans_stats(&population.humans);
         println!("propate_stats: {:?}", propagate_stats);
         println!("population: {:?}", stats);
         assert_eq!(propagate_stats, [8, 1, 0, 0]);
-        assert_eq!(stats.normal, 8);
+        assert_eq!(stats.normal, 0);
-        assert_eq!(stats.infected, 0);
+        assert_eq!(stats.infected, 8);
         assert_eq!(stats.immune, 1);
         assert_eq!(stats.dead, 0);
 
+        // cure every one
+        propagate_stats = population.propagate();
+        stats = humans_stats(&population.humans);
+        println!("propate_stats: {:?}", propagate_stats);
+        println!("population: {:?}", stats);
+        assert_eq!(propagate_stats, [0, 8, 1, 0]);
+        assert_eq!(stats.normal, 0);
+        assert_eq!(stats.infected, 0);
+        assert_eq!(stats.immune, 9);
+        assert_eq!(stats.dead, 0);
+
+        // then
         for _x in 0..100 {
             propagate_stats = population.propagate();
             stats = humans_stats(&population.humans);
             println!("propate_stats: {:?}", propagate_stats);
             println!("population: {:?}", stats);
-            assert_eq!(propagate_stats, [8, 0, 1, 0]);
+            assert_eq!(propagate_stats, [0, 0, 9, 0]);
-            assert_eq!(stats.normal, 8);
+            assert_eq!(stats.normal, 0);
             assert_eq!(stats.infected, 0);
-            assert_eq!(stats.immune, 1);
+            assert_eq!(stats.immune, 9);
             assert_eq!(stats.dead, 0);
         }
     }
 
+    #[parameterized(infection_start = {0, 50, 100})]
+    fn propagate_harmless(infection_start: i32) {
+        let disease: Disease = Disease::new(0, 0, 0, String::from("Harmless"));
+        let (width, height) = (100, 100);
+        let mut population: Population = Population::new(infection_start, 0, 0, width, height, disease);
+        let stats_before: Stats;
+        let stats_after: Stats;
+
+        stats_before = humans_stats(&population.humans);
+        let should_be_infected = population.size as i32 * infection_start / 100;
+        let tolerance = (should_be_infected as f32 * 0.20) as i32;
+        println!("{:?}", stats_before);
+        assert!(stats_before.infected <= should_be_infected + tolerance, "{} infected, should be less than {}", stats_before.infected, should_be_infected + tolerance);
+        assert!(stats_before.infected >= should_be_infected - tolerance, "{} infected, should be more than {}", stats_before.infected, should_be_infected - tolerance);
+
+        population.propagate();
+
+        stats_after = humans_stats(&population.humans);
+        assert_eq!(stats_before.infected, stats_after.infected, "no one should have been infected");
+    }
+
     #[parameterized(infection_rate = {0, 100, 0}, death_rate = {0, 0, 100}, infected_expected = {0, 1, 0})]
     fn propagate_test(infection_rate: i32, death_rate: i32, infected_expected: i32) {
         let disease: Disease;
@@ -598,9 +713,8 @@ mod tests {
         disease = Disease::new(infection_rate, 0, death_rate, String::from("Test"));
         population = Population::new(start_infected, 0, 0, width, height, disease);
 
-        population.generate();
         stats = humans_stats(&population.humans);
-        println!("Population after generate: {:?}", stats);
+        println!("Population after init: {:?}", stats);
 
         // total * proba - 20% < infected < total * proba + 20%
         let infected_at_start_proba = width * height * start_infected / 100;