45 #ifndef __IFPACK2_FILU_DEF_HPP__ 46 #define __IFPACK2_FILU_DEF_HPP__ 48 #include "Ifpack2_Details_Filu_decl.hpp" 50 #include <impl/Kokkos_Timer.hpp> 51 #include <shylu_fastilu.hpp> 58 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
60 Filu(Teuchos::RCP<const TRowMatrix> A) :
61 FastILU_Base<Scalar, LocalOrdinal, GlobalOrdinal, Node>(A) {}
63 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
67 return localPrec_->getNFact();
70 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
74 return localPrec_->getNTrisol();
77 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
81 localPrec_->checkILU();
84 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
88 localPrec_->checkIC();
91 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
95 auto nRows = this->mat_->getNodeNumRows();
96 auto& p = this->params_;
97 localPrec_ = Teuchos::rcp(
new LocalFILU(this->localRowPtrs_, this->localColInds_, this->localValues_, nRows,
98 p.nFact, p.nTrisol, p.level, p.omega,
99 p.shift, p.guessFlag ? 1 : 0, p.blockSize));
100 localPrec_->initialize();
101 this->initTime_ = localPrec_->getInitializeTime();
104 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
109 localPrec_->setValues(this->localValues_);
110 localPrec_->compute();
111 this->computeTime_ = localPrec_->getComputeTime();
114 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
118 localPrec_->apply(x, y);
120 this->applyTime_ += localPrec_->getApplyTime();
123 template<
typename Scalar,
typename LocalOrdinal,
typename GlobalOrdinal,
typename Node>
130 #define IFPACK2_DETAILS_FILU_INSTANT(S, L, G, N) \ 131 template class Ifpack2::Details::Filu<S, L, G, N>; void applyLocalPrec(ScalarArray x, ScalarArray y) const
Apply the local preconditioner with 1-D views of the local parts of X and Y (one vector only) ...
Definition: Ifpack2_Details_Filu_def.hpp:116
int getSweeps() const
Get the sweeps ("nFact") from localPrec_.
Definition: Ifpack2_Details_Filu_def.hpp:65
void initLocalPrec()
Construct the underlying preconditioner (localPrec_) using given params and then call localPrec_->ini...
Definition: Ifpack2_Details_Filu_def.hpp:93
std::string getName() const
Get the name of the underlying preconditioner ("Filu", "Fildl" or "Fic")
Definition: Ifpack2_Details_Filu_def.hpp:125
Ifpack2 implementation details.
Filu(Teuchos::RCP< const TRowMatrix > mat_)
Constructor.
Definition: Ifpack2_Details_Filu_def.hpp:60
The base class of the Ifpack2 FastILU wrappers (Filu, Fildl and Fic)
Definition: Ifpack2_Details_FastILU_Base_decl.hpp:66
void checkLocalIC() const
Verify and print debug info about the internal IC preconditioner.
Definition: Ifpack2_Details_Filu_def.hpp:86
void checkLocalILU() const
Verify and print debug info about the internal ILU preconditioner.
Definition: Ifpack2_Details_Filu_def.hpp:79
Kokkos::View< Scalar *, execution_space > ScalarArray
Array of Scalar on device.
Definition: Ifpack2_Details_FastILU_Base_decl.hpp:88
Preconditioners and smoothers for Tpetra sparse matrices.
Definition: Ifpack2_AdditiveSchwarz_decl.hpp:73
int getNTrisol() const
Get the number of triangular solves ("nTrisol") from localPrec_.
Definition: Ifpack2_Details_Filu_def.hpp:72
Provides functions for retrieving local CRS arrays (row pointers, column indices, and values) from Tp...
void computeLocalPrec()
Get values array from the matrix and then call compute() on the underlying preconditioner.
Definition: Ifpack2_Details_Filu_def.hpp:106