PETSc/LinearSolverOptionsParser.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "LinearSolverOptionsParser.h"
 
#include "BaseLib/ConfigTree.h"
#include "MathLib/LinAlg/LinearSolverOptions.h"
 
namespace MathLib
{
std::tuple<std::string, std::string>
LinearSolverOptionsParser<PETScLinearSolver>::parseNameAndOptions(
    std::string solver_prefix, BaseLib::ConfigTree const* const config) const
{
    // Insert options into petsc database. Default options are given in the
    // string below.
    std::string petsc_options =
        "-ksp_type cg -pc_type bjacobi -ksp_rtol 1e-16 -ksp_max_it 10000";
 
    if (config)
    {
        ignoreOtherLinearSolvers(*config, "petsc");

        if (auto const subtree = config->getConfigSubtreeOptional("petsc"))
        {
            if (auto const parameters =
                subtree->getConfigParameterOptional<std::string>("parameters"))
            {
                petsc_options = *parameters;
            }
 
            if (auto const pre =
                subtree->getConfigParameterOptional<std::string>("prefix"))
            {
                if (!pre->empty())
                    solver_prefix = *pre + "_";
            }
        }
    }
    return {solver_prefix, petsc_options};
}
}  // namespace MathLib