42 #ifndef BELOS_DETAILS_LINEARSOLVER_HPP 43 #define BELOS_DETAILS_LINEARSOLVER_HPP 49 #include "Trilinos_Details_LinearSolver.hpp" 75 template<
class MV,
class OP,
class ScalarType,
class NormType>
77 public Trilinos::Details::LinearSolver<MV, OP, NormType>
109 }
else if (A !=
problem_->getOperator ()) {
117 return Teuchos::null;
125 TEUCHOS_TEST_FOR_EXCEPTION
127 std::runtime_error,
"Belos::Details::LinearSolver::solve: " 128 "The matrix A in the linear system to solve has not yet been set. " 129 "Please call setMatrix() with a nonnull input before calling solve().");
130 Teuchos::RCP<MV> X_ptr = Teuchos::rcpFromRef (X);
131 Teuchos::RCP<const MV> B_ptr = Teuchos::rcpFromRef (B);
152 if (!
solver_.is_null () && ! params.is_null ()) {
153 solver_->setParameters (params);
160 TEUCHOS_TEST_FOR_EXCEPTION
162 std::runtime_error,
"Belos::Details::LinearSolver::symbolic: " 163 "The matrix A in the linear system to solve has not yet been set. " 164 "Please call setMatrix() with a nonnull input before calling this method.");
174 TEUCHOS_TEST_FOR_EXCEPTION
176 std::runtime_error,
"Belos::Details::LinearSolver::numeric: " 177 "The matrix A in the linear system to solve has not yet been set. " 178 "Please call setMatrix() with a nonnull input before calling this method.");
void solve(MV &X, const MV &B)
Solve the linear system AX=B for X.
std::string solverName_
The name of the Belos solver to create.
Teuchos::RCP< Teuchos::ParameterList > params_
The Belos solver's list of parameters.
Belos' implementation of Trilinos::Details::LinearSolver.
virtual ~LinearSolver()
Destructor (virtual for memory safety).
void setMatrix(const Teuchos::RCP< const OP > &A)
Set the solver's matrix.
void symbolic()
Precompute for matrix structure changes.
void numeric()
Precompute for matrix values' changes.
Belos::LinearProblem< ScalarType, MV, OP > problem_type
Belos::LinearProblem specialization used by this class.
A linear system to solve, and its associated information.
LinearSolver(const std::string &solverName)
Constructor.
The Belos::SolverManager is a templated virtual base class that defines the basic interface that any ...
void setParameters(const Teuchos::RCP< Teuchos::ParameterList > ¶ms)
Set the solver's parameters.
Teuchos::RCP< solver_type > solver_
The Belos solver (SolverManager instance).
Teuchos::RCP< const OP > getMatrix() const
Get the solver's matrix.
Teuchos::RCP< problem_type > problem_
The LinearProblem instance to give to the Belos solver.
typename ::Belos::Impl::SolverFactorySelector< SC, MV, OP >::type SolverFactory
Belos::SolverManager< ScalarType, MV, OP > solver_type
Belos' own solver type.