From 5b3bedc6fbfed63854cc6a7d47454580b2734c08 Mon Sep 17 00:00:00 2001
From: westcott <westcott>
Date: Mon, 10 Jan 2011 18:12:21 +0000
Subject: [PATCH] pca command

---
 linearalgebra.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

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));
-- 
2.39.2