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_BORDEREDSOLVER_NESTED_H
00043 #define LOCA_BORDEREDSOLVER_NESTED_H
00044
00045 #include "LOCA_BorderedSolver_AbstractStrategy.H"
00046
00047
00048 namespace LOCA {
00049 class GlobalData;
00050 namespace Parameter {
00051 class SublistParser;
00052 }
00053 namespace BorderedSystem {
00054 class AbstractGroup;
00055 }
00056 }
00057
00058 namespace LOCA {
00059
00060 namespace BorderedSolver {
00061
00063
00093 class Nested :
00094 public LOCA::BorderedSolver::AbstractStrategy {
00095
00096 public:
00097
00099
00104 Nested(
00105 const Teuchos::RCP<LOCA::GlobalData>& global_data,
00106 const Teuchos::RCP<LOCA::Parameter::SublistParser>& topParams,
00107 const Teuchos::RCP<Teuchos::ParameterList>& solverParams);
00108
00110 virtual ~Nested();
00111
00113
00118 virtual void setMatrixBlocks(
00119 const Teuchos::RCP<const LOCA::BorderedSolver::AbstractOperator>& op,
00120 const Teuchos::RCP<const NOX::Abstract::MultiVector>& blockA,
00121 const Teuchos::RCP<const LOCA::MultiContinuation::ConstraintInterface>& blockB,
00122 const Teuchos::RCP<const NOX::Abstract::MultiVector::DenseMatrix>& blockC);
00123
00125
00129 virtual NOX::Abstract::Group::ReturnType
00130 initForSolve();
00131
00133
00137 virtual NOX::Abstract::Group::ReturnType
00138 initForTransposeSolve();
00139
00164 virtual NOX::Abstract::Group::ReturnType
00165 apply(const NOX::Abstract::MultiVector& X,
00166 const NOX::Abstract::MultiVector::DenseMatrix& Y,
00167 NOX::Abstract::MultiVector& U,
00168 NOX::Abstract::MultiVector::DenseMatrix& V) const;
00169
00194 virtual NOX::Abstract::Group::ReturnType
00195 applyTranspose(const NOX::Abstract::MultiVector& X,
00196 const NOX::Abstract::MultiVector::DenseMatrix& Y,
00197 NOX::Abstract::MultiVector& U,
00198 NOX::Abstract::MultiVector::DenseMatrix& V) const;
00199
00208 virtual NOX::Abstract::Group::ReturnType
00209 applyInverse(Teuchos::ParameterList& params,
00210 const NOX::Abstract::MultiVector* F,
00211 const NOX::Abstract::MultiVector::DenseMatrix* G,
00212 NOX::Abstract::MultiVector& X,
00213 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00214
00224 virtual NOX::Abstract::Group::ReturnType
00225 applyInverseTranspose(Teuchos::ParameterList& params,
00226 const NOX::Abstract::MultiVector* F,
00227 const NOX::Abstract::MultiVector::DenseMatrix* G,
00228 NOX::Abstract::MultiVector& X,
00229 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00230
00231 private:
00232
00234 Nested(const Nested&);
00235
00237 Nested& operator = (const Nested&);
00238
00239 protected:
00240
00242 Teuchos::RCP<LOCA::GlobalData> globalData;
00243
00245 Teuchos::RCP<Teuchos::ParameterList> solverParams;
00246
00248 Teuchos::RCP<LOCA::BorderedSolver::AbstractStrategy> solver;
00249
00251 Teuchos::RCP<const LOCA::BorderedSystem::AbstractGroup> grp;
00252
00254 Teuchos::RCP<const NOX::Abstract::Group> unbordered_grp;
00255
00257 int myWidth;
00258
00260 int underlyingWidth;
00261
00263 int numConstraints;
00264
00265 };
00266 }
00267 }
00268
00269 #endif