42#ifndef ANASAZI_RTR_SOLMGR_HPP
43#define ANASAZI_RTR_SOLMGR_HPP
69#include "Teuchos_TimeMonitor.hpp"
70#include "Teuchos_FancyOStream.hpp"
83template<
class ScalarType,
class MV,
class OP>
89 typedef Teuchos::ScalarTraits<ScalarType> SCT;
90 typedef typename Teuchos::ScalarTraits<ScalarType>::magnitudeType MagnitudeType;
91 typedef Teuchos::ScalarTraits<MagnitudeType> MT;
117 Teuchos::ParameterList &pl );
136 Teuchos::Array<Teuchos::RCP<Teuchos::Time> >
getTimers()
const {
137 return Teuchos::tuple(_timerSolve);
163 Teuchos::RCP<Eigenproblem<ScalarType,MV,OP> > problem_;
168 MagnitudeType convtol_;
176 Teuchos::RCP<Teuchos::Time> _timerSolve;
177 Teuchos::RCP<OutputManager<ScalarType> > printer_;
178 Teuchos::ParameterList& pl_;
183template<
class ScalarType,
class MV,
class OP>
186 Teuchos::ParameterList &pl ) :
190 convtol_(MT::prec()),
194#ifdef ANASAZI_TEUCHOS_TIME_MONITOR
195 _timerSolve(Teuchos::TimeMonitor::getNewTimer(
"Anasazi: RTRSolMgr::solve()")),
199 TEUCHOS_TEST_FOR_EXCEPTION(problem_ == Teuchos::null, std::invalid_argument,
"Problem not given to solver manager.");
200 TEUCHOS_TEST_FOR_EXCEPTION(!problem_->isProblemSet(), std::invalid_argument,
"Problem not set.");
201 TEUCHOS_TEST_FOR_EXCEPTION(!problem_->isHermitian(), std::invalid_argument,
"Problem not symmetric.");
202 TEUCHOS_TEST_FOR_EXCEPTION(problem_->getInitVec() == Teuchos::null,std::invalid_argument,
"Problem does not contain initial vectors to clone from.");
206 whch_ = pl_.get(
"Which",
"SR");
207 TEUCHOS_TEST_FOR_EXCEPTION(whch_ !=
"SR" && whch_ !=
"LR",
208 std::invalid_argument,
"Anasazi::RTRSolMgr: Invalid sorting string. RTR solvers compute only LR or SR.");
211 convtol_ = pl_.get(
"Convergence Tolerance",convtol_);
212 relconvtol_ = pl_.get(
"Relative Convergence Tolerance",relconvtol_);
213 strtmp = pl_.get(
"Convergence Norm",std::string(
"2"));
215 convNorm_ = RES_2NORM;
217 else if (strtmp ==
"M") {
218 convNorm_ = RES_ORTH;
221 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::invalid_argument,
222 "Anasazi::RTRSolMgr: Invalid Convergence Norm.");
227 maxIters_ = pl_.get(
"Maximum Iterations",maxIters_);
230 skinny_ = pl_.get(
"Skinny Solver",skinny_);
233 numICGS_ = pl_.get(
"Num ICGS",2);
236 blkSize_ = pl_.get(
"Block Size", problem_->getNEV());
240 int osProc = pl.get(
"Output Processor", 0);
243 Teuchos::RCP<Teuchos::FancyOStream> osp;
245 if (pl.isParameter(
"Output Stream")) {
246 osp = Teuchos::getParameter<Teuchos::RCP<Teuchos::FancyOStream> >(pl,
"Output Stream");
253 if (pl_.isParameter(
"Verbosity")) {
254 if (Teuchos::isParameterType<int>(pl_,
"Verbosity")) {
255 verbosity = pl_.get(
"Verbosity", verbosity);
257 verbosity = (int)Teuchos::getParameter<Anasazi::MsgType>(pl_,
"Verbosity");
266template<
class ScalarType,
class MV,
class OP>
274 const int nev = problem_->getNEV();
280 Teuchos::RCP<Teuchos::FancyOStream>
281 out = Teuchos::getFancyOStream(Teuchos::rcpFromRef(printer_->stream(
Debug)));
282 out->setShowAllFrontMatter(
false).setShowProcRank(
true);
283 *out <<
"Entering Anasazi::RTRSolMgr::solve()\n";
293 Teuchos::RCP<StatusTest<ScalarType,MV,OP> > maxtest;
294 Teuchos::RCP<StatusTest<ScalarType,MV,OP> > ordertest;
295 Teuchos::RCP<StatusTest<ScalarType,MV,OP> > combotest;
296 Teuchos::RCP<StatusTest<ScalarType,MV,OP> > convtest;
302 maxtest = Teuchos::null;
310 Teuchos::Array<Teuchos::RCP<StatusTest<ScalarType,MV,OP> > > alltests;
311 alltests.push_back(ordertest);
312 if (maxtest != Teuchos::null) alltests.push_back(maxtest);
316 Teuchos::RCP<StatusTestOutput<ScalarType,MV,OP> > outputtest;
317 if ( printer_->isVerbosity(
Debug) ) {
326 Teuchos::RCP<ICGSOrthoManager<ScalarType,MV,OP> > ortho
332 bool leftMost =
true;
333 if (whch_ ==
"LR" || whch_ ==
"LM") {
336 pl_.set<
bool>(
"Leftmost",leftMost);
337 Teuchos::RCP<RTRBase<ScalarType,MV,OP> > rtr_solver;
338 if (skinny_ ==
false) {
347 Teuchos::RCP< const MV > probauxvecs = problem_->getAuxVecs();
348 if (probauxvecs != Teuchos::null) {
349 rtr_solver->setAuxVecs( Teuchos::tuple< Teuchos::RCP<const MV> >(probauxvecs) );
352 TEUCHOS_TEST_FOR_EXCEPTION(rtr_solver->getBlockSize() < problem_->getNEV(),std::logic_error,
353 "Anasazi::RTRSolMgr requires block size >= requested number of eigenvalues.");
356 Teuchos::RCP<MV> foundvecs;
357 std::vector<MagnitudeType> foundvals;
361#ifdef ANASAZI_TEUCHOS_TIME_MONITOR
362 Teuchos::TimeMonitor slvtimer(*_timerSolve);
364 rtr_solver->iterate();
365 numIters_ = rtr_solver->getNumIters();
367 catch (
const std::exception &e) {
369 printer_->stream(
Anasazi::Errors) <<
"Exception: " << e.what() << endl;
372 problem_->setSolution(sol);
377 if (convtest->getStatus() ==
Passed || (maxtest != Teuchos::null && maxtest->getStatus() ==
Passed))
379 int num = convtest->howMany();
381 std::vector<int> ind = convtest->whichVecs();
383 foundvecs = MVT::CloneCopy(*rtr_solver->getRitzVectors(),ind);
385 foundvals.resize(num);
386 std::vector<Value<ScalarType> > all = rtr_solver->getRitzValues();
387 for (
int i=0; i<num; i++) {
388 foundvals[i] = all[ind[i]].realpart;
394 TEUCHOS_TEST_FOR_EXCEPTION(
true,std::logic_error,
"Anasazi::RTRSolMgr::solve(): solver returned without satisfy status test.");
400 sol.
Evecs = foundvecs;
403 for (
int i=0; i<sol.
numVecs; i++) {
404 sol.
Evals[i].realpart = foundvals[i];
407 sol.
index.resize(numfound,0);
410 rtr_solver->currentStatus(printer_->stream(
FinalSummary));
413#ifdef ANASAZI_TEUCHOS_TIME_MONITOR
415 Teuchos::TimeMonitor::summarize( printer_->stream(
TimingDetails ) );
420 problem_->setSolution(sol);
421 printer_->stream(
Debug) <<
"Returning " << sol.
numVecs <<
" eigenpairs to eigenproblem." << endl;
Basic implementation of the Anasazi::SortManager class.
Anasazi header file which uses auto-configuration information to include necessary C++ headers.
Abstract base class which defines the interface required by an eigensolver and status test class to c...
Basic implementation of the Anasazi::OrthoManager class.
Abstract class definition for Anasazi Output Managers.
Abstract class definition for Anasazi output stream.
Pure virtual base class which describes the basic interface for a solver manager.
Class which provides internal utilities for the Anasazi solvers.
Status test for forming logical combinations of other status tests.
Status test for testing the number of iterations.
Special StatusTest for printing status tests.
A status test for testing the norm of the eigenvectors residuals.
A status test for testing the norm of the eigenvectors residuals along with a set of auxiliary eigenv...
Types and exceptions used within Anasazi solvers and interfaces.
An implementation of the Anasazi::SortManager that performs a collection of common sorting techniques...
This class defines the interface required by an eigensolver and status test class to compute solution...
An implementation of the Anasazi::GenOrthoManager that performs orthogonalization using iterated clas...
Traits class which defines basic operations on multivectors.
Virtual base class which defines basic traits for the operator type.
Output managers remove the need for the eigensolver to know any information about the required output...
The Anasazi::RTRSolMgr provides a simple solver manager over the RTR eigensolver. For more informatio...
ReturnType solve()
This method performs possibly repeated calls to the underlying eigensolver's iterate() routine until ...
RTRSolMgr(const Teuchos::RCP< Eigenproblem< ScalarType, MV, OP > > &problem, Teuchos::ParameterList &pl)
Basic constructor for RTRSolMgr.
int getNumIters() const
Get the iteration count for the most recent call to solve.
Teuchos::Array< Teuchos::RCP< Teuchos::Time > > getTimers() const
Return the timers for this object.
const Eigenproblem< ScalarType, MV, OP > & getProblem() const
Return the eigenvalue problem.
virtual ~RTRSolMgr()
Destructor.
The Anasazi::SolverManager is a templated virtual base class that defines the basic interface that an...
Status test for forming logical combinations of other status tests.
A status test for testing the number of iterations.
A special StatusTest for printing other status tests.
A status test for testing the norm of the eigenvectors residuals.
A status test for testing the norm of the eigenvectors residuals along with a set of auxiliary eigenv...
Namespace Anasazi contains the classes, structs, enums and utilities used by the Anasazi package.
ReturnType
Enumerated type used to pass back information from a solver manager.
ResType
Enumerated type used to specify which residual norm used by residual norm status tests.
Struct for storing an eigenproblem solution.
Teuchos::RCP< MV > Evecs
The computed eigenvectors.
std::vector< int > index
An index into Evecs to allow compressed storage of eigenvectors for real, non-Hermitian problems.
int numVecs
The number of computed eigenpairs.
Teuchos::RCP< MV > Espace
An orthonormal basis for the computed eigenspace.
std::vector< Value< ScalarType > > Evals
The computed eigenvalues.
Output managers remove the need for the eigensolver to know any information about the required output...