42 #ifndef BELOS_OPERATOR_HPP 43 #define BELOS_OPERATOR_HPP 52 #ifdef HAVE_BELOS_EXPERIMENTAL 54 #endif // HAVE_BELOS_EXPERIMENTAL 80 template <
class ScalarType>
162 template<
class ScalarType>
173 Op.
Apply (x, y, trans);
184 #ifdef HAVE_BELOS_EXPERIMENTAL 199 template<
class Scalar>
200 class OperatorInnerSolver :
public Operator<Scalar> {
202 typedef Scalar scalar_type;
203 typedef MultiVec<Scalar> multivector_type;
204 typedef Operator<Scalar> operator_type;
205 typedef InnerSolver<scalar_type, multivector_type, operator_type> inner_solver_type;
214 virtual ~OperatorInnerSolver() {}
251 Apply (
const multivector_type& X,
255 using Teuchos::rcpFromRef;
258 "Belos::OperatorInnerSolver only supports applying the" 259 " operator itself, not its transpose or conjugate " 261 solver_->solve (rcpFromRef (Y), rcpFromRef (X));
268 virtual bool HasApplyTranspose()
const {
274 OperatorInnerSolver ();
291 template <
class Scalar>
292 class InnerSolverTraits<Scalar, MultiVec<Scalar>, Operator<Scalar> > {
297 typedef InnerSolver<scalar_type, multivector_type, operator_type>
inner_solver_type;
308 using Teuchos::rcp_implicit_cast;
325 using Teuchos::rcp_dynamic_cast;
326 RCP<wrapper_type> wrapper = rcp_dynamic_cast<
wrapper_type> (op,
true);
327 return wrapper->getInnerSolver();
330 #endif // HAVE_BELOS_EXPERIMENTAL virtual void Apply(const MultiVec< ScalarType > &x, MultiVec< ScalarType > &y, ETrans trans=NOTRANS) const =0
Apply the operator to x, putting the result in y.
static Teuchos::RCP< inner_solver_type > getInnerSolver(const Teuchos::RCP< operator_type > &op)
Return the given wrapper's inner solver object.
static bool HasApplyTranspose(const Operator< ScalarType > &Op)
Specialization of HasApplyTranspose() for Operator objects.
static Teuchos::RCP< OP > makeInnerSolverOperator(const Teuchos::RCP< InnerSolver< Scalar, MV, OP > > &solver)
Wrap the given inner solver in a wrapper_type.
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
virtual ~Operator()
Virtual destructor, for memory safety of derived classes.
Class which defines basic traits for the operator type.
ETrans
Whether to apply the (conjugate) transpose of an operator.
static void Apply(const Operator< ScalarType > &Op, const MultiVec< ScalarType > &x, MultiVec< ScalarType > &y, ETrans trans=NOTRANS)
Specialization of Apply() for Operator and MultiVec objects.
Alternative run-time polymorphic interface for operators.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
UndefinedWrapperType< Scalar, MV, OP > wrapper_type
InnerSolver< scalar_type, multivector_type, operator_type > inner_solver_type
Operator()
Default constructor (does nothing).
virtual bool HasApplyTranspose() const
Whether this operator implements applying the transpose.
Interface for multivectors used by Belos' linear solvers.
Class which defines basic traits for the operator type.
Belos header file which uses auto-configuration information to include necessary C++ headers...
Interface for multivectors used by Belos' linear solvers.