46 #ifndef MUELU_NULLSPACEFACTORY_DEF_HPP
47 #define MUELU_NULLSPACEFACTORY_DEF_HPP
60 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
64 validParamList->
set< std::string >(
"Fine level nullspace",
"Nullspace",
"Variable name which is used to store null space multi vector on the finest level (default=\"Nullspace\"). For block matrices also \"Nullspace1\" to \"Nullspace9\" are accepted to describe the null space vectors for the (i,i) block (i=1..9).");
66 validParamList->
set<
RCP<const FactoryBase> >(
"A", Teuchos::null,
"Generating factory of the fine level matrix (only needed if default null space is generated)");
67 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace", Teuchos::null,
"Generating factory of the fine level null space");
71 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace1", Teuchos::null,
"Generating factory of the fine level null space associated with the (1,1) block in your n x n block matrix.");
72 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace2", Teuchos::null,
"Generating factory of the fine level null space associated with the (2,2) block in your n x n block matrix.");
73 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace3", Teuchos::null,
"Generating factory of the fine level null space associated with the (3,3) block in your n x n block matrix.");
74 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace4", Teuchos::null,
"Generating factory of the fine level null space associated with the (4,4) block in your n x n block matrix.");
75 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace5", Teuchos::null,
"Generating factory of the fine level null space associated with the (5,5) block in your n x n block matrix.");
76 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace6", Teuchos::null,
"Generating factory of the fine level null space associated with the (6,6) block in your n x n block matrix.");
77 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace7", Teuchos::null,
"Generating factory of the fine level null space associated with the (7,7) block in your n x n block matrix.");
78 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace8", Teuchos::null,
"Generating factory of the fine level null space associated with the (8,8) block in your n x n block matrix.");
79 validParamList->
set<
RCP<const FactoryBase> >(
"Nullspace9", Teuchos::null,
"Generating factory of the fine level null space associated with the (9,9) block in your n x n block matrix.");
81 return validParamList;
84 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
88 std::string nspName = pL.
get<std::string>(
"Fine level nullspace");
94 Input(currentLevel,
"A");
102 TEUCHOS_TEST_FOR_EXCEPTION(GetFactory(nspName)==Teuchos::null,
Exceptions::RuntimeError,
"MueLu::NullspaceFactory::DeclareInput(): You must declare an existing factory which produces the variable \"Nullspace\" in the NullspaceFactory (e.g. a TentativePFactory).");
103 currentLevel.
DeclareInput(
"Nullspace", GetFactory(nspName).
get(),
this);
107 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
115 std::string nspName = pL.
get<std::string>(
"Fine level nullspace");
123 GetOStream(
Runtime1) <<
"Use user-given nullspace " << nspName <<
": nullspace dimension=" << nullspace->getNumVectors() <<
" nullspace length=" << nullspace->getGlobalLength() << std::endl;
126 RCP<Matrix> A = Get< RCP<Matrix> >(currentLevel,
"A");
129 LocalOrdinal numPDEs = 1;
130 if(A->IsView(
"stridedMaps")==
true) {
133 numPDEs = Teuchos::rcp_dynamic_cast<
const StridedMap>(A->getRowMap())->getFixedBlockSize();
134 oldView = A->SwitchToView(oldView);
137 GetOStream(
Runtime1) <<
"Generating canonical nullspace: dimension = " << numPDEs << std::endl;
138 nullspace = MultiVectorFactory::Build(A->getDomainMap(), numPDEs);
141 if(bnsp.is_null() ==
true) {
142 for (
int i=0; i<numPDEs; ++i) {
144 int numBlocks = nsValues.
size() / numPDEs;
145 for (
int j=0; j< numBlocks; ++j) {
146 nsValues[j*numPDEs + i] = 1.0;
150 fillNullspaceVector(bnsp,numPDEs);
162 Set(currentLevel,
"Nullspace", nullspace);
166 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
170 for(
size_t r = 0; r < bmap->getNumMaps(); r++) {
174 for (
int i=0; i<numPDEs; ++i) {
176 int numBlocks = nsValues.
size() / numPDEs;
177 for (
int j=0; j< numBlocks; ++j) {
178 nsValues[j*numPDEs + i] = 1.0;
183 fillNullspaceVector(bpart,numPDEs);
190 #endif // MUELU_NULLSPACEFACTORY_DEF_HPP
Exception indicating invalid cast attempted.
T & Get(const std::string &ename, const FactoryBase *factory=NoFactory::get())
Get data without decrementing associated storage counter (i.e., read-only access). Usage: Level->Get< RCP<Matrix> >("A", factory) if factory == NULL => use default factory.
RCP< const ParameterList > GetValidParameterList() const
Define valid parameters for internal factory parameters.
T & get(const std::string &name, T def_value)
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
Timer to be used in factories. Similar to Monitor but with additional timers.
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
static const NoFactory * get()
void DeclareInput(Level ¤tLevel) const
Specifies the data that this class needs, and the factories that generate that data.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Class that holds all level-specific information.
int GetLevelID() const
Return level number.
Exception throws to report errors in the internal logical of the program.
Description of what is happening (more verbose)
void DeclareInput(const std::string &ename, const FactoryBase *factory, const FactoryBase *requestedBy=NoFactory::get())
Callback from FactoryBase::CallDeclareInput() and FactoryBase::DeclareInput()
void fillNullspaceVector(const RCP< BlockedMultiVector > &nsp, LocalOrdinal numPDEs) const
Helper function to recursively fill BlockedMultiVector with default null space vectors.
bool IsAvailable(const std::string &ename, const FactoryBase *factory=NoFactory::get()) const
Test whether a need's value has been saved.
void Build(Level ¤tLevel) const
Build an object with this factory.