X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=linearalgebra.cpp;h=192701af4738ffebc9190a54e5a4da1095679bfc;hb=5b3bedc6fbfed63854cc6a7d47454580b2734c08;hp=5de5a71dc35c1ab8888ff28ab8d7e92107f8cc72;hpb=37eac2026d91179acda0494c4dcca22f176551b9;p=mothur.git

diff --git a/linearalgebra.cpp b/linearalgebra.cpp
index 5de5a71..192701a 100644
--- a/linearalgebra.cpp
+++ b/linearalgebra.cpp
@@ -27,7 +27,7 @@ vector<vector<double> > LinearAlgebra::matrix_mult(vector<vector<double> > first
 		
 		product.resize(first_rows);
 		for(int i=0;i<first_rows;i++){
-			product[i].resize(first_cols);
+			product[i].resize(second_cols);
 		}
 		
 		for(int i=0;i<first_rows;i++){
@@ -172,7 +172,7 @@ int LinearAlgebra::qtli(vector<double>& d, vector<double>& e, vector<vector<doub
 					if(fabs(e[myM])+dd == dd) break;
 				}
 				if(myM != l){
-					if(iter++ == 30) cerr << "Too many iterations in tqli\n";
+					if(iter++ == 3000) cerr << "Too many iterations in tqli\n";
 					g = (d[l+1]-d[l]) / (2.0 * e[l]);
 					r = pythag(g, 1.0);
 					g = d[myM] - d[l] + e[l] / (g + SIGN(r,g));