1 #ifndef __Teko_MLLinearOp_hpp__
2 #define __Teko_MLLinearOp_hpp__
5 #include "Teko_BlockImplicitLinearOp.hpp"
9 class MultiLevelPreconditioner;
15 class MappingStrategy;
16 class EpetraOperatorWrapper;
19 class MLLinearOp :
public BlockImplicitLinearOp {
21 MLLinearOp(
const Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> &mlPrecOp);
24 virtual VectorSpace range()
const {
return productRange_; }
27 virtual VectorSpace domain()
const {
return productDomain_; }
41 virtual void implicitApply(
const BlockedMultiVector &x, BlockedMultiVector &y,
42 const double alpha = 1.0,
const double beta = 0.0)
const;
44 virtual void describe(Teuchos::FancyOStream &out_arg,
45 const Teuchos::EVerbosityLevel verbLevel)
const;
47 Teuchos::RCP<const ML_Epetra::MultiLevelPreconditioner> getMLPreconditioner()
const;
48 Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> getMLPreconditioner();
51 void extractConversionInformation(ML_Epetra::MultiLevelPreconditioner &mlPrec);
53 Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> >
55 Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> >
58 Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> mlPrecOp_;
59 Teuchos::RCP<Epetra::EpetraOperatorWrapper> Amat_;
60 Teuchos::RCP<const Epetra::MappingStrategy>
63 mutable Teuchos::RCP<Epetra_MultiVector> eX_, eY_;
68 MLLinearOp(
const MLLinearOp &);
71 Teuchos::RCP<const ML_Epetra::MultiLevelPreconditioner> getMLPreconditioner(
72 const Teko::LinearOp &lo);
virtual void implicitApply(const Thyra::EOpTransp M_trans, const BlockedMultiVector &x, BlockedMultiVector &y, const double alpha=1.0, const double beta=0.0) const
Perform a matrix vector multiply with this implicitly defined blocked operator.