Hi all, I have been trying to solve generalized eigenvalue problem using matrices of size 10K. Sparsity of the matrix is 6%. I am using the following command
./solver -f1 hamold.petsc -f2 ovlbaby.petsc -st_ksp_type preonly -st_pc_type jacobi -st_pc_factor_mat_solver_package mumps -eps_interval -2,0 -eps_nev 1000 * solver is the program * f1 and f2 are the input file for both matrices in petsc binary (mpiaij) I am getting the following output: Generalized eigenproblem stored in file. Reading REAL matrices from binary files... TYPE OF MATRIX A: mpiaij TYPE OF MATRIX B: mpiaij Solving for Eigen values... Solved! 1: -9771.8339 0 2: -9559.8347 0 3: -9408.5603 0 4: -9387.423 0 5: -9235.9137 0 6: -9102.5334 0 7: -9098.1307 0 8: -8970.3594 0 9: -8854.4964 0 10: -8850.3629 0 11: -8736.6619 0 12: -8637.1749 0 13: -8628.214 0 14: -8524.2494 0 15: -8440.0801 0 16: -8424.1789 0 17: -8327.5389 0 18: -8257.7763 0 19: -8233.9564 0 20: -8143.1251 0 21: -8086.9865 0 22: -8054.7899 0 23: -7968.7355 0 24: -7925.5421 0 25: -7884.7777 0 26: -7802.7577 0 27: -7771.913 0 28: -7722.537 0 29: -7643.9943 0 30: -7624.9684 0 ......................... ......................... 541: -24.947288 0 542: -24.945875 0 543: -24.94017 0 First column is the eigenvalues. My concern is, * Eigenvalues are not right * I defined the interval but still it's giving me eigenvalues outside of that interval Please help me out. M Hassan
