From 058e7d3d12e3abc0088af98e94e2e6caad703c96 Mon Sep 17 00:00:00 2001 From: depristo Date: Wed, 7 Apr 2010 00:21:35 +0000 Subject: [PATCH] Bug fix for Gregory git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@3127 348d0f76-0448-11de-a6fe-93d51630548a --- .../sting/gatk/walkers/genotyper/DiploidGenotypePriors.java | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) 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 +}