CreateSaturationDependentSwelling.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "BaseLib/ConfigTree.h"
#include "ParameterLib/CoordinateSystem.h"
#include "SaturationDependentSwelling.h"
 
namespace MaterialPropertyLib
{
std::unique_ptr<SaturationDependentSwelling> createSaturationDependentSwelling(
    BaseLib::ConfigTree const& config,
    ParameterLib::CoordinateSystem const* const local_coordinate_system)
{
    config.checkConfigParameter("type", "SaturationDependentSwelling");
 
    // Second access for storage.
    auto property_name = config.peekConfigParameter<std::string>("name");
 
    DBUG("Create SaturationDependentSwelling solid phase property {:s}.",
         property_name);
 
    auto const swelling_pressures =
        config.getConfigParameter<std::vector<double>>("swelling_pressures");
 
    if (swelling_pressures.size() != 3)
    {
        OGS_FATAL(
            "The number of swelling pressures must be three, but {:d} were "
            "given.",
            swelling_pressures.size());
    }
 
    auto const exponents =
        config.getConfigParameter<std::vector<double>>("exponents");
 
    if (exponents.size() != 3)
    {
        OGS_FATAL("The number of exponents must be three, but {:d} were given.",
                  exponents.size());
    }
 
    if (swelling_pressures.size() != exponents.size())
    {
        OGS_FATAL(
            "The number of swelling pressures and exponents must be equal, but "
            "they are {:d} and {:d}, respectively.",
            swelling_pressures.size(), exponents.size());
    }
 
    auto const lower_saturation_limit =
        config.getConfigParameter<double>("lower_saturation_limit");
 
    auto const upper_saturation_limit =
        config.getConfigParameter<double>("upper_saturation_limit");
 
    return std::make_unique<SaturationDependentSwelling>(
        std::move(property_name),
        std::array<double, 3>{swelling_pressures[0], swelling_pressures[1],
                              swelling_pressures[2]},
        std::array<double, 3>{exponents[0], exponents[1], exponents[2]},
        lower_saturation_limit, upper_saturation_limit,
        local_coordinate_system);
}
}  // namespace MaterialPropertyLib