MCPcopy Create free account
hub / github.com/FEniCS/dolfinx / get_eigenpair

Method get_eigenpair

cpp/dolfinx/la/slepc.cpp:140–156  ·  view source on GitHub ↗

-----------------------------------------------------------------------------

Source from the content-addressed store, hash-verified

138}
139//-----------------------------------------------------------------------------
140void SLEPcEigenSolver::get_eigenpair(PetscScalar& lr, PetscScalar& lc, Vec r,
141 Vec c, int i) const
142{
143 assert(_eps);
144 PetscInt ii = static_cast<PetscInt>(i);
145
146 // Get number of computed eigenvectors/values
147 PetscInt num_computed_eigenvalues;
148 EPSGetConverged(_eps, &num_computed_eigenvalues);
149 if (ii < num_computed_eigenvalues)
150 EPSGetEigenpair(_eps, ii, &lr, &lc, r, c);
151 else
152 {
153 throw std::runtime_error("Requested eigenpair (" + std::to_string(i)
154 + ") has not been computed");
155 }
156}
157//-----------------------------------------------------------------------------
158std::int64_t SLEPcEigenSolver::get_number_converged() const
159{

Callers

nothing calls this directly

Calls

no outgoing calls

Tested by

no test coverage detected