phynexis-deps/arpack-ng
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix solver parameters for pyarpackSparse{LDLT/LU}

Browse Source
This commit is contained in:
Fabien Péan
2023-09-14 12:12:59 +02:00
committed by Franck HOUSSEN
parent 447b4387e9
commit 93c048fae3
2 changed files with 9 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files
EXAMPLES/PYARPACK/pyarpackSparseLDLT.py.in
+6 -2
View File
@@ -32,7 +32,9 @@ arpackSlv.nbEV = 1
arpackSlv.nbCV = 2*arpackSlv.nbEV + 1
arpackSlv.mag = 'LM'
arpackSlv.maxIt = 200
arpackSlv.slvPvtThd = 1.e-6
arpackSlv.slvOffset = 0.
arpackSlv.slvScale = 1.
arpackSlv.symPb = True
# Solve eigen problem.
@@ -92,8 +94,10 @@ arpackSlv.nbEV = 2
arpackSlv.nbCV = 2*arpackSlv.nbEV + 1
arpackSlv.mag = 'LM'
arpackSlv.maxIt = 200
arpackSlv.slvPvtThd = 1.e-6
arpackSlv.slvOffset = 0.
arpackSlv.slvScale = 1.
arpackSlv.sigmaReal = 1
arpackSlv.symPb = True
# Solve eigen problem.
EXAMPLES/PYARPACK/pyarpackSparseLU.py.in
+3 -5
View File
@@ -33,8 +33,7 @@ arpackSlv.nbEV = 1
arpackSlv.nbCV = 2*arpackSlv.nbEV + 1
arpackSlv.mag = 'LM'
arpackSlv.maxIt = 200
arpackSlv.slvOffset = 0.
arpackSlv.slvScale = 1.
arpackSlv.slvPvtThd = 1.e-6
arpackSlv.symPb = False
# Solve eigen problem.
@@ -98,8 +97,7 @@ arpackSlv.nbEV = 2
arpackSlv.nbCV = 2*arpackSlv.nbEV + 1
arpackSlv.mag = 'LM'
arpackSlv.maxIt = 200
arpackSlv.slvOffset = 0.
arpackSlv.slvScale = 1.
arpackSlv.slvPvtThd = 1.e-6
arpackSlv.sigmaReal = 1
arpackSlv.sigmaImag = 1
arpackSlv.symPb = False
@@ -122,4 +120,4 @@ for val, vec in zip(arpackSlv.val, arpackSlv.vec):
print("eigen value:", val)
print("eigen vector:")
for v in range(n):
print(vec[v])
print(vec[v])