46 #ifndef MUELU_GENERALGEOMETRICPFACTORY_DECL_HPP 47 #define MUELU_GENERALGEOMETRICPFACTORY_DECL_HPP 49 #include <Teuchos_SerialDenseVector.hpp> 51 #include <Xpetra_MultiVector.hpp> 52 #include <Xpetra_Matrix_fwd.hpp> 55 #include "MueLu_PFactory.hpp" 62 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
115 template <
class Scalar =
double,
class LocalOrdinal =
int,
class GlobalOrdinal = LocalOrdinal,
class Node = KokkosClassic::DefaultNode::DefaultNodeType>
117 #undef MUELU_GENERALGEOMETRICPFACTORY_SHORT 202 RCP<const Map> fineCoordsMap, RCP<GeometricData> myGeometry,
203 RCP<NodesIDs> ghostedCoarseNodes,
204 Array<Array<GO> >& lCoarseNodesGIDs)
const;
207 RCP<const Map> fineCoordsMap, RCP<GeometricData> myGeometry,
208 RCP<NodesIDs> ghostedCoarseNodes,
209 Array<Array<GO> >& lCoarseNodesGIDs)
const;
212 const RCP<Xpetra::MultiVector<double,LO,GO,NO> >& fCoords,
213 const LO nnzP,
const LO dofsPerNode,
214 RCP<const Map>& stridedDomainMapP,
215 RCP<Matrix> & Amat, RCP<Matrix>& P,
216 RCP<Xpetra::MultiVector<double,LO,GO,NO> >& cCoords,
217 RCP<NodesIDs> ghostedCoarseNodes, Array<Array<GO> > coarseNodesGIDs,
218 int interpolationOrder)
const;
220 void ComputeStencil(
const LO numDimension,
const Array<GO> currentNodeIndices,
221 const Array<GO> coarseNodeIndices,
const LO rate[3],
222 const Array<Array<double> > coord,
const int interpolationOrder,
223 std::vector<double>& stencil)
const;
226 const Array<GO> currentNodeIndices,
227 const Array<GO> coarseNodeIndices,
228 const LO rate[3], std::vector<double>& stencil)
const;
231 std::vector<double>& stencil)
const;
233 const Teuchos::SerialDenseVector<LO,double> parameters,
234 double functions[4][8])
const;
237 const typename Teuchos::Array<LocalOrdinal>::iterator& first1,
238 const typename Teuchos::Array<LocalOrdinal>::iterator& last1,
239 const typename Teuchos::Array<LocalOrdinal>::iterator& first2,
240 const typename Teuchos::Array<LocalOrdinal>::iterator& last2)
const;
243 const typename Teuchos::Array<LocalOrdinal>::iterator& first1,
244 const typename Teuchos::Array<LocalOrdinal>::iterator& last1,
245 const typename Teuchos::Array<LocalOrdinal>::iterator& first2,
246 const typename Teuchos::Array<LocalOrdinal>::iterator& last2)
const;
249 const Array<LO> coarseNodeFineIndices,
250 const RCP<GeometricData> myGeo,
const LO myRankIndex,
const LO pi,
251 const LO pj,
const LO pk,
252 const typename std::vector<std::vector<GO> >::iterator blockStart,
253 const typename std::vector<std::vector<GO> >::iterator blockEnd,
254 GO& myGID, LO& myPID, LO& myLID)
const;
260 #define MUELU_GENERALGEOMETRICPFACTORY_SHORT 261 #endif // MUELU_GENERALGEOMETRICPFACTORY_DECL_HPP Array< GO > gFineNodesPerDir
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 ComputeLinearInterpolationStencil(const LO numDimension, const Array< Array< double > > coord, std::vector< double > &stencil) const
std::vector< GO > colInds
virtual ~GeneralGeometricPFactory()
Destructor.
RCP< const ParameterList > GetValidParameterList() const
Return a const parameter list of valid parameters that setParameterList() will accept.
GeneralGeometricPFactory()
Constructor.
std::vector< std::vector< GO > > meshData
void MakeGeneralGeometricP(RCP< GeometricData > myGeo, const RCP< Xpetra::MultiVector< double, LO, GO, NO > > &fCoords, const LO nnzP, const LO dofsPerNode, RCP< const Map > &stridedDomainMapP, RCP< Matrix > &Amat, RCP< Matrix > &P, RCP< Xpetra::MultiVector< double, LO, GO, NO > > &cCoords, RCP< NodesIDs > ghostedCoarseNodes, Array< Array< GO > > coarseNodesGIDs, int interpolationOrder) const
Array< GO > gCoarseNodesPerDir
Namespace for MueLu classes and methods.
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 BuildP(Level &fineLevel, Level &coarseLevel) const
Abstract Build method.
Array< LO > ghostedCoarseNodesPerDir
Array< GO > startGhostedCoarseNode
void MeshLayoutInterface(const int interpolationOrder, const LO blkSize, RCP< const Map > fineCoordsMap, RCP< GeometricData > myGeometry, RCP< NodesIDs > ghostedCoarseNodes, Array< Array< GO > > &lCoarseNodesGIDs) const
Prolongator factory performing geometric coarsening.
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
Class that holds all level-specific information.
void ComputeStencil(const LO numDimension, const Array< GO > currentNodeIndices, const Array< GO > coarseNodeIndices, const LO rate[3], const Array< Array< double > > coord, const int interpolationOrder, std::vector< double > &stencil) const
void GetInterpolationFunctions(const LO numDimension, const Teuchos::SerialDenseVector< LO, double > parameters, double functions[4][8]) const
void ComputeConstantInterpolationStencil(const LO numDimension, const Array< GO > currentNodeIndices, const Array< GO > coarseNodeIndices, const LO rate[3], std::vector< double > &stencil) const
Array< LO > lCoarseNodesPerDir
void GetCoarsePoints(const int interpolationOrder, const LO blkSize, RCP< const Map > fineCoordsMap, RCP< GeometricData > myGeometry, RCP< NodesIDs > ghostedCoarseNodes, Array< Array< GO > > &lCoarseNodesGIDs) const
void Build(Level &fineLevel, Level &coarseLevel) const
Build an object with this factory.
Factory that provides an interface for a concrete implementation of a prolongation operator...
void DeclareInput(Level &fineLevel, Level &coarseLevel) const
Input.
Array< LO > lFineNodesPerDir