1#ifndef __Teko_MLLinearOp_hpp__
2#define __Teko_MLLinearOp_hpp__
5#include "Teko_BlockImplicitLinearOp.hpp"
9class MultiLevelPreconditioner;
16class EpetraOperatorWrapper;
19class 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();
52 void extractConversionInformation(ML_Epetra::MultiLevelPreconditioner & mlPrec);
54 Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> > productRange_;
55 Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> > productDomain_;
57 Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> mlPrecOp_;
58 Teuchos::RCP<Epetra::EpetraOperatorWrapper> Amat_;
59 Teuchos::RCP<const Epetra::MappingStrategy> mappingStrategy_;
61 mutable Teuchos::RCP<Epetra_MultiVector> eX_, eY_;
66 MLLinearOp(
const MLLinearOp &);
69Teuchos::RCP<const ML_Epetra::MultiLevelPreconditioner> getMLPreconditioner(
const Teko::LinearOp & lo);