43 #ifndef IFPACK2_DETAILS_ONELEVELFACTORY_DEF_HPP 44 #define IFPACK2_DETAILS_ONELEVELFACTORY_DEF_HPP 46 #include "Ifpack2_Factory.hpp" 48 #include "Ifpack2_Chebyshev.hpp" 49 #include "Ifpack2_Details_DenseSolver.hpp" 50 #include "Ifpack2_Diagonal.hpp" 51 #include "Ifpack2_IdentitySolver.hpp" 52 #include "Ifpack2_ILUT.hpp" 53 #include "Ifpack2_Relaxation.hpp" 54 #include "Ifpack2_RILUK.hpp" 55 #include "Ifpack2_Experimental_RBILUK.hpp" 56 #include "Ifpack2_BlockRelaxation.hpp" 57 #include "Ifpack2_BandedContainer.hpp" 58 #include "Ifpack2_DenseContainer.hpp" 59 #include "Ifpack2_SparseContainer.hpp" 60 #include "Ifpack2_TriDiContainer.hpp" 61 #include "Ifpack2_LocalSparseTriangularSolver.hpp" 62 #include "Ifpack2_Hiptmair.hpp" 64 #ifdef HAVE_IFPACK2_SHYLU_NODEFASTILU 65 #include "Ifpack2_Details_Fic.hpp" 66 #include "Ifpack2_Details_Fildl.hpp" 67 #include "Ifpack2_Details_Filu.hpp" 68 #endif // HAVE_IFPACK2_SHYLU_NODEFASTILU 70 #ifdef HAVE_IFPACK2_AMESOS2 71 # include "Ifpack2_Details_Amesos2Wrapper.hpp" 72 #endif // HAVE_IFPACK2_AMESOS2 74 #ifdef HAVE_IFPACK2_HYPRE 75 # include "Ifpack2_Hypre.hpp" 76 #endif // HAVE_IFPACK2_HYPRE 82 template<
class MatrixType>
83 Teuchos::RCP<typename OneLevelFactory<MatrixType>::prec_type>
85 const Teuchos::RCP<const row_matrix_type>& matrix)
const 93 std::string precTypeUpper = canonicalize(precType);
95 if (precTypeUpper ==
"CHEBYSHEV") {
98 prec = rcp (new ::Ifpack2::Chebyshev<row_matrix_type> (matrix));
100 else if (precTypeUpper ==
"DENSE" || precTypeUpper ==
"LAPACK") {
103 else if (precTypeUpper ==
"AMESOS2") {
104 #ifdef HAVE_IFPACK2_AMESOS2 107 TEUCHOS_TEST_FOR_EXCEPTION(
108 true, std::invalid_argument,
"Ifpack2::Details::OneLevelFactory: " 109 "You may not ask for the preconditioner \"AMESOS2\" unless " 110 "you have built Trilinos with the Amesos2 package enabled.");
111 #endif // HAVE_IFPACK2_AMESOS2 113 else if (precTypeUpper ==
"DIAGONAL") {
114 prec = rcp (
new Diagonal<row_matrix_type> (matrix));
116 else if (precTypeUpper ==
"ILUT") {
119 else if (precTypeUpper ==
"RELAXATION") {
122 else if (precTypeUpper ==
"RILUK") {
125 else if (precTypeUpper ==
"RBILUK") {
128 else if (precTypeUpper ==
"FAST_IC" || precTypeUpper ==
"FAST_ILU" || precTypeUpper ==
"FAST_ILDL") {
129 #ifdef HAVE_IFPACK2_SHYLU_NODEFASTILU 131 if(precTypeUpper ==
"FAST_IC")
133 else if(precTypeUpper ==
"FAST_ILU")
135 else if(precTypeUpper ==
"FAST_ILDL")
140 throw std::invalid_argument(
"The Ifpack2 FastIC, FastILU and FastILDL preconditioners require the FastILU subpackage of ShyLU to be enabled\n" 141 "To enable FastILU, set the CMake option Trilinos_ENABLE_ShyLU_NodeFastILU=ON");
145 else if (precTypeUpper ==
"KRYLOV") {
146 TEUCHOS_TEST_FOR_EXCEPTION
147 (
true, std::invalid_argument,
"The \"KRYLOV\" preconditioner option has " 148 "been deprecated and removed. If you want a Krylov solver, use the " 151 else if (precTypeUpper ==
"BLOCK_RELAXATION" ||
152 precTypeUpper ==
"BLOCK RELAXATION" ||
153 precTypeUpper ==
"BLOCKRELAXATION" ||
154 precTypeUpper ==
"DENSE_BLOCK_RELAXATION" ||
155 precTypeUpper ==
"DENSE BLOCK RELAXATION" ||
156 precTypeUpper ==
"DENSEBLOCKRELAXATION" ) {
161 Teuchos::ParameterList params;
162 params.set (
"relaxation: container",
"Dense");
163 prec->setParameters (params);
165 else if (precTypeUpper ==
"SPARSE_BLOCK_RELAXATION" ||
166 precTypeUpper ==
"SPARSE BLOCK RELAXATION" ||
167 precTypeUpper ==
"SPARSEBLOCKRELAXATION" ) {
174 #ifdef HAVE_IFPACK2_AMESOS2 176 Teuchos::ParameterList params;
177 params.set (
"relaxation: container",
"SparseAmesos2");
178 prec->setParameters (params);
180 TEUCHOS_TEST_FOR_EXCEPTION
181 (
true, std::invalid_argument,
"Ifpack2::Details::OneLevelFactory: " 182 "\"SPARSE BLOCK RELAXATION\" requires building Trilinos with Amesos2 enabled.");
185 else if (precTypeUpper ==
"TRIDI_RELAXATION" ||
186 precTypeUpper ==
"TRIDI RELAXATION" ||
187 precTypeUpper ==
"TRIDIRELAXATION" ||
188 precTypeUpper ==
"TRIDIAGONAL_RELAXATION" ||
189 precTypeUpper ==
"TRIDIAGONAL RELAXATION" ||
190 precTypeUpper ==
"TRIDIAGONALRELAXATION") {
192 Teuchos::ParameterList params;
193 params.set (
"relaxation: container",
"TriDi");
194 prec->setParameters (params);
196 else if (precTypeUpper ==
"BANDED_RELAXATION" ||
197 precTypeUpper ==
"BANDED RELAXATION" ||
198 precTypeUpper ==
"BANDEDRELAXATION") {
200 Teuchos::ParameterList params;
201 params.set (
"relaxation: container",
"Banded");
202 prec->setParameters (params);
204 else if (precTypeUpper ==
"IDENTITY" || precTypeUpper ==
"IDENTITY_SOLVER") {
208 else if (precTypeUpper ==
"LOCAL SPARSE TRIANGULAR SOLVER" ||
209 precTypeUpper ==
"LOCAL_SPARSE_TRIANGULAR_SOLVER" ||
210 precTypeUpper ==
"LOCALSPARSETRIANGULARSOLVER" ||
211 precTypeUpper ==
"SPARSE TRIANGULAR SOLVER" ||
212 precTypeUpper ==
"SPARSE_TRIANGULAR_SOLVER" ||
213 precTypeUpper ==
"SPARSETRIANGULARSOLVER") {
216 else if(precTypeUpper ==
"HIPTMAIR") {
219 #ifdef HAVE_IFPACK2_HYPRE 220 else if (precTypeUpper ==
"HYPRE") {
221 prec = rcp (
new Hypre<row_matrix_type> (matrix));
225 TEUCHOS_TEST_FOR_EXCEPTION(
226 true, std::invalid_argument,
"Ifpack2::Details::OneLevelFactory::create: " 227 "Invalid preconditioner type \"" << precType <<
"\".");
230 TEUCHOS_TEST_FOR_EXCEPTION(
231 prec.is_null (), std::logic_error,
"Ifpack2::Details::OneLevelFactory::" 232 "create: Return value is null right before return. This should never " 233 "happen. Please report this bug to the Ifpack2 developers.");
237 template<
class MatrixType>
242 std::string precTypeUpper = canonicalize(precType);
243 std::vector<std::string> supportedNames = {
244 "CHEBYSHEV",
"DENSE",
"LAPACK",
245 #ifdef HAVE_IFPACK2_AMESOS2 248 "DIAGONAL",
"ILUT",
"RELAXATION",
"RILUK",
"RBILUK",
249 #ifdef HAVE_IFPACK2_SHYLU_NODEFASTILU 250 "FAST_IC",
"FAST_ILU",
"FAST_ILDL",
252 "BLOCK_RELAXATION",
"BLOCK RELAXATION",
"BLOCKRELAXATION",
"DENSE_BLOCK_RELAXATION",
"DENSE BLOCK RELAXATION",
"DENSEBLOCKRELAXATION",
253 #ifdef HAVE_IFPACK2_AMESOS2 254 "SPARSE_BLOCK_RELAXATION",
"SPARSE BLOCK RELAXATION",
"SPARSEBLOCKRELAXATION",
256 #ifdef HAVE_IFPACK2_HYPRE 259 "TRIDI_RELAXATION",
"TRIDI RELAXATION",
"TRIDIRELAXATION",
"TRIDIAGONAL_RELAXATION",
"TRIDIAGONAL RELAXATION",
"TRIDIAGONALRELAXATION",
260 "BANDED_RELAXATION",
"BANDED RELAXATION",
"BANDEDRELAXATION",
261 "IDENTITY",
"IDENTITY_SOLVER",
262 "LOCAL SPARSE TRIANGULAR SOLVER",
"LOCAL_SPARSE_TRIANGULAR_SOLVER",
"LOCALSPARSETRIANGULARSOLVER",
"SPARSE TRIANGULAR SOLVER",
"SPARSE_TRIANGULAR_SOLVER",
"SPARSETRIANGULARSOLVER",
267 auto it = std::find(std::begin(supportedNames), std::end(supportedNames), precTypeUpper);
268 return it != std::end(supportedNames);
274 #define IFPACK2_DETAILS_ONELEVELFACTORY_INSTANT(S,LO,GO,N) \ 275 template class Ifpack2::Details::OneLevelFactory< Tpetra::RowMatrix<S, LO, GO, N> >; 277 #endif // IFPACK2_DETAILS_ONELEVELFACTORY_DEF_HPP ILU(k) factorization of a given Tpetra::RowMatrix.
Definition: Ifpack2_RILUK_decl.hpp:245
"Preconditioner" that uses LAPACK's dense LU.
Definition: Ifpack2_Details_DenseSolver_decl.hpp:75
The Ifpack2 wrapper to the ILDL preconditioner of ShyLU FastILU.
Definition: Ifpack2_Details_Fildl_decl.hpp:62
The Ifpack2 wrapper to the ILU preconditioner of ShyLU FastILU.
Definition: Ifpack2_Details_Filu_decl.hpp:62
ILUT (incomplete LU factorization with threshold) of a Tpetra sparse matrix.
Definition: Ifpack2_ILUT_decl.hpp:91
Block relaxation preconditioners (or smoothers) for Tpetra::RowMatrix and Tpetra::CrsMatrix sparse ma...
Definition: Ifpack2_BlockRelaxation_decl.hpp:82
Ifpack2 implementation details.
"Identity" preconditioner.
Definition: Ifpack2_IdentitySolver_decl.hpp:60
Wrapper for Hiptmair smoothers.
Definition: Ifpack2_Hiptmair_decl.hpp:71
Wrapper class for direct solvers in Amesos2.
Definition: Ifpack2_Details_Amesos2Wrapper_decl.hpp:102
"Preconditioner" that solves local sparse triangular systems.
Definition: Ifpack2_LocalSparseTriangularSolver_decl.hpp:79
File for utility functions.
The Ifpack2 wrapper to the incomplete Chebyshev preconditioner of ShyLU FastILU.
Definition: Ifpack2_Details_Fic_decl.hpp:62
"Factory" for creating single-level preconditioners.
Definition: Ifpack2_Details_OneLevelFactory_decl.hpp:124
Relaxation preconditioners for Tpetra::RowMatrix and Tpetra::CrsMatrix sparse matrices.
Definition: Ifpack2_Relaxation_decl.hpp:237
Preconditioners and smoothers for Tpetra sparse matrices.
Definition: Ifpack2_AdditiveSchwarz_decl.hpp:73
ILU(k) factorization of a given Tpetra::BlockCrsMatrix.
Definition: Ifpack2_Experimental_RBILUK_decl.hpp:128
Teuchos::RCP< prec_type > create(const std::string &precType, const Teuchos::RCP< const row_matrix_type > &matrix) const
Create an instance of Preconditioner given the string name of the preconditioner type.
Definition: Ifpack2_Details_OneLevelFactory_def.hpp:84