From: pschloss Date: Thu, 12 Aug 2010 13:04:22 +0000 (+0000) Subject: fixed int bug in simpson.cpp X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=commitdiff_plain;h=afed323beddfc099c365104b9071600aea9731bd fixed int bug in simpson.cpp --- diff --git a/simpson.cpp b/simpson.cpp index fa5b435..bf3365d 100644 --- a/simpson.cpp +++ b/simpson.cpp @@ -28,13 +28,13 @@ EstOutput Simpson::getValues(SAbundVector* rank){ if(sobs != 0){ double simnum=0.0000; - for(int i=1;i<=maxRank;i++){ + for(unsigned long int i=1;i<=maxRank;i++){ simnum += (double)(rank->get(i)*i*(i-1)); } simpson = simnum / (sampled*(sampled-1)); - for(int i=1;i<=maxRank;i++){ + for(unsigned long int i=1;i<=maxRank;i++){ double piI = (double) i / (double)sampled; firstTerm += rank->get(i) * pow(piI, 3); secondTerm += rank->get(i) * pow(piI, 2);