[Bf-blender-cvs] SVN commit: /data/svn/bf-blender [39092] branches/soc-2011-avocado/blender/ intern/autoseam/EigenSolverArpack.cpp: Tuning of arpack eigen calculation parameters.

[shuvro sarker]

Revision: 39092
          http://projects.blender.org/scm/viewvc.php?view=rev&root=bf-blender&revision=39092
Author:   shuvro
Date:     2011-08-06 05:00:40 +0000 (Sat, 06 Aug 2011)
Log Message:
-----------
Tuning of arpack eigen calculation parameters.

Modified Paths:
--------------
    branches/soc-2011-avocado/blender/intern/autoseam/EigenSolverArpack.cpp

Modified: branches/soc-2011-avocado/blender/intern/autoseam/EigenSolverArpack.cpp
===================================================================
--- branches/soc-2011-avocado/blender/intern/autoseam/EigenSolverArpack.cpp	2011-08-06 04:41:19 UTC (rev 39091)
+++ branches/soc-2011-avocado/blender/intern/autoseam/EigenSolverArpack.cpp	2011-08-06 05:00:40 UTC (rev 39092)
@@ -39,7 +39,8 @@
 
 
 //#define EIGEN_DEBUG_ARPACK 1
-#define NUM_EIGEN_VAL 100
+//#define NUM_EIGEN_VAL 100
+#define NUM_EIGEN_VAL 40
 #define minimum(a,b)  (a <= b) ? a : b;
 
 extern "C"{
@@ -137,7 +138,11 @@
 	nev = minimum(matrix_dimension - 1, NUM_EIGEN_VAL);
 	//num_eigen_values = nev;
 	num_eigen_vectors = nev;
-	ncv = minimum(nev + 1, n);
+	//ncv = minimum(nev + 1, n);
+	ncv = nev*2.5;
+	if(ncv > matrix_dimension) { ncv = matrix_dimension ; }
+	if(nev > matrix_dimension) { nev = matrix_dimension ; }
+	if(nev + 2 > ncv) {  nev = ncv - 2 ;   } 
 	
 	
 	*(unsigned char *)bmat = 'I';
@@ -175,7 +180,7 @@
 	
 	
 	lworkl = ncv * (ncv + 8);
-	tol = 0.;
+	tol = 0.00000075;
 	info = 0;
 	ido = 0;