46 #ifndef MUELU_FACTORYFACTORY_DECL_HPP 47 #define MUELU_FACTORYFACTORY_DECL_HPP 52 #include <Teuchos_ParameterEntry.hpp> 53 #include <Teuchos_Array.hpp> 61 #include "MueLu_FactoryManager.hpp" 70 #include "MueLu_AggregationExportFactory.hpp" 71 #include "MueLu_AmalgamationFactory.hpp" 72 #ifdef HAVE_MUELU_EXPERIMENTAL 73 #include "MueLu_BlockedCoarseMapFactory.hpp" 74 #include "MueLu_BlockedDirectSolver.hpp" 75 #include "MueLu_BlockedGaussSeidelSmoother.hpp" 76 #include "MueLu_BlockedPFactory.hpp" 77 #include "MueLu_BlockedRAPFactory.hpp" 78 #include "MueLu_BraessSarazinSmoother.hpp" 80 #include "MueLu_BrickAggregationFactory.hpp" 81 #include "MueLu_CloneRepartitionInterface.hpp" 82 #include "MueLu_CoalesceDropFactory.hpp" 83 #include "MueLu_CoarseMapFactory.hpp" 84 #include "MueLu_CoarseningVisualizationFactory.hpp" 85 #include "MueLu_ConstraintFactory.hpp" 86 #include "MueLu_CoupledAggregationFactory.hpp" 87 #include "MueLu_CoordinatesTransferFactory.hpp" 88 #include "MueLu_DirectSolver.hpp" 89 #include "MueLu_DropNegativeEntriesFactory.hpp" 90 #include "MueLu_EminPFactory.hpp" 91 #include "MueLu_FilteredAFactory.hpp" 92 #include "MueLu_GenericRFactory.hpp" 93 #ifdef HAVE_MUELU_EXPERIMENTAL 94 #include "MueLu_IndefBlockedDiagonalSmoother.hpp" 96 #include "MueLu_IsorropiaInterface.hpp" 97 #include "MueLu_LineDetectionFactory.hpp" 98 #include "MueLu_RepartitionInterface.hpp" 99 #include "MueLu_MapTransferFactory.hpp" 100 #include "MueLu_MatrixAnalysisFactory.hpp" 101 #include "MueLu_MultiVectorTransferFactory.hpp" 102 #include "MueLu_NullspaceFactory.hpp" 103 #include "MueLu_NullspacePresmoothFactory.hpp" 104 #include "MueLu_PatternFactory.hpp" 105 #include "MueLu_PgPFactory.hpp" 106 #ifdef HAVE_MUELU_EXPERIMENTAL 107 #include "MueLu_RebalanceBlockInterpolationFactory.hpp" 108 #include "MueLu_RebalanceBlockRestrictionFactory.hpp" 109 #include "MueLu_RebalanceBlockAcFactory.hpp" 111 #include "MueLu_RebalanceTransferFactory.hpp" 112 #include "MueLu_RepartitionFactory.hpp" 113 #include "MueLu_RepartitionHeuristicFactory.hpp" 114 #include "MueLu_RAPFactory.hpp" 115 #include "MueLu_RebalanceAcFactory.hpp" 116 #include "MueLu_ReorderBlockAFactory.hpp" 117 #include "MueLu_SaPFactory.hpp" 118 #include "MueLu_SegregatedAFactory.hpp" 119 #ifdef HAVE_MUELU_EXPERIMENTAL 120 #include "MueLu_SchurComplementFactory.hpp" 121 #include "MueLu_SimpleSmoother.hpp" 123 #include "MueLu_SmootherFactory.hpp" 124 #include "MueLu_SubBlockAFactory.hpp" 125 #ifdef HAVE_MUELU_EXPERIMENTAL 126 #ifdef HAVE_MUELU_TEKO 127 #include "MueLu_TekoSmoother.hpp" 130 #include "MueLu_TentativePFactory.hpp" 131 #include "MueLu_ToggleCoordinatesTransferFactory.hpp" 132 #include "MueLu_TogglePFactory.hpp" 133 #include "MueLu_TrilinosSmoother.hpp" 134 #include "MueLu_TransPFactory.hpp" 135 #include "MueLu_UncoupledAggregationFactory.hpp" 136 #include "MueLu_UserAggregationFactory.hpp" 137 #include "MueLu_UserPFactory.hpp" 138 #include "MueLu_SemiCoarsenPFactory.hpp" 139 #ifdef HAVE_MUELU_EXPERIMENTAL 140 #include "MueLu_UzawaSmoother.hpp" 142 #include "MueLu_ZoltanInterface.hpp" 143 #include "MueLu_Zoltan2Interface.hpp" 145 #ifdef HAVE_MUELU_MATLAB 147 #include "../matlab/src/MueLu_SingleLevelMatlabFactory_decl.hpp" 148 #include "../matlab/src/MueLu_SingleLevelMatlabFactory_def.hpp" 149 #include "../matlab/src/MueLu_TwoLevelMatlabFactory_decl.hpp" 150 #include "../matlab/src/MueLu_TwoLevelMatlabFactory_def.hpp" 151 #include "../matlab/src/MueLu_MatlabSmoother_decl.hpp" 152 #include "../matlab/src/MueLu_MatlabSmoother_def.hpp" 163 template <
class Scalar =
double,
class LocalOrdinal =
int,
class GlobalOrdinal = LocalOrdinal,
class Node = KokkosClassic::DefaultNode::DefaultNodeType>
165 #undef MUELU_FACTORYFACTORY_SHORT 168 typedef std::map<std::string, RCP<const FactoryBase> >
FactoryMap;
186 std::string factoryName;
187 Teuchos::ParameterList paramList;
188 if (!param.isList()) {
189 factoryName = Teuchos::getValue<std::string>(param);
191 paramList = Teuchos::getValue<Teuchos::ParameterList>(param);
192 factoryName = paramList.get<std::string>(
"factory");
196 if (factoryName ==
"AggregationExportFactory")
return Build2<AggregationExportFactory> (paramList, factoryMapIn, factoryManagersIn);
197 if (factoryName ==
"AmalgamationFactory")
return Build2<AmalgamationFactory> (paramList, factoryMapIn, factoryManagersIn);
198 #ifdef HAVE_MUELU_EXPERIMENTAL 199 if (factoryName ==
"BlockedCoarseMapFactory")
return Build2<BlockedCoarseMapFactory> (paramList, factoryMapIn, factoryManagersIn);
200 if (factoryName ==
"BlockedRAPFactory")
return BuildRAPFactory<BlockedRAPFactory> (paramList, factoryMapIn, factoryManagersIn);
202 if (factoryName ==
"BrickAggregationFactory")
return Build2<BrickAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
203 if (factoryName ==
"CloneRepartitionInterface")
return Build2<CloneRepartitionInterface> (paramList, factoryMapIn, factoryManagersIn);
204 if (factoryName ==
"CoarseMapFactory")
return Build2<CoarseMapFactory> (paramList, factoryMapIn, factoryManagersIn);
205 if (factoryName ==
"CoarseningVisualizationFactory")
return Build2<CoarseningVisualizationFactory>(paramList, factoryMapIn, factoryManagersIn);
206 if (factoryName ==
"CoalesceDropFactory")
return Build2<CoalesceDropFactory> (paramList, factoryMapIn, factoryManagersIn);
207 if (factoryName ==
"ConstraintFactory")
return Build2<ConstraintFactory> (paramList, factoryMapIn, factoryManagersIn);
209 if (factoryName ==
"CoordinatesTransferFactory")
return Build2<CoordinatesTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
210 if (factoryName ==
"DirectSolver")
return BuildDirectSolver (paramList, factoryMapIn, factoryManagersIn);
211 if (factoryName ==
"DropNegativeEntriesFactory")
return Build2<DropNegativeEntriesFactory> (paramList, factoryMapIn, factoryManagersIn);
212 if (factoryName ==
"EminPFactory")
return Build2<EminPFactory> (paramList, factoryMapIn, factoryManagersIn);
213 if (factoryName ==
"FilteredAFactory")
return Build2<FilteredAFactory> (paramList, factoryMapIn, factoryManagersIn);
214 if (factoryName ==
"GenericRFactory")
return Build2<GenericRFactory> (paramList, factoryMapIn, factoryManagersIn);
215 if (factoryName ==
"LineDetectionFactory")
return Build2<LineDetectionFactory> (paramList, factoryMapIn, factoryManagersIn);
216 if (factoryName ==
"MapTransferFactory")
return Build2<MapTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
217 if (factoryName ==
"MatrixAnalysisFactory")
return Build2<MatrixAnalysisFactory> (paramList, factoryMapIn, factoryManagersIn);
218 if (factoryName ==
"MultiVectorTransferFactory")
return Build2<MultiVectorTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
220 if (factoryName ==
"NoSmoother")
return rcp(
new SmootherFactory(Teuchos::null));
221 if (factoryName ==
"NullspaceFactory")
return Build2<NullspaceFactory> (paramList, factoryMapIn, factoryManagersIn);
222 if (factoryName ==
"NullspacePresmoothFactory")
return Build2<NullspacePresmoothFactory> (paramList, factoryMapIn, factoryManagersIn);
223 if (factoryName ==
"PatternFactory")
return Build2<PatternFactory> (paramList, factoryMapIn, factoryManagersIn);
224 if (factoryName ==
"PgPFactory")
return Build2<PgPFactory> (paramList, factoryMapIn, factoryManagersIn);
225 if (factoryName ==
"SaPFactory")
return Build2<SaPFactory> (paramList, factoryMapIn, factoryManagersIn);
226 if (factoryName ==
"RAPFactory")
return BuildRAPFactory<RAPFactory> (paramList, factoryMapIn, factoryManagersIn);
227 if (factoryName ==
"RebalanceAcFactory")
return Build2<RebalanceAcFactory> (paramList, factoryMapIn, factoryManagersIn);
228 if (factoryName ==
"RebalanceTransferFactory")
return Build2<RebalanceTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
229 if (factoryName ==
"ReorderBlockAFactory")
return Build2<ReorderBlockAFactory> (paramList, factoryMapIn, factoryManagersIn);
230 if (factoryName ==
"SegregatedAFactory")
return Build2<SegregatedAFactory> (paramList, factoryMapIn, factoryManagersIn);
231 if (factoryName ==
"SubBlockAFactory")
return Build2<SubBlockAFactory> (paramList, factoryMapIn, factoryManagersIn);
232 if (factoryName ==
"TentativePFactory")
return Build2<TentativePFactory> (paramList, factoryMapIn, factoryManagersIn);
234 if (factoryName ==
"TogglePFactory")
return BuildTogglePFactory<TogglePFactory> (paramList, factoryMapIn, factoryManagersIn);
235 if (factoryName ==
"TransPFactory")
return Build2<TransPFactory> (paramList, factoryMapIn, factoryManagersIn);
236 if (factoryName ==
"TrilinosSmoother")
return BuildTrilinosSmoother (paramList, factoryMapIn, factoryManagersIn);
237 if (factoryName ==
"UncoupledAggregationFactory")
return Build2<UncoupledAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
238 if (factoryName ==
"UserAggregationFactory")
return Build2<UserAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
239 if (factoryName ==
"UserPFactory")
return Build2<UserPFactory> (paramList, factoryMapIn, factoryManagersIn);
240 if (factoryName ==
"SemiCoarsenPFactory")
return Build2<SemiCoarsenPFactory> (paramList, factoryMapIn, factoryManagersIn);
241 if (factoryName ==
"RepartitionInterface")
return Build2<RepartitionInterface> (paramList, factoryMapIn, factoryManagersIn);
243 if (factoryName ==
"ZoltanInterface") {
244 #if defined(HAVE_MUELU_ZOLTAN) && defined(HAVE_MPI) 245 return Build2<ZoltanInterface>(paramList, factoryMapIn, factoryManagersIn);
247 TEUCHOS_TEST_FOR_EXCEPTION(
true,
Exceptions::RuntimeError,
"MueLu::FactoryFactory:BuildFactory(): Cannot create a ZoltanInterface object: Zoltan is disabled: HAVE_MUELU_ZOLTAN && HAVE_MPI == false.");
248 #endif // HAVE_MUELU_ZOLTAN && HAVE_MPI 250 if (factoryName ==
"Zoltan2Interface") {
251 #if defined(HAVE_MUELU_ZOLTAN2) && defined(HAVE_MPI) 252 return Build2<Zoltan2Interface>(paramList, factoryMapIn, factoryManagersIn);
254 TEUCHOS_TEST_FOR_EXCEPTION(
true,
Exceptions::RuntimeError,
"MueLu::FactoryFactory:BuildFactory(): Cannot create a Zoltan2Interface object: Zoltan2 is disabled: HAVE_MUELU_ZOLTAN2 && HAVE_MPI == false.");
255 #endif // HAVE_MUELU_ZOLTAN2 && HAVE_MPI 257 if (factoryName ==
"IsorropiaInterface") {
258 #if defined(HAVE_MUELU_ISORROPIA) && defined(HAVE_MPI) 259 return Build2<IsorropiaInterface>(paramList, factoryMapIn, factoryManagersIn);
261 TEUCHOS_TEST_FOR_EXCEPTION(
true,
Exceptions::RuntimeError,
"MueLu::FactoryFactory:BuildFactory(): Cannot create a IsorropiaInterface object: Isorropia is disabled: HAVE_MUELU_ISORROPIA && HAVE_MPI == false.");
262 #endif // HAVE_MUELU_ZOLTAN2 && HAVE_MPI 265 if (factoryName ==
"RepartitionFactory") {
267 return Build2<RepartitionFactory>(paramList, factoryMapIn, factoryManagersIn);
269 TEUCHOS_TEST_FOR_EXCEPTION(
true,
Exceptions::RuntimeError,
"MueLu::FactoryFactory:BuildFactory(): Cannot create a RepartitionFactory object: HAVE_MPI == false.");
272 if (factoryName ==
"RepartitionHeuristicFactory") {
274 return Build2<RepartitionHeuristicFactory>(paramList, factoryMapIn, factoryManagersIn);
276 TEUCHOS_TEST_FOR_EXCEPTION(
true,
Exceptions::RuntimeError,
"MueLu::FactoryFactory:BuildFactory(): Cannot create a RepartitionHeuristicFactory object: HAVE_MPI == false.");
280 #ifdef HAVE_MUELU_EXPERIMENTAL 281 if (factoryName ==
"BlockedDirectSolver")
return BuildBlockedDirectSolver(paramList, factoryMapIn, factoryManagersIn);
282 if (factoryName ==
"BlockedGaussSeidelSmoother")
return BuildBlockedSmoother<BlockedGaussSeidelSmoother>(paramList, factoryMapIn, factoryManagersIn);
283 if (factoryName ==
"BlockedPFactory")
return BuildBlockedFactory<BlockedPFactory>(paramList, factoryMapIn, factoryManagersIn);
284 if (factoryName ==
"BraessSarazinSmoother")
return BuildBlockedSmoother<BraessSarazinSmoother>(paramList, factoryMapIn, factoryManagersIn);
285 if (factoryName ==
"IndefiniteBlockDiagonalSmoother")
return BuildBlockedSmoother<IndefBlockedDiagonalSmoother>(paramList, factoryMapIn, factoryManagersIn);
286 if (factoryName ==
"SimpleSmoother")
return BuildBlockedSmoother<SimpleSmoother>(paramList, factoryMapIn, factoryManagersIn);
287 if (factoryName ==
"SchurComplementFactory")
return Build2<SchurComplementFactory> (paramList, factoryMapIn, factoryManagersIn);
288 if (factoryName ==
"RebalanceBlockRestrictionFactory")
return BuildBlockedFactory<RebalanceBlockRestrictionFactory>(paramList, factoryMapIn, factoryManagersIn);
289 if (factoryName ==
"RebalanceBlockAcFactory")
return BuildBlockedFactory<RebalanceBlockAcFactory>(paramList, factoryMapIn, factoryManagersIn);
290 if (factoryName ==
"RebalanceBlockInterpolationFactory")
return BuildBlockedFactory<RebalanceBlockInterpolationFactory>(paramList, factoryMapIn, factoryManagersIn);
291 #ifdef HAVE_MUELU_TEKO 292 if (factoryName ==
"TekoSmoother")
return BuildTekoSmoother(paramList, factoryMapIn, factoryManagersIn);
294 if (factoryName ==
"UzawaSmoother")
return BuildBlockedSmoother<UzawaSmoother>(paramList, factoryMapIn, factoryManagersIn);
298 #ifdef HAVE_MUELU_MATLAB 299 if (factoryName ==
"TwoLevelMatlabFactory")
return Build2<TwoLevelMatlabFactory> (paramList, factoryMapIn, factoryManagersIn);
300 if (factoryName ==
"SingleLevelMatlabFactory")
return Build2<SingleLevelMatlabFactory> (paramList, factoryMapIn, factoryManagersIn);
301 if (factoryName ==
"MatlabSmoother")
return BuildMatlabSmoother (paramList, factoryMapIn, factoryManagersIn);
305 if (factoryMapIn.find(factoryName) != factoryMapIn.end()) {
307 "MueLu::FactoryFactory: Error during the parsing of: " << std::endl << paramList << std::endl
308 <<
"'" << factoryName <<
"' is not a factory name but an existing instance of a factory." << std::endl
309 <<
"Extra parameters cannot be specified after the creation of the object." << std::endl << std::endl
310 <<
"Correct syntaxes includes:" << std::endl
311 <<
" <Parameter name=\"...\" type=\"string\" value=\"" << factoryName <<
"\"/>" << std::endl
313 <<
" <ParameterList name=\"...\"><Parameter name=\"factory\" type=\"string\" value=\"" << factoryName <<
"\"/></ParameterList>" << std::endl
316 return factoryMapIn.find(factoryName)->second;
319 TEUCHOS_TEST_FOR_EXCEPTION(
true,
Exceptions::RuntimeError,
"MueLu::FactoryFactory: unknown factory name : " << factoryName);
321 return Teuchos::null;
334 #define arraysize(ar) (sizeof(ar) / sizeof(ar[0])) 338 RCP<T> factory = rcp(
new T());
340 const char* strarray[] = {
"A",
"P",
"R",
"Graph",
"UnAmalgamationInfo",
"Aggregates",
"Nullspace",
"TransferFactory",
"DofsPerNode"};
341 std::vector<std::string> v(strarray, strarray +
arraysize(strarray));
342 for (
size_t i = 0; i < v.size(); ++i)
343 if (paramList.isParameter(v[i]))
344 factory->SetFactory(v[i],
BuildFactory(paramList.getEntry(v[i]), factoryMapIn, factoryManagersIn));
351 RCP<T> factory = rcp(
new T());
353 ParameterList paramListWithFactories;
356 RCP<const ParameterList> validParamList = factory->GetValidParameterList();
357 for (ParameterList::ConstIterator param = validParamList->begin(); param != validParamList->end(); ++param) {
358 const std::string& pName = validParamList->name(param);
360 if (!paramList.isParameter(pName)) {
365 if (validParamList->isType< RCP<const FactoryBase> >(pName)) {
367 RCP<const FactoryBase> generatingFact =
BuildFactory(paramList.getEntry(pName), factoryMapIn, factoryManagersIn);
368 paramListWithFactories.set(pName, generatingFact);
369 }
else if (validParamList->isType<RCP<const ParameterList> >(pName)) {
370 if (pName ==
"ParameterList") {
375 RCP<const ParameterList> subList = Teuchos::sublist(rcp(
new ParameterList(paramList)), pName);
376 paramListWithFactories.set(pName, subList);
379 paramListWithFactories.setEntry(pName, paramList.getEntry(pName));
384 factory->SetParameterList(paramListWithFactories);
392 if (paramList.isSublist(
"TransferFactories") ==
false) {
393 factory = Build2<T>(paramList, factoryMapIn, factoryManagersIn);
396 RCP<Teuchos::ParameterList> paramListNonConst = rcp(
new Teuchos::ParameterList(paramList));
397 RCP<const Teuchos::ParameterList> transferFactories = rcp(
new Teuchos::ParameterList(*sublist(paramListNonConst,
"TransferFactories")));
399 paramListNonConst->remove(
"TransferFactories");
401 factory = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
403 for (Teuchos::ParameterList::ConstIterator param = transferFactories->begin(); param != transferFactories->end(); ++param) {
404 RCP<const FactoryBase> p =
BuildFactory(transferFactories->entry(param), factoryMapIn, factoryManagersIn);
405 factory->AddTransferFactory(p);
415 if (paramList.isSublist(
"TransferFactories") ==
false) {
417 factory = Build2<T>(paramList, factoryMapIn, factoryManagersIn);
420 RCP<Teuchos::ParameterList> paramListNonConst = rcp(
new Teuchos::ParameterList(paramList));
421 RCP<const Teuchos::ParameterList> transferFactories = rcp(
new Teuchos::ParameterList(*sublist(paramListNonConst,
"TransferFactories")));
423 paramListNonConst->remove(
"TransferFactories");
426 factory = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
430 int numProlongatorFactories = 0;
431 int numPtentFactories = 0;
432 int numCoarseNspFactories = 0;
433 for (Teuchos::ParameterList::ConstIterator param = transferFactories->begin(); param != transferFactories->end(); ++param) {
434 size_t foundNsp = transferFactories->name(param).find(
"Nullspace");
435 if (foundNsp != std::string::npos && foundNsp == 0 && transferFactories->name(param).length()==10) {
436 numCoarseNspFactories++;
439 size_t foundPtent = transferFactories->name(param).find(
"Ptent");
440 if (foundPtent != std::string::npos && foundPtent == 0 && transferFactories->name(param).length()==6) {
444 size_t foundP = transferFactories->name(param).find(
"P");
445 if (foundP != std::string::npos && foundP == 0 && transferFactories->name(param).length()==2) {
446 numProlongatorFactories++;
450 TEUCHOS_TEST_FOR_EXCEPTION(numProlongatorFactories!=numCoarseNspFactories,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleP: The user has to provide the same number of prolongator and coarse nullspace factories!");
451 TEUCHOS_TEST_FOR_EXCEPTION(numPtentFactories!=numCoarseNspFactories,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleP: The user has to provide the same number of ptent and coarse nullspace factories!");
452 TEUCHOS_TEST_FOR_EXCEPTION(numProlongatorFactories < 2,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleP: The TogglePFactory needs at least two different prolongation operators. The factories have to be provided using the names P%i and Nullspace %i, where %i denotes a number between 1 and 9.");
455 std::vector<Teuchos::ParameterEntry> prolongatorFactoryNames(numProlongatorFactories);
456 std::vector<Teuchos::ParameterEntry> coarseNspFactoryNames(numProlongatorFactories);
457 std::vector<Teuchos::ParameterEntry> ptentFactoryNames(numProlongatorFactories);
459 for (Teuchos::ParameterList::ConstIterator param = transferFactories->begin(); param != transferFactories->end(); ++param) {
460 size_t foundNsp = transferFactories->name(param).find(
"Nullspace");
461 if (foundNsp != std::string::npos && foundNsp == 0 && transferFactories->name(param).length()==10) {
462 int number = atoi(&(transferFactories->name(param).at(9)));
463 TEUCHOS_TEST_FOR_EXCEPTION(number < 1 || number > numProlongatorFactories,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleP: Please use the format Nullspace%i with %i an integer between 1 and the maximum number of prolongation operators in TogglePFactory!");
464 coarseNspFactoryNames[number-1] = transferFactories->entry(param);
467 size_t foundPtent = transferFactories->name(param).find(
"Ptent");
468 if (foundPtent != std::string::npos && foundPtent == 0 && transferFactories->name(param).length()==6) {
469 int number = atoi(&(transferFactories->name(param).at(5)));
470 TEUCHOS_TEST_FOR_EXCEPTION(number < 1 || number > numPtentFactories,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleP: Please use the format Ptent%i with %i an integer between 1 and the maximum number of prolongation operators in TogglePFactory!");
471 ptentFactoryNames[number-1] = transferFactories->entry(param);
474 size_t foundP = transferFactories->name(param).find(
"P");
475 if (foundP != std::string::npos && foundP == 0 && transferFactories->name(param).length()==2) {
476 int number = atoi(&(transferFactories->name(param).at(1)));
477 TEUCHOS_TEST_FOR_EXCEPTION(number < 1 || number > numProlongatorFactories,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleP: Please use the format P%i with %i an integer between 1 and the maximum number of prolongation operators in TogglePFactory!");
478 prolongatorFactoryNames[number-1] = transferFactories->entry(param);
484 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = prolongatorFactoryNames.begin(); it != prolongatorFactoryNames.end(); ++it) {
485 RCP<const FactoryBase> p =
BuildFactory(*it, factoryMapIn, factoryManagersIn);
486 factory->AddProlongatorFactory(p);
490 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = ptentFactoryNames.begin(); it != ptentFactoryNames.end(); ++it) {
491 RCP<const FactoryBase> p =
BuildFactory(*it, factoryMapIn, factoryManagersIn);
492 factory->AddPtentFactory(p);
496 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = coarseNspFactoryNames.begin(); it != coarseNspFactoryNames.end(); ++it) {
497 RCP<const FactoryBase> p =
BuildFactory(*it, factoryMapIn, factoryManagersIn);
498 factory->AddCoarseNullspaceFactory(p);
505 RCP<ToggleCoordinatesTransferFactory> factory;
506 TEUCHOS_TEST_FOR_EXCEPTION(paramList.isSublist(
"TransferFactories") ==
false,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleCoordinatesTransferFactory: the ToggleCoordinatesTransferFactory needs a sublist 'TransferFactories' containing information about the subfactories for coordinate transfer!");
508 RCP<Teuchos::ParameterList> paramListNonConst = rcp(
new Teuchos::ParameterList(paramList));
509 RCP<const Teuchos::ParameterList> transferFactories = rcp(
new Teuchos::ParameterList(*sublist(paramListNonConst,
"TransferFactories")));
510 paramListNonConst->remove(
"TransferFactories");
513 factory = Build2<ToggleCoordinatesTransferFactory>(*paramListNonConst, factoryMapIn, factoryManagersIn);
517 int numCoordTransferFactories = 0;
518 for (Teuchos::ParameterList::ConstIterator param = transferFactories->begin(); param != transferFactories->end(); ++param) {
519 size_t foundCoordinates = transferFactories->name(param).find(
"Coordinates");
520 if (foundCoordinates != std::string::npos && foundCoordinates == 0 && transferFactories->name(param).length()==12) {
521 numCoordTransferFactories++;
525 TEUCHOS_TEST_FOR_EXCEPTION(numCoordTransferFactories != 2,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleCoordinatesTransfer: The ToggleCoordinatesTransferFactory needs two (different) coordinate transfer factories. The factories have to be provided using the names Coordinates%i, where %i denotes a number between 1 and 9.");
528 std::vector<Teuchos::ParameterEntry> coarseCoordsFactoryNames(numCoordTransferFactories);
530 for (Teuchos::ParameterList::ConstIterator param = transferFactories->begin(); param != transferFactories->end(); ++param) {
531 size_t foundCoords = transferFactories->name(param).find(
"Coordinates");
532 if (foundCoords != std::string::npos && foundCoords == 0 && transferFactories->name(param).length()==12) {
533 int number = atoi(&(transferFactories->name(param).at(11)));
534 TEUCHOS_TEST_FOR_EXCEPTION(number < 1 || number > numCoordTransferFactories,
Exceptions::RuntimeError,
"FactoryFactory::BuildToggleCoordinatesTransfer: Please use the format Coordinates%i with %i an integer between 1 and the maximum number of coordinate transfer factories in ToggleCoordinatesTransferFactory!");
535 coarseCoordsFactoryNames[number-1] = transferFactories->entry(param);
541 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = coarseCoordsFactoryNames.begin(); it != coarseCoordsFactoryNames.end(); ++it) {
542 RCP<const FactoryBase> p =
BuildFactory(*it, factoryMapIn, factoryManagersIn);
543 factory->AddCoordTransferFactory(p);
551 RCP<CoupledAggregationFactory> factory = Build<CoupledAggregationFactory>(paramList, factoryMapIn, factoryManagersIn);
553 if (paramList.isParameter(
"aggregation: ordering"))
554 factory->SetOrdering(paramList.get<std::string>(
"aggregation: ordering"));
556 if (paramList.isParameter(
"aggregation: max selected neighbors"))
557 factory->SetMaxNeighAlreadySelected(paramList.get<
int>(
"aggregation: max selected neighbors"));
559 if (paramList.isParameter(
"Phase3AggCreation"))
560 factory->SetPhase3AggCreation(paramList.get<
double>(
"Phase3AggCreation"));
562 if(paramList.isParameter(
"aggregation: min agg size"))
563 factory->SetMinNodesPerAggregate(paramList.get<
int>(
"aggregation: min agg size"));
579 if (paramList.begin() == paramList.end())
582 TEUCHOS_TEST_FOR_EXCEPTION(paramList.get<std::string>(
"factory") !=
"TrilinosSmoother",
Exceptions::RuntimeError,
"");
587 std::string type=
"";
if(paramList.isParameter(
"type")) type = paramList.get<std::string>(
"type");
588 int overlap=0;
if(paramList.isParameter(
"overlap")) overlap = paramList.get<
int> (
"overlap");
590 Teuchos::ParameterList params;
if(paramList.isParameter(
"ParameterList")) params = paramList.get<Teuchos::ParameterList>(
"ParameterList");
599 Teuchos::RCP<TrilinosSmoother> trilSmoo = Teuchos::rcp(
new TrilinosSmoother(type, params, overlap));
601 if (paramList.isParameter(
"LineDetection_Layers")) {
602 RCP<const FactoryBase> generatingFact =
BuildFactory(paramList.getEntry(
"LineDetection_Layers"), factoryMapIn, factoryManagersIn);
603 trilSmoo->SetFactory(
"LineDetection_Layers", generatingFact);
605 if (paramList.isParameter(
"LineDetection_VertLineIds")) {
606 RCP<const FactoryBase> generatingFact =
BuildFactory(paramList.getEntry(
"LineDetection_Layers"), factoryMapIn, factoryManagersIn);
607 trilSmoo->SetFactory(
"LineDetection_Layers", generatingFact);
609 if (paramList.isParameter(
"CoarseNumZLayers")) {
610 RCP<const FactoryBase> generatingFact =
BuildFactory(paramList.getEntry(
"CoarseNumZLayers"), factoryMapIn, factoryManagersIn);
611 trilSmoo->SetFactory(
"CoarseNumZLayers", generatingFact);
614 RCP<SmootherFactory> smooFact = rcp(
new SmootherFactory(Teuchos::null));
615 Teuchos::ParameterList smooFactParams;
617 smooFact->SetParameterList(smooFactParams);
618 smooFact->SetSmootherPrototypes(trilSmoo);
622 #ifdef HAVE_MUELU_MATLAB 635 if (paramList.begin() == paramList.end())
638 TEUCHOS_TEST_FOR_EXCEPTION(paramList.get<std::string>(
"factory") !=
"MatlabSmoother",
Exceptions::RuntimeError,
"");
644 Teuchos::RCP<MatlabSmoother> matSmoo = Teuchos::rcp(
new MatlabSmoother(paramList));
651 if (paramList.begin() == paramList.end())
656 std::string type;
if(paramList.isParameter(
"type")) type = paramList.get<std::string>(
"type");
658 Teuchos::ParameterList params;
if(paramList.isParameter(
"ParameterList")) params = paramList.get<Teuchos::ParameterList>(
"ParameterList");
663 #ifdef HAVE_MUELU_EXPERIMENTAL 667 RCP<ParameterList> paramListNonConst = rcp(
new ParameterList(paramList));
670 std::vector<RCP<FactoryManager> > facManagers;
674 bool blockExists =
true;
675 while (blockExists ==
true) {
676 std::stringstream ss;
677 ss <<
"block" << blockid;
679 if(paramList.isSublist(ss.str()) ==
true) {
681 RCP<const ParameterList> b = rcp(
new ParameterList(*sublist(paramListNonConst, ss.str())));
683 RCP<FactoryManager> M = Teuchos::null;
685 if (b->isParameter(
"group")) {
687 std::string facManagerName = b->get< std::string >(
"group");
688 TEUCHOS_TEST_FOR_EXCEPTION(factoryManagersIn.count(facManagerName) != 1,
Exceptions::RuntimeError,
"Factory manager has not been found. Please check the spelling of the factory managers in your xml file.");
689 RCP<FactoryManagerBase> Mb = factoryManagersIn.find(facManagerName)->second;
691 TEUCHOS_TEST_FOR_EXCEPTION(M==Teuchos::null,
Exceptions::RuntimeError,
"Failed to cast FactoryManagerBase object to FactoryManager.");
695 for (ParameterList::ConstIterator param = b->begin(); param != b->end(); ++param) {
696 RCP<const FactoryBase> p =
BuildFactory(b->entry(param), factoryMapIn, factoryManagersIn);
697 M->SetFactory(b->name(param),p);
702 M->SetIgnoreUserData(
true);
703 facManagers.push_back(M);
704 paramListNonConst->remove(ss.str());
714 RCP<T> bs = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
722 for (
int i = 0; i<Teuchos::as<int>(facManagers.size()); i++) {
723 bs->AddFactoryManager(facManagers[i],i);
729 #ifdef HAVE_MUELU_TEKO 732 RCP<ParameterList> paramListNonConst = rcp(
new ParameterList(paramList));
733 RCP<ParameterList> tekoParams = rcp(
new ParameterList(paramListNonConst->sublist(
"Inverse Factory Library")));
734 paramListNonConst->remove(
"Inverse Factory Library");
737 RCP<TekoSmoother> bs = Build2<TekoSmoother>(*paramListNonConst, factoryMapIn, factoryManagersIn);
746 bs->SetTekoParameters(tekoParams);
770 RCP<T> pfac = Teuchos::null;
773 RCP<ParameterList> paramListNonConst = rcp(
new ParameterList(paramList));
776 std::vector<RCP<FactoryManager> > facManagers;
780 bool blockExists =
true;
781 while (blockExists ==
true) {
782 std::stringstream ss;
783 ss <<
"block" << blockid;
785 if(paramList.isSublist(ss.str()) ==
true) {
787 RCP<const ParameterList> b = rcp(
new ParameterList(*sublist(paramListNonConst, ss.str())));
789 RCP<FactoryManager> M = Teuchos::null;
791 if (b->isParameter(
"group")) {
793 std::string facManagerName = b->get< std::string >(
"group");
794 TEUCHOS_TEST_FOR_EXCEPTION(factoryManagersIn.count(facManagerName) != 1,
Exceptions::RuntimeError,
"Factory manager has not been found. Please check the spelling of the factory managers in your xml file.");
795 RCP<FactoryManagerBase> Mb = factoryManagersIn.find(facManagerName)->second;
797 TEUCHOS_TEST_FOR_EXCEPTION(M==Teuchos::null,
Exceptions::RuntimeError,
"Failed to cast FactoryManagerBase object to FactoryManager.");
801 for (ParameterList::ConstIterator param = b->begin(); param != b->end(); ++param) {
802 RCP<const FactoryBase> p =
BuildFactory(b->entry(param), factoryMapIn, factoryManagersIn);
803 M->SetFactory(b->name(param),p);
808 M->SetIgnoreUserData(
true);
809 facManagers.push_back(M);
810 paramListNonConst->remove(ss.str());
820 pfac = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
823 for(
size_t i = 0; i<facManagers.size(); i++) {
824 pfac->AddFactoryManager(facManagers[i]);
833 #define MUELU_FACTORYFACTORY_SHORT 834 #endif // MUELU_FACTORYFACTORY_DECL_HPP RCP< FactoryBase > BuildMatlabSmoother(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
MatlabSmoother.
Generic Smoother Factory for generating the smoothers of the MG hierarchy.
This class specifies the default factory that should generate some data on a Level if the data does n...
RCP< FactoryBase > BuildBlockedSmoother(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
RCP< T > Build(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
Factory that can generate other factories from.
Class that encapsulates external library smoothers.
virtual RCP< const FactoryBase > BuildFactory(const Teuchos::ParameterEntry ¶m, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
RCP< T > BuildBlockedFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
RCP< ToggleCoordinatesTransferFactory > BuildToggleCoordinatesTransferFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
RCP< T > Build2(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
Namespace for MueLu classes and methods.
Class that encapsulates direct solvers. Autoselection of AmesosSmoother or Amesos2Smoother according ...
std::map< std::string, RCP< FactoryManagerBase > > FactoryManagerMap
std::map< std::string, RCP< const FactoryBase > > FactoryMap
RCP< FactoryBase > BuildBlockedDirectSolver(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
direct solver for nxn blocked matrices
RCP< T > BuildRAPFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
Base class for MueLu classes.
Class that encapsulates Matlab smoothers.
RCP< FactoryBase > BuildTekoSmoother(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
RCP< FactoryBase > BuildTrilinosSmoother(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
TrilinosSmoother.
RCP< FactoryBase > BuildCoupledAggregationFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
CoupledAggregationFactory.
Exception throws to report errors in the internal logical of the program.
RCP< T > BuildTogglePFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
RCP< FactoryBase > BuildDirectSolver(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
static const RCP< const NoFactory > getRCP()
Static Get() functions.