fix duplication

src/libaevol/7/DnaMutator.cpp

@@ -79,12 +79,10 @@ std::unique_ptr<MutationEvent> DnaMutator::generate_large_duplication_event(Jump
   }
 
   Dna_7::size_type pos_1, pos_2, pos_3;
-  pos_1 = prng.random(genome_size);
-  pos_2 = prng.random(genome_size);
-  // TODO: why are complete duplications forbidden ?
-  while (pos_2 == pos_1) {
-      pos_2 = prng.random(genome_size);
-  }
+  pos_1 = prng.random(genome_size + (exp_setup->linear_chrsm() ? 1 : 0));
+  do { // TODO: why are WGD forbidden ?
+    pos_2 = prng.random(genome_size + (exp_setup->linear_chrsm() ? 1 : 0));
+  } while (pos_2 == pos_1);
   pos_3 = prng.random(exp_setup->linear_chrsm() ? genome_size + 1 : genome_size);
 
   Dna_7::size_type seqlen = 0;