MueLu_LineDetectionFactory_decl.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//        MueLu: A package for multigrid based preconditioning
//
// Copyright 2012 NTESS and the MueLu contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef MUELU_LINEDETECTIONFACTORY_DECL_HPP
#define MUELU_LINEDETECTIONFACTORY_DECL_HPP
 
// same as in SemiCoarsenPFactory (TODO rework this)
#define VERTICAL 1
#define HORIZONTAL 2
#define GRID_SUPPLIED -1
 
#include "MueLu_ConfigDefs.hpp"
#include "MueLu_LineDetectionFactory_fwd.hpp"
 
#include "MueLu_Level_fwd.hpp"
#include "MueLu_SingleLevelFactoryBase.hpp"
 
namespace MueLu {

 
template <class Scalar        = DefaultScalar,
          class LocalOrdinal  = DefaultLocalOrdinal,
          class GlobalOrdinal = DefaultGlobalOrdinal,
          class Node          = DefaultNode>
class LineDetectionFactory : public SingleLevelFactoryBase {
#undef MUELU_LINEDETECTIONFACTORY_SHORT
#include "MueLu_UseShortNames.hpp"
 
 public:
  using coordinate_type       = typename Teuchos::ScalarTraits<SC>::coordinateType;
  using CoordinateMultiVector = typename Xpetra::MultiVector<coordinate_type, LO, GO, NO>;


 
  LineDetectionFactory()
    : Zorientation_(VERTICAL) {}

  virtual ~LineDetectionFactory() {}
 
  RCP<const ParameterList> GetValidParameterList() const;



 
  void DeclareInput(Level& currentLevel) const;




  void Build(Level& currentLevel) const;

 
 private:
  void sort_coordinates(LO numCoords, LO* OrigLoc,
                        coordinate_type* xvals,
                        coordinate_type* yvals,
                        coordinate_type* zvals,
                        coordinate_type* xtemp,
                        coordinate_type* ytemp,
                        coordinate_type* ztemp,
                        bool flipXY = false) const;
 
  LO ML_compute_line_info(LO LayerId[], LO VertLineId[],
                          LO Ndof, LO DofsPerNode,
                          LO MeshNumbering, LO NumNodesPerVertLine,
                          coordinate_type* xvals, coordinate_type* yvals, coordinate_type* zvals,
                          const Teuchos::Comm<int>& comm) const;
 
  void ML_az_dsort2(coordinate_type dlist[], LO N, LO list2[]) const;

  mutable LO Zorientation_;
 
};  // class LineDetectionFactory
 
}  // namespace MueLu
 
#define MUELU_LINEDETECTIONFACTORY_SHORT
#endif  // MUELU_LINEDETECTIONFACTORY_DECL_HPP