Ifpack2_Details_Fic_decl.hpp

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


#ifndef __IFPACK2_FIC_DECL_HPP__
#define __IFPACK2_FIC_DECL_HPP__

#include <Ifpack2_Details_FastILU_Base.hpp>

// forward-declare the local preconditioner type
template <typename LocalOrdinal, typename Scalar, typename execution_space>
class FastICPrec;

namespace Ifpack2 {
namespace Details {

template <typename Scalar, typename LocalOrdinal, typename GlobalOrdinal, typename Node>
class Fic : public FastILU_Base<Scalar, LocalOrdinal, GlobalOrdinal, Node> {
 public:
  typedef FastILU_Base<Scalar, LocalOrdinal, GlobalOrdinal, Node> Base;
  typedef typename Base::TRowMatrix TRowMatrix;
  typedef typename Base::ImplScalar ImplScalar;
  typedef typename Base::ImplScalarArray ImplScalarArray;
  typedef FastICPrec<LocalOrdinal, ImplScalar, typename Base::execution_space> LocalFIC;

  Fic(Teuchos::RCP<const TRowMatrix> mat_);

  int getSweeps() const;

  std::string getSpTrsvType() const;

  int getNTrisol() const;

  void checkLocalIC() const;

 protected:
  mutable Teuchos::RCP<LocalFIC> localPrec_;

  void initLocalPrec();
  // compute() takes A's local values
  void computeLocalPrec();
  void applyLocalPrec(ImplScalarArray x, ImplScalarArray y) const;
  std::string getName() const;
};

}  // namespace Details
}  // namespace Ifpack2

#endif