CreateTimeDecayDirichletBoundaryCondition.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateTimeDecayDirichletBoundaryCondition.h"
 
#include <numeric>
 
#include "BaseLib/ConfigTree.h"
#include "BoundaryCondition.h"
#include "MeshLib/Mesh.h"
#include "NumLib/DOF/LocalToGlobalIndexMap.h"
#include "ParameterLib/Parameter.h"
#include "ParameterLib/Utils.h"
#include "TimeDecayDirichletBoundaryCondition.h"
 
namespace ProcessLib
{
TimeDecayDirichletBoundaryConditionConfig
parseTimeDecayDirichletBoundaryConditionConfig(
    BaseLib::ConfigTree const& config)
{
    config.checkConfigParameter("type", "TimeDecayDirichlet");
 
    auto const parameter_name =
        config.getConfigParameter<std::string>("time_decay_parameter");
    DBUG("Using parameter {:s}", parameter_name);
 
    auto const lower_limit =
        config.getConfigParameter<double>("lower_limit");
 
    return {parameter_name, lower_limit};
}
 
std::unique_ptr<BoundaryCondition> createTimeDecayDirichletBoundaryCondition(
    TimeDecayDirichletBoundaryConditionConfig const& config,
    int const variable_id, int const component_id, MeshLib::Mesh const& bc_mesh,
    NumLib::LocalToGlobalIndexMap const& dof_table_bulk,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters)
{
    DBUG("Create TimeDecayDirichlet.");
    auto const& time_decay_parameter = ParameterLib::findParameter<double>(
        config.parameter_name, parameters, 1, &bc_mesh);
 
    // In case of partitioned mesh the boundary could be empty, i.e. there
    // is no boundary condition.
#ifdef USE_PETSC
    // This can be extracted to createBoundaryCondition() but then the config
    // parameters are not read and will cause an error.
    // TODO (naumov): Add a function to ConfigTree for skipping the tags of the
    // subtree and move the code up in createBoundaryCondition().
    if (bc_mesh.getDimension() == 0 && bc_mesh.getNumberOfNodes() == 0 &&
        bc_mesh.getNumberOfElements() == 0)
    {
        return nullptr;
    }
#endif  // USE_PETSC
 
    return std::make_unique<TimeDecayDirichletBoundaryCondition>(
        variable_id, component_id, bc_mesh, dof_table_bulk,
        time_decay_parameter, config.lower_limit);
}
 
}  // namespace ProcessLib