DeactivatedSubdomain.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 <Eigen/Core>
#include <memory>
#include <optional>
#include <string>
#include <unordered_set>
#include <vector>
 
#include "MathLib/InterpolationAlgorithms/PiecewiseLinearInterpolation.h"
#include "MeshLib/Mesh.h"
#include "processlib_export.h"
 
namespace MeshLib
{
class Element;
}  // namespace MeshLib
 
namespace ParameterLib
{
template <typename T>
struct Parameter;
}
 
namespace ProcessLib
{
bool isTimeInSupportInterval(
    MathLib::PiecewiseLinearInterpolation const& deactivation_curve,
    double const t);
 
struct DeactivatedSubdomainMesh
{
    MeshLib::Mesh mesh;
    std::unordered_set<std::size_t> bulk_element_ids;

    std::vector<std::size_t> inner_nodes;
    std::vector<std::size_t> outer_nodes;

    std::vector<std::vector<std::size_t>> outer_nodes_elements;
};
 
namespace detail
{
struct Ball
{
    Eigen::Vector3d center;
    double radius;
};
}  // namespace detail
struct DeactivatedSubdomain
{
    bool isDeactivated(MeshLib::Element const& element,
                       double const time) const;
 
    MathLib::PiecewiseLinearInterpolation time_interval;

    std::optional<std::pair<Eigen::Vector3d, Eigen::Vector3d>> line_segment;

    //  problems.
    std::optional<detail::Ball> ball;
 
    DeactivatedSubdomainMesh deactivated_subdomain_mesh;

    ParameterLib::Parameter<double> const* boundary_value_parameter = nullptr;
 
    PROCESSLIB_EXPORT static const std::string zero_parameter_name;
};
}  // namespace ProcessLib