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_BORDERING_H
00043 #define LOCA_BORDEREDSOLVER_BORDERING_H
00044
00045 #include "LOCA_BorderedSolver_AbstractStrategy.H"
00046
00047
00048 namespace LOCA {
00049 class GlobalData;
00050 namespace Parameter {
00051 class SublistParser;
00052 }
00053 }
00054
00055 namespace LOCA {
00056
00057 namespace BorderedSolver {
00058
00060
00095 class Bordering : public LOCA::BorderedSolver::AbstractStrategy {
00096
00097 public:
00098
00100
00106 Bordering(
00107 const Teuchos::RCP<LOCA::GlobalData>& global_data,
00108 const Teuchos::RCP<LOCA::Parameter::SublistParser>& topParams,
00109 const Teuchos::RCP<Teuchos::ParameterList>& solverParams);
00110
00112 virtual ~Bordering();
00113
00115
00120 virtual void setMatrixBlocks(
00121 const Teuchos::RCP<const LOCA::BorderedSolver::AbstractOperator>& op,
00122 const Teuchos::RCP<const NOX::Abstract::MultiVector>& blockA,
00123 const Teuchos::RCP<const LOCA::MultiContinuation::ConstraintInterface>& blockB,
00124 const Teuchos::RCP<const NOX::Abstract::MultiVector::DenseMatrix>& blockC);
00125
00127
00131 virtual NOX::Abstract::Group::ReturnType
00132 initForSolve();
00133
00135
00139 virtual NOX::Abstract::Group::ReturnType
00140 initForTransposeSolve();
00141
00166 virtual NOX::Abstract::Group::ReturnType
00167 apply(const NOX::Abstract::MultiVector& X,
00168 const NOX::Abstract::MultiVector::DenseMatrix& Y,
00169 NOX::Abstract::MultiVector& U,
00170 NOX::Abstract::MultiVector::DenseMatrix& V) const;
00171
00196 virtual NOX::Abstract::Group::ReturnType
00197 applyTranspose(const NOX::Abstract::MultiVector& X,
00198 const NOX::Abstract::MultiVector::DenseMatrix& Y,
00199 NOX::Abstract::MultiVector& U,
00200 NOX::Abstract::MultiVector::DenseMatrix& V) const;
00201
00210 virtual NOX::Abstract::Group::ReturnType
00211 applyInverse(Teuchos::ParameterList& params,
00212 const NOX::Abstract::MultiVector* F,
00213 const NOX::Abstract::MultiVector::DenseMatrix* G,
00214 NOX::Abstract::MultiVector& X,
00215 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00216
00228 virtual NOX::Abstract::Group::ReturnType
00229 applyInverseTranspose(Teuchos::ParameterList& params,
00230 const NOX::Abstract::MultiVector* F,
00231 const NOX::Abstract::MultiVector::DenseMatrix* G,
00232 NOX::Abstract::MultiVector& X,
00233 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00234
00235 protected:
00236
00240 NOX::Abstract::Group::ReturnType
00241 solveFZero(Teuchos::ParameterList& params,
00242 const NOX::Abstract::MultiVector* AA,
00243 const LOCA::MultiContinuation::ConstraintInterface* BB,
00244 const NOX::Abstract::MultiVector::DenseMatrix* CC,
00245 const NOX::Abstract::MultiVector::DenseMatrix* G,
00246 NOX::Abstract::MultiVector& X,
00247 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00248
00252 NOX::Abstract::Group::ReturnType
00253 solveContiguous(Teuchos::ParameterList& params,
00254 const NOX::Abstract::MultiVector* AA,
00255 const LOCA::MultiContinuation::ConstraintInterface* BB,
00256 const NOX::Abstract::MultiVector::DenseMatrix* CC,
00257 vector<int>& indexF,
00258 vector<int>& indexA,
00259 const NOX::Abstract::MultiVector* F,
00260 const NOX::Abstract::MultiVector::DenseMatrix* G,
00261 NOX::Abstract::MultiVector& X,
00262 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00263
00267 NOX::Abstract::Group::ReturnType
00268 solveFZeroTrans(Teuchos::ParameterList& params,
00269 const NOX::Abstract::MultiVector* AA,
00270 const NOX::Abstract::MultiVector* BB,
00271 const NOX::Abstract::MultiVector::DenseMatrix* CC,
00272 const NOX::Abstract::MultiVector::DenseMatrix* G,
00273 NOX::Abstract::MultiVector& X,
00274 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00275
00280 NOX::Abstract::Group::ReturnType
00281 solveContiguousTrans(
00282 Teuchos::ParameterList& params,
00283 const NOX::Abstract::MultiVector* AA,
00284 const NOX::Abstract::MultiVector* BB,
00285 const NOX::Abstract::MultiVector::DenseMatrix* CC,
00286 vector<int>& indexF,
00287 vector<int>& indexA,
00288 const NOX::Abstract::MultiVector* F,
00289 const NOX::Abstract::MultiVector::DenseMatrix* G,
00290 NOX::Abstract::MultiVector& X,
00291 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00292
00293 private:
00294
00296 Bordering(const Bordering&);
00297
00299 Bordering& operator = (const Bordering&);
00300
00301 protected:
00302
00304 Teuchos::RCP<LOCA::GlobalData> globalData;
00305
00307 Teuchos::RCP<Teuchos::ParameterList> solverParams;
00308
00310 Teuchos::RCP<const LOCA::BorderedSolver::AbstractOperator> op;
00311
00313 Teuchos::RCP<const NOX::Abstract::MultiVector> A;
00314
00316 Teuchos::RCP<const LOCA::MultiContinuation::ConstraintInterface> B;
00317
00319 Teuchos::RCP<const NOX::Abstract::MultiVector::DenseMatrix> C;
00320
00322 bool isZeroA;
00323
00325 bool isZeroB;
00326
00328 bool isZeroC;
00329
00331 mutable bool isZeroF;
00332
00334 mutable bool isZeroG;
00335
00336 };
00337 }
00338 }
00339
00340 #endif