[petsc-users] Help with SLEPc eigenvectors convergence.
Eric Chamberland
Eric.Chamberland at giref.ulaval.ca
Wed Mar 6 13:57:43 CST 2024
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240306/64bcbfb5/attachment-0001.html>
-------------- next part --------------
1 EPS nconv=8 Values (Errors) 1.72355e-12 (0.00000000e+00) 4.83987e-12 (0.00000000e+00) 7.16714e-12 (0.00000000e+00) 199.286 (1.47893576e-39) 415.653 (1.99313185e-32) 570.22 (3.09568215e-29) 1294.72 (7.71230444e-19) 1295.67 (7.95775669e-19) 3769.79 (1.62599075e-08) 3771.16 (7.52750700e-09) 4323.24 (8.61648126e-09) 5141.1 (1.25820965e-08) 6023.68 (1.91065876e-05) 11729.6 (1.18367157e-02) 11845.2 (1.49674912e-02) 14148.3 (4.04353216e-02) 15071.8 (3.24127647e-02) 18157.1 (1.10833352e-01) 25128.7 (5.14116230e-02) 32499.7 (1.88695529e-01) 48189.7 (1.41699930e-01) 86602.1 (2.54292847e-01) 150776. (5.17558147e-01) 356904. (9.25768638e-01) 1.79343e+06 (1.80518276e+00)
2 EPS nconv=12 Values (Errors) 1.72355e-12 (0.00000000e+00) 4.83987e-12 (0.00000000e+00) 7.16714e-12 (0.00000000e+00) 199.286 (1.47893576e-39) 415.653 (1.99313185e-32) 570.22 (3.09568215e-29) 1294.72 (7.71230444e-19) 1295.67 (7.95775669e-19) 3769.79 (1.36933620e-19) 3771.16 (6.36276545e-20) 4323.24 (2.95572007e-19) 5141.1 (2.67431497e-18) 6023.68 (2.28621509e-14) 11727.3 (8.57510515e-08) 11841.5 (1.21364008e-07) 14111.7 (1.44584686e-05) 15046.9 (1.24205257e-03) 17784.4 (2.90225977e-03) 24962.4 (8.04532077e-02) 29636.1 (7.20333497e-02) 41881.1 (1.57409624e-01) 56572.2 (3.00061077e-01) 111415. (5.91435311e-01) 277740. (9.37846805e-01) 2.06838e+06 (2.35079865e+00)
EPS Object: SolveurEPSGen (options_eps_gen) 1 MPI process
type: krylovschur
50% of basis vectors kept after restart
using the locking variant
problem type: generalized symmetric eigenvalue problem
selected portion of the spectrum: closest to target: 0. (in magnitude)
postprocessing eigenvectors with purification
computing all residuals (for tracking convergence)
number of eigenvalues (nev): 10
number of column vectors (ncv): 25
maximum dimension of projected problem (mpd): 25
maximum number of iterations: 10000
tolerance: 1e-14
convergence test: relative to the eigenvalue
BV Object: (options_eps_gen) 1 MPI process
type: svec
26 columns of global length 150
vector orthogonalization method: modified Gram-Schmidt
orthogonalization refinement: if needed (eta: 0.7071)
block orthogonalization method: GS
non-standard inner product
tolerance for definite inner product: 2.22045e-15
inner product matrix:
Mat Object: (MatB) 1 MPI process
type: seqaij
rows=150, cols=150
total: nonzeros=4932, allocated nonzeros=4932
total number of mallocs used during MatSetValues calls=0
using I-node routines: found 50 nodes, limit used is 5
doing matmult as a single matrix-matrix product
DS Object: (options_eps_gen) 1 MPI process
type: hep
solving the problem with: Implicit QR method (_steqr)
ST Object: (options_eps_gen) 1 MPI process
type: sinvert
shift: 0.
number of matrices: 2
nonzero pattern of the matrices: UNKNOWN
KSP Object: (options_eps_genst_) 1 MPI process
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (options_eps_genst_) 1 MPI process
type: lu
out-of-place factorization
tolerance for zero pivot 2.22045e-14
matrix ordering: external
factor fill ratio given 0., needed 0.
Factored matrix follows:
Mat Object: (options_eps_genst_) 1 MPI process
type: mumps
rows=150, cols=150
package used to perform factorization: mumps
total: nonzeros=8028, allocated nonzeros=8028
MUMPS run parameters:
Use -options_eps_genst_ksp_view ::ascii_info_detail to display information for all processes
RINFOG(1) (global estimated flops for the elimination after analysis): 244937.
RINFOG(2) (global estimated flops for the assembly after factorization): 7740.
RINFOG(3) (global estimated flops for the elimination after factorization): 244937.
(RINFOG(12) RINFOG(13))*2^INFOG(34) (determinant): (0.,0.)*(2^0)
INFOG(3) (estimated real workspace for factors on all processors after analysis): 8028
INFOG(4) (estimated integer workspace for factors on all processors after analysis): 1236
INFOG(5) (estimated maximum front size in the complete tree): 54
INFOG(6) (number of nodes in the complete tree): 15
INFOG(7) (ordering option effectively used after analysis): 2
INFOG(8) (structural symmetry in percent of the permuted matrix after analysis): 100
INFOG(9) (total real/complex workspace to store the matrix factors after factorization): 8028
INFOG(10) (total integer space store the matrix factors after factorization): 1236
INFOG(11) (order of largest frontal matrix after factorization): 54
INFOG(12) (number of off-diagonal pivots): 0
INFOG(13) (number of delayed pivots after factorization): 0
INFOG(14) (number of memory compress after factorization): 0
INFOG(15) (number of steps of iterative refinement after solution): 0
INFOG(16) (estimated size (in MB) of all MUMPS internal data for factorization after analysis: value on the most memory consuming processor): 0
INFOG(17) (estimated size of all MUMPS internal data for factorization after analysis: sum over all processors): 0
INFOG(18) (size of all MUMPS internal data allocated during factorization: value on the most memory consuming processor): 0
INFOG(19) (size of all MUMPS internal data allocated during factorization: sum over all processors): 0
INFOG(20) (estimated number of entries in the factors): 8028
INFOG(21) (size in MB of memory effectively used during factorization - value on the most memory consuming processor): 0
INFOG(22) (size in MB of memory effectively used during factorization - sum over all processors): 0
INFOG(23) (after analysis: value of ICNTL(6) effectively used): 0
INFOG(24) (after analysis: value of ICNTL(12) effectively used): 1
INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
INFOG(28) (after factorization: number of null pivots encountered): 0
INFOG(29) (after factorization: effective number of entries in the factors (sum over all processors)): 8028
INFOG(30, 31) (after solution: size in Mbytes of memory used during solution phase): 0, 0
INFOG(32) (after analysis: type of analysis done): 1
INFOG(33) (value used for ICNTL(8)): 7
INFOG(34) (exponent of the determinant if determinant is requested): 0
INFOG(35) (after factorization: number of entries taking into account BLR factor compression - sum over all processors): 8028
INFOG(36) (after analysis: estimated size of all MUMPS internal data for running BLR in-core - value on the most memory consuming processor): 0
INFOG(37) (after analysis: estimated size of all MUMPS internal data for running BLR in-core - sum over all processors): 0
INFOG(38) (after analysis: estimated size of all MUMPS internal data for running BLR out-of-core - value on the most memory consuming processor): 0
INFOG(39) (after analysis: estimated size of all MUMPS internal data for running BLR out-of-core - sum over all processors): 0
linear system matrix = precond matrix:
Mat Object: (MatA) 1 MPI process
type: seqaij
rows=150, cols=150
total: nonzeros=4932, allocated nonzeros=4932
total number of mallocs used during MatSetValues calls=0
using I-node routines: found 50 nodes, limit used is 5
options_eps_gen Linear eigensolve converged (12 eigenpairs) due to CONVERGED_TOL; iterations 2
Problem: some of the first 10 relative errors are higher than the tolerance
-------------- next part --------------
A non-text attachment was scrubbed...
Name: matrice0_gen.m
Type: text/x-matlab
Size: 154373 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240306/64bcbfb5/attachment-0002.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: matrice1_gen.m
Type: text/x-matlab
Size: 153890 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240306/64bcbfb5/attachment-0003.bin>
More information about the petsc-users
mailing list