identity.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_TENSOR_TRAIN_OPERATOR_IDENTITY_H
#define pRC_TENSOR_TRAIN_OPERATOR_IDENTITY_H
 
#include <prc/core/basic/functions/identity.hpp>
#include <prc/core/functors/eval.hpp>
#include <prc/core/tensor/functions/direct_sum.hpp>
#include <prc/tensor_train/operator/views/enumerate.hpp>
 
namespace pRC
{
    template<class T>
    struct Identity<T, If<TensorTrain::IsOperatorView<T>>>
        : Identity<ResultOf<Eval, T>>
    {
    };
 
    template<class T>
    struct Identity<T, If<TensorTrain::IsOperator<T>>>
    {
        constexpr auto operator()()
        {
            auto const f = []<Index C>()
            {
                using Core = typename T::template Cores<C>;
                constexpr auto CRL = Core::size(0);
                constexpr auto CM = Core::size(1);
                constexpr auto CN = Core::size(2);
                constexpr auto CRR = Core::size(3);
 
                return permute<2, 0, 1, 3>(exclude<DirectSum, 1, 2>(
                    reshape<1, CM, CN, 1>(identity<
                        typename Core::template ChangeSizes<CM, CN>>()),
                    zero<typename Core::template ChangeSizes<CRL - 1, CM, CN,
                        CRR - 1>>()));
            };
 
            using F = RemoveConstReference<decltype(f)>;
            using M = typename T::M;
            using N = typename T::N;
            using Ranks = typename T::Ranks;
 
            return TensorTrain::OperatorViews::Enumerate<typename T::Type, M, N,
                Ranks, F>(f);
        }
 
        template<class X, If<IsConstructible<typename T::Type, X>> = 0>
        constexpr auto operator()(X &&value)
        {
            auto const f =
                [value = typename T::Type(forward<X>(value))]<Index C>()
            {
                using Core = typename T::template Cores<C>;
                constexpr auto CRL = Core::size(0);
                constexpr auto CM = Core::size(1);
                constexpr auto CN = Core::size(2);
                constexpr auto CRR = Core::size(3);
 
                if constexpr(C == 0)
                {
                    return permute<2, 0, 1, 3>(exclude<DirectSum, 1, 2>(
                        reshape<1, CM, CN, 1>(identity<
                            typename Core::template ChangeSizes<CM, CN>>(
                            value)),
                        zero<typename Core::template ChangeSizes<CRL - 1, CM,
                            CN, CRR - 1>>()));
                }
                else
                {
                    return permute<2, 0, 1, 3>(exclude<DirectSum, 1, 2>(
                        reshape<1, CM, CN, 1>(identity<
                            typename Core::template ChangeSizes<CM, CN>>()),
                        zero<typename Core::template ChangeSizes<CRL - 1, CM,
                            CN, CRR - 1>>()));
                }
            };
 
            using F = RemoveConstReference<decltype(f)>;
            using M = typename T::M;
            using N = typename T::N;
            using Ranks = typename T::Ranks;
 
            return TensorTrain::OperatorViews::Enumerate<typename T::Type, M, N,
                Ranks, F>(f);
        }
    };
}
#endif // pRC_TENSOR_TRAIN_OPERATOR_IDENTITY_H