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_SIMPLE_OPERATOR_SOURCE_H
00030 #define MEROS_SIMPLE_OPERATOR_SOURCE_H
00031
00032 #include "Thyra_LinearOpSourceBase.hpp"
00033 #include "Teuchos_ConstNonconstObjectContainer.hpp"
00034 #include "Thyra_VectorImpl.hpp"
00035 #include "Thyra_VectorSpaceImpl.hpp"
00036 #include "Thyra_LinearOperatorDecl.hpp"
00037 #include "Epetra_RowMatrix.h"
00038
00039
00040 namespace Meros
00041 {
00042 using namespace Thyra;
00043
00048
00049 class SIMPLEOperatorSource : virtual public LinearOpSourceBase<double>
00050 {
00051 public:
00052
00055
00058 SIMPLEOperatorSource();
00059
00060
00063 SIMPLEOperatorSource(ConstLinearOperator<double> op);
00064
00065
00068 SIMPLEOperatorSource(Epetra_RowMatrix* S00,
00069 Epetra_RowMatrix* S01,
00070 Epetra_RowMatrix* S10,
00071 Epetra_RowMatrix* S11);
00072
00073
00074
00077 void initialize(ConstLinearOperator<double> op);
00078
00085 void uninitialize();
00086
00087
00089 bool isOpConst() const;
00090
00092 RCP<const LinearOpBase<double> > getOp() const;
00093
00095 RCP<LinearOpBase<double> > getNonconstOp() ;
00096
00098 ConstLinearOperator<double> getDinvOp() const;
00099
00101
00102
00103
00104
00105 private:
00106
00107 Teuchos::ConstNonconstObjectContainer<LinearOpBase<double> > op_;
00108 Teuchos::ConstNonconstObjectContainer<LinearOpBase<double> > dinv_;
00109
00110
00111 };
00112
00113
00114 }
00115
00116 #endif // MEROS_SIMPLE_OPERATOR_SOURCE_H