46 #ifndef MUELU_HIERARCHYUTILS_DEF_HPP 47 #define MUELU_HIERARCHYUTILS_DEF_HPP 49 #include "Teuchos_ScalarTraits.hpp" 51 #include <Xpetra_Matrix.hpp> 52 #include <Xpetra_Operator.hpp> 57 #include "MueLu_SmootherFactory.hpp" 58 #include "MueLu_FactoryManager.hpp" 61 #ifdef HAVE_MUELU_INTREPID2 70 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
72 typedef typename Xpetra::MultiVector<typename Teuchos::ScalarTraits<Scalar>::coordinateType,
75 for (ParameterList::ConstIterator it = paramList.begin(); it != paramList.end(); it++) {
76 const std::string& levelName = it->first;
78 if (paramList.isSublist(levelName) && levelName.find(
"level ") == 0 && levelName.size() > 6) {
79 int levelID = strtol(levelName.substr(6).c_str(), 0, 0);
86 RCP<Level> level = H.
GetLevel(levelID);
89 TEUCHOS_TEST_FOR_EXCEPTION(M.is_null(),
Exceptions::InvalidArgument,
"MueLu::Utils::AddNonSerializableDataToHierarchy: cannot get FactoryManager");
92 const ParameterList& levelList = paramList.sublist(levelName);
93 for (ParameterList::ConstIterator it2 = levelList.begin(); it2 != levelList.end(); it2++) {
94 const std::string& name = it2->first;
95 TEUCHOS_TEST_FOR_EXCEPTION(name !=
"A" && name !=
"P" && name !=
"R" && name !=
"K" && name !=
"M" && name !=
"Mdiag" &&
96 name !=
"Nullspace" && name !=
"Coordinates" && name !=
"pcoarsen: element to node map" &&
97 name !=
"Node Comm" && name !=
"DualNodeID2PrimalNodeID" &&
99 std::string(
"MueLu::Utils::AddNonSerializableDataToHierarchy: parameter list contains unknown data type(") + name +
")");
102 level->Set(name, Teuchos::getValue<RCP<Matrix > > (it2->second),
NoFactory::get());
107 else if(name ==
"P" || name ==
"R" || name ==
"K" || name ==
"M") {
109 level->Set(name, Teuchos::getValue<RCP<Matrix > > (it2->second),
NoFactory::get());
111 else if (name ==
"Mdiag")
114 level->Set(name, Teuchos::getValue<RCP<Vector > > (it2->second),
NoFactory::get());
116 else if (name ==
"Nullspace")
119 level->Set(name, Teuchos::getValue<RCP<MultiVector > >(it2->second),
NoFactory::get());
123 else if(name ==
"Coordinates")
126 level->Set(name, Teuchos::getValue<RCP<realvaluedmultivector_type> >(it2->second),
NoFactory::get());
129 else if(name ==
"Node Comm")
134 else if(name ==
"DualNodeID2PrimalNodeID")
137 level->Set(name, Teuchos::getValue<RCP<std::map<LO, LO>>>(it2->second),
NoFactory::get());
139 #ifdef HAVE_MUELU_INTREPID2 140 else if (name ==
"pcoarsen: element to node map")
143 level->Set(name, Teuchos::getValue<RCP<Kokkos::DynRankView<LocalOrdinal,typename Node::device_type> > >(it2->second),
NoFactory::get());
147 #ifdef HAVE_MUELU_MATLAB 150 size_t typeNameStart = name.find_first_not_of(
' ');
151 size_t typeNameEnd = name.find(
' ', typeNameStart);
152 std::string typeName = name.substr(typeNameStart, typeNameEnd - typeNameStart);
153 std::transform(typeName.begin(), typeName.end(), typeName.begin(),
::tolower);
155 if(typeName ==
"matrix")
156 level->Set(name, Teuchos::getValue<RCP<Matrix> >(it2->second),
NoFactory::get());
157 else if(typeName ==
"multivector")
158 level->Set(name, Teuchos::getValue<RCP<MultiVector> >(it2->second),
NoFactory::get());
159 else if(typeName ==
"map")
160 level->Set(name, Teuchos::getValue<RCP<Xpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > >(it2->second),
NoFactory::get());
161 else if(typeName ==
"ordinalvector")
162 level->Set(name, Teuchos::getValue<RCP<Xpetra::Vector<LocalOrdinal, LocalOrdinal, GlobalOrdinal, Node> > >(it2->second),
NoFactory::get());
163 else if(typeName ==
"scalar")
164 level->Set(name, Teuchos::getValue<Scalar>(it2->second),
NoFactory::get());
165 else if(typeName ==
"double")
166 level->Set(name, Teuchos::getValue<double>(it2->second),
NoFactory::get());
167 else if(typeName ==
"complex")
168 level->Set(name, Teuchos::getValue<std::complex<double> >(it2->second),
NoFactory::get());
169 else if(typeName ==
"int")
170 level->Set(name, Teuchos::getValue<int>(it2->second),
NoFactory::get());
171 else if(typeName ==
"string")
172 level->Set(name, Teuchos::getValue<std::string>(it2->second),
NoFactory::get());
176 throw std::runtime_error(
"Invalid non-serializable data on list");
180 }
else if (paramList.isSublist(levelName) && levelName.find(
"user data") != std::string::npos) {
183 RCP<Level> level = H.
GetLevel(levelID);
186 TEUCHOS_TEST_FOR_EXCEPTION(M.is_null(),
Exceptions::InvalidArgument,
"MueLu::Utils::AddNonSerializableDataToHierarchy: cannot get FactoryManager");
189 const ParameterList& userList = paramList.sublist(levelName);
190 for (ParameterList::ConstIterator it2 = userList.begin(); it2 != userList.end(); it2++) {
191 const std::string& name = it2->first;
192 TEUCHOS_TEST_FOR_EXCEPTION(name !=
"P" && name !=
"R" && name !=
"K" && name !=
"M" && name !=
"Mdiag" &&
193 name !=
"Nullspace" && name !=
"Coordinates" && name !=
"pcoarsen: element to node map" &&
194 name !=
"Node Comm" && name !=
"DualNodeID2PrimalNodeID" &&
195 name !=
"output stream" &&
197 std::string(
"MueLu::Utils::AddNonSerializableDataToHierarchy: user data parameter list contains unknown data type (") + name +
")");
198 if( name ==
"P" || name ==
"R" || name ==
"K" || name ==
"M") {
200 level->Set(name, Teuchos::getValue<RCP<Matrix > > (it2->second),
NoFactory::get());
201 }
else if (name ==
"Mdiag") {
203 level->Set(name, Teuchos::getValue<RCP<Vector > >(it2->second),
NoFactory::get());
204 }
else if (name ==
"Nullspace") {
206 level->Set(name, Teuchos::getValue<RCP<MultiVector > >(it2->second),
NoFactory::get());
209 }
else if(name ==
"Coordinates") {
211 level->Set(name, Teuchos::getValue<RCP<realvaluedmultivector_type> >(it2->second),
NoFactory::get());
212 level->print(std::cout, MueLu::VERB_EXTREME);
214 else if(name ==
"Node Comm") {
218 else if(name ==
"DualNodeID2PrimalNodeID")
221 level->Set(name, Teuchos::getValue<RCP<std::map<LO, LO>>>(it2->second),
NoFactory::get());
223 #ifdef HAVE_MUELU_INTREPID2 224 else if (name ==
"pcoarsen: element to node map")
227 level->Set(name, Teuchos::getValue<RCP<Kokkos::DynRankView<LocalOrdinal,typename Node::device_type> > >(it2->second),
NoFactory::get());
230 else if (name ==
"output stream")
232 H.
SetMueLuOStream(Teuchos::getValue<RCP<Teuchos::FancyOStream> >(it2->second));
236 size_t typeNameStart = name.find_first_not_of(
' ');
237 size_t typeNameEnd = name.find(
' ', typeNameStart);
238 std::string typeName = name.substr(typeNameStart, typeNameEnd - typeNameStart);
239 size_t varNameStart = name.find_first_not_of(
' ', typeNameEnd);
240 std::string varName = name.substr(varNameStart, name.size());
241 std::transform(typeName.begin(), typeName.end(), typeName.begin(),
::tolower);
243 if(typeName ==
"matrix")
244 level->Set(varName, Teuchos::getValue<RCP<Matrix> >(it2->second),
NoFactory::get());
245 else if(typeName ==
"multivector")
246 level->Set(varName, Teuchos::getValue<RCP<MultiVector> >(it2->second),
NoFactory::get());
247 else if(typeName ==
"vector")
248 level->Set(varName, Teuchos::getValue<RCP<Vector> >(it2->second),
NoFactory::get());
249 else if(typeName ==
"map")
250 level->Set(varName, Teuchos::getValue<RCP<Xpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > >(it2->second),
NoFactory::get());
251 else if(typeName ==
"ordinalvector")
252 level->Set(varName, Teuchos::getValue<RCP<Xpetra::Vector<LocalOrdinal, LocalOrdinal, GlobalOrdinal, Node> > >(it2->second),
NoFactory::get());
253 else if(typeName ==
"scalar")
254 level->Set(varName, Teuchos::getValue<Scalar>(it2->second),
NoFactory::get());
255 else if(typeName ==
"double")
256 level->Set(varName, Teuchos::getValue<double>(it2->second),
NoFactory::get());
257 else if(typeName ==
"complex")
258 level->Set(varName, Teuchos::getValue<std::complex<double> >(it2->second),
NoFactory::get());
259 else if(typeName ==
"int")
260 level->Set(varName, Teuchos::getValue<int>(it2->second),
NoFactory::get());
261 else if(typeName ==
"string")
262 level->Set(varName, Teuchos::getValue<std::string>(it2->second),
NoFactory::get());
263 else if(typeName ==
"array<go>")
264 level->Set(varName, Teuchos::getValue<Array<GlobalOrdinal> > (it2->second),
NoFactory::get());
265 else if(typeName ==
"array<lo>")
266 level->Set(varName, Teuchos::getValue<Array<LocalOrdinal> >(it2->second),
NoFactory::get());
267 else if(typeName ==
"arrayrcp<lo>")
268 level->Set(varName, Teuchos::getValue<ArrayRCP<LocalOrdinal> >(it2->second),
NoFactory::get());
269 else if(typeName ==
"arrayrcp<go>")
270 level->Set(varName, Teuchos::getValue<ArrayRCP<GlobalOrdinal> >(it2->second),
NoFactory::get());
272 throw std::runtime_error(
"Invalid non-serializable data on list");
280 #define MUELU_HIERARCHY_UTILS_SHORT 281 #endif // MUELU_HIERARCHYHELPERS_DEF_HPP This class specifies the default factory that should generate some data on a Level if the data does n...
MueLu::DefaultLocalOrdinal LocalOrdinal
RCP< Level > & GetLevel(const int levelID=0)
Retrieve a certain level from hierarchy.
static void AddNonSerializableDataToHierarchy(HierarchyManager &HM, Hierarchy &H, const ParameterList &nonSerialList)
Add non-serializable data to Hierarchy.
void AddNewLevel()
Add a new level at the end of the hierarchy.
User data are always kept. This flag is set automatically when Level::Set("data", data) is used...
std::string tolower(const std::string &str)
Namespace for MueLu classes and methods.
bool IsParamMuemexVariable(const std::string &name)
RCP< FactoryManagerBase > GetFactoryManager(int levelID) const
static const NoFactory * get()
MueLu::DefaultGlobalOrdinal GlobalOrdinal
static void SetMueLuOStream(const Teuchos::RCP< Teuchos::FancyOStream > &mueluOStream)
bool IsParamValidVariable(const std::string &name)
Provides methods to build a multigrid hierarchy and apply multigrid cycles.
Exception throws to report invalid user entry.
static const RCP< const NoFactory > getRCP()
Static Get() functions.