# HG changeset patch
# User Sean Farley <sean.michael.farley@gmail.com>
# Date 1354650287 21600
# Node ID 5687b9076dd8350facc03263810765de6079924b
# Parent  87072fd48e1569d5db45aa55d5fccca7c020f01a
update: rename eigen_symmetric to eigsh

diff --git a/binfile_parcellation.py b/binfile_parcellation.py
--- a/binfile_parcellation.py
+++ b/binfile_parcellation.py
@@ -32,11 +32,11 @@
 # distributed with this script to be installed in a directory that is
 # accessible through PythonPath (the current directory will do for python_ncult_lib.py).
 import time as time
 from numpy import *
 from scipy.sparse import csc_matrix
-from scipy.sparse.linalg.eigen.arpack import eigen_symmetric
+from scipy.sparse.linalg.eigen.arpack import eigsh
 from scipy.linalg import norm, svd
 from python_ncut_lib import *
 
 # binfile_parcellate( infile, outfile, K ):
 #
diff --git a/group_mean_binfile_parcellation.py b/group_mean_binfile_parcellation.py
--- a/group_mean_binfile_parcellation.py
+++ b/group_mean_binfile_parcellation.py
@@ -34,11 +34,11 @@
 # distributed with this script to be installed in a directory that is
 # accessible through PythonPath (the current directory will do for python_ncult_lib.py).
 import time as time
 from numpy import *
 from scipy.sparse import csc_matrix
-from scipy.sparse.linalg.eigen.arpack import eigen_symmetric
+from scipy.sparse.linalg.eigen.arpack import eigsh
 from scipy.linalg import norm, svd
 from python_ncut_lib import *
 
 
 # group_mean_binfile_parcellate( infiles, outfile, K ):
diff --git a/make_local_connectivity_scorr.py b/make_local_connectivity_scorr.py
--- a/make_local_connectivity_scorr.py
+++ b/make_local_connectivity_scorr.py
@@ -35,11 +35,11 @@
 # accessible through PythonPath 
 from nifti import *
 from numpy import array
 from scipy import *
 from scipy.sparse import *
-from scipy.sparse.linalg.eigen.arpack import eigen_symmetric
+from scipy.sparse.linalg.eigen.arpack import eigsh
 from scipy.linalg import norm, svd
 from numpy import corrcoef
 
 # simple function to translate 1D vector coordinates to 3D matrix coordinates, for a 3D matrix of 
 # size sz
diff --git a/make_local_connectivity_tcorr.py b/make_local_connectivity_tcorr.py
--- a/make_local_connectivity_tcorr.py
+++ b/make_local_connectivity_tcorr.py
@@ -35,11 +35,11 @@
 # accessible through PythonPath 
 from nifti import *
 from numpy import array
 from scipy import *
 from scipy.sparse import *
-from scipy.sparse.linalg.eigen.arpack import eigen_symmetric
+from scipy.sparse.linalg.eigen.arpack import eigsh
 from scipy.linalg import norm, svd
 from numpy import corrcoef
 
 # simple function to translate 1D vector coordinates to 3D matrix coordinates, for a 3D matrix of 
 # size sz
diff --git a/python_ncut_lib.py b/python_ncut_lib.py
--- a/python_ncut_lib.py
+++ b/python_ncut_lib.py
@@ -41,11 +41,11 @@
 from numpy import array,reshape,shape,matrix,ones,zeros,sqrt,sort,arange
 from numpy import nonzero,fromfile,tile,append,prod,double,argsort,sign
 from numpy import kron,multiply,divide,abs,reshape,asarray
 from scipy import rand
 from scipy.sparse import csc_matrix, spdiags
-from scipy.sparse.linalg.eigen.arpack import eigen_symmetric
+from scipy.sparse.linalg.eigen.arpack import eigsh
 from scipy.linalg import norm, svd, LinAlgError
 
 # exception hander for singular value decomposition
 class SVDError(Exception):
 	def __init__(self,value):
@@ -97,11 +97,11 @@
 	W=W+spdiags(dr,[0],m,m,"csc")
 	Dinvsqrt=spdiags((1.0/sqrt(d+eps)),[0],m,m,"csc")
 	P=Dinvsqrt*(W*Dinvsqrt);
 	
 	# perform the eigen decomposition
-	eigen_val,eigen_vec=eigen_symmetric(P,nbEigenValues,maxiter=maxiterations,tol=eigsErrorTolerence,which='LA')
+	eigen_val,eigen_vec=eigsh(P,nbEigenValues,maxiter=maxiterations,tol=eigsErrorTolerence,which='LA')
 	
 	# sort the eigen_vals so that the first
 	# is the largest
 	i=argsort(-eigen_val)
 	eigen_val=eigen_val[i]