TimeDecayDirichletBoundaryCondition.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <memory>
 
#include "MathLib/LinAlg/GlobalMatrixVectorTypes.h"
#include "NumLib/DOF/LocalToGlobalIndexMap.h"
#include "ProcessLib/BoundaryConditionAndSourceTerm/BoundaryCondition.h"
 
namespace ParameterLib
{
template <typename T>
struct Parameter;
}  // namespace ParameterLib
 
namespace MeshLib
{
class Mesh;
}
 
namespace ProcessLib
{
class TimeDecayDirichletBoundaryCondition final : public BoundaryCondition
{
public:
    explicit TimeDecayDirichletBoundaryCondition(
        int const variable_id, int const component_id,
        MeshLib::Mesh const& boundary_mesh,
        NumLib::LocalToGlobalIndexMap const& dof_table_bulk,
        ParameterLib::Parameter<double> const& time_decay_parameter,
        double const lower_limit);
 
    void getEssentialBCValues(
        const double t, GlobalVector const& x,
        NumLib::IndexValueVector<GlobalIndexType>& bc_values) const override;
 
private:
    int const variable_id_;
    int const component_id_;
 
    MeshLib::Mesh const& boundary_mesh_;
 
    std::unique_ptr<NumLib::LocalToGlobalIndexMap> const dof_table_boundary_;

    ParameterLib::Parameter<double> const& time_decay_parameter_;

    double const lower_limit_;

    mutable std::vector<double> initial_x_values_;
 
    mutable bool initial_values_are_set_ = false;
 
    void setInitialValues(GlobalVector const& x) const;
};
 
}  // namespace ProcessLib