46 #ifndef MUELU_FACTORYFACTORY_DECL_HPP 47 #define MUELU_FACTORYFACTORY_DECL_HPP 61 #include "MueLu_FactoryManager.hpp" 70 #include "MueLu_AggregationExportFactory.hpp" 71 #include "MueLu_AmalgamationFactory.hpp" 72 #include "MueLu_BlackBoxPFactory.hpp" 73 #include "MueLu_BlockedCoarseMapFactory.hpp" 74 #include "MueLu_BlockedCoordinatesTransferFactory.hpp" 75 #include "MueLu_BlockedDirectSolver.hpp" 76 #include "MueLu_BlockedGaussSeidelSmoother.hpp" 77 #include "MueLu_BlockedJacobiSmoother.hpp" 78 #include "MueLu_BlockedPFactory.hpp" 79 #include "MueLu_BlockedRAPFactory.hpp" 80 #include "MueLu_BraessSarazinSmoother.hpp" 81 #include "MueLu_BrickAggregationFactory.hpp" 82 #include "MueLu_CloneRepartitionInterface.hpp" 83 #include "MueLu_CoalesceDropFactory.hpp" 84 #include "MueLu_CoarseMapFactory.hpp" 85 #include "MueLu_CoarseningVisualizationFactory.hpp" 86 #include "MueLu_ConstraintFactory.hpp" 87 #include "MueLu_CoupledAggregationFactory.hpp" 88 #include "MueLu_CoordinatesTransferFactory.hpp" 89 #include "MueLu_DirectSolver.hpp" 90 #include "MueLu_DropNegativeEntriesFactory.hpp" 91 #include "MueLu_EminPFactory.hpp" 92 #include "MueLu_FilteredAFactory.hpp" 93 #include "MueLu_FineLevelInputDataFactory.hpp" 94 #include "MueLu_GeneralGeometricPFactory.hpp" 95 #include "MueLu_GenericRFactory.hpp" 96 #include "MueLu_GeometricInterpolationPFactory.hpp" 97 #include "MueLu_IndefBlockedDiagonalSmoother.hpp" 98 #include "MueLu_IsorropiaInterface.hpp" 99 #include "MueLu_LineDetectionFactory.hpp" 100 #include "MueLu_RepartitionInterface.hpp" 101 #include "MueLu_RepartitionBlockDiagonalFactory.hpp" 102 #include "MueLu_MapTransferFactory.hpp" 103 #include "MueLu_MatrixAnalysisFactory.hpp" 104 #include "MueLu_MultiVectorTransferFactory.hpp" 105 #include "MueLu_NullspaceFactory.hpp" 106 #include "MueLu_NullspacePresmoothFactory.hpp" 107 #include "MueLu_PatternFactory.hpp" 108 #include "MueLu_PgPFactory.hpp" 109 #include "MueLu_RebalanceBlockInterpolationFactory.hpp" 110 #include "MueLu_RebalanceBlockRestrictionFactory.hpp" 111 #include "MueLu_RebalanceBlockAcFactory.hpp" 112 #include "MueLu_RebalanceTransferFactory.hpp" 113 #include "MueLu_RepartitionFactory.hpp" 114 #include "MueLu_RepartitionHeuristicFactory.hpp" 115 #include "MueLu_RAPFactory.hpp" 116 #include "MueLu_RAPShiftFactory.hpp" 117 #include "MueLu_RebalanceAcFactory.hpp" 118 #include "MueLu_ReorderBlockAFactory.hpp" 119 #include "MueLu_SaPFactory.hpp" 120 #include "MueLu_SegregatedAFactory.hpp" 121 #include "MueLu_SemiCoarsenPFactory.hpp" 122 #include "MueLu_SchurComplementFactory.hpp" 123 #include "MueLu_SimpleSmoother.hpp" 124 #include "MueLu_SmootherFactory.hpp" 125 #include "MueLu_StructuredAggregationFactory.hpp" 126 #include "MueLu_StructuredLineDetectionFactory.hpp" 127 #include "MueLu_SubBlockAFactory.hpp" 128 #ifdef HAVE_MUELU_TEKO 129 #include "MueLu_TekoSmoother.hpp" 131 #include "MueLu_TentativePFactory.hpp" 132 #include "MueLu_ToggleCoordinatesTransferFactory.hpp" 133 #include "MueLu_TogglePFactory.hpp" 134 #include "MueLu_TrilinosSmoother.hpp" 135 #include "MueLu_TransPFactory.hpp" 136 #include "MueLu_UncoupledAggregationFactory.hpp" 137 #include "MueLu_HybridAggregationFactory.hpp" 138 #include "MueLu_UnsmooshFactory.hpp" 139 #include "MueLu_UserAggregationFactory.hpp" 140 #include "MueLu_UserPFactory.hpp" 141 #include "MueLu_UzawaSmoother.hpp" 142 #include "MueLu_VariableDofLaplacianFactory.hpp" 143 #include "MueLu_ZoltanInterface.hpp" 144 #include "MueLu_Zoltan2Interface.hpp" 146 #ifdef HAVE_MUELU_KOKKOS_REFACTOR 147 #include "MueLu_CoalesceDropFactory_kokkos.hpp" 148 #include "MueLu_CoarseMapFactory_kokkos.hpp" 149 #include "MueLu_CoordinatesTransferFactory_kokkos.hpp" 150 #include "MueLu_NullspaceFactory_kokkos.hpp" 151 #include "MueLu_SaPFactory_kokkos.hpp" 152 #include "MueLu_TentativePFactory_kokkos.hpp" 153 #include "MueLu_UncoupledAggregationFactory_kokkos.hpp" 156 #ifdef HAVE_MUELU_MATLAB 158 #include "../matlab/src/MueLu_SingleLevelMatlabFactory_decl.hpp" 159 #include "../matlab/src/MueLu_SingleLevelMatlabFactory_def.hpp" 160 #include "../matlab/src/MueLu_TwoLevelMatlabFactory_decl.hpp" 161 #include "../matlab/src/MueLu_TwoLevelMatlabFactory_def.hpp" 162 #include "../matlab/src/MueLu_MatlabSmoother_decl.hpp" 163 #include "../matlab/src/MueLu_MatlabSmoother_def.hpp" 166 #ifdef HAVE_MUELU_INTREPID2 167 #include "MueLu_IntrepidPCoarsenFactory.hpp" 178 template <
class Scalar =
double,
class LocalOrdinal =
int,
class GlobalOrdinal = LocalOrdinal,
class Node = KokkosClassic::DefaultNode::DefaultNodeType>
180 #undef MUELU_FACTORYFACTORY_SHORT 183 typedef std::map<std::string, RCP<const FactoryBase> >
FactoryMap;
207 std::string factoryName;
210 factoryName = Teuchos::getValue<std::string>(param);
212 paramList = Teuchos::getValue<Teuchos::ParameterList>(param);
213 factoryName = paramList.
get<std::string>(
"factory");
217 if (factoryName ==
"AggregationExportFactory")
return Build2<AggregationExportFactory> (paramList, factoryMapIn, factoryManagersIn);
218 if (factoryName ==
"AmalgamationFactory")
return Build2<AmalgamationFactory> (paramList, factoryMapIn, factoryManagersIn);
219 if (factoryName ==
"BlockedCoarseMapFactory")
return Build2<BlockedCoarseMapFactory> (paramList, factoryMapIn, factoryManagersIn);
220 if (factoryName ==
"BlockedCoordinatesTransferFactory")
return Build2<BlockedCoordinatesTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
221 if (factoryName ==
"BlockedRAPFactory")
return BuildRAPFactory<BlockedRAPFactory> (paramList, factoryMapIn, factoryManagersIn);
222 if (factoryName ==
"BrickAggregationFactory")
return Build2<BrickAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
223 if (factoryName ==
"CloneRepartitionInterface")
return Build2<CloneRepartitionInterface> (paramList, factoryMapIn, factoryManagersIn);
224 if (factoryName ==
"CoarseMapFactory")
return Build2<CoarseMapFactory> (paramList, factoryMapIn, factoryManagersIn);
225 if (factoryName ==
"CoarseningVisualizationFactory")
return Build2<CoarseningVisualizationFactory> (paramList, factoryMapIn, factoryManagersIn);
226 if (factoryName ==
"CoalesceDropFactory")
return Build2<CoalesceDropFactory> (paramList, factoryMapIn, factoryManagersIn);
227 if (factoryName ==
"ConstraintFactory")
return Build2<ConstraintFactory> (paramList, factoryMapIn, factoryManagersIn);
229 if (factoryName ==
"CoordinatesTransferFactory")
return Build2<CoordinatesTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
230 if (factoryName ==
"DirectSolver")
return BuildDirectSolver (paramList, factoryMapIn, factoryManagersIn);
231 if (factoryName ==
"DropNegativeEntriesFactory")
return Build2<DropNegativeEntriesFactory> (paramList, factoryMapIn, factoryManagersIn);
232 if (factoryName ==
"EminPFactory")
return Build2<EminPFactory> (paramList, factoryMapIn, factoryManagersIn);
233 if (factoryName ==
"FilteredAFactory")
return Build2<FilteredAFactory> (paramList, factoryMapIn, factoryManagersIn);
234 if (factoryName ==
"FineLevelInputDataFactory")
return Build2<FineLevelInputDataFactory> (paramList, factoryMapIn, factoryManagersIn);
235 if (factoryName ==
"GeneralGeometricPFactory")
return Build2<GeneralGeometricPFactory> (paramList, factoryMapIn, factoryManagersIn);
236 if (factoryName ==
"GenericRFactory")
return Build2<GenericRFactory> (paramList, factoryMapIn, factoryManagersIn);
237 if (factoryName ==
"GeometricInterpolationPFactory")
return Build2<GeometricInterpolationPFactory> (paramList, factoryMapIn, factoryManagersIn);
238 if (factoryName ==
"LineDetectionFactory")
return Build2<LineDetectionFactory> (paramList, factoryMapIn, factoryManagersIn);
239 if (factoryName ==
"MapTransferFactory")
return Build2<MapTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
240 if (factoryName ==
"MatrixAnalysisFactory")
return Build2<MatrixAnalysisFactory> (paramList, factoryMapIn, factoryManagersIn);
241 if (factoryName ==
"MultiVectorTransferFactory")
return Build2<MultiVectorTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
244 if (factoryName ==
"NullspaceFactory")
return Build2<NullspaceFactory> (paramList, factoryMapIn, factoryManagersIn);
245 if (factoryName ==
"NullspacePresmoothFactory")
return Build2<NullspacePresmoothFactory> (paramList, factoryMapIn, factoryManagersIn);
246 if (factoryName ==
"PatternFactory")
return Build2<PatternFactory> (paramList, factoryMapIn, factoryManagersIn);
247 if (factoryName ==
"PgPFactory")
return Build2<PgPFactory> (paramList, factoryMapIn, factoryManagersIn);
248 if (factoryName ==
"SaPFactory")
return Build2<SaPFactory> (paramList, factoryMapIn, factoryManagersIn);
249 if (factoryName ==
"RAPFactory")
return BuildRAPFactory<RAPFactory> (paramList, factoryMapIn, factoryManagersIn);
250 if (factoryName ==
"RAPShiftFactory")
return BuildRAPFactory<RAPShiftFactory> (paramList, factoryMapIn, factoryManagersIn);
251 if (factoryName ==
"RebalanceAcFactory")
return Build2<RebalanceAcFactory> (paramList, factoryMapIn, factoryManagersIn);
252 if (factoryName ==
"RebalanceTransferFactory")
return Build2<RebalanceTransferFactory> (paramList, factoryMapIn, factoryManagersIn);
253 if (factoryName ==
"ReorderBlockAFactory")
return Build2<ReorderBlockAFactory> (paramList, factoryMapIn, factoryManagersIn);
254 if (factoryName ==
"RepartitionInterface")
return Build2<RepartitionInterface> (paramList, factoryMapIn, factoryManagersIn);
255 if (factoryName ==
"SegregatedAFactory")
return Build2<SegregatedAFactory> (paramList, factoryMapIn, factoryManagersIn);
256 if (factoryName ==
"SemiCoarsenPFactory")
return Build2<SemiCoarsenPFactory> (paramList, factoryMapIn, factoryManagersIn);
257 if (factoryName ==
"StructuredAggregationFactory")
return Build2<StructuredAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
258 if (factoryName ==
"StructuredLineDetectionFactory")
return Build2<StructuredLineDetectionFactory> (paramList, factoryMapIn, factoryManagersIn);
259 if (factoryName ==
"SubBlockAFactory")
return Build2<SubBlockAFactory> (paramList, factoryMapIn, factoryManagersIn);
260 if (factoryName ==
"TentativePFactory")
return Build2<TentativePFactory> (paramList, factoryMapIn, factoryManagersIn);
262 if (factoryName ==
"TogglePFactory")
return BuildTogglePFactory<TogglePFactory> (paramList, factoryMapIn, factoryManagersIn);
263 if (factoryName ==
"TransPFactory")
return Build2<TransPFactory> (paramList, factoryMapIn, factoryManagersIn);
264 if (factoryName ==
"TrilinosSmoother")
return BuildTrilinosSmoother (paramList, factoryMapIn, factoryManagersIn);
265 if (factoryName ==
"UncoupledAggregationFactory")
return Build2<UncoupledAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
266 if (factoryName ==
"HybridAggregationFactory")
return Build2<HybridAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
267 if (factoryName ==
"UnsmooshFactory")
return Build2<UnsmooshFactory> (paramList, factoryMapIn, factoryManagersIn);
268 if (factoryName ==
"UserAggregationFactory")
return Build2<UserAggregationFactory> (paramList, factoryMapIn, factoryManagersIn);
269 if (factoryName ==
"UserPFactory")
return Build2<UserPFactory> (paramList, factoryMapIn, factoryManagersIn);
270 if (factoryName ==
"VariableDofLaplacianFactory")
return Build2<VariableDofLaplacianFactory> (paramList, factoryMapIn, factoryManagersIn);
271 #ifdef HAVE_MUELU_KOKKOS_REFACTOR 272 if (factoryName ==
"CoalesceDropFactory_kokkos")
return Build2<CoalesceDropFactory_kokkos> (paramList, factoryMapIn, factoryManagersIn);
273 if (factoryName ==
"CoarseMapFactory_kokkos")
return Build2<CoarseMapFactory_kokkos> (paramList, factoryMapIn, factoryManagersIn);
274 if (factoryName ==
"CoordinatesTransferFactory_kokkos")
return Build2<CoordinatesTransferFactory_kokkos> (paramList, factoryMapIn, factoryManagersIn);
275 if (factoryName ==
"NullspaceFactory_kokkos")
return Build2<NullspaceFactory_kokkos> (paramList, factoryMapIn, factoryManagersIn);
276 if (factoryName ==
"SaPFactory_kokkos")
return Build2<SaPFactory_kokkos> (paramList, factoryMapIn, factoryManagersIn);
277 if (factoryName ==
"TentativePFactory_kokkos")
return Build2<TentativePFactory_kokkos> (paramList, factoryMapIn, factoryManagersIn);
278 if (factoryName ==
"UncoupledAggregationFactory_kokkos")
return Build2<UncoupledAggregationFactory_kokkos> (paramList, factoryMapIn, factoryManagersIn);
281 if (factoryName ==
"ZoltanInterface") {
282 #if defined(HAVE_MUELU_ZOLTAN) && defined(HAVE_MPI) 283 return Build2<ZoltanInterface>(paramList, factoryMapIn, factoryManagersIn);
286 #endif // HAVE_MUELU_ZOLTAN && HAVE_MPI 288 if (factoryName ==
"Zoltan2Interface") {
289 #if defined(HAVE_MUELU_ZOLTAN2) && defined(HAVE_MPI) 290 return Build2<Zoltan2Interface>(paramList, factoryMapIn, factoryManagersIn);
293 #endif // HAVE_MUELU_ZOLTAN2 && HAVE_MPI 295 if (factoryName ==
"IsorropiaInterface") {
296 #if defined(HAVE_MUELU_ISORROPIA) && defined(HAVE_MPI) 297 return Build2<IsorropiaInterface>(paramList, factoryMapIn, factoryManagersIn);
300 #endif // HAVE_MUELU_ZOLTAN2 && HAVE_MPI 303 if (factoryName ==
"RepartitionFactory") {
305 return Build2<RepartitionFactory>(paramList, factoryMapIn, factoryManagersIn);
310 if (factoryName ==
"RepartitionHeuristicFactory") {
312 return Build2<RepartitionHeuristicFactory>(paramList, factoryMapIn, factoryManagersIn);
318 if (factoryName ==
"BlockedDirectSolver")
return BuildBlockedDirectSolver(paramList, factoryMapIn, factoryManagersIn);
319 if (factoryName ==
"BlockedGaussSeidelSmoother")
return BuildBlockedSmoother<BlockedGaussSeidelSmoother>(paramList, factoryMapIn, factoryManagersIn);
320 if (factoryName ==
"BlockedJacobiSmoother")
return BuildBlockedSmoother<BlockedJacobiSmoother>(paramList, factoryMapIn, factoryManagersIn);
321 if (factoryName ==
"BlockedPFactory")
return BuildBlockedFactory<BlockedPFactory>(paramList, factoryMapIn, factoryManagersIn);
322 if (factoryName ==
"BraessSarazinSmoother")
return BuildBlockedSmoother<BraessSarazinSmoother>(paramList, factoryMapIn, factoryManagersIn);
323 if (factoryName ==
"IndefiniteBlockDiagonalSmoother")
return BuildBlockedSmoother<IndefBlockedDiagonalSmoother>(paramList, factoryMapIn, factoryManagersIn);
324 if (factoryName ==
"SimpleSmoother")
return BuildBlockedSmoother<SimpleSmoother>(paramList, factoryMapIn, factoryManagersIn);
325 if (factoryName ==
"SchurComplementFactory")
return Build2<SchurComplementFactory> (paramList, factoryMapIn, factoryManagersIn);
326 if (factoryName ==
"RebalanceBlockRestrictionFactory")
return BuildBlockedFactory<RebalanceBlockRestrictionFactory>(paramList, factoryMapIn, factoryManagersIn);
327 if (factoryName ==
"RebalanceBlockAcFactory")
return BuildBlockedFactory<RebalanceBlockAcFactory>(paramList, factoryMapIn, factoryManagersIn);
328 if (factoryName ==
"RebalanceBlockInterpolationFactory")
return BuildBlockedFactory<RebalanceBlockInterpolationFactory>(paramList, factoryMapIn, factoryManagersIn);
330 if (factoryName ==
"RepartitionBlockDiagonalFactory")
return Build2<RepartitionBlockDiagonalFactory> (paramList, factoryMapIn, factoryManagersIn);
332 #ifdef HAVE_MUELU_TEKO 333 if (factoryName ==
"TekoSmoother")
return BuildTekoSmoother(paramList, factoryMapIn, factoryManagersIn);
335 if (factoryName ==
"UzawaSmoother")
return BuildBlockedSmoother<UzawaSmoother>(paramList, factoryMapIn, factoryManagersIn);
338 #ifdef HAVE_MUELU_MATLAB 339 if (factoryName ==
"TwoLevelMatlabFactory")
return Build2<TwoLevelMatlabFactory> (paramList, factoryMapIn, factoryManagersIn);
340 if (factoryName ==
"SingleLevelMatlabFactory")
return Build2<SingleLevelMatlabFactory> (paramList, factoryMapIn, factoryManagersIn);
341 if (factoryName ==
"MatlabSmoother")
return BuildMatlabSmoother (paramList, factoryMapIn, factoryManagersIn);
344 #ifdef HAVE_MUELU_INTREPID2 345 if (factoryName ==
"IntrepidPCoarsenFactory")
return Build2<IntrepidPCoarsenFactory> (paramList, factoryMapIn, factoryManagersIn);
349 if (factoryMapIn.find(factoryName) != factoryMapIn.end()) {
351 "MueLu::FactoryFactory: Error during the parsing of: " << std::endl << paramList << std::endl
352 <<
"'" << factoryName <<
"' is not a factory name but an existing instance of a factory." << std::endl
353 <<
"Extra parameters cannot be specified after the creation of the object." << std::endl << std::endl
354 <<
"Correct syntaxes includes:" << std::endl
355 <<
" <Parameter name=\"...\" type=\"string\" value=\"" << factoryName <<
"\"/>" << std::endl
357 <<
" <ParameterList name=\"...\"><Parameter name=\"factory\" type=\"string\" value=\"" << factoryName <<
"\"/></ParameterList>" << std::endl
360 return factoryMapIn.find(factoryName)->second;
378 #define arraysize(ar) (sizeof(ar) / sizeof(ar[0])) 384 const char* strarray[] = {
"A",
"P",
"R",
"Graph",
"UnAmalgamationInfo",
"Aggregates",
"Nullspace",
"TransferFactory",
"DofsPerNode"};
385 std::vector<std::string> v(strarray, strarray +
arraysize(strarray));
386 for (
size_t i = 0; i < v.size(); ++i)
388 factory->SetFactory(v[i],
BuildFactory(paramList.
getEntry(v[i]), factoryMapIn, factoryManagersIn));
403 const std::string& pName = validParamList->name(param);
413 paramListWithFactories.
set(pName, generatingFact);
415 if (pName ==
"ParameterList") {
421 paramListWithFactories.
set(pName, subList);
429 factory->SetParameterList(paramListWithFactories);
437 if (paramList.
isSublist(
"TransferFactories") ==
false) {
438 factory = Build2<T>(paramList, factoryMapIn, factoryManagersIn);
444 paramListNonConst->
remove(
"TransferFactories");
446 factory = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
450 factory->AddTransferFactory(p);
460 if (paramList.
isSublist(
"TransferFactories") ==
false) {
462 factory = Build2<T>(paramList, factoryMapIn, factoryManagersIn);
468 paramListNonConst->
remove(
"TransferFactories");
471 factory = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
475 int numProlongatorFactories = 0;
476 int numPtentFactories = 0;
477 int numCoarseNspFactories = 0;
479 size_t foundNsp = transferFactories->
name(param).find(
"Nullspace");
480 if (foundNsp != std::string::npos && foundNsp == 0 && transferFactories->
name(param).length()==10) {
481 numCoarseNspFactories++;
484 size_t foundPtent = transferFactories->
name(param).find(
"Ptent");
485 if (foundPtent != std::string::npos && foundPtent == 0 && transferFactories->
name(param).length()==6) {
489 size_t foundP = transferFactories->
name(param).find(
"P");
490 if (foundP != std::string::npos && foundP == 0 && transferFactories->
name(param).length()==2) {
491 numProlongatorFactories++;
497 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.");
500 std::vector<Teuchos::ParameterEntry> prolongatorFactoryNames(numProlongatorFactories);
501 std::vector<Teuchos::ParameterEntry> coarseNspFactoryNames(numProlongatorFactories);
502 std::vector<Teuchos::ParameterEntry> ptentFactoryNames(numProlongatorFactories);
505 size_t foundNsp = transferFactories->
name(param).find(
"Nullspace");
506 if (foundNsp != std::string::npos && foundNsp == 0 && transferFactories->
name(param).length()==10) {
507 int number = atoi(&(transferFactories->
name(param).at(9)));
508 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!");
509 coarseNspFactoryNames[number-1] = transferFactories->
entry(param);
512 size_t foundPtent = transferFactories->
name(param).find(
"Ptent");
513 if (foundPtent != std::string::npos && foundPtent == 0 && transferFactories->
name(param).length()==6) {
514 int number = atoi(&(transferFactories->
name(param).at(5)));
515 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!");
516 ptentFactoryNames[number-1] = transferFactories->
entry(param);
519 size_t foundP = transferFactories->
name(param).find(
"P");
520 if (foundP != std::string::npos && foundP == 0 && transferFactories->
name(param).length()==2) {
521 int number = atoi(&(transferFactories->
name(param).at(1)));
522 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!");
523 prolongatorFactoryNames[number-1] = transferFactories->
entry(param);
529 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = prolongatorFactoryNames.begin(); it != prolongatorFactoryNames.end(); ++it) {
531 factory->AddProlongatorFactory(p);
535 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = ptentFactoryNames.begin(); it != ptentFactoryNames.end(); ++it) {
537 factory->AddPtentFactory(p);
541 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = coarseNspFactoryNames.begin(); it != coarseNspFactoryNames.end(); ++it) {
543 factory->AddCoarseNullspaceFactory(p);
551 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!");
555 paramListNonConst->remove(
"TransferFactories");
558 factory = Build2<ToggleCoordinatesTransferFactory>(*paramListNonConst, factoryMapIn, factoryManagersIn);
562 int numCoordTransferFactories = 0;
564 size_t foundCoordinates = transferFactories->
name(param).find(
"Coordinates");
565 if (foundCoordinates != std::string::npos && foundCoordinates == 0 && transferFactories->
name(param).length()==12) {
566 numCoordTransferFactories++;
570 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.");
573 std::vector<Teuchos::ParameterEntry> coarseCoordsFactoryNames(numCoordTransferFactories);
576 size_t foundCoords = transferFactories->
name(param).find(
"Coordinates");
577 if (foundCoords != std::string::npos && foundCoords == 0 && transferFactories->
name(param).length()==12) {
578 int number = atoi(&(transferFactories->
name(param).at(11)));
579 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!");
580 coarseCoordsFactoryNames[number-1] = transferFactories->
entry(param);
586 for (std::vector<Teuchos::ParameterEntry>::const_iterator it = coarseCoordsFactoryNames.begin(); it != coarseCoordsFactoryNames.end(); ++it) {
588 factory->AddCoordTransferFactory(p);
598 if (paramList.isParameter(
"aggregation: ordering"))
599 factory->SetOrdering(paramList.get<std::string>(
"aggregation: ordering"));
601 if (paramList.isParameter(
"aggregation: max selected neighbors"))
602 factory->SetMaxNeighAlreadySelected(paramList.get<
int>(
"aggregation: max selected neighbors"));
604 if (paramList.isParameter(
"Phase3AggCreation"))
605 factory->SetPhase3AggCreation(paramList.get<
double>(
"Phase3AggCreation"));
607 if(paramList.isParameter(
"aggregation: min agg size"))
608 factory->SetMinNodesPerAggregate(paramList.get<
int>(
"aggregation: min agg size"));
624 if (paramList.
begin() == paramList.
end())
632 std::string type=
"";
if(paramList.
isParameter(
"type")) type = paramList.
get<std::string>(
"type");
633 int overlap=0;
if(paramList.
isParameter(
"overlap")) overlap = paramList.
get<
int> (
"overlap");
646 if (paramList.
isParameter(
"LineDetection_Layers")) {
648 trilSmoo->SetFactory(
"LineDetection_Layers", generatingFact);
650 if (paramList.
isParameter(
"LineDetection_VertLineIds")) {
652 trilSmoo->SetFactory(
"LineDetection_Layers", generatingFact);
656 trilSmoo->SetFactory(
"CoarseNumZLayers", generatingFact);
662 smooFact->SetParameterList(smooFactParams);
663 smooFact->SetSmootherPrototypes(trilSmoo);
667 #ifdef HAVE_MUELU_MATLAB 680 if (paramList.
begin() == paramList.
end())
696 if (paramList.
begin() == paramList.
end())
701 std::string type;
if(paramList.
isParameter(
"type")) type = paramList.
get<std::string>(
"type");
714 std::vector<RCP<FactoryManager> > facManagers;
718 bool blockExists =
true;
719 while (blockExists ==
true) {
720 std::stringstream ss;
721 ss <<
"block" << blockid;
723 if(paramList.
isSublist(ss.str()) ==
true) {
729 if (b->isParameter(
"group")) {
731 std::string facManagerName = b->
get< std::string >(
"group");
741 M->SetFactory(b->name(param),p);
746 M->SetIgnoreUserData(
true);
747 facManagers.push_back(M);
748 paramListNonConst->remove(ss.str());
758 RCP<T> bs = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
766 for (
int i = 0; i<Teuchos::as<int>(facManagers.size()); i++) {
767 bs->AddFactoryManager(facManagers[i],i);
773 #ifdef HAVE_MUELU_TEKO 778 paramListNonConst->remove(
"Inverse Factory Library");
781 RCP<TekoSmoother> bs = Build2<TekoSmoother>(*paramListNonConst, factoryMapIn, factoryManagersIn);
790 bs->SetTekoParameters(tekoParams);
814 RCP<T> pfac = Teuchos::null;
820 std::vector<RCP<FactoryManager> > facManagers;
824 bool blockExists =
true;
825 while (blockExists ==
true) {
826 std::stringstream ss;
827 ss <<
"block" << blockid;
829 if(paramList.
isSublist(ss.str()) ==
true) {
835 if (b->isParameter(
"group")) {
837 std::string facManagerName = b->
get< std::string >(
"group");
847 M->SetFactory(b->name(param),p);
852 M->SetIgnoreUserData(
true);
853 facManagers.push_back(M);
854 paramListNonConst->remove(ss.str());
864 pfac = Build2<T>(*paramListNonConst, factoryMapIn, factoryManagersIn);
867 for(
size_t i = 0; i<facManagers.size(); i++) {
868 pfac->AddFactoryManager(facManagers[i]);
876 #define MUELU_FACTORYFACTORY_SHORT 877 #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.
ConstIterator begin() const
T & get(const std::string &name, T def_value)
Class that encapsulates external library smoothers.
virtual RCP< const FactoryBase > BuildFactory(const Teuchos::ParameterEntry ¶m, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
: Interpret Factory parameter list and build new factory
RCP< T > BuildBlockedFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
const ParameterEntry & entry(ConstIterator i) const
RCP< ToggleCoordinatesTransferFactory > BuildToggleCoordinatesTransferFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
bool isSublist(const std::string &name) 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 ...
ParameterList & setEntry(const std::string &name, const ParameterEntry &entry)
ConstIterator end() const
std::map< std::string, RCP< FactoryManagerBase > > FactoryManagerMap
bool remove(std::string const &name, bool throwIfNotExists=true)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
std::map< std::string, RCP< const FactoryBase > > FactoryMap
params_t::ConstIterator ConstIterator
RCP< FactoryBase > BuildBlockedDirectSolver(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
direct solver for nxn blocked matrices
const std::string & name() const
RCP< T > BuildRAPFactory(const Teuchos::ParameterList ¶mList, const FactoryMap &factoryMapIn, const FactoryManagerMap &factoryManagersIn) const
Base class for MueLu classes.
Class that encapsulates Matlab smoothers.
#define TEUCHOS_UNREACHABLE_RETURN(dummyReturnVal)
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.
bool isParameter(const std::string &name) const
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.
ParameterEntry & getEntry(const std::string &name)
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.