42#ifndef BELOS_DETAILS_LINEARSOLVER_HPP
43#define BELOS_DETAILS_LINEARSOLVER_HPP
49#include "Trilinos_Details_LinearSolver.hpp"
75template<
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.");
Belos' implementation of Trilinos::Details::LinearSolver.
void setParameters(const Teuchos::RCP< Teuchos::ParameterList > ¶ms)
Set the solver's parameters.
void symbolic()
Precompute for matrix structure changes.
void setMatrix(const Teuchos::RCP< const OP > &A)
Set the solver's matrix.
Teuchos::RCP< problem_type > problem_
The LinearProblem instance to give to the Belos solver.
Teuchos::RCP< const OP > getMatrix() const
Get the solver's matrix.
Belos::SolverManager< ScalarType, MV, OP > solver_type
Belos' own solver type.
LinearSolver(const std::string &solverName)
Constructor.
virtual ~LinearSolver()
Destructor (virtual for memory safety).
std::string solverName_
The name of the Belos solver to create.
Belos::LinearProblem< ScalarType, MV, OP > problem_type
Belos::LinearProblem specialization used by this class.
void numeric()
Precompute for matrix values' changes.
Teuchos::RCP< solver_type > solver_
The Belos solver (SolverManager instance).
void solve(MV &X, const MV &B)
Solve the linear system AX=B for X.
Teuchos::RCP< Teuchos::ParameterList > params_
The Belos solver's list of parameters.
A linear system to solve, and its associated information.
The Belos::SolverManager is a templated virtual base class that defines the basic interface that any ...
typename ::Belos::Impl::SolverFactorySelector< SC, MV, OP >::type SolverFactory