Skip to content

Commit

Permalink
Updated the reference to eigen, substituted for eigs.
Browse files Browse the repository at this point in the history
Checked that the results are the same when use_arpack=False, when I look at the
propagation constants for example 2 step_index.
  • Loading branch information
dbs660 committed Sep 22, 2015
1 parent 3e8077a commit a7d5551
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions Polymode/TriDiagonalSolver.py
Original file line number Diff line number Diff line change
@@ -39,7 +39,7 @@
use_arpack=True
if use_arpack:
#ARPACK interface - this keeps changing in scipy
from scipy.sparse.linalg.eigen.arpack import eigen
from scipy.sparse.linalg.eigen.arpack import eigs
from scipy.sparse.linalg.interface import LinearOperator
else:
from .mathlink.eigensolver import eigs
@@ -356,7 +356,7 @@ def calculate_one(self, evsearch, searchnum):

#Solve linear eigenproblem
if use_arpack:
evals, revecs = eigen(self.si, k=searchnum, \
evals, revecs = eigs(self.si, k=searchnum, \
which='LM', return_eigenvectors=True)
evals = self.si.eigenvalue_transform(evals)
else:

0 comments on commit a7d5551

Please sign in to comment.