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
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042 #ifndef LOCA_EPETRA_TRANSPOSELINEARSYSTEM_LEFTPRECONDITIONING_H
00043 #define LOCA_EPETRA_TRANSPOSELINEARSYSTEM_LEFTPRECONDITIONING_H
00044
00045 #include "LOCA_Epetra_TransposeLinearSystem_AbstractStrategy.H"
00046
00047
00048 namespace Teuchos {
00049 class ParameterList;
00050 }
00051 namespace NOX {
00052 namespace Epetra {
00053 class LinearSystem;
00054 }
00055 }
00056 namespace LOCA {
00057 class GlobalData;
00058 }
00059
00060 namespace LOCA {
00061
00062 namespace Epetra {
00063
00064 namespace TransposeLinearSystem {
00065
00071
00072
00073
00074 class LeftPreconditioning : public virtual AbstractStrategy {
00075
00076 public:
00077
00079 LeftPreconditioning(
00080 const Teuchos::RCP<LOCA::GlobalData>& global_data,
00081 const Teuchos::RCP<Teuchos::ParameterList>& solverParams,
00082 const Teuchos::RCP<NOX::Epetra::LinearSystem>& linsys);
00083
00085 virtual ~LeftPreconditioning();
00086
00098 virtual bool
00099 applyJacobianTransposeInverse(Teuchos::ParameterList ¶ms,
00100 const NOX::Epetra::Vector &input,
00101 NOX::Epetra::Vector &result);
00102
00103
00104
00106
00111 virtual bool
00112 createJacobianTranspose();
00113
00123 virtual bool
00124 createTransposePreconditioner(const NOX::Epetra::Vector& x,
00125 Teuchos::ParameterList& p);
00126
00128 virtual Teuchos::RCP<Epetra_Operator>
00129 getJacobianTransposeOperator();
00130
00132 virtual Teuchos::RCP<Epetra_Operator>
00133 getTransposePreconditioner();
00134
00136 virtual void setJacobianTransposeOperator(
00137 const Teuchos::RCP<Epetra_Operator>& new_jac_trans);
00138
00140 virtual void setTransposePreconditioner(
00141 const Teuchos::RCP<Epetra_Operator>& new_prec_trans);
00142
00143 protected:
00144
00146 Teuchos::RCP<LOCA::GlobalData> globalData;
00147
00149 Teuchos::RCP<NOX::Epetra::LinearSystem> linsys;
00150
00152 Teuchos::RCP<Epetra_Operator> jac;
00153
00155 Teuchos::RCP<Epetra_Operator> prec;
00156
00157 };
00158
00159 }
00160
00161 }
00162
00163 }
00164
00165
00166 #endif