Tensor.h

Go to the documentation of this file.

/*
 * \file
 * \copyright
 * Copyright (c) 2012-2024, OpenGeoSys Community (http://www.opengeosys.org)
 *            Distributed under a Modified BSD License.
 *              See accompanying file LICENSE.txt or
 *              http://www.opengeosys.org/project/license
 */
 
#pragma once
 
#include <Eigen/Core>
 
#include "MaterialLib/MPL/Property.h"
 
namespace MaterialPropertyLib
{
constexpr int tensorSize(int dim)
{
    if (dim == 1)
    {
        return 3;  // Diagonal entries.
    }
    if (dim == 2)
    {
        return 5;  // 2x2 matrix and the 3rd diagonal entry.
    }
    if (dim == 3)
    {
        return 9;  // Full 3x3 matrix.
    }
    OGS_FATAL("Tensor size for dimension {} is not defined.", dim);
}
 
template <int Dim>
using Tensor = Eigen::Matrix<double, tensorSize(Dim), 1>;
 
}  // namespace MaterialPropertyLib