direct_sum.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_CORE_TENSOR_FUNCTIONS_DIRECT_SUM_H
#define pRC_CORE_TENSOR_FUNCTIONS_DIRECT_SUM_H
 
#include <prc/core/functors/view.hpp>
#include <prc/core/tensor/declarations.hpp>
#include <prc/core/tensor/views/direct_sum.hpp>
 
namespace pRC
{
    template<class XA, class XB, IsTensorish RA = RemoveReference<XA>,
        IsTensorish RB = RemoveReference<XB>>
        requires(RA::Dimension == RB::Dimension) &&
        HasCommon<typename RA::Type, typename RB::Type>
    static inline constexpr auto directSum(XA &&a, XB &&b)
    {
        if constexpr(IsInvocable<View, XA> && IsInvocable<View, XB>)
        {
            using VA = RemoveReference<decltype(view(forward<XA>(a)))>;
            using VB = RemoveReference<decltype(view(forward<XB>(b)))>;
            using T = Common<typename RA::Type, typename RB::Type>;
            return TensorViews::DirectSum<T,
                decltype(typename RA::Sizes() + typename RB::Sizes()), VA, VB>(
                view(forward<XA>(a)), view(forward<XB>(b)));
        }
        else
        {
            return eval(directSum(a, b));
        }
    }
}
#endif // pRC_CORE_TENSOR_FUNCTIONS_DIRECT_SUM_H