Tensor.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 "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