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_AZTEC_SOLVE_STRATEGY_H
00030 #define MEROS_AZTEC_SOLVE_STRATEGY_H
00031
00032 #include "Meros_LinearSolveStrategy.hpp"
00033 #include "Meros_PreconditioningStrategy.hpp"
00034 #include "Meros_Preconditioner.hpp"
00035 #include "AztecOO.h"
00036 #include "Thyra_AztecOOLinearOpWithSolveFactory.hpp"
00037 #include "Thyra_AztecOOLinearOpWithSolve.hpp"
00038
00039
00040 namespace Meros
00041 {
00042 using namespace Teuchos;
00043 using namespace Thyra;
00044
00048 class AztecSolveStrategy
00049 : public AztecOOLinearOpWithSolveFactory,
00050 public Handleable<LinearOpWithSolveFactoryBase<double> >
00051 {
00052 public:
00054 AztecSolveStrategy(const ParameterList& params)
00055 : AztecOOLinearOpWithSolveFactory()
00056 {
00057 RCP<ParameterList> p = rcp(new ParameterList(params));
00058 this->setParameterList(p);
00059 }
00060
00062 AztecSolveStrategy(const ParameterList& params,
00063 const PreconditioningStrategy<double>& prec)
00064 : AztecOOLinearOpWithSolveFactory()
00065 {
00066 RCP<ParameterList> p = rcp(new ParameterList(params));
00067 this->setParameterList(p);
00068 this->setPreconditionerFactory(prec.ptr(), "dummy name");
00069 }
00070
00071 TEUCHOS_GET_RCP(LinearOpWithSolveFactoryBase<double>);
00072
00074 static ParameterList getParameters()
00075 {
00076 AztecOOLinearOpWithSolveFactory dummy;
00077 return *dummy.getValidParameters();
00078 }
00079
00080
00081
00082 };
00083
00084
00085 }
00086
00087 #endif // MEROS_PCD_OPERATOR_SOURCE_H