FormEigenVector.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "FormEigenVector.h"
 
#include "MaterialLib/MPL/PropertyType.h"
#include "MathLib/FormattingUtils.h"
 
namespace MaterialPropertyLib
{
template <int GlobalDim>
struct FormEigenVector
{
    Eigen::Matrix<double, GlobalDim, 1> operator()(double const value) const
    {
        if constexpr (GlobalDim == 1)
        {
            return Eigen::Matrix<double, 1, 1>{value};
        }
        if constexpr (GlobalDim == 2)
        {
            return Eigen::Matrix<double, 2, 1>{value, value};
        }
        if constexpr (GlobalDim == 3)
        {
            return Eigen::Matrix<double, 3, 1>{value, value, value};
        }
        OGS_FATAL("Cannot convert a scalar {} to a {:d}d vector.", value,
                  GlobalDim);
    }
 
    Eigen::Matrix<double, GlobalDim, 1> operator()(
        Eigen::Vector2d const& values) const
    {
        if constexpr (GlobalDim == 2)
        {
            return values;
        }
        OGS_FATAL(
            "Cannot convert a 2d vector with values [{}] to a {:d}d vector.",
            values, GlobalDim);
    }
 
    Eigen::Matrix<double, GlobalDim, 1> operator()(
        Eigen::Vector3d const& values) const
    {
        if constexpr (GlobalDim == 3)
        {
            return values;
        }
        OGS_FATAL(
            "Cannot convert a 3d vector with values [{}] to a {:d}d vector.",
            values, GlobalDim);
    }
 
    Eigen::Matrix<double, GlobalDim, 1> operator()(
        Eigen::Matrix<double, 2, 2> const& values) const
    {
        OGS_FATAL(
            "Cannot convert a 2d tensor with values [{}] to a {:d}d Vector.",
            values, GlobalDim);
    }
    Eigen::Matrix<double, GlobalDim, 1> operator()(
        Eigen::Matrix<double, 3, 3> const& values) const
    {
        OGS_FATAL(
            "Cannot convert a 3d tensor with values [{}] to a {:d}d Vector.",
            values, GlobalDim);
    }
 
    Eigen::Matrix<double, GlobalDim, 1> operator()(
        Eigen::Matrix<double, 4, 1> const& values) const
    {
        OGS_FATAL(
            "Cannot convert a 4d vector with values [{}] to a {:d}d vector.",
            values, GlobalDim);
    }
 
    Eigen::Matrix<double, GlobalDim, 1> operator()(
        Eigen::Matrix<double, 6, 1> const& values) const
    {
        OGS_FATAL(
            "Cannot convert a 6d vector with values [{}] to a {:d}d vector.",
            values, GlobalDim);
    }
 
    Eigen::Matrix<double, GlobalDim, 1> operator()(
        Eigen::MatrixXd const& values) const
    {
        OGS_FATAL(
            "Cannot convert a dynamic {}x{} Eigen matrix with values [{}] to a "
            "{:d}d vector ",
            values.rows(), values.cols(), values, GlobalDim);
    }
};
 
template <int GlobalDim>
Eigen::Matrix<double, GlobalDim, 1> formEigenVector(
    MaterialPropertyLib::PropertyDataType const& values)
{
    return std::visit(FormEigenVector<GlobalDim>(), values);
}
 
template Eigen::Matrix<double, 1, 1> formEigenVector<1>(
    MaterialPropertyLib::PropertyDataType const& values);
 
template Eigen::Matrix<double, 2, 1> formEigenVector<2>(
    MaterialPropertyLib::PropertyDataType const& values);
 
template Eigen::Matrix<double, 3, 1> formEigenVector<3>(
    MaterialPropertyLib::PropertyDataType const& values);
 
}  // namespace MaterialPropertyLib