Fixed arpack wrapper

shogun-toolbox · Jan 26, 2012 · 73a2d9b · 73a2d9b
1 parent 4347344
commit 73a2d9b
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/src/shogun/mathematics/arpack.cpp b/src/shogun/mathematics/arpack.cpp
@@ -258,7 +258,7 @@ void arpack_dsxupd(double* matrix, double* rhs, bool is_rhs_diag, int n, int nev
 		{
 			if (mode==1)
 			{
-				if (cov) 
+				if (!cov) 
 				{
 					// compute (workd+ipntr[1]-1) = A*(workd+ipntr[0]-1)
 					cblas_dsymv(CblasColMajor,CblasUpper,
@@ -268,12 +268,12 @@ void arpack_dsxupd(double* matrix, double* rhs, bool is_rhs_diag, int n, int nev
 				}
 				else
 				{
-					cblas_dgemv(CblasColMajor,CblasNoTrans,
-					            n,n,1.0,matrix,n,
+					cblas_dsymv(CblasColMajor,CblasUpper,
+					            n,1.0,matrix,n,
 					            (workd+ipntr[0]-1),1,
 					            0.0,(workd+ipntr[1]-1),1);
-					cblas_dgemv(CblasColMajor,CblasTrans,
-					            n,n,1.0,matrix,n,
+					cblas_dsymv(CblasColMajor,CblasUpper,
+					            n,1.0,matrix,n,
 					            (workd+ipntr[1]-1),1,
 					            0.0,(workd+ipntr[0]-1),1);
 					cblas_dcopy(n,workd+ipntr[0]-1,1,workd+ipntr[1]-1,1);