From ad13352c18b72ea8e12e95732cb988ea45404dfc Mon Sep 17 00:00:00 2001
From: Nicolas Marsic <marsic@temf.tu-darmstadt.de>
Date: Fri, 2 Jun 2017 18:55:43 +0200
Subject: [PATCH] fixing a wrong comment in beyn.py + better default setings
---
bin/beyn.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/bin/beyn.py b/bin/beyn.py
index 208f20c..0e1a34a 100644
--- a/bin/beyn.py
+++ b/bin/beyn.py
@@ -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;
--
GitLab