Skip to content
Snippets Groups Projects
Commit ad13352c authored by Nicolas Marsic's avatar Nicolas Marsic
Browse files

fixing a wrong comment in beyn.py + better default setings

parent 8e6b33ea
No related branches found
No related tags found
No related merge requests found
......@@ -13,7 +13,7 @@ import numpy as np
import numpy.matlib
def simple(operator, origin, radius,
nNode=100, maxIt=10, lStart=1, lStep=1, rankTol=1e-4, verbose=True):
nNode=100, maxIt=10, lStart=4, lStep=3, rankTol=1e-4, verbose=True):
"""Solves an eigenvalue problem using Beyn's algorithm (simple version)
Keyword arguments:
......@@ -87,7 +87,7 @@ def simple(operator, origin, radius,
# Compute B
B = V0.H * A1 * W0 * S0Inv
# Eigenvalues & eigenvectors (by projecting QHat onto V0)
# Eigenvalues & eigenvectors
if(verbose): print "Solving linear EVP..."
myLambda, QHat = numpy.linalg.eig(B)
Q = V0 * QHat;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment