106 typedef typename Xpetra::MultiVector<typename Teuchos::ScalarTraits<Scalar>::coordinateType,
109 for (ParameterList::ConstIterator nonSerialEntry = nonSerialList.begin(); nonSerialEntry != nonSerialList.end(); nonSerialEntry++) {
110 const std::string& levelName = nonSerialEntry->first;
112 if (nonSerialList.isSublist(levelName) && levelName.find(
"level ") == 0 && levelName.size() > 6) {
113 int levelID = strtol(levelName.substr(6).c_str(), 0, 0);
120 RCP<Level> level = H.
GetLevel(levelID);
122 RCP<FactoryManager> M = Teuchos::rcp_dynamic_cast<FactoryManager>(HM.
GetFactoryManager(levelID));
123 TEUCHOS_TEST_FOR_EXCEPTION(M.is_null(),
Exceptions::InvalidArgument,
"MueLu::Utils::AddNonSerializableDataToHierarchy: cannot get FactoryManager");
126 const ParameterList& levelList = nonSerialList.sublist(levelName);
127 for (ParameterList::ConstIterator levelListEntry = levelList.begin(); levelListEntry != levelList.end(); levelListEntry++) {
128 const std::string& name = levelListEntry->first;
129 TEUCHOS_TEST_FOR_EXCEPTION(name !=
"A" && name !=
"P" && name !=
"R" && name !=
"K" && name !=
"M" && name !=
"Mdiag" &&
130 name !=
"D0" && name !=
"M1" && name !=
"Ms" && name !=
"M0inv" &&
131 name !=
"Pnodal" && name !=
"NodeMatrix" && name !=
"NodeAggMatrix" &&
132 name !=
"Nullspace" && name !=
"Coordinates" && name !=
"pcoarsen: element to node map" &&
133 name !=
"Node Comm" && name !=
"DualNodeID2PrimalNodeID" && name !=
"Primal interface DOF map" &&
135 std::string(
"MueLu::Utils::AddNonSerializableDataToHierarchy: parameter list contains unknown data type(") + name +
")");
138 RCP<const Teuchos::Comm<int> > comm;
139 if (!level->GetComm().is_null())
140 comm = level->GetComm();
141 else if (level->IsAvailable(
"A")) {
143 level->Get(
"A", mat);
144 comm = mat->getMap()->getComm();
147 if (!level0->GetComm().is_null())
148 comm = level0->GetComm();
151 level0->Get(
"A", mat);
152 comm = mat->getMap()->getComm();
155 Xpetra::UnderlyingLib lib = level->lib();
159 if (levelListEntry->second.isType<std::string>())
161 mat = Xpetra::IO<Scalar,LocalOrdinal,GlobalOrdinal,Node>::Read(Teuchos::getValue<std::string>(levelListEntry->second), lib, comm);
163 mat = Teuchos::getValue<RCP<Matrix > > (levelListEntry->second);
169 else if(name ==
"P" || name ==
"R" || name ==
"K" || name ==
"M" ) {
170 if (levelListEntry->second.isType<RCP<Operator> >()) {
172 mat = Teuchos::getValue<RCP<Operator> > (levelListEntry->second);
174 RCP<const FactoryBase> fact = M->GetFactory(name);
176 level->Set(name, mat, fact.get());
182 if (levelListEntry->second.isType<std::string>())
184 mat = Xpetra::IO<Scalar,LocalOrdinal,GlobalOrdinal,Node>::Read(Teuchos::getValue<std::string>(levelListEntry->second), lib, comm);
186 mat = Teuchos::getValue<RCP<Matrix > > (levelListEntry->second);
188 RCP<const FactoryBase> fact = M->GetFactory(name);
190 level->Set(name, mat, fact.get());
196 else if (name ==
"D0" || name ==
"M1" || name ==
"Ms" || name ==
"M0inv" || name ==
"Pnodal" || name ==
"NodeMatrix" || name ==
"NodeAggMatrix") {
198 if (levelListEntry->second.isType<RCP<Operator> >())
199 level->Set(name, Teuchos::getValue<RCP<Operator> > (levelListEntry->second),
NoFactory::get());
201 level->Set(name, Teuchos::getValue<RCP<Matrix> > (levelListEntry->second),
NoFactory::get());
203 else if (name ==
"Mdiag")
206 level->Set(name, Teuchos::getValue<RCP<Vector > > (levelListEntry->second),
NoFactory::get());
208 else if (name ==
"Nullspace")
210 RCP<MultiVector> vec;
211 if (levelListEntry->second.isType<std::string>()) {
212 TEUCHOS_ASSERT(level->IsAvailable(
"A"));
214 level->Get(
"A", mat);
215 auto map = mat->getMap();
216 vec = Xpetra::IO<Scalar,LocalOrdinal,GlobalOrdinal,Node>::ReadMultiVector(Teuchos::getValue<std::string>(levelListEntry->second), map);
218 vec = Teuchos::getValue<RCP<MultiVector> > (levelListEntry->second);
224 else if(name ==
"Coordinates")
226 RCP<realvaluedmultivector_type> vec;
227 if (levelListEntry->second.isType<std::string>()) {
228 TEUCHOS_ASSERT(level->IsAvailable(
"A"));
230 level->Get(
"A", mat);
231 size_t blkSize = mat->GetFixedBlockSize();
232 RCP<const Map> nodeMap = mat->getRowMap();
235 RCP<const Map> dofMap = mat->getRowMap();
236 GO indexBase = dofMap->getIndexBase();
237 size_t numLocalDOFs = dofMap->getLocalNumElements();
239 "HierarchyUtils: block size (" << blkSize <<
") is incompatible with the number of local dofs in a row map (" << numLocalDOFs);
240 ArrayView<const GO> GIDs = dofMap->getLocalElementList();
242 Array<GO> nodeGIDs(numLocalDOFs/blkSize);
243 for (
size_t i = 0; i < numLocalDOFs; i += blkSize)
244 nodeGIDs[i/blkSize] = (GIDs[i] - indexBase)/blkSize + indexBase;
246 Xpetra::global_size_t INVALID = Teuchos::OrdinalTraits<Xpetra::global_size_t>::invalid();
247 nodeMap = MapFactory::Build(dofMap->lib(), INVALID, nodeGIDs(), indexBase, dofMap->getComm());
249 vec = Xpetra::IO<typename Teuchos::ScalarTraits<Scalar>::coordinateType,
LocalOrdinal,
GlobalOrdinal,
Node>::ReadMultiVector(Teuchos::getValue<std::string>(levelListEntry->second), nodeMap);
251 vec = Teuchos::getValue<RCP<realvaluedmultivector_type > > (levelListEntry->second);
256 else if(name ==
"Node Comm")
261 else if(name ==
"DualNodeID2PrimalNodeID")
264 level->Set(name, Teuchos::getValue<RCP<std::map<LO, LO>>>(levelListEntry->second),
NoFactory::get());
266 else if(name ==
"Primal interface DOF map")
269 level->Set(name, Teuchos::getValue<RCP<const Map>>(levelListEntry->second),
NoFactory::get());
271#ifdef HAVE_MUELU_INTREPID2
272 else if (name ==
"pcoarsen: element to node map")
275 level->Set(name, Teuchos::getValue<RCP<Kokkos::DynRankView<LocalOrdinal,typename Node::device_type> > >(levelListEntry->second),
NoFactory::get());
279#ifdef HAVE_MUELU_MATLAB
282 size_t typeNameStart = name.find_first_not_of(
' ');
283 size_t typeNameEnd = name.find(
' ', typeNameStart);
284 std::string typeName = name.substr(typeNameStart, typeNameEnd - typeNameStart);
285 std::transform(typeName.begin(), typeName.end(), typeName.begin(), ::tolower);
287 if(typeName ==
"matrix")
288 level->Set(name, Teuchos::getValue<RCP<Matrix> >(levelListEntry->second),
NoFactory::get());
289 else if(typeName ==
"multivector")
290 level->Set(name, Teuchos::getValue<RCP<MultiVector> >(levelListEntry->second),
NoFactory::get());
291 else if(typeName ==
"map")
292 level->Set(name, Teuchos::getValue<RCP<Xpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > >(levelListEntry->second),
NoFactory::get());
293 else if(typeName ==
"ordinalvector")
294 level->Set(name, Teuchos::getValue<RCP<Xpetra::Vector<LocalOrdinal, LocalOrdinal, GlobalOrdinal, Node> > >(levelListEntry->second),
NoFactory::get());
295 else if(typeName ==
"scalar")
296 level->Set(name, Teuchos::getValue<Scalar>(levelListEntry->second),
NoFactory::get());
297 else if(typeName ==
"double")
298 level->Set(name, Teuchos::getValue<double>(levelListEntry->second),
NoFactory::get());
299 else if(typeName ==
"complex")
300 level->Set(name, Teuchos::getValue<std::complex<double> >(levelListEntry->second),
NoFactory::get());
301 else if(typeName ==
"int")
302 level->Set(name, Teuchos::getValue<int>(levelListEntry->second),
NoFactory::get());
303 else if(typeName ==
"string")
304 level->Set(name, Teuchos::getValue<std::string>(levelListEntry->second),
NoFactory::get());
308 throw std::runtime_error(
"Invalid non-serializable data on list");
312 }
else if (nonSerialList.isSublist(levelName) && levelName.find(
"user data") != std::string::npos) {
315 RCP<Level> level = H.
GetLevel(levelID);
317 RCP<FactoryManager> M = Teuchos::rcp_dynamic_cast<FactoryManager>(HM.
GetFactoryManager(levelID));
318 TEUCHOS_TEST_FOR_EXCEPTION(M.is_null(),
Exceptions::InvalidArgument,
"MueLu::Utils::AddNonSerializableDataToHierarchy: cannot get FactoryManager");
321 const ParameterList& userList = nonSerialList.sublist(levelName);
322 for (ParameterList::ConstIterator userListEntry = userList.begin(); userListEntry != userList.end(); userListEntry++) {
323 const std::string& name = userListEntry->first;
324 TEUCHOS_TEST_FOR_EXCEPTION(name !=
"P" && name !=
"R" && name !=
"K" && name !=
"M" && name !=
"Mdiag" &&
325 name !=
"D0" && name !=
"M1" && name !=
"Ms" && name !=
"M0inv" &&
326 name !=
"Nullspace" && name !=
"Coordinates" && name !=
"pcoarsen: element to node map" &&
327 name !=
"Node Comm" && name !=
"DualNodeID2PrimalNodeID" && name !=
"Primal interface DOF map" &&
328 name !=
"output stream" &&
330 std::string(
"MueLu::Utils::AddNonSerializableDataToHierarchy: user data parameter list contains unknown data type (") + name +
")");
331 if( name ==
"P" || name ==
"R" || name ==
"K" || name ==
"M" || name ==
"D0" || name ==
"M1" || name ==
"Ms" || name ==
"M0inv" ) {
333 level->Set(name, Teuchos::getValue<RCP<Matrix > > (userListEntry->second),
NoFactory::get());
334 }
else if (name ==
"Mdiag") {
336 level->Set(name, Teuchos::getValue<RCP<Vector > >(userListEntry->second),
NoFactory::get());
337 }
else if (name ==
"Nullspace") {
339 level->Set(name, Teuchos::getValue<RCP<MultiVector > >(userListEntry->second),
NoFactory::get());
342 }
else if(name ==
"Coordinates") {
344 level->Set(name, Teuchos::getValue<RCP<realvaluedmultivector_type> >(userListEntry->second),
NoFactory::get());
346 else if(name ==
"Node Comm") {
350 else if(name ==
"DualNodeID2PrimalNodeID")
353 level->Set(name, Teuchos::getValue<RCP<std::map<LO, LO>>>(userListEntry->second),
NoFactory::get());
355 else if(name ==
"Primal interface DOF map")
358 level->Set(name, Teuchos::getValue<RCP<const Map>>(userListEntry->second),
NoFactory::get());
360#ifdef HAVE_MUELU_INTREPID2
361 else if (name ==
"pcoarsen: element to node map")
364 level->Set(name, Teuchos::getValue<RCP<Kokkos::DynRankView<LocalOrdinal,typename Node::device_type> > >(userListEntry->second),
NoFactory::get());
367 else if (name ==
"output stream")
369 H.
SetMueLuOStream(Teuchos::getValue<RCP<Teuchos::FancyOStream> >(userListEntry->second));
373 size_t typeNameStart = name.find_first_not_of(
' ');
374 size_t typeNameEnd = name.find(
' ', typeNameStart);
375 std::string typeName = name.substr(typeNameStart, typeNameEnd - typeNameStart);
376 size_t varNameStart = name.find_first_not_of(
' ', typeNameEnd);
377 std::string varName = name.substr(varNameStart, name.size());
378 std::transform(typeName.begin(), typeName.end(), typeName.begin(), ::tolower);
380 if(typeName ==
"matrix")
381 level->Set(varName, Teuchos::getValue<RCP<Matrix> >(userListEntry->second),
NoFactory::get());
382 else if(typeName ==
"multivector")
383 level->Set(varName, Teuchos::getValue<RCP<MultiVector> >(userListEntry->second),
NoFactory::get());
384 else if(typeName ==
"vector")
385 level->Set(varName, Teuchos::getValue<RCP<Vector> >(userListEntry->second),
NoFactory::get());
386 else if(typeName ==
"map")
387 level->Set(varName, Teuchos::getValue<RCP<Xpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > >(userListEntry->second),
NoFactory::get());
388 else if(typeName ==
"ordinalvector")
389 level->Set(varName, Teuchos::getValue<RCP<Xpetra::Vector<LocalOrdinal, LocalOrdinal, GlobalOrdinal, Node> > >(userListEntry->second),
NoFactory::get());
390 else if(typeName ==
"scalar")
391 level->Set(varName, Teuchos::getValue<Scalar>(userListEntry->second),
NoFactory::get());
392 else if(typeName ==
"double")
393 level->Set(varName, Teuchos::getValue<double>(userListEntry->second),
NoFactory::get());
394 else if(typeName ==
"complex")
395 level->Set(varName, Teuchos::getValue<std::complex<double> >(userListEntry->second),
NoFactory::get());
396 else if(typeName ==
"int")
397 level->Set(varName, Teuchos::getValue<int>(userListEntry->second),
NoFactory::get());
398 else if(typeName ==
"string")
399 level->Set(varName, Teuchos::getValue<std::string>(userListEntry->second),
NoFactory::get());
400 else if(typeName ==
"array<go>")
401 level->Set(varName, Teuchos::getValue<Array<GlobalOrdinal> > (userListEntry->second),
NoFactory::get());
402 else if(typeName ==
"array<lo>")
403 level->Set(varName, Teuchos::getValue<Array<LocalOrdinal> >(userListEntry->second),
NoFactory::get());
404 else if(typeName ==
"arrayrcp<lo>")
405 level->Set(varName, Teuchos::getValue<ArrayRCP<LocalOrdinal> >(userListEntry->second),
NoFactory::get());
406 else if(typeName ==
"arrayrcp<go>")
407 level->Set(varName, Teuchos::getValue<ArrayRCP<GlobalOrdinal> >(userListEntry->second),
NoFactory::get());
409 throw std::runtime_error(
"Invalid non-serializable data on list");