[petsc-users] algorithm to get pure real eigen valule for general eigenvalue problem(Non-hermitian type)

Jose E. Roman jroman at dsic.upv.es
Tue Apr 16 01:54:17 CDT 2013


El 15/04/2013, a las 23:18, Zhang Wei escribió:

> hi
> Thanks very much for your help!
> Indeed, I am dealing with linearised Navier Stokes.  Actually the algorithm is already tested in a FEM code by directly using Arpack. I tried slepc with wrapped Arpack. I noticed that the shit between each iteration is controlled by the number of converged eigenvalue and a predefined portion,which is the same as restart parameter "keep" in krylovschur solver, but can't be changed. While in the FEM code, the shift value is exact one. But I can't find different in the implementation. Could it be the problem?
> 
> Yours Sincerely
> ------------------------
> Wei Zhang
> Ph.D
> Hydrodynamic Group
> Dept. of Shipping and Marine Technology
> Chalmers University of Technology
> Sweden
> Phone:+46-31 772 2703

Which shift? Do you mean doing shift-and-invert? This has nothing to do with the restart parameter.

Jose



More information about the petsc-users mailing list