47 #include "Teko_GaussSeidelPreconditionerFactory.hpp"
49 #include "Teko_BlockUpperTriInverseOp.hpp"
50 #include "Teko_BlockLowerTriInverseOp.hpp"
58 const LinearOp& invD0,
59 const LinearOp& invD1)
63 TriSolveType solveType,
const RCP<const BlockInvDiagonalStrategy>& strategy)
64 : invOpsStrategy_(strategy), solveType_(solveType) {}
67 : solveType_(GS_UseLowerTriangle) {}
71 int rows = blockRowCount(blo);
72 int cols = blockColCount(blo);
74 TEUCHOS_ASSERT(rows == cols);
77 std::vector<LinearOp> invDiag;
79 TEUCHOS_ASSERT(rows == (
int)invDiag.size());
81 if (solveType_ == GS_UseUpperTriangle) {
83 BlockedLinearOp U = getUpperTriBlocks(blo);
85 return createBlockUpperTriInverseOp(U, invDiag,
"Gauss Seidel");
86 }
else if (solveType_ == GS_UseLowerTriangle) {
88 BlockedLinearOp L = getLowerTriBlocks(blo);
90 return createBlockLowerTriInverseOp(L, invDiag,
"Gauss Seidel");
93 TEUCHOS_ASSERT(
false);
98 const Teuchos::ParameterList& pl) {
99 Teko_DEBUG_SCOPE(
"GaussSeidelPreconditionerFactory::initializeFromParameterList", 10);
100 Teko_DEBUG_MSG_BEGIN(9);
101 DEBUG_STREAM <<
"Parameter list: " << std::endl;
102 pl.print(DEBUG_STREAM);
103 Teko_DEBUG_MSG_END();
105 const std::string inverse_type =
"Inverse Type";
106 const std::string preconditioner_type =
"Preconditioner Type";
107 std::vector<RCP<InverseFactory> > inverses;
108 std::vector<RCP<InverseFactory> > preconditioners;
113 std::string invStr =
"";
114 #if defined(Teko_ENABLE_Amesos)
116 #elif defined(Teko_ENABLE_Amesos2)
119 std::string precStr =
"None";
120 if (pl.isParameter(inverse_type)) invStr = pl.get<std::string>(inverse_type);
121 if (pl.isParameter(preconditioner_type)) precStr = pl.get<std::string>(preconditioner_type);
122 if (pl.isParameter(
"Use Upper Triangle"))
123 solveType_ = pl.get<
bool>(
"Use Upper Triangle") ? GS_UseUpperTriangle : GS_UseLowerTriangle;
125 Teko_DEBUG_MSG(
"GSPrecFact: Building default inverse \"" << invStr <<
"\"", 5);
126 RCP<InverseFactory> defaultInverse = invLib->getInverseFactory(invStr);
127 RCP<InverseFactory> defaultPrec;
128 if (precStr !=
"None") defaultPrec = invLib->getInverseFactory(precStr);
131 Teuchos::ParameterList::ConstIterator itr;
132 for (itr = pl.begin(); itr != pl.end(); ++itr) {
133 std::string fieldName = itr->first;
134 Teko_DEBUG_MSG(
"GSPrecFact: checking fieldName = \"" << fieldName <<
"\"", 9);
137 if (fieldName.compare(0, inverse_type.length(), inverse_type) == 0 &&
138 fieldName != inverse_type) {
140 std::string inverse, type;
143 std::stringstream ss(fieldName);
144 ss >> inverse >> type >> position;
147 Teko_DEBUG_MSG(
"Gauss-Seidel \"Inverse Type\" must be a (strictly) positive integer", 1);
151 std::string invStr2 = pl.get<std::string>(fieldName);
152 Teko_DEBUG_MSG(
"GSPrecFact: Building inverse " << position <<
" \"" << invStr2 <<
"\"", 5);
153 if (position > (
int)inverses.size()) {
154 inverses.resize(position, defaultInverse);
155 inverses[position - 1] = invLib->getInverseFactory(invStr2);
157 inverses[position - 1] = invLib->getInverseFactory(invStr2);
158 }
else if (fieldName.compare(0, preconditioner_type.length(), preconditioner_type) == 0 &&
159 fieldName != preconditioner_type) {
161 std::string preconditioner, type;
164 std::stringstream ss(fieldName);
165 ss >> preconditioner >> type >> position;
168 Teko_DEBUG_MSG(
"Gauss-Seidel \"Preconditioner Type\" must be a (strictly) positive integer",
173 std::string precStr2 = pl.get<std::string>(fieldName);
175 "GSPrecFact: Building preconditioner " << position <<
" \"" << precStr2 <<
"\"", 5);
176 if (position > (
int)preconditioners.size()) {
177 preconditioners.resize(position, defaultPrec);
178 preconditioners[position - 1] = invLib->getInverseFactory(precStr2);
180 preconditioners[position - 1] = invLib->getInverseFactory(precStr2);
185 if (inverses.size() == 0) inverses.push_back(defaultInverse);
An implementation of a state object for block preconditioners.
Teuchos::RCP< const BlockInvDiagonalStrategy > invOpsStrategy_
some members
LinearOp buildPreconditionerOperator(BlockedLinearOp &blo, BlockPreconditionerState &state) const
Create the Gauss-Seidel preconditioner operator.
GaussSeidelPreconditionerFactory()
Teuchos::RCP< const InverseLibrary > getInverseLibrary() const
Get the inverse library used by this preconditioner factory.
virtual void initializeFromParameterList(const Teuchos::ParameterList &pl)
Initialize from a parameter list.