29 #ifndef ANASAZI_RTR_SOLMGR_HPP 30 #define ANASAZI_RTR_SOLMGR_HPP 56 #include <Teuchos_FancyOStream.hpp> 69 template<
class ScalarType,
class MV,
class OP>
120 return Teuchos::tuple(_timerSolve);
151 MagnitudeType convtol_;
165 template<
class ScalarType,
class MV,
class OP>
172 convtol_(
MT::prec()),
176 #ifdef ANASAZI_TEUCHOS_TIME_MONITOR
177 _timerSolve(
Teuchos::TimeMonitor::getNewTimer(
"Anasazi: RTRSolMgr::solve()")),
184 TEUCHOS_TEST_FOR_EXCEPTION(problem_->getInitVec() == Teuchos::null,std::invalid_argument,
"Problem does not contain initial vectors to clone from.");
188 whch_ = pl_.
get(
"Which",
"SR");
190 std::invalid_argument,
"Anasazi::RTRSolMgr: Invalid sorting string. RTR solvers compute only LR or SR.");
193 convtol_ = pl_.
get(
"Convergence Tolerance",convtol_);
194 relconvtol_ = pl_.
get(
"Relative Convergence Tolerance",relconvtol_);
195 strtmp = pl_.
get(
"Convergence Norm",std::string(
"2"));
197 convNorm_ = RES_2NORM;
199 else if (strtmp ==
"M") {
200 convNorm_ = RES_ORTH;
204 "Anasazi::RTRSolMgr: Invalid Convergence Norm.");
209 maxIters_ = pl_.
get(
"Maximum Iterations",maxIters_);
212 skinny_ = pl_.
get(
"Skinny Solver",skinny_);
215 numICGS_ = pl_.
get(
"Num ICGS",2);
218 std::string fntemplate =
"";
219 bool allProcs =
false;
221 if (Teuchos::isParameterType<bool>(pl_,
"Output on all processors")) {
222 allProcs = pl_.
get(
"Output on all processors",allProcs);
224 allProcs = ( Teuchos::getParameter<int>(pl_,
"Output on all processors") != 0 );
227 fntemplate = pl_.
get(
"Output filename template",fntemplate);
232 MPI_Initialized(&mpiStarted);
233 if (mpiStarted) MPI_Comm_rank(MPI_COMM_WORLD, &MyPID);
238 if (fntemplate !=
"") {
239 std::ostringstream MyPIDstr;
243 while ( (pos = fntemplate.find(
"%d",start)) != -1 ) {
244 fntemplate.replace(pos,2,MyPIDstr.str());
249 if (fntemplate !=
"") {
250 osp =
Teuchos::rcp(
new std::ofstream(fntemplate.c_str(),std::ios::out | std::ios::app) );
252 osp = Teuchos::rcpFromRef(std::cout);
253 std::cout <<
"Anasazi::RTRSolMgr::constructor(): Could not open file for write: " << fntemplate << std::endl;
257 osp = Teuchos::rcpFromRef(std::cout);
262 if (Teuchos::isParameterType<int>(pl_,
"Verbosity")) {
263 verbosity = pl_.
get(
"Verbosity", verbosity);
265 verbosity = (int)Teuchos::getParameter<Anasazi::MsgType>(pl_,
"Verbosity");
280 template<
class ScalarType,
class MV,
class OP>
288 const int nev = problem_->getNEV();
295 out = Teuchos::getFancyOStream(Teuchos::rcpFromRef(printer_->stream(
Debug)));
297 *out <<
"Entering Anasazi::RTRSolMgr::solve()\n";
316 maxtest = Teuchos::null;
326 if (maxtest != Teuchos::null) alltests.
push_back(maxtest);
331 if ( printer_->isVerbosity(
Debug) ) {
346 bool leftMost =
true;
347 if (whch_ ==
"LR" || whch_ ==
"LM") {
350 pl_.set<
bool>(
"Leftmost",leftMost);
352 if (skinny_ ==
false) {
362 if (probauxvecs != Teuchos::null) {
367 "Anasazi::RTRSolMgr requires block size >= requested number of eigenvalues.");
371 std::vector<MagnitudeType> foundvals;
375 #ifdef ANASAZI_TEUCHOS_TIME_MONITOR 378 rtr_solver->iterate();
379 numIters_ = rtr_solver->getNumIters();
381 catch (
const std::exception &e) {
383 printer_->stream(
Anasazi::Errors) <<
"Exception: " << e.what() << endl;
386 problem_->setSolution(sol);
391 if (convtest->getStatus() ==
Passed || (maxtest != Teuchos::null && maxtest->getStatus() ==
Passed))
393 int num = convtest->howMany();
395 std::vector<int> ind = convtest->whichVecs();
397 foundvecs = MVT::CloneCopy(*rtr_solver->getRitzVectors(),ind);
399 foundvals.resize(num);
400 std::vector<Value<ScalarType> > all = rtr_solver->getRitzValues();
401 for (
int i=0; i<num; i++) {
402 foundvals[i] = all[ind[i]].realpart;
408 TEUCHOS_TEST_FOR_EXCEPTION(
true,std::logic_error,
"Anasazi::RTRSolMgr::solve(): solver returned without satisfy status test.");
414 sol.
Evecs = foundvecs;
417 for (
int i=0; i<sol.
numVecs; i++) {
418 sol.
Evals[i].realpart = foundvals[i];
421 sol.
index.resize(numfound,0);
424 rtr_solver->currentStatus(printer_->stream(
FinalSummary));
427 #ifdef ANASAZI_TEUCHOS_TIME_MONITOR 434 problem_->setSolution(sol);
435 printer_->stream(
Debug) <<
"Returning " << sol.
numVecs <<
" eigenpairs to eigenproblem." << endl;
Pure virtual base class which describes the basic interface for a solver manager. ...
Teuchos::Array< Teuchos::RCP< Teuchos::Time > > getTimers() const
Return the timers for this object.
std::vector< Value< ScalarType > > Evals
The computed eigenvalues.
ResType
Enumerated type used to specify which residual norm used by residual norm status tests.
A special StatusTest for printing other status tests.
This class defines the interface required by an eigensolver and status test class to compute solution...
RTRSolMgr(const Teuchos::RCP< Eigenproblem< ScalarType, MV, OP > > &problem, Teuchos::ParameterList &pl)
Basic constructor for RTRSolMgr.
An implementation of the Anasazi::SortManager that performs a collection of common sorting techniques...
basic_FancyOStream & setShowProcRank(const bool showProcRank)
T & get(const std::string &name, T def_value)
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Virtual base class which defines basic traits for the operator type.
Teuchos::RCP< MV > Evecs
The computed eigenvectors.
An implementation of the Anasazi::GenOrthoManager that performs orthogonalization using iterated clas...
Status test for forming logical combinations of other status tests.
The Anasazi::SolverManager is a templated virtual base class that defines the basic interface that an...
int getNumIters() const
Get the iteration count for the most recent call to solve.
The Anasazi::RTRSolMgr provides a simple solver manager over the RTR eigensolver. For more informatio...
virtual ~RTRSolMgr()
Destructor.
Basic implementation of the Anasazi::SortManager class.
Namespace Anasazi contains the classes, structs, enums and utilities used by the Anasazi package...
A status test for testing the norm of the eigenvectors residuals along with a set of auxiliary eigenv...
const Eigenproblem< ScalarType, MV, OP > & getProblem() const
Return the eigenvalue problem.
int numVecs
The number of computed eigenpairs.
Basic output manager for sending information of select verbosity levels to the appropriate output str...
ReturnType solve()
This method performs possibly repeated calls to the underlying eigensolver's iterate() routine until ...
Teuchos::RCP< MV > Espace
An orthonormal basis for the computed eigenspace.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Anasazi's basic output manager for sending information of select verbosity levels to the appropriate ...
Abstract base class which defines the interface required by an eigensolver and status test class to c...
static void summarize(Ptr< const Comm< int > > comm, std::ostream &out=std::cout, const bool alwaysWriteLocal=false, const bool writeGlobalStats=true, const bool writeZeroTimers=true, const ECounterSetOp setOp=Intersection, const std::string &filter="", const bool ignoreZeroTimers=false)
ReturnType
Enumerated type used to pass back information from a solver manager.
A status test for testing the norm of the eigenvectors residuals.
Traits class which defines basic operations on multivectors.
std::vector< int > index
An index into Evecs to allow compressed storage of eigenvectors for real, non-Hermitian problems...
basic_FancyOStream & setShowAllFrontMatter(const bool showAllFrontMatter)
Anasazi header file which uses auto-configuration information to include necessary C++ headers...
Struct for storing an eigenproblem solution.
A status test for testing the number of iterations.
Status test for testing the number of iterations.
void push_back(const value_type &x)
Special StatusTest for printing status tests.
A status test for testing the norm of the eigenvectors residuals along with a set of auxiliary eigenv...
Status test for forming logical combinations of other status tests.
Types and exceptions used within Anasazi solvers and interfaces.
bool isParameter(const std::string &name) const
A status test for testing the norm of the eigenvectors residuals.
Basic implementation of the Anasazi::OrthoManager class.
Class which provides internal utilities for the Anasazi solvers.