working on pam

[mothur.git] / simpson.cpp

diff --git a/simpson.cpp b/simpson.cpp
index bf3365ddeadd942a53733e85cbb1eeb30cdc0bb7..edc35bd6c04f42c3d687dce377df6b6442e36a18 100644 (file)
--- a/simpson.cpp
+++ b/simpson.cpp
@@ -28,13 +28,13 @@ EstOutput Simpson::getValues(SAbundVector* rank){
                if(sobs != 0){
                        double simnum=0.0000;
                
-                       for(unsigned long int i=1;i<=maxRank;i++){
+                       for(unsigned long long i=1;i<=maxRank;i++){
                                simnum += (double)(rank->get(i)*i*(i-1));
                        }
                        
                        simpson = simnum / (sampled*(sampled-1));
                
-                       for(unsigned long int i=1;i<=maxRank;i++){
+                       for(unsigned long long 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);