SymmTensor.hpp

/*
  Copyright 2025 Equinor ASA

  This file is part of the Open Porous Media Project (OPM).

  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.

  OPM is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef OPM_UTIL_SYMM_TENSOR_HPP
#define OPM_UTIL_SYMM_TENSOR_HPP

namespace Dune {
template<class K, int Size> class FieldVector;
}

#include <opm/common/utility/VoigtArray.hpp>

#include <utility>

namespace Opm {

template<class T>
class SymmTensor : public VoigtContainer<T>
{
public:
    using field_type = T;

    SymmTensor() = default;
    SymmTensor(std::initializer_list<T> value)
        : VoigtContainer<T>(std::move(value))
    {}

    void operator+=(const T data);
    void operator+=(const SymmTensor<T>& data);

    void operator*=(const T data);

    SymmTensor<T>& operator=(const T value);

    void reset();

    T trace() const;
#if HAVE_DUNE_COMMON
    T traction(const Dune::FieldVector<T,3>& normal) const;
#endif
};

template<class T1, class T2>
SymmTensor<T1> operator*(const T2 value, SymmTensor<T1> t1);
template<class T1, class T2>
SymmTensor<T1> operator*(SymmTensor<T1> t1, const T2 value);
template<class T>
SymmTensor<T> operator+(SymmTensor<T> t1, const SymmTensor<T>& t2);

} // namespace Opm

#endif // OPM_UTIL_SYMM_TENSOR_HPP