diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/DiploidGenotypePriors.java b/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/DiploidGenotypePriors.java index 4a8483b40..be8ca8316 100755 --- a/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/DiploidGenotypePriors.java +++ b/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/DiploidGenotypePriors.java @@ -216,8 +216,8 @@ public class DiploidGenotypePriors { double pHet = heterozygosity2HetProbability(heterozyosity); double pHomVar = heterozygosity2HomVarProbability(heterozyosity); - if ((pHomRef + pHet + pHomVar) != 1) { - throw new RuntimeException(String.format("Prior probabilities don't sum to one => %f, %f, %f", pHomRef, pHet, pHomVar)); + if (MathUtils.compareDoubles(pHomRef + pHet + pHomVar, 1.0) != 0) { + throw new RuntimeException(String.format("BUG: Prior probabilities don't sum to one => %f, %f, %f", pHomRef, pHet, pHomVar)); } double[] priors = new double[DiploidGenotype.values().length]; @@ -288,4 +288,4 @@ public class DiploidGenotypePriors { flatPriors[g.ordinal()] = Math.log10(1.0 / DiploidGenotype.values().length); } } -} \ No newline at end of file +}