00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifndef MEROS_PRECONDITIONER_H
00030 #define MEROS_PRECONDITIONER_H
00031
00032 #include "Thyra_VectorImpl.hpp"
00033 #include "Thyra_VectorSpaceImpl.hpp"
00034 #include "Thyra_LinearOperatorImpl.hpp"
00035 #include "Thyra_PreconditionerBase.hpp"
00036 #include "Teuchos_RCP.hpp"
00037
00038
00039 namespace Meros
00040 {
00041
00042 using namespace Teuchos;
00043 using namespace Thyra;
00044
00048 template <class Scalar>
00049 class Preconditioner
00050 {
00051 public:
00053 Preconditioner() : p_() {;}
00054
00056 Preconditioner(const RCP<PreconditionerBase<Scalar> >& p)
00057 : p_(p) {;}
00058
00060 Preconditioner(Handleable<PreconditionerBase<Scalar> >* p)
00061 : p_(p->getRcp()) {;}
00062
00064 LinearOperator<Scalar> getRight() const
00065 {
00066 return p_->getNonconstRightPrecOp();
00067 }
00068
00070 bool hasRight() const
00071 {
00072 return p_->getNonconstRightPrecOp().get() != 0;
00073 }
00074
00076 LinearOperator<Scalar> getLeft() const {return p_->getNonconstLeftPrecOp();}
00077
00079 bool hasLeft() const
00080 {
00081 return p_->getNonconstLeftPrecOp().get() != 0;
00082 }
00083
00085 LinearOperator<Scalar> getGenericPrec() const {return p_->getNonconstGenericPrecOp();}
00086
00087
00089 bool hasGeneric() const
00090 {
00091 return p_->getNonconstUnspecifiedPrecOp().get() != 0;
00092 }
00093
00095 const RCP<PreconditionerBase<Scalar> >& ptr() const {return p_;}
00097 RCP<PreconditionerBase<Scalar> > ptr() {return p_;}
00098
00099 private:
00100 RCP<PreconditionerBase<Scalar> > p_;
00101 };
00102
00103
00104 }
00105
00106 #endif // MEROS_PCD_OPERATOR_SOURCE_H