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_UPPERTRIANGULARBLOCKELIMINATION_H
00043 #define LOCA_BORDEREDSOLVER_UPPERTRIANGULARBLOCKELIMINATION_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
00090 class UpperTriangularBlockElimination {
00091
00092 public:
00093
00095
00098 UpperTriangularBlockElimination(
00099 const Teuchos::RCP<LOCA::GlobalData>& global_data);
00100
00102 virtual ~UpperTriangularBlockElimination();
00103
00110 NOX::Abstract::Group::ReturnType
00111 solve(Teuchos::ParameterList& params,
00112 const LOCA::BorderedSolver::AbstractOperator& op,
00113 const NOX::Abstract::MultiVector* A,
00114 const NOX::Abstract::MultiVector::DenseMatrix& C,
00115 const NOX::Abstract::MultiVector* F,
00116 const NOX::Abstract::MultiVector::DenseMatrix* G,
00117 NOX::Abstract::MultiVector& X,
00118 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00119
00127 NOX::Abstract::Group::ReturnType
00128 solveTranspose(Teuchos::ParameterList& params,
00129 const LOCA::BorderedSolver::AbstractOperator& op,
00130 const NOX::Abstract::MultiVector* A,
00131 const NOX::Abstract::MultiVector::DenseMatrix& C,
00132 const NOX::Abstract::MultiVector* F,
00133 const NOX::Abstract::MultiVector::DenseMatrix* G,
00134 NOX::Abstract::MultiVector& X,
00135 NOX::Abstract::MultiVector::DenseMatrix& Y) const;
00136
00137 private:
00138
00140 UpperTriangularBlockElimination(const UpperTriangularBlockElimination&);
00141
00143 UpperTriangularBlockElimination&
00144 operator = (const UpperTriangularBlockElimination&);
00145
00146 protected:
00147
00149 Teuchos::RCP<LOCA::GlobalData> globalData;
00150
00151 };
00152 }
00153 }
00154
00155 #endif