jacobi.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 jacobi_preconditioner_h
#define jacobi_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 JacobiPreconditioner
    : public dealii::ParameterAcceptor,
      public dealii::TrilinosWrappers::PreconditionJacobi
  {
  public:
    JacobiPreconditioner(const std::string  &name         = "",
                         const double       &omega        = 1,
                         const double       &min_diagonal = 0,
                         const unsigned int &n_sweeps     = 1);
 
    template <typename Matrix>
    void
    initialize_preconditioner(const Matrix &matrix);
 
    using dealii::TrilinosWrappers::PreconditionJacobi::initialize;
 
  private:
    void
    add_parameters();
 
    double omega;
 
    double min_diagonal;
 
    unsigned int n_sweeps;
  };
 
} // namespace ParsedLAC
 
#endif // DEAL_II_WITH_TRILINOS
 
#endif