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