amg.h

Go to the documentation of this file.

//-----------------------------------------------------------
//
//    Copyright (C) 2015 by the deal2lkit authors
//
//    This file is part of the deal2lkit library.
//
//    The deal2lkit library 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 2.1 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 deal2lkit distribution.
//
//-----------------------------------------------------------
 
#ifndef amg_preconditioner_h
#define amg_preconditioner_h
 
#include <deal.II/base/config.h>
 
#ifdef DEAL_II_WITH_TRILINOS
 
#  include <deal.II/base/parameter_acceptor.h>
 
#  include <deal.II/lac/trilinos_precondition.h>
 
namespace ParsedLAC
{
  class AMGPreconditioner : public dealii::ParameterAcceptor,
                            public dealii::TrilinosWrappers::PreconditionAMG
  {
  public:
    AMGPreconditioner(const std::string  &name                  = "",
                      const bool         &elliptic              = true,
                      const bool         &higher_order_elements = false,
                      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(const Matrix &matrix);
 
    using dealii::TrilinosWrappers::PreconditionAMG::initialize;
 
    void
    set_constant_modes(const std::vector<std::vector<bool>> &constant_modes);
 
  private:
    void
    add_parameters();
 
    bool elliptic;
 
    bool higher_order_elements;
 
    unsigned int n_cycles;
 
    bool w_cycle;
 
    double aggregation_threshold;
 
    std::string var_const_modes;
 
    unsigned int smoother_sweeps;
 
    unsigned int smoother_overlap;
 
    bool output_details;
 
    std::string smoother_type;
 
    std::string coarse_type;
 
    std::vector<std::vector<bool>> constant_modes;
  };
 
 
} // namespace ParsedLAC
 
#endif // DEAL_II_WITH_TRILINOS
 
#endif