9 #ifndef MUELU_MATLABSMOOTHER_DEF_HPP
10 #define MUELU_MATLABSMOOTHER_DEF_HPP
14 #if defined(HAVE_MUELU_MATLAB)
19 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
21 SetParameterList(paramList);
24 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
28 setupFunction_ = pL.
get(
"Setup Function",
"");
29 solveFunction_ = pL.
get(
"Solve Function",
"");
30 solveDataSize_ = pL.
get(
"Number of Solver Args", 0);
33 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
36 this->Input(currentLevel,
"A");
38 needsSetup_ = pL.
get<
string>(
"Needs");
40 for (
size_t i = 0; i < needsList.size(); i++) {
42 this->Input(currentLevel, needsList[i]);
46 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
50 if (this->IsSetup() ==
true)
51 this->GetOStream(
Warnings0) <<
"MueLu::MatlabSmoother::Setup(): Setup() has already been called";
52 vector<RCP<MuemexArg>> InputArgs = processNeeds<Scalar, LocalOrdinal, GlobalOrdinal, Node>(
this, needsSetup_, currentLevel);
53 A_ = Factory::Get<RCP<Matrix>>(currentLevel,
"A");
56 InputArgs.insert(InputArgs.begin(), AmatArg);
58 if (!setupFunction_.length())
59 throw runtime_error(
"Invalid matlab function name");
60 solveData_ =
callMatlab(setupFunction_, solveDataSize_, InputArgs);
61 this->GetOStream(
Statistics1) << description() << endl;
65 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
68 "MueLu::MatlabSmoother::Apply(): Setup() has not been called");
69 using namespace Teuchos;
71 if (InitialGuessIsZero)
74 vector<RCP<MuemexArg>> InputArgs;
78 MultiVector* BPtrNonConst = (MultiVector*)&B;
82 InputArgs.push_back(XData);
83 InputArgs.push_back(BData);
84 for (
size_t i = 0; i < solveData_.size(); i++)
85 InputArgs.push_back(solveData_[i]);
86 if (!solveFunction_.length())
throw std::runtime_error(
"Invalid matlab function name");
87 vector<Teuchos::RCP<MuemexArg>> mexOutput =
callMatlab(solveFunction_, 1, InputArgs);
89 X = *(mydata->getData());
92 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
95 smoother->SetParameterList(this->GetParameterList());
99 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
101 std::ostringstream out;
103 out <<
"Matlab Smoother(" << setupFunction_ <<
"/" << solveFunction_ <<
")";
110 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
115 out <<
"Matlab Smoother(" << setupFunction_ <<
"/" << solveFunction_ <<
")";
118 out0 <<
"Parameter list: " << std::endl;
120 out << this->GetParameterList();
123 if (verbLevel &
Debug) {
149 #endif // HAVE_MUELU_MATLAB
150 #endif // MUELU_MATLABSMOOTHER_DEF_HPP
Important warning messages (one line)
void SetParameterList(const Teuchos::ParameterList ¶mList)
Set parameters from a parameter list and return with default values.
void DeclareInput(Level ¤tLevel) const
Input.
std::vector< std::string > tokenizeList(const std::string ¶ms)
void Apply(MultiVector &X, const MultiVector &B, bool InitialGuessIsZero=false) const
Apply the preconditioner.
RCP< SmootherPrototype > Copy() const
std::vector< RCP< MuemexArg > > callMatlab(std::string function, int numOutputs, std::vector< RCP< MuemexArg > > args)
T & get(const std::string &name, T def_value)
Timer to be used in factories. Similar to Monitor but with additional timers.
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
void Setup(Level ¤tLevel)
Set up the smoother.
Print additional debugging information.
bool IsParamMuemexVariable(const std::string &name)
virtual void SetParameterList(const Teuchos::ParameterList ¶mList)
Set parameters from a parameter list and return with default values.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Class that holds all level-specific information.
void print(Teuchos::FancyOStream &out, const VerbLevel verbLevel=Default) const
Print the object with some verbosity level to an FancyOStream object.
bool IsSetup() const
Get the state of a smoother prototype.
#define MUELU_DESCRIBE
Helper macro for implementing Describable::describe() for BaseClass objects.
friend class MatlabSmoother
Constructor.
Class that encapsulates Matlab smoothers.
std::string description() const
Return a simple one-line description of this object.
Print class parameters (more parameters, more verbose)
Exception throws to report errors in the internal logical of the program.
virtual std::string description() const
Return a simple one-line description of this object.
std::string toString(const T &t)