linearsolverparameters.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
// SPDX-License-Identifier: GPL-3.0-or-later
//
#ifndef DUMUX_LINEAR_SOLVER_PARAMETERS_HH
#define DUMUX_LINEAR_SOLVER_PARAMETERS_HH
 
#include <string>
#include <array>
#include <vector>
 
#include <dune/common/parametertree.hh>
#include <dune/common/std/type_traits.hh>
#include <dumux/common/parameters.hh>
 
namespace Dumux::Detail::LinearSolverParameters {
template <typename T>
using GVDetector = typename T::GridView;
 
template <typename T>
constexpr bool hasGridView = Dune::Std::is_detected<GVDetector, T>::value;
} // end namespace Dumux::Detail::LinearSolverParameters
 
namespace Dumux {
 
template<class LinearSolverTraits>
class LinearSolverParameters
{
public:
 
    static const std::vector<std::array<std::string, 2>> dumuxToIstlSolverParams;
 
    static Dune::ParameterTree createParameterTree(const std::string& paramGroup = "")
    {
        Dune::ParameterTree params;
        setDefaultParameters(params, paramGroup);
        fillValuesForIstlKeys(params, paramGroup);
        return params;
    }
 
    static void setDefaultParameters(Dune::ParameterTree& params, const std::string& paramGroup = "")
    {
        params["restart"] = "10";
        params["maxit"] = "250";
        params["reduction"] = "1e-13";
        params["verbose"] = "0";
        params["preconditioner.iterations"] = "1";
        params["preconditioner.relaxation"] = "1.0";
        params["preconditioner.verbosity"] = "0";
        params["preconditioner.defaultAggregationSizeMode"] = "isotropic";
        if constexpr (Detail::LinearSolverParameters::hasGridView<LinearSolverTraits>)
            params["preconditioner.defaultAggregationDimension"] = std::to_string(LinearSolverTraits::GridView::dimension);
        else
            params["preconditioner.defaultAggregationDimension"] = "3";
        params["preconditioner.maxLevel"] = "100";
        params["ParameterGroup"] = paramGroup;
        params["preconditioner.ParameterGroup"] = paramGroup;
    }
 
    static void fillValuesForIstlKeys(Dune::ParameterTree& params, const std::string& paramGroup = "")
    {
        const auto linearSolverGroups = getParamSubGroups("LinearSolver", paramGroup);
        if (linearSolverGroups.empty()) // no linear solver parameters were specified
            return;
 
        for (const auto& [dumuxKey, istlKey] : dumuxToIstlSolverParams)
        {
            for (const auto& group : linearSolverGroups)
            {
                const auto fullDumuxKey = group + "." + dumuxKey;
                const auto value = getParam<std::string>(fullDumuxKey, "");
                if (!value.empty())
                {
                    params[istlKey] = value;
                    break; // skip groups with smaller depth in the tree
                }
            }
        }
    }
};
 
template<class LinearSolverTraits>
const std::vector<std::array<std::string, 2>>
LinearSolverParameters<LinearSolverTraits>::dumuxToIstlSolverParams =
{
    // solver params
    {"Verbosity", "verbose"},
    {"MaxIterations", "maxit"},
    {"ResidualReduction", "reduction"},
    {"Type", "type"},
    {"GMResRestart", "restart"}, // cycles before restarting
    {"Restart", "restart"}, // cycles before restarting
    {"MaxOrthogonalizationVectors", "mmax"},
 
    // preconditioner params
    {"Preconditioner.Verbosity", "preconditioner.verbosity"},
    {"Preconditioner.Type", "preconditioner.type"},
    {"Preconditioner.Iterations", "preconditioner.iterations"},
    {"Preconditioner.Relaxation", "preconditioner.relaxation"},
    {"Preconditioner.ILUOrder", "preconditioner.n"},
    {"Preconditioner.ILUResort", "preconditioner.resort"},
    {"Preconditioner.AmgSmootherRelaxation", "preconditioner.smootherRelaxation"},
    {"Preconditioner.AmgSmootherIterations", "preconditioner.smootherIterations"},
    {"Preconditioner.AmgMaxLevel", "preconditioner.maxLevel"},
    {"Preconditioner.AmgCoarsenTarget", "preconditioner.coarsenTarget"},
    {"Preconditioner.AmgMinCoarseningRate", "preconditioner.minCoarseningRate"},
    {"Preconditioner.AmgAccumulationMode", "preconditioner.accumulationMode"},
    {"Preconditioner.AmgProlongationDampingFactor", "preconditioner.prolongationDampingFactor"},
    {"Preconditioner.AmgAlpha", "preconditioner.alpha"},
    {"Preconditioner.AmgBeta", "preconditioner.beta"},
    {"Preconditioner.AmgAdditive", "preconditioner.additive"},
    {"Preconditioner.AmgGamma", "preconditioner.gamma"},
    {"Preconditioner.AmgPreSmoothingSteps", "preconditioner.preSteps"},
    {"Preconditioner.AmgPostSmoothingSteps", "preconditioner.postSteps"},
    {"Preconditioner.AmgCriterionSymmetric", "preconditioner.criterionSymmetric"},
    {"Preconditioner.AmgStrengthMeasure", "preconditioner.strengthMeasure"},
    {"Preconditioner.AmgDiagonalRowIndex", "preconditioner.diagonalRowIndex"},
    {"Preconditioner.AmgDefaultAggregationSizeMode", "preconditioner.defaultAggregationSizeMode"},
    {"Preconditioner.AmgDefaultAggregationDimension", "preconditioner.defaultAggregationDimension"},
    {"Preconditioner.AmgMaxAggregateDistance", "preconditioner.maxAggregateDistance"},
    {"Preconditioner.AmgMinAggregateSize", "preconditioner.minAggregateSize"},
    {"Preconditioner.AmgMaxAggregateSize", "preconditioner.maxAggregateSize"}
};
 
} // end namespace Dumux
 
#endif