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 #ifndef LOCA_EPETRA_ANASAZIOPERATOR_FLOQUET_H
00034 #define LOCA_EPETRA_ANASAZIOPERATOR_FLOQUET_H
00035
00036 #include "Teuchos_RCP.hpp"
00037
00038 #include "LOCA_AnasaziOperator_AbstractStrategy.H"
00039 #include "LOCA_TimeDependent_AbstractGroup.H"
00040 #include "LOCA_Epetra_Interface_xyzt.H"
00041 #include "NOX_Epetra_Group.H"
00042
00043
00044 namespace LOCA {
00045 class GlobalData;
00046 namespace Parameter {
00047 class SublistParser;
00048 }
00049 }
00050
00051 namespace LOCA {
00052
00053 namespace Epetra {
00054 namespace AnasaziOperator {
00055
00085 class Floquet : public LOCA::AnasaziOperator::AbstractStrategy {
00086
00087 public:
00088
00090
00094 Floquet(
00095 const Teuchos::RCP<LOCA::GlobalData>& global_data,
00096 const Teuchos::RCP<LOCA::Parameter::SublistParser>& topParams,
00097 const Teuchos::RCP<Teuchos::ParameterList>& eigenParams,
00098 const Teuchos::RCP<Teuchos::ParameterList>& solverParams,
00099 const Teuchos::RCP<NOX::Abstract::Group>& grp);
00100
00102 virtual ~Floquet();
00103
00107 virtual const string& label() const;
00108
00118 virtual void
00119 apply(const NOX::Abstract::MultiVector& input,
00120 NOX::Abstract::MultiVector& output) const;
00121
00128 virtual void
00129 transformEigenvalue(double& ev_r, double& ev_i) const;
00130
00137 virtual NOX::Abstract::Group::ReturnType
00138 rayleighQuotient(const NOX::Abstract::Vector& evec_r,
00139 const NOX::Abstract::Vector& evec_i,
00140 double& rq_r, double& rq_i) const;
00141
00142 protected:
00143
00145 Teuchos::RCP<LOCA::GlobalData> globalData;
00146
00148 string myLabel;
00149
00151 Teuchos::RCP<Teuchos::ParameterList> eigenParams;
00152
00154 Teuchos::RCP<Teuchos::ParameterList> solverParams;
00155
00157 Teuchos::RCP<NOX::Abstract::Group> grp;
00158
00160 Teuchos::RCP<LOCA::Epetra::Interface::xyzt> xyztInterface;
00161
00162 };
00163 }
00164 }
00165 }
00166
00167 #endif