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

[Jose E. Roman]

El 15/04/2013, a las 11:28, Zhang Wei escribió:

> Hi
> Sorry I forgot to say it is matrix free case. So if I can export matrix like that? Or any another way to do so? 
> 
> Yours Sincerely
> ------------------------
> Wei Zhang
> Ph.D
> Hydrodynamic Group
> Dept. of Shipping and Marine Technology
> Chalmers University of Technology
> Sweden
> Phone:+46-31 772 2703

You could run the SLEPc program with "-eps_type lapack", then a dense matrix will be created using matrix-vector products (you will see it if you set the -mat_view_matlab flag). This is very rudimentary and only viable for very small dimension. A better alternative would be to call MatComputeExplicitOperator() from within your code, use a binary viewer to save the resulting matrix and load it in Matlab with PetscBinaryRead.m.

Jose