Hi
I'm trying to use slepc wrapper to solve the generalized eigen value problem
(schrodinger equation). I tried to use Jacobi Davidson method. The piece of
code is as follows
* SolverControl solver_control (dof_handler.n_dofs(), 1e-10);
SLEPcWrappers::SolverDavidson eigensolver (solver_control);
eigensolver.set_which_eigenpairs (EPS_SMALLEST_REAL);
eigensolver.solve (system_matrix_psi, mass_matrix_psi, eigenvalues,
eigenfunctions, eigenfunctions.size());
Upsi=eigenfunctions[0];
*However the problem is i need to specify EPSProblemType to be EPS_GHEP
since only this option helps to get an M orthonormal eigen vectors. The
issue is i currently don't know of an option in the slepc wrappers to
specify this input. It will be of great help if you can help me figuring
this issue out.
Regards
Ram
<http://sitemaker.umich.edu/ram/home/>
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
