ilu.cc

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.
//
// ---------------------------------------------------------------------
 
 
 
#include <parsed_lac/ilu.h>
 
#ifdef DEAL_II_WITH_TRILINOS
 
using namespace dealii;
 
namespace ParsedLAC
{
  ILUPreconditioner::ILUPreconditioner(const std::string  &name,
                                       const unsigned int &ilu_fill,
                                       const double       &ilu_atol,
                                       const double       &ilu_rtol,
                                       const unsigned int &overlap)
    : ParameterAcceptor(name)
    , PreconditionILU()
    , ilu_fill(ilu_fill)
    , ilu_atol(ilu_atol)
    , ilu_rtol(ilu_rtol)
    , overlap(overlap)
  {
    add_parameters();
  }
 
  void
  ILUPreconditioner::add_parameters()
  {
    add_parameter("Fill-in", ilu_fill, "Additional fill-in.");
 
    add_parameter("ILU atol",
                  ilu_atol,
                  "The amount of perturbation to add to diagonal entries.");
 
    add_parameter("ILU rtol", ilu_rtol, "Scaling factor for diagonal entries.");
 
    add_parameter("Overlap", overlap, "Overlap between processors.");
  }
 
  template <typename Matrix>
  void
  ILUPreconditioner::initialize_preconditioner(const Matrix &matrix)
  {
    TrilinosWrappers::PreconditionILU::AdditionalData data;
 
    data.ilu_fill = ilu_fill;
    data.ilu_atol = ilu_atol;
    data.ilu_rtol = ilu_rtol;
    data.overlap  = overlap;
    this->initialize(matrix, data);
  }
} // namespace ParsedLAC
 
template void
ParsedLAC::ILUPreconditioner::initialize_preconditioner<
  dealii::TrilinosWrappers::SparseMatrix>(
  const dealii::TrilinosWrappers::SparseMatrix &);
 
#endif