Ifpack2_AdditiveSchwarz_decl.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//       Ifpack2: Templated Object-Oriented Algebraic Preconditioner Package
//
// Copyright 2009 NTESS and the Ifpack2 contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER


#ifndef IFPACK2_ADDITIVESCHWARZ_DECL_HPP
#define IFPACK2_ADDITIVESCHWARZ_DECL_HPP

#include "Ifpack2_ConfigDefs.hpp"
#include "Ifpack2_Preconditioner.hpp"
#include "Ifpack2_Details_CanChangeMatrix.hpp"
#include "Ifpack2_Details_NestedPreconditioner.hpp"
#include "Ifpack2_OverlappingRowMatrix.hpp"
#include "Tpetra_Map.hpp"
#include "Tpetra_MultiVector.hpp"
#include "Tpetra_RowMatrix.hpp"
#include <memory>
#include <type_traits>

namespace Trilinos {
namespace Details {
template <class MV, class OP, class NormType>
class LinearSolver;  // forward declaration
}  // namespace Details
}  // namespace Trilinos

namespace Ifpack2 {

template <class MatrixType,
          class LocalInverseType =
              Preconditioner<typename MatrixType::scalar_type,
                             typename MatrixType::local_ordinal_type,
                             typename MatrixType::global_ordinal_type,
                             typename MatrixType::node_type>>
class AdditiveSchwarz : virtual public Preconditioner<typename MatrixType::scalar_type,
                                                      typename MatrixType::local_ordinal_type,
                                                      typename MatrixType::global_ordinal_type,
                                                      typename MatrixType::node_type>,
                        virtual public Details::CanChangeMatrix<Tpetra::RowMatrix<typename MatrixType::scalar_type,
                                                                                  typename MatrixType::local_ordinal_type,
                                                                                  typename MatrixType::global_ordinal_type,
                                                                                  typename MatrixType::node_type>>,
                        virtual public Details::NestedPreconditioner<Preconditioner<typename MatrixType::scalar_type,
                                                                                    typename MatrixType::local_ordinal_type,
                                                                                    typename MatrixType::global_ordinal_type,
                                                                                    typename MatrixType::node_type>> {
 public:
  static_assert(std::is_same<LocalInverseType,
                             Preconditioner<typename MatrixType::scalar_type,
                                            typename MatrixType::local_ordinal_type,
                                            typename MatrixType::global_ordinal_type,
                                            typename MatrixType::node_type>>::value,
                "Ifpack2::AdditiveSchwarz: You are not allowed to use nondefault values for the LocalInverseType template parameter.  Please stop specifying this explicitly.  The default template parameter is perfectly fine.");

  static_assert(std::is_same<MatrixType,
                             Tpetra::RowMatrix<typename MatrixType::scalar_type,
                                               typename MatrixType::local_ordinal_type,
                                               typename MatrixType::global_ordinal_type,
                                               typename MatrixType::node_type>>::value,
                "Ifpack2::AdditiveSchwarz: Please use MatrixType = Tpetra::RowMatrix instead of MatrixType = Tpetra::CrsMatrix.  Don't worry, AdditiveSchwarz's constructor can take either type of matrix; it does a dynamic cast if necessary inside.  Restricting the set of allowed types here will improve build times and reduce library and executable sizes.");



  using scalar_type = typename MatrixType::scalar_type;

  using local_ordinal_type = typename MatrixType::local_ordinal_type;

  using global_ordinal_type = typename MatrixType::global_ordinal_type;

  using node_type = typename MatrixType::node_type;

  using magnitude_type =
      typename Teuchos::ScalarTraits<scalar_type>::magnitudeType;

  using row_matrix_type =
      Tpetra::RowMatrix<scalar_type, local_ordinal_type,
                        global_ordinal_type, node_type>;

  using crs_matrix_type =
      Tpetra::CrsMatrix<scalar_type, local_ordinal_type,
                        global_ordinal_type, node_type>;

  using coord_type =
      Tpetra::MultiVector<magnitude_type, local_ordinal_type,
                          global_ordinal_type, node_type>;

  // \name Constructors and destructor

  AdditiveSchwarz(const Teuchos::RCP<const row_matrix_type>& A);

  AdditiveSchwarz(const Teuchos::RCP<const row_matrix_type>& A,
                  const Teuchos::RCP<const coord_type>& coordinates);

  AdditiveSchwarz(const Teuchos::RCP<const row_matrix_type>& A,
                  const int overlapLevel);

  virtual ~AdditiveSchwarz() = default;



  virtual Teuchos::RCP<const Tpetra::Map<local_ordinal_type, global_ordinal_type, node_type>> getDomainMap() const;

  virtual Teuchos::RCP<const Tpetra::Map<local_ordinal_type, global_ordinal_type, node_type>> getRangeMap() const;

  virtual void
  apply(const Tpetra::MultiVector<scalar_type, local_ordinal_type, global_ordinal_type, node_type>& X,
        Tpetra::MultiVector<scalar_type, local_ordinal_type, global_ordinal_type, node_type>& Y,
        Teuchos::ETransp mode = Teuchos::NO_TRANS,
        scalar_type alpha     = Teuchos::ScalarTraits<scalar_type>::one(),
        scalar_type beta      = Teuchos::ScalarTraits<scalar_type>::zero()) const;



  virtual void
  setInnerPreconditioner(const Teuchos::RCP<Preconditioner<scalar_type,
                                                           local_ordinal_type,
                                                           global_ordinal_type,
                                                           node_type>>& innerPrec);



  virtual void
  setMatrix(const Teuchos::RCP<const row_matrix_type>& A);

  virtual Teuchos::RCP<const row_matrix_type> getMatrix() const;

  void
  setCoord(const Teuchos::RCP<const coord_type>& Coordinates);

  Teuchos::RCP<const coord_type> getCoord() const;

  virtual void setParameters(const Teuchos::ParameterList& plist);

  void
  setParameterList(const Teuchos::RCP<Teuchos::ParameterList>& plist);

  bool supportsZeroStartingSolution() { return true; }

  void setZeroStartingSolution(bool zeroStartingSolution) { ZeroStartingSolution_ = zeroStartingSolution; };

  Teuchos::RCP<const Teuchos::ParameterList> getValidParameters() const;

  virtual void initialize();

  virtual bool isInitialized() const;

  virtual void compute();

  virtual bool isComputed() const;

  virtual int getNumInitialize() const;

  virtual int getNumCompute() const;

  virtual int getNumApply() const;

  virtual double getInitializeTime() const;

  virtual double getComputeTime() const;

  virtual double getApplyTime() const;



  std::string description() const;

  void describe(Teuchos::FancyOStream& out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default) const;


  virtual std::ostream& print(std::ostream& os) const;

  virtual int getOverlapLevel() const;

 private:
  typedef Tpetra::Map<local_ordinal_type,
                      global_ordinal_type,
                      node_type>
      map_type;
  typedef Tpetra::Import<local_ordinal_type,
                         global_ordinal_type,
                         node_type>
      import_type;
  typedef Tpetra::MultiVector<scalar_type,
                              local_ordinal_type,
                              global_ordinal_type,
                              node_type>
      MV;
  typedef Tpetra::Operator<scalar_type,
                           local_ordinal_type,
                           global_ordinal_type,
                           node_type>
      OP;
  typedef Preconditioner<scalar_type,
                         local_ordinal_type,
                         global_ordinal_type,
                         node_type>
      prec_type;

  typedef Trilinos::Details::LinearSolver<MV, OP, typename MV::mag_type> inner_solver_type;

  typedef Kokkos::DualView<local_ordinal_type*, typename crs_matrix_type::device_type> perm_dualview_type;

  AdditiveSchwarz(const AdditiveSchwarz& RHS);

  void setup();

  void localApply(MV& OverlappingX, MV& OverlappingY) const;

  bool hasInnerPrecName() const;

  std::string innerPrecName() const;

  void removeInnerPrecName();

  std::pair<Teuchos::ParameterList, bool> innerPrecParams() const;

  void removeInnerPrecParams();

  static std::string defaultInnerPrecName();

  Teuchos::RCP<const row_matrix_type> Matrix_;

  Teuchos::RCP<OverlappingRowMatrix<row_matrix_type>> OverlappingMatrix_;

  Teuchos::RCP<row_matrix_type> ReorderedLocalizedMatrix_;

  Teuchos::RCP<row_matrix_type> innerMatrix_;

  Teuchos::RCP<const coord_type> Coordinates_ = Teuchos::null;

  bool IsInitialized_ = false;
  bool IsComputed_ = false;
  bool IsOverlapping_ = false;
  int OverlapLevel_ = 0;

  Teuchos::ParameterList List_;

  mutable Teuchos::RCP<const Teuchos::ParameterList> validParams_;

  Tpetra::CombineMode CombineMode_ = Tpetra::ZERO;
  bool AvgOverlap_                 = false;
  bool UseReordering_ = false;
  std::string ReorderingAlgorithm_ = "none";
  bool FilterSingletons_ = false;
  Teuchos::RCP<row_matrix_type> SingletonMatrix_;
  int NumIterations_ = 1;
  bool ZeroStartingSolution_ = true;
  scalar_type UpdateDamping_ = Teuchos::ScalarTraits<scalar_type>::one();

  int NumInitialize_ = 0;
  int NumCompute_ = 0;
  mutable int NumApply_ = 0;
  double InitializeTime_ = 0.0;
  double ComputeTime_ = 0.0;
  mutable double ApplyTime_ = 0.0;
  double InitializeFlops_ = 0.0;
  double ComputeFlops_ = 0.0;
  mutable double ApplyFlops_ = 0.0;
  Teuchos::RCP<inner_solver_type> Inverse_;
  Teuchos::RCP<const map_type> localMap_;
  mutable std::unique_ptr<MV> overlapping_B_;
  mutable std::unique_ptr<MV> overlapping_Y_;
  mutable std::unique_ptr<MV> reduced_reordered_B_;
  mutable std::unique_ptr<MV> reduced_reordered_Y_;
  mutable std::unique_ptr<MV> reduced_B_;
  mutable std::unique_ptr<MV> reduced_Y_;
  mutable std::unique_ptr<MV> reordered_B_;
  mutable std::unique_ptr<MV> reordered_Y_;
  mutable std::unique_ptr<MV> num_overlap_copies_;
  mutable std::unique_ptr<MV> R_;
  mutable std::unique_ptr<MV> C_;
  perm_dualview_type perm_coors;

  mutable Teuchos::RCP<const import_type> DistributedImporter_;
};  // class AdditiveSchwarz

}  // namespace Ifpack2

#endif  // IFPACK2_ADDITIVESCHWARZ_DECL_HPP