47#ifndef PACKAGES_MUELU_ADAPTERS_XPETRA_MUELU_CREATEXPETRAPRECONDITIONER_HPP_
48#define PACKAGES_MUELU_ADAPTERS_XPETRA_MUELU_CREATEXPETRAPRECONDITIONER_HPP_
53#include <Teuchos_XMLParameterListHelpers.hpp>
54#include <Xpetra_CrsMatrix.hpp>
55#include <Xpetra_MultiVector.hpp>
60#include <MueLu_Hierarchy.hpp>
62#include <MueLu_MLParameterListInterpreter.hpp>
63#include <MueLu_ParameterListInterpreter.hpp>
64#include <MueLu_Utilities.hpp>
65#include <MueLu_HierarchyUtils.hpp>
79 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
80 Teuchos::RCP<MueLu::Hierarchy<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
82 const Teuchos::ParameterList& inParamList) {
94 bool hasParamList = inParamList.numParams();
96 RCP<HierarchyManager> mueLuFactory;
99 Teuchos::ParameterList nonSerialList,paramList;
103 if (hasParamList && paramList.isParameter(
"hierarchy label")) {
104 label = paramList.get<std::string>(
"hierarchy label");
105 paramList.remove(
"hierarchy label");
107 label = op->getObjectLabel();
109 std::string timerName;
111 timerName =
"MueLu setup time (" + label +
")";
113 timerName =
"MueLu setup time";
114 RCP<Teuchos::Time> tm = Teuchos::TimeMonitor::getNewTimer(timerName);
117 std::string syntaxStr =
"parameterlist: syntax";
118 if (hasParamList && paramList.isParameter(syntaxStr) && paramList.get<std::string>(syntaxStr) ==
"ml") {
119 paramList.remove(syntaxStr);
126 RCP<Hierarchy> H = mueLuFactory->CreateHierarchy(label);
127 H->setlib(op->getDomainMap()->lib());
130 H->GetLevel(0)->Set(
"A", op);
131 H->SetProcRankVerbose(op->getDomainMap()->getComm()->getRank());
136 mueLuFactory->SetupHierarchy(*H);
139 tm->incrementNumCalls();
142 const bool alwaysWriteLocal =
true;
143 const bool writeGlobalStats =
true;
144 const bool writeZeroTimers =
false;
145 const bool ignoreZeroTimers =
true;
146 const std::string filter = timerName;
147 Teuchos::TimeMonitor::summarize(op->getRowMap()->getComm().ptr(), H->GetOStream(
Statistics0), alwaysWriteLocal, writeGlobalStats,
148 writeZeroTimers, Teuchos::Union, filter, ignoreZeroTimers);
164 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
165 Teuchos::RCP<MueLu::Hierarchy<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
167 const std::string& xmlFileName) {
168 Teuchos::ParameterList paramList;
169 Teuchos::updateParametersFromXmlFileAndBroadcast(xmlFileName, Teuchos::Ptr<Teuchos::ParameterList>(¶mList), *op->getDomainMap()->getComm());
170 return CreateXpetraPreconditioner<Scalar, LocalOrdinal, GlobalOrdinal, Node>(op, paramList);
180 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
181 Teuchos::RCP<MueLu::Hierarchy<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
183 Teuchos::ParameterList paramList;
184 return CreateXpetraPreconditioner<Scalar, LocalOrdinal, GlobalOrdinal, Node>(op, paramList);
194 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
197 std::string label = H->GetLevel(0)->getObjectLabel();
199 std::string timerName;
201 timerName =
"MueLu setup time (" + label +
")";
203 timerName =
"MueLu setup time";
204 RCP<Teuchos::Time> tm = Teuchos::TimeMonitor::getNewTimer(timerName);
212 typedef Xpetra::Matrix<SC,LO,GO,NO> Matrix;
213 typedef Xpetra::Operator<SC,LO,GO,NO> Operator;
216 "MueLu::ReuseXpetraPreconditioner: Hierarchy has no levels in it");
218 "MueLu::ReuseXpetraPreconditioner: Hierarchy has no fine level operator");
219 RCP<Level> level0 = H->GetLevel(0);
221 RCP<Operator> O0 = level0->Get<RCP<Operator> >(
"A");
222 RCP<Matrix> A0 = Teuchos::rcp_dynamic_cast<Matrix>(O0);
228 A->SetFixedBlockSize(A0->GetFixedBlockSize());
235 tm->incrementNumCalls();
238 const bool alwaysWriteLocal =
true;
239 const bool writeGlobalStats =
true;
240 const bool writeZeroTimers =
false;
241 const bool ignoreZeroTimers =
true;
242 const std::string filter = timerName;
243 Teuchos::TimeMonitor::summarize(A->getRowMap()->getComm().ptr(), H->GetOStream(
Statistics0), alwaysWriteLocal, writeGlobalStats,
244 writeZeroTimers, Teuchos::Union, filter, ignoreZeroTimers);
254#define XPETRA_CREATEXPETRAPRECONDITIONER_SHORT
MueLu::DefaultLocalOrdinal LocalOrdinal
MueLu::DefaultScalar Scalar
MueLu::DefaultGlobalOrdinal GlobalOrdinal
Exception throws to report errors in the internal logical of the program.
static void AddNonSerializableDataToHierarchy(HierarchyManager &HM, Hierarchy &H, const ParameterList &nonSerialList)
Add non-serializable data to Hierarchy.
Provides methods to build a multigrid hierarchy and apply multigrid cycles.
Class that accepts ML-style parameters and builds a MueLu preconditioner. This interpreter uses the s...
Namespace for MueLu classes and methods.
long ExtractNonSerializableData(const Teuchos::ParameterList &inList, Teuchos::ParameterList &serialList, Teuchos::ParameterList &nonSerialList)
Extract non-serializable data from level-specific sublists and move it to a separate parameter list.
@ Statistics0
Print statistics that do not involve significant additional computation.
void ReuseXpetraPreconditioner(const Teuchos::RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &A, Teuchos::RCP< MueLu::Hierarchy< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &H)
Helper function to reuse an existing MueLu preconditioner.
Teuchos::RCP< MueLu::Hierarchy< Scalar, LocalOrdinal, GlobalOrdinal, Node > > CreateXpetraPreconditioner(Teuchos::RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > op, const Teuchos::ParameterList &inParamList)
Helper function to create a MueLu preconditioner that can be used by Xpetra.Given an Xpetra::Matrix,...