#include "slepcds.h" PetscErrorCode DSVectors(DS ds,DSMatType mat,PetscInt *j,PetscReal *rnorm)Logically Collective
ds | - the direct solver context | |
mat | - the matrix, used to indicate which vectors are required |
j | - (optional) index of vector to be computed |
rnorm | - (optional) computed residual norm |
If NULL is passed in argument j then all vectors are computed, otherwise j indicates which vector must be computed. In real non-symmetric problems, on exit the index j will be incremented when a complex conjugate pair is found.
This function can be invoked after the dense problem has been solved, to get the residual norm estimate of the associated Ritz pair. In that case, the relevant information is returned in rnorm.
For computing eigenvectors, LAPACK's _trevc is used so the matrix must be in (quasi-)triangular form, or call DSSolve() first.