43#ifndef BELOS_STATUS_TEST_USER_OUTPUT_HPP
44#define BELOS_STATUS_TEST_USER_OUTPUT_HPP
69template <
class ScalarType,
class MV,
class OP>
144 TEUCHOS_TEST_FOR_EXCEPTION(
iterTest_ == Teuchos::null,
StatusTestError,
"StatusTestUserOutput::checkStatus(): iteration test pointer is null.");
145 TEUCHOS_TEST_FOR_EXCEPTION(
resTestVec_.size() == 0,
StatusTestError,
"StatusTestUserOutput::checkStatus(): residual test pointer is null.");
200 Teuchos::RCP<StatusTestCombo_t> comboTest = Teuchos::rcp_dynamic_cast<StatusTestCombo_t>(test);
201 TEUCHOS_TEST_FOR_EXCEPTION(comboTest == Teuchos::null,
StatusTestError,
"StatusTestUserOutput::setChild: The parameter \"test\" must be a Belos::StatusTestCombo.");
202 std::vector<Teuchos::RCP<StatusTest<ScalarType,MV,OP> > > tmpVec = comboTest->getStatusTests();
205 int numTests = tmpVec.size();
208 for (
int i=0; i<numTests; ++i) {
211 Teuchos::RCP<StatusTestMaxIters_t> tmpItrTest = Teuchos::rcp_dynamic_cast<StatusTestMaxIters_t>(tmpVec[i]);
212 if (tmpItrTest != Teuchos::null) {
220 Teuchos::RCP<StatusTestResNorm_t> tmpResTest = Teuchos::rcp_dynamic_cast<StatusTestResNorm_t>(tmpVec[i]);
222 if (tmpResTest != Teuchos::null) {
233 Teuchos::RCP<StatusTestCombo_t> tmpComboTest = Teuchos::rcp_dynamic_cast<StatusTestCombo_t>(tmpVec[i]);
234 TEUCHOS_TEST_FOR_EXCEPTION(tmpComboTest == Teuchos::null,
StatusTestError,
"StatusTestUserOutput(): test must be Belos::StatusTest[MaxIters|ResNorm|Combo].");
235 tmpVec = tmpComboTest->getStatusTests();
240 typename std::map<std::string,Teuchos::RCP<StatusTest<ScalarType,MV,OP> > >::iterator it;
241 for (
size_t j=0; j<tmpVec.size(); ++j) {
242 tmpResTest = Teuchos::rcp_dynamic_cast<StatusTestResNorm_t>(tmpVec[j]);
244 if(tmpResTest == Teuchos::null)
continue;
248 if(tmpVec[j] == it->second) { bFound =
true;
break; }
270 Teuchos::RCP<StatusTest<ScalarType,MV,OP> >
getChild()
const {
310 void print(std::ostream& os,
int indent = 0)
const {
311 std::string ind(indent,
' ');
312 std::string starLine(55,
'*');
313 std::string starFront(5,
'*');
315 std::ios_base::fmtflags osFlags(os.flags());
317 os.setf(std::ios::scientific, std::ios::floatfield);
322 os << std::endl << ind << starLine << std::endl;
323 os << ind << starFront <<
" Belos Iterative Solver: " <<
solverDesc_ << std::endl;
325 os << ind << starFront <<
" Preconditioner: " <<
precondDesc_ << std::endl;
326 os << ind << starFront <<
" Maximum Iterations: " <<
iterTest_->getMaxIters() << std::endl;
327 os << ind << starFront <<
" Block Size: " <<
blockSize_ << std::endl;
328 os << ind << starFront <<
" Status tests: " << std::endl;
329 test_->print(os,indent + 3);
330 os << ind << starLine << std::endl;
331 os.setf(std::ios_base::right, std::ios_base::adjustfield);
337 os << ind << indheader;
349 os.setf(std::ios_base::right, std::ios_base::adjustfield);
361 Teuchos::RCP<StatusTestResNorm_t> tempResTest = Teuchos::rcp_dynamic_cast<StatusTestResNorm_t>(
resTestVec_[j]);
362 if(tempResTest != Teuchos::null)
363 os << std::setw(15) << (*tempResTest->getTestValue())[
currIdx_[i]];
366 os << std::setw(15) <<
"Passed";
368 os << std::setw(15) <<
"Failed";
369 else os << std::setw(15) <<
"Undefined";
372 os << std::setw(15) <<
"---";
388 Teuchos::RCP<StatusTest<ScalarType,MV,OP> >
test_;
391 Teuchos::RCP<std::map<std::string,Teuchos::RCP<StatusTest<ScalarType,MV,OP> > > >
taggedTests_;
394 Teuchos::RCP<StatusTestMaxIters<ScalarType,MV,OP> >
iterTest_;
397 std::vector<Teuchos::RCP<StatusTest<ScalarType,MV,OP> > >
resTestVec_;
Belos header file which uses auto-configuration information to include necessary C++ headers.
Pure virtual base class which describes the basic interface to the linear solver iteration.
Belos::StatusTest for logically combining several status tests.
Belos::StatusTest class for specifying a maximum number of iterations.
Virtual base class for StatusTest that printing status tests.
Belos::StatusTest abstract class for specifying a residual norm stopping criteria.
Pure virtual base class for defining the status testing capabilities of Belos.
Collection of types and exceptions used within the Belos solvers.
virtual int getBlockSize() const =0
Get the blocksize to be used by the iterative solver in solving this linear problem.
virtual const LinearProblem< ScalarType, MV, OP > & getProblem() const =0
Get a constant reference to the linear problem.
A linear system to solve, and its associated information.
int getLSNumber() const
The number of linear systems that have been set.
const std::vector< int > & getLSIndex() const
(Zero-based) indices of the linear system(s) currently being solved.
Teuchos::RCP< const MV > getRHS() const
A pointer to the right-hand side B.
Traits class which defines basic operations on multivectors.
static int GetNumberVecs(const MV &mv)
Obtain the number of vectors in mv.
Belos's basic output manager for sending information of select verbosity levels to the appropriate ou...
A class for extending the status testing capabilities of Belos via logical combinations.
Exception thrown to signal error in a status test during Belos::StatusTest::checkStatus().
A Belos::StatusTest class for specifying a maximum number of iterations.
A virtual base class for StatusTest that print other status tests.
An abstract class of StatusTest for stopping criteria using residual norms.
A special StatusTest for printing other status tests in a simple format.
void setSolverDesc(const std::string &solverDesc)
Set a short solver description for output clarity.
MultiVecTraits< ScalarType, MV > MVT
void setOutputFrequency(int mod)
Set how often the child test is printed.
Teuchos::RCP< StatusTest< ScalarType, MV, OP > > getChild() const
Get child test.
StatusType checkStatus(Iteration< ScalarType, MV, OP > *solver)
StatusTestUserOutput(const Teuchos::RCP< OutputManager< ScalarType > > &printer, Teuchos::RCP< StatusTest< ScalarType, MV, OP > > test, Teuchos::RCP< std::map< std::string, Teuchos::RCP< StatusTest< ScalarType, MV, OP > > > > taggedTests, int mod=1, int printStates=Passed)
Constructor.
void setPrecondDesc(const std::string &precondDesc)
Set a short preconditioner description for output clarity.
Teuchos::RCP< std::map< std::string, Teuchos::RCP< StatusTest< ScalarType, MV, OP > > > > taggedTests_
Teuchos::RCP< StatusTestMaxIters< ScalarType, MV, OP > > iterTest_
void resetNumCalls()
Informs the outputting status test that it should reset the number of calls to zero.
void print(std::ostream &os, int indent=0) const
Output formatted description of stopping test to output stream.
StatusType getStatus() const
Return the result of the most recent checkStatus call, or undefined if it has not been run.
std::vector< std::string > resTestNamesVec_
Name tags for status tests.
Belos::StatusTestResNorm< ScalarType, MV, OP > StatusTestResNorm_t
void reset()
Informs the status test that it should reset its internal configuration to the uninitialized state.
Teuchos::RCP< StatusTest< ScalarType, MV, OP > > test_
Belos::StatusTestMaxIters< ScalarType, MV, OP > StatusTestMaxIters_t
Teuchos::RCP< OutputManager< ScalarType > > printer_
void setChild(Teuchos::RCP< StatusTest< ScalarType, MV, OP > > test)
Set child test, which must be a combination of a Belos::StatusTestMaxIters AND a combination of Belos...
std::vector< int > currIdx_
void setOutputManager(const Teuchos::RCP< OutputManager< ScalarType > > &printer)
Set the output manager.
virtual ~StatusTestUserOutput()
Destructor.
std::vector< Teuchos::RCP< StatusTest< ScalarType, MV, OP > > > resTestVec_
Vector of residual status tests.
Belos::StatusTestCombo< ScalarType, MV, OP > StatusTestCombo_t
A pure virtual class for defining the status tests for the Belos iterative solvers.
StatusType
Whether the StatusTest wants iteration to stop.