Simplify num_genotypes expression

luntergroup · Sep 29, 2020 · b83ce11 · b83ce11
1 parent b82d319
commit b83ce11
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/core/types/genotype.cpp b/src/core/types/genotype.cpp
@@ -21,7 +21,7 @@ Genotype<Haplotype> remap(const Genotype<Haplotype>& genotype, const GenomicRegi
 
 std::size_t num_genotypes(const unsigned num_elements, const unsigned ploidy)
 {
-    return boost::math::binomial_coefficient<double>(num_elements + ploidy - 1, num_elements - 1);
+    return boost::math::binomial_coefficient<double>(num_elements + ploidy - 1, ploidy);
 }
 
 boost::optional<std::size_t> num_genotypes_noexcept(const unsigned num_elements, const unsigned ploidy) noexcept