#include "slepcnep.h" PetscErrorCode NEPSetInitialSpace(NEP nep,PetscInt n,Vec is[])Collective
nep | - the nonlinear eigensolver context | |
n | - number of vectors | |
is | - set of basis vectors of the initial space |
These vectors do not persist from one NEPSolve() call to the other, so the initial space should be set every time.
The vectors do not need to be mutually orthonormal, since they are explicitly orthonormalized internally.
Common usage of this function is when the user can provide a rough approximation of the wanted eigenspace. Then, convergence may be faster.