46#ifndef MUELU_SHIFTEDLAPLACIAN_DECL_HPP
47#define MUELU_SHIFTEDLAPLACIAN_DECL_HPP
50#include <Xpetra_Matrix_fwd.hpp>
51#include <Xpetra_VectorFactory_fwd.hpp>
52#include <Xpetra_MultiVectorFactory_fwd.hpp>
53#include <Xpetra_TpetraMultiVector.hpp>
59#if defined(HAVE_MUELU_IFPACK2) and defined(HAVE_MUELU_TPETRA)
77#include <MueLu_ShiftedLaplacianOperator.hpp>
86#ifdef HAVE_MUELU_TPETRA_INST_INT_INT
87#include <BelosConfigDefs.hpp>
88#include <BelosLinearProblem.hpp>
89#include <BelosSolverFactory.hpp>
90#include <BelosTpetraAdapter.hpp>
109#undef MUELU_SHIFTEDLAPLACIAN_SHORT
112 typedef Tpetra::Vector<SC,LO,GO,NO>
TVEC;
113 typedef Tpetra::MultiVector<SC,LO,GO,NO>
TMV;
114 typedef Tpetra::Operator<SC,LO,GO,NO>
OP;
115#ifdef HAVE_MUELU_TPETRA_INST_INT_INT
116 typedef Belos::LinearProblem<SC,TMV,OP> LinearProblem;
117 typedef Belos::SolverManager<SC,TMV,OP> SolverManager;
118 typedef Belos::SolverFactory<SC,TMV,OP> SolverFactory;
178 void setParameters(Teuchos::RCP< Teuchos::ParameterList > paramList);
186 void setstiff(RCP< Tpetra::CrsMatrix<SC,LO,GO,NO> >& TpetraK);
188 void setmass(RCP< Tpetra::CrsMatrix<SC,LO,GO,NO> >& TpetraM);
189 void setcoords(RCP<MultiVector>& Coords);
215 int solve(
const RCP<TMV> B, RCP<TMV>& X);
217 RCP<MultiVector>& X);
219 RCP<Tpetra::MultiVector<SC,LO,GO,NO> >& X);
221 typename Teuchos::ScalarTraits<Scalar>::magnitudeType
GetResidual();
299 RCP< MueLu::ShiftedLaplacianOperator<SC,LO,GO,NO> >
MueLuOp_;
302#ifdef HAVE_MUELU_TPETRA_INST_INT_INT
304 RCP<LinearProblem> LinearProblem_;
305 RCP<SolverManager> SolverManager_;
306 RCP<SolverFactory> SolverFactory_;
307 RCP<Teuchos::ParameterList> BelosList_;
314#define MUELU_SHIFTEDLAPLACIAN_SHORT
MueLu::DefaultLocalOrdinal LocalOrdinal
MueLu::DefaultScalar Scalar
MueLu::DefaultGlobalOrdinal GlobalOrdinal
Base class for MueLu classes.
Shifted Laplacian Helmholtz solver.
void setPreconditioningMatrix(RCP< Matrix > &P)
RCP< TransPFactory > TransPfact_
RCP< MultiVector > NullSpace_
void setmass(RCP< Matrix > &M)
Tpetra::MultiVector< SC, LO, GO, NO > TMV
double ilu_diagpivotthresh_
std::string ilu_normtype_
int krylov_preconditioner_
RCP< CoarseMapFactory > CoarseMapfact_
std::string ilu_milutype_
RCP< GenericRFactory > Rfact_
void resetLinearProblem()
Tpetra::Vector< SC, LO, GO, NO > TVEC
RCP< MultiVector > Coords_
void setNullSpace(RCP< MultiVector > NullSpace)
RCP< RAPFactory > Acfact_
RCP< MueLu::ShiftedLaplacianOperator< SC, LO, GO, NO > > MueLuOp_
RCP< SmootherPrototype > coarsestSmooProto_
std::vector< SC > levelshifts_
Tpetra::Operator< SC, LO, GO, NO > OP
RCP< RAPShiftFactory > Acshift_
void setcoords(RCP< MultiVector > &Coords)
Teuchos::ScalarTraits< Scalar >::magnitudeType GetResidual()
RCP< TentativePFactory > TentPfact_
RCP< Tpetra::CrsMatrix< SC, LO, GO, NO > > TpetraA_
Teuchos::ParameterList precList_
RCP< UncoupledAggregationFactory > UCaggfact_
RCP< SmootherPrototype > smooProto_
RCP< SmootherFactory > smooFact_
int solve(const RCP< TMV > B, RCP< TMV > &X)
void setLevelShifts(std::vector< Scalar > levelshifts)
void setParameters(Teuchos::RCP< Teuchos::ParameterList > paramList)
std::string schwarz_ordermethod_
RCP< FactoryManager > Manager_
RCP< CoupledAggregationFactory > Aggfact_
virtual ~ShiftedLaplacian()
void setstiff(RCP< Matrix > &K)
void setProblemMatrix(RCP< Matrix > &A)
void multigrid_apply(const RCP< MultiVector > B, RCP< MultiVector > &X)
Teuchos::ParameterList coarsestSmooList_
RCP< Hierarchy > Hierarchy_
ShiftedLaplacian()
Constructors.
RCP< AmalgamationFactory > Amalgfact_
std::string ilu_drop_rule_
RCP< CoalesceDropFactory > Dropfact_
RCP< SmootherFactory > coarsestSmooFact_
RCP< PgPFactory > PgPfact_
Tpetra::CombineMode schwarz_combinemode_
Namespace for MueLu classes and methods.
KokkosClassic::DefaultNode::DefaultNodeType DefaultNode
Tpetra::Details::DefaultTypes::scalar_type DefaultScalar