[petsc-users] Incorrect eigenvalues

Hassan Md Mahmudulla mhassan at miners.utep.edu
Mon Jul 18 14:01:28 CDT 2016


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


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160718/2af9bd31/attachment.html>


More information about the petsc-users mailing list