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_ANASAZIOPERATOR_JACOBIANINVERSE_H
00043 #define LOCA_ANASAZIOPERATOR_JACOBIANINVERSE_H
00044
00045 #include "Teuchos_RCP.hpp"
00046
00047 #include "LOCA_AnasaziOperator_AbstractStrategy.H"
00048
00049
00050 namespace LOCA {
00051 class GlobalData;
00052 namespace Parameter {
00053 class SublistParser;
00054 }
00055 }
00056
00057 namespace LOCA {
00058
00059 namespace AnasaziOperator {
00060
00062
00066 class JacobianInverse : public LOCA::AnasaziOperator::AbstractStrategy {
00067
00068 public:
00069
00071 JacobianInverse(
00072 const Teuchos::RCP<LOCA::GlobalData>& global_data,
00073 const Teuchos::RCP<LOCA::Parameter::SublistParser>& topParams,
00074 const Teuchos::RCP<Teuchos::ParameterList>& eigenParams,
00075 const Teuchos::RCP<Teuchos::ParameterList>& solverParams,
00076 const Teuchos::RCP<NOX::Abstract::Group>& grp);
00077
00079 virtual ~JacobianInverse();
00080
00084 virtual const string& label() const;
00085
00091 virtual void
00092 apply(const NOX::Abstract::MultiVector& input,
00093 NOX::Abstract::MultiVector& output) const;
00094
00101 virtual void
00102 transformEigenvalue(double& ev_r, double& ev_i) const;
00103
00110 virtual NOX::Abstract::Group::ReturnType
00111 rayleighQuotient(const NOX::Abstract::Vector& evec_r,
00112 const NOX::Abstract::Vector& evec_i,
00113 double& rq_r, double& rq_i) const;
00114
00115 protected:
00116
00118 Teuchos::RCP<LOCA::GlobalData> globalData;
00119
00121 string myLabel;
00122
00124 Teuchos::RCP<Teuchos::ParameterList> eigenParams;
00125
00127 Teuchos::RCP<Teuchos::ParameterList> solverParams;
00128
00130 Teuchos::RCP<NOX::Abstract::Group> grp;
00131
00133 mutable Teuchos::RCP<NOX::Abstract::Vector> tmp_r;
00134
00136 mutable Teuchos::RCP<NOX::Abstract::Vector> tmp_i;
00137
00138 };
00139 }
00140 }
00141
00142 #endif