MueLu  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node > Class Template Reference

Prolongator factory performing geometric coarsening. More...

#include <MueLu_GeneralGeometricPFactory_decl.hpp>

Inheritance diagram for MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >:
MueLu::PFactory MueLu::TwoLevelFactoryBase MueLu::Factory MueLu::FactoryBase MueLu::FactoryAcceptor MueLu::ParameterListAcceptorImpl MueLu::BaseClass MueLu::ParameterListAcceptor MueLu::VerboseObject MueLu::Describable Teuchos::VerboseObject< VerboseObject > Teuchos::Describable Teuchos::VerboseObjectBase Teuchos::LabeledObject

Classes

struct  GeometricData
 
struct  NodeID
 
struct  NodesIDs
 

Public Member Functions

RCP< const ParameterListGetValidParameterList () const
 Return a const parameter list of valid parameters that setParameterList() will accept. More...
 
- Public Member Functions inherited from MueLu::PFactory
 PFactory ()
 Constructor. More...
 
virtual ~PFactory ()
 Destructor. More...
 
void setRestrictionMode (bool bRestrictionMode=false)
 
bool isRestrictionModeSet ()
 returns restrictionMode flag More...
 
- Public Member Functions inherited from MueLu::TwoLevelFactoryBase
 TwoLevelFactoryBase ()
 Constructor. More...
 
virtual ~TwoLevelFactoryBase ()
 Destructor. More...
 
virtual void CallDeclareInput (Level &requestedLevel) const
 
virtual void CallBuild (Level &requestedLevel) const
 
- Public Member Functions inherited from MueLu::Factory
void EnableMultipleCallCheck () const
 Constructor. More...
 
void DisableMultipleCallCheck () const
 
void ResetDebugData () const
 
- Public Member Functions inherited from MueLu::FactoryBase
 FactoryBase ()
 Constructor. More...
 
virtual ~FactoryBase ()
 Destructor. More...
 
int GetID () const
 return unique factory id More...
 
- Public Member Functions inherited from MueLu::BaseClass
virtual ~BaseClass ()
 Destructor. More...
 
- Public Member Functions inherited from MueLu::VerboseObject
VerbLevel GetVerbLevel () const
 Get the verbosity level. More...
 
void SetVerbLevel (const VerbLevel verbLevel)
 Set the verbosity level of this object. More...
 
int GetProcRankVerbose () const
 Get proc rank used for printing. Do not use this information for any other purpose. More...
 
int SetProcRankVerbose (int procRank) const
 Set proc rank used for printing. More...
 
bool IsPrint (MsgType type, int thisProcRankOnly=-1) const
 Find out whether we need to print out information for a specific message type. More...
 
Teuchos::FancyOStreamGetOStream (MsgType type, int thisProcRankOnly=0) const
 Get an output stream for outputting the input message type. More...
 
Teuchos::FancyOStreamGetBlackHole () const
 
 VerboseObject ()
 
virtual ~VerboseObject ()
 Destructor. More...
 
- Public Member Functions inherited from Teuchos::VerboseObject< VerboseObject >
 VerboseObject (const EVerbosityLevel verbLevel=VERB_DEFAULT, const RCP< FancyOStream > &oStream=Teuchos::null)
 
virtual const VerboseObjectsetVerbLevel (const EVerbosityLevel verbLevel) const
 
virtual const VerboseObjectsetOverridingVerbLevel (const EVerbosityLevel verbLevel) const
 
virtual EVerbosityLevel getVerbLevel () const
 
TEUCHOSPARAMETERLIST_LIB_DLL_EXPORT
RCP< const ParameterList
getValidVerboseObjectSublist ()
 
TEUCHOSPARAMETERLIST_LIB_DLL_EXPORT
void 
setupVerboseObjectSublist (ParameterList *paramList)
 
TEUCHOSPARAMETERLIST_LIB_DLL_EXPORT
void 
readVerboseObjectSublist (ParameterList *paramList, RCP< FancyOStream > *oStream, EVerbosityLevel *verbLevel)
 
void readVerboseObjectSublist (ParameterList *paramList, VerboseObject< ObjectType > *verboseObject)
 
- Public Member Functions inherited from Teuchos::VerboseObjectBase
virtual ~VerboseObjectBase ()
 
 VerboseObjectBase (const RCP< FancyOStream > &oStream=Teuchos::null)
 
virtual const VerboseObjectBasesetOStream (const RCP< FancyOStream > &oStream) const
 
virtual const VerboseObjectBasesetOverridingOStream (const RCP< FancyOStream > &oStream) const
 
virtual VerboseObjectBasesetLinePrefix (const std::string &linePrefix)
 
virtual RCP< FancyOStreamgetOStream () const
 
virtual RCP< FancyOStreamgetOverridingOStream () const
 
virtual std::string getLinePrefix () const
 
virtual OSTab getOSTab (const int tabs=1, const std::string &linePrefix="") const
 
- Public Member Functions inherited from MueLu::Describable
virtual ~Describable ()
 Destructor. More...
 
virtual std::string ShortClassName () const
 Return the class name of the object, without template parameters and without namespace. More...
 
virtual void describe (Teuchos::FancyOStream &out_arg, const VerbLevel verbLevel=Default) const
 
virtual std::string description () const
 Return a simple one-line description of this object. More...
 
void describe (Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
 Print the object with some verbosity level to an FancyOStream object. More...
 
- Public Member Functions inherited from Teuchos::Describable
void describe (std::ostream &out, const EVerbosityLevel verbLevel=verbLevel_default) const
 
 LabeledObject ()
 
virtual ~LabeledObject ()
 
virtual void setObjectLabel (const std::string &objectLabel)
 
virtual std::string getObjectLabel () const
 
DescribableStreamManipulatorState describe (const Describable &describable, const EVerbosityLevel verbLevel=Describable::verbLevel_default)
 
std::ostream & operator<< (std::ostream &os, const DescribableStreamManipulatorState &d)
 
- Public Member Functions inherited from MueLu::FactoryAcceptor
virtual ~FactoryAcceptor ()
 
virtual void SetFactory (const std::string &varName, const RCP< const FactoryBase > &factory)=0
 Configuration. More...
 
virtual const RCP< const
FactoryBase
GetFactory (const std::string &varName) const =0
 
- Public Member Functions inherited from MueLu::ParameterListAcceptorImpl
 ParameterListAcceptorImpl ()
 
virtual ~ParameterListAcceptorImpl ()=default
 
virtual void SetParameterList (const Teuchos::ParameterList &paramList)
 Set parameters from a parameter list and return with default values. More...
 
virtual const
Teuchos::ParameterList
GetParameterList () const
 
void SetParameter (const std::string &name, const ParameterEntry &entry)
 Set a parameter directly as a ParameterEntry. More...
 
const ParameterEntryGetParameter (const std::string &name) const
 Retrieves a const entry with the name name. More...
 
virtual void GetDocumentation (std::ostream &os) const
 
- Public Member Functions inherited from MueLu::ParameterListAcceptor
 ParameterListAcceptor ()
 
virtual ~ParameterListAcceptor ()=default
 

Private Member Functions

void MeshLayoutInterface (const int interpolationOrder, const LO blkSize, RCP< const Map > fineCoordsMap, RCP< GeometricData > myGeometry, RCP< NodesIDs > ghostedCoarseNodes, Array< Array< GO > > &lCoarseNodesGIDs) const
 
void GetCoarsePoints (const int interpolationOrder, const LO blkSize, RCP< const Map > fineCoordsMap, RCP< GeometricData > myGeometry, RCP< NodesIDs > ghostedCoarseNodes, Array< Array< GO > > &lCoarseNodesGIDs) const
 
void MakeGeneralGeometricP (RCP< GeometricData > myGeo, const RCP< Xpetra::MultiVector< typename Teuchos::ScalarTraits< Scalar >::coordinateType, LO, GO, NO > > &fCoords, const LO nnzP, const LO dofsPerNode, RCP< const Map > &stridedDomainMapP, RCP< Matrix > &Amat, RCP< Matrix > &P, RCP< Xpetra::MultiVector< typename Teuchos::ScalarTraits< Scalar >::coordinateType, LO, GO, NO > > &cCoords, RCP< NodesIDs > ghostedCoarseNodes, Array< Array< GO > > coarseNodesGIDs, int interpolationOrder) const
 
void ComputeStencil (const LO numDimension, const Array< GO > currentNodeIndices, const Array< GO > coarseNodeIndices, const LO rate[3], const Array< Array< typename Teuchos::ScalarTraits< Scalar >::coordinateType > > coord, const int interpolationOrder, std::vector< double > &stencil) const
 
void ComputeConstantInterpolationStencil (const LO numDimension, const Array< GO > currentNodeIndices, const Array< GO > coarseNodeIndices, const LO rate[3], std::vector< double > &stencil) const
 
void ComputeLinearInterpolationStencil (const LO numDimension, const Array< Array< typename Teuchos::ScalarTraits< Scalar >::coordinateType > > coord, std::vector< double > &stencil) const
 
void GetInterpolationFunctions (const LO numDimension, const Teuchos::SerialDenseVector< LO, double > parameters, double functions[4][8]) const
 
void sh_sort_permute (const typename Teuchos::Array< LocalOrdinal >::iterator &first1, const typename Teuchos::Array< LocalOrdinal >::iterator &last1, const typename Teuchos::Array< LocalOrdinal >::iterator &first2, const typename Teuchos::Array< LocalOrdinal >::iterator &last2) const
 
void sh_sort2 (const typename Teuchos::Array< LocalOrdinal >::iterator &first1, const typename Teuchos::Array< LocalOrdinal >::iterator &last1, const typename Teuchos::Array< LocalOrdinal >::iterator &first2, const typename Teuchos::Array< LocalOrdinal >::iterator &last2) const
 
void GetGIDLocalLexicographic (const GO i, const GO j, const GO k, const Array< LO > coarseNodeFineIndices, const RCP< GeometricData > myGeo, const LO myRankIndex, const LO pi, const LO pj, const LO pk, const typename std::vector< std::vector< GO > >::iterator blockStart, const typename std::vector< std::vector< GO > >::iterator blockEnd, GO &myGID, LO &myPID, LO &myLID) const
 

Friends

class MueLuTests::GeneralGeometricPFactoryTester< Scalar, LocalOrdinal, GlobalOrdinal, Node >
 

Constructors/Destructors.

 GeneralGeometricPFactory ()
 Constructor. More...
 
virtual ~GeneralGeometricPFactory ()
 Destructor. More...
 
void DeclareInput (Level &fineLevel, Level &coarseLevel) const
 Input. More...
 

Build methods.

void Build (Level &fineLevel, Level &coarseLevel) const
 Build an object with this factory. More...
 
void BuildP (Level &fineLevel, Level &coarseLevel) const
 Abstract Build method. More...
 

Additional Inherited Members

- Static Public Member Functions inherited from MueLu::Factory
static void EnableMultipleCheckGlobally ()
 
static void DisableMultipleCheckGlobally ()
 
- Static Public Member Functions inherited from MueLu::VerboseObject
static void SetMueLuOStream (const Teuchos::RCP< Teuchos::FancyOStream > &mueluOStream)
 
static void SetMueLuOFileStream (const std::string &filename)
 
static Teuchos::RCP
< Teuchos::FancyOStream
GetMueLuOStream ()
 
static void SetDefaultVerbLevel (const VerbLevel defaultVerbLevel)
 Set the default (global) verbosity level. More...
 
static VerbLevel GetDefaultVerbLevel ()
 Get the default (global) verbosity level. More...
 
- Static Public Member Functions inherited from Teuchos::VerboseObject< VerboseObject >
static void setDefaultVerbLevel (const EVerbosityLevel defaultVerbLevel)
 
static EVerbosityLevel getDefaultVerbLevel ()
 
- Static Public Member Functions inherited from Teuchos::VerboseObjectBase
static void setDefaultOStream (const RCP< FancyOStream > &defaultOStream)
 
static RCP< FancyOStreamgetDefaultOStream ()
 
- Static Public Attributes inherited from Teuchos::Describable
static const EVerbosityLevel verbLevel_default
 
- Protected Member Functions inherited from Teuchos::VerboseObject< VerboseObject >
void initializeVerboseObject (const EVerbosityLevel verbLevel=VERB_DEFAULT, const RCP< FancyOStream > &oStream=Teuchos::null)
 
- Protected Member Functions inherited from Teuchos::VerboseObjectBase
void initializeVerboseObjectBase (const RCP< FancyOStream > &oStream=Teuchos::null)
 
virtual void informUpdatedVerbosityState () const
 
- Protected Attributes inherited from MueLu::PFactory
bool restrictionMode_
 

Detailed Description

template<class Scalar = DefaultScalar, class LocalOrdinal = DefaultLocalOrdinal, class GlobalOrdinal = DefaultGlobalOrdinal, class Node = DefaultNode>
class MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >

Prolongator factory performing geometric coarsening.

The geometric algorithm assumes the underlying mesh is reasonably structured. Any rate of coarsening can be applied, and the rate is automatically decrease along an edge if the number of element is not divisible by the coarsening rate. The coarsening rate is allowed to be different in all direction which means that semi-coarsening can be achieved within this algorithm in 1 or 2 directions. The main difficulty is to obtain the number of elements/nodes in each directions to identify coarse nodes and fine nodes.

Input/output of GeneralGeometricPFactory

User parameters of SemiCoarsenPFactory

Parameter type default master.xml validated requested description
Coarsen string -1 A string that specify the coarsening rate, if it is a single character, it will
indicate a unique coarsening rate in each direction, if it is longer, it will be
processed as a vector with 3 entries, one for each spatial direction
A Factory null * * Generating factory of the matrix A used during the prolongator smoothing process
Nullspace Factory null * * Generating factory of the nullspace. The GeneralGeometricPFactory provides
a coarse version of the given Nullspace.
Coordinates Factory NoFactory * * Generating factory for coorindates. The coordinates are expected to be provided
on the finest level using the NoFactory mechanism. The coordinates are used to
compute the coarsening stencil and coarse coordinates are generated for the next
level.

The * in the master.xml column denotes that the parameter is defined in the master.xml file.
The * in the validated column means that the parameter is declared in the list of valid input parameters (see GeneralGeometricCoarsenPFactory::GetValidParameters).
The * in the requested column states that the data is requested as input with all dependencies (see GeneralGeometricCoarsenPFactory::DeclareInput).

Variables provided by GeneralGeometricPFactory

After GeneralGeometricPFactory::Build the following data is available (if requested)

Parameter generated by description
P GeneralGeometricPFactory Prolongator
Nullspace GeneralGeometricPFactory Coarse nullspace (the fine level nullspace information is coarsened using P to generate a coarse version
of the nullspace. No scaling is applied.
coarseCoordinates NoFactory Coarse coordinates that will be used on the next fine level to compute the coarsening stencils

Definition at line 119 of file MueLu_GeneralGeometricPFactory_decl.hpp.

Constructor & Destructor Documentation

template<class Scalar = DefaultScalar, class LocalOrdinal = DefaultLocalOrdinal, class GlobalOrdinal = DefaultGlobalOrdinal, class Node = DefaultNode>
MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::GeneralGeometricPFactory ( )
inline

Constructor.

Definition at line 130 of file MueLu_GeneralGeometricPFactory_decl.hpp.

template<class Scalar = DefaultScalar, class LocalOrdinal = DefaultLocalOrdinal, class GlobalOrdinal = DefaultGlobalOrdinal, class Node = DefaultNode>
virtual MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::~GeneralGeometricPFactory ( )
inlinevirtual

Destructor.

Definition at line 133 of file MueLu_GeneralGeometricPFactory_decl.hpp.

Member Function Documentation

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
RCP< const ParameterList > MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::GetValidParameterList ( ) const
virtual

Return a const parameter list of valid parameters that setParameterList() will accept.

Also define the default values of parameters according to the input parameter list.

Implements MueLu::ParameterListAcceptor.

Definition at line 73 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::DeclareInput ( Level fineLevel,
Level coarseLevel 
) const
virtual

Input.

Implements MueLu::TwoLevelFactoryBase.

Definition at line 105 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::Build ( Level fineLevel,
Level coarseLevel 
) const
virtual

Build an object with this factory.

Implements MueLu::TwoLevelFactoryBase.

Definition at line 137 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::BuildP ( Level fineLevel,
Level coarseLevel 
) const
virtual

Abstract Build method.

Implements MueLu::PFactory.

Definition at line 143 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::MeshLayoutInterface ( const int  interpolationOrder,
const LO  blkSize,
RCP< const Map >  fineCoordsMap,
RCP< GeometricData myGeometry,
RCP< NodesIDs ghostedCoarseNodes,
Array< Array< GO > > &  lCoarseNodesGIDs 
) const
private

Definition at line 324 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::GetCoarsePoints ( const int  interpolationOrder,
const LO  blkSize,
RCP< const Map >  fineCoordsMap,
RCP< GeometricData myGeometry,
RCP< NodesIDs ghostedCoarseNodes,
Array< Array< GO > > &  lCoarseNodesGIDs 
) const
private

Definition at line 559 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::MakeGeneralGeometricP ( RCP< GeometricData myGeo,
const RCP< Xpetra::MultiVector< typename Teuchos::ScalarTraits< Scalar >::coordinateType, LO, GO, NO > > &  fCoords,
const LO  nnzP,
const LO  dofsPerNode,
RCP< const Map > &  stridedDomainMapP,
RCP< Matrix > &  Amat,
RCP< Matrix > &  P,
RCP< Xpetra::MultiVector< typename Teuchos::ScalarTraits< Scalar >::coordinateType, LO, GO, NO > > &  cCoords,
RCP< NodesIDs ghostedCoarseNodes,
Array< Array< GO > >  coarseNodesGIDs,
int  interpolationOrder 
) const
private

Definition at line 755 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::ComputeStencil ( const LO  numDimension,
const Array< GO currentNodeIndices,
const Array< GO coarseNodeIndices,
const LO  rate[3],
const Array< Array< typename Teuchos::ScalarTraits< Scalar >::coordinateType > >  coord,
const int  interpolationOrder,
std::vector< double > &  stencil 
) const
private

Definition at line 1566 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::ComputeConstantInterpolationStencil ( const LO  numDimension,
const Array< GO currentNodeIndices,
const Array< GO coarseNodeIndices,
const LO  rate[3],
std::vector< double > &  stencil 
) const
private

Definition at line 1586 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::ComputeLinearInterpolationStencil ( const LO  numDimension,
const Array< Array< typename Teuchos::ScalarTraits< Scalar >::coordinateType > >  coord,
std::vector< double > &  stencil 
) const
private

Definition at line 1609 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::GetInterpolationFunctions ( const LO  numDimension,
const Teuchos::SerialDenseVector< LO, double >  parameters,
double  functions[4][8] 
) const
private

Definition at line 1697 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::sh_sort_permute ( const typename Teuchos::Array< LocalOrdinal >::iterator &  first1,
const typename Teuchos::Array< LocalOrdinal >::iterator &  last1,
const typename Teuchos::Array< LocalOrdinal >::iterator &  first2,
const typename Teuchos::Array< LocalOrdinal >::iterator &  last2 
) const
private

Definition at line 1754 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::sh_sort2 ( const typename Teuchos::Array< LocalOrdinal >::iterator &  first1,
const typename Teuchos::Array< LocalOrdinal >::iterator &  last1,
const typename Teuchos::Array< LocalOrdinal >::iterator &  first2,
const typename Teuchos::Array< LocalOrdinal >::iterator &  last2 
) const
private

Definition at line 1777 of file MueLu_GeneralGeometricPFactory_def.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >
void MueLu::GeneralGeometricPFactory< Scalar, LocalOrdinal, GlobalOrdinal, Node >::GetGIDLocalLexicographic ( const GO  i,
const GO  j,
const GO  k,
const Array< LO coarseNodeFineIndices,
const RCP< GeometricData myGeo,
const LO  myRankIndex,
const LO  pi,
const LO  pj,
const LO  pk,
const typename std::vector< std::vector< GO > >::iterator  blockStart,
const typename std::vector< std::vector< GO > >::iterator  blockEnd,
GO myGID,
LO myPID,
LO myLID 
) const
private

Definition at line 1802 of file MueLu_GeneralGeometricPFactory_def.hpp.

Friends And Related Function Documentation

template<class Scalar = DefaultScalar, class LocalOrdinal = DefaultLocalOrdinal, class GlobalOrdinal = DefaultGlobalOrdinal, class Node = DefaultNode>
friend class MueLuTests::GeneralGeometricPFactoryTester< Scalar, LocalOrdinal, GlobalOrdinal, Node >
friend

Definition at line 124 of file MueLu_GeneralGeometricPFactory_decl.hpp.


The documentation for this class was generated from the following files: