mods in testing 1.16.0

[mothur.git] / linearalgebra.cpp

diff --git a/linearalgebra.cpp b/linearalgebra.cpp
index 27d35ac2cd7c40d4502aa469ffa62f2d01ff31ec..a2ee221b9d79c0997f25f6c9baf09cdd0b752bf2 100644 (file)
--- a/linearalgebra.cpp
+++ b/linearalgebra.cpp
@@ -372,6 +372,9 @@ double LinearAlgebra::calcPearson(vector< vector<double> >& euclidDists, vector<
                double denom = (sqrt(denomTerm1) * sqrt(denomTerm2));
                double r = numerator / denom;
                
+               //divide by zero error
+               if (isnan(r) || isinf(r)) { r = 0.0; }
+               
                return r;
                
        }