Teko_MLLinearOp.hpp

#ifndef __Teko_MLLinearOp_hpp__
#define __Teko_MLLinearOp_hpp__
 
#include "Teko_Utilities.hpp"
#include "Teko_BlockImplicitLinearOp.hpp"
 
// forward declarations!
namespace ML_Epetra {
class MultiLevelPreconditioner;
}
 
namespace Teko {
 
namespace Epetra {
class MappingStrategy;
class EpetraOperatorWrapper;
}
 
class MLLinearOp : public BlockImplicitLinearOp {
public:
   MLLinearOp(const Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> & mlPrecOp);

   virtual VectorSpace range() const { return productRange_; }

   virtual VectorSpace domain() const { return productDomain_; }

   virtual void implicitApply(const BlockedMultiVector & x, BlockedMultiVector & y,
              const double alpha = 1.0, const double beta = 0.0) const;
 
   virtual void describe(Teuchos::FancyOStream & out_arg,
                         const Teuchos::EVerbosityLevel verbLevel) const;
 
   Teuchos::RCP<const ML_Epetra::MultiLevelPreconditioner> getMLPreconditioner() const;
   Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> getMLPreconditioner();
 
protected:
 
   void extractConversionInformation(ML_Epetra::MultiLevelPreconditioner & mlPrec);
 
   Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> > productRange_; 
   Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> > productDomain_; 
 
   Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> mlPrecOp_;
   Teuchos::RCP<Epetra::EpetraOperatorWrapper> Amat_;
   Teuchos::RCP<const Epetra::MappingStrategy> mappingStrategy_; 
 
   mutable Teuchos::RCP<Epetra_MultiVector> eX_, eY_; // storage to avoid repeated reallocation
 
private:
   // hide me!
   MLLinearOp();
   MLLinearOp(const MLLinearOp &);
};
 
Teuchos::RCP<const ML_Epetra::MultiLevelPreconditioner> getMLPreconditioner(const Teko::LinearOp & lo);
 
}
 
#endif