fixed logic error in generateMutations

FaridKaveh · Mar 14, 2024 · 0acfe86 · 0acfe86
1 parent 4d3e766
commit 0acfe86
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/MoranModel/MoranModel.cpp b/MoranModel/MoranModel.cpp
@@ -67,13 +67,13 @@ void MoranProcess::generateMuts(){
     }
 
     mutations.insert(mutations.end(), population*events, 0); 
-    std::uniform_int_distribution<> pick_line (0,population); 
+    std::uniform_int_distribution<> pick_line (0,population-1); 
     int line; 
 
     for (int i = 0; i < events; ++i) { 
         while (allocations.at(i) > 0){
         line = pick_line(engine);
-        ++mutations.at(i+line);
+        ++mutations.at(population*i+line);
 
         --allocations.at(i);
         }