amg_muelu.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2022 by Luca Heltai
//
// This file is part of the FSI-suite platform, based on the deal.II library.
//
// The FSI-suite platform is free software; you can use it, redistribute it,
// and/or modify it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 3.0 of the License,
// or (at your option) any later version. The full text of the license can be
// found in the file LICENSE at the top level of the FSI-suite platform
// distribution.
//
// ---------------------------------------------------------------------
 
#ifndef amg_muelu_preconditioner_h
#define amg_muelu_preconditioner_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/parameter_acceptor.h>
 
#if defined(DEAL_II_WITH_TRILINOS) && defined(DEAL_II_TRILINOS_WITH_MUELU)
 
#  include <deal.II/lac/trilinos_precondition.h>
 
namespace ParsedLAC
{
  class AMGMueLuPreconditioner
    : public dealii::ParameterAcceptor,
      public dealii::TrilinosWrappers::PreconditionAMGMueLu
  {
  public:
    AMGMueLuPreconditioner(const std::string  &name                  = "",
                           const bool         &elliptic              = true,
                           const unsigned int &n_cycles              = 1,
                           const bool         &w_cycle               = false,
                           const double       &aggregation_threshold = 1e-4,
                           const unsigned int &smoother_sweeps       = 2,
                           const unsigned int &smoother_overlap      = 0,
                           const bool         &output_details        = false,
                           const std::string  &smoother_type = "Chebyshev",
                           const std::string  &coarse_type   = "Amesos-KLU");
 
    template <typename Matrix>
    void
    initialize_preconditioner(const Matrix &matrix);
 
    using dealii::TrilinosWrappers::PreconditionAMGMueLu::initialize;
 
  private:
    void
    add_parameters();
 
    bool elliptic;
 
    unsigned int n_cycles;
 
    bool w_cycle;
 
    double aggregation_threshold;
 
    unsigned int smoother_sweeps;
 
    unsigned int smoother_overlap;
 
    bool output_details;
 
    std::string smoother_type;
 
    std::string coarse_type;
  };
 
 
} // namespace ParsedLAC
 
#endif // DEAL_II_WITH_TRILINOS
 
#endif