identity.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_CORE_TENSOR_IDENTITY_H
#define pRC_CORE_TENSOR_IDENTITY_H
 
#include <prc/core/basic/functions/unit.hpp>
#include <prc/core/basic/identity.hpp>
#include <prc/core/functors/eval.hpp>
#include <prc/core/tensor/operator/functions/diagonal.hpp>
 
namespace pRC
{
    template<IsTensor T>
        requires(isEven(T::Dimension))
    struct Identity<T>
    {
        constexpr auto operator()() const
        {
            return diagonal(unit<T>());
        }
 
        template<class X>
            requires IsConstructible<typename T::Type, X>
        constexpr auto operator()(X &&value) const
        {
            return diagonal(unit<T>(forward<X>(value)));
        }
    };
 
    template<IsTensorView T>
    struct Identity<T> : Identity<ResultOf<Eval, T>>
    {
    };
}
#endif // pRC_CORE_TENSOR_IDENTITY_H