block.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_CORE_TENSOR_VIEWS_BLOCK_H
#define pRC_CORE_TENSOR_VIEWS_BLOCK_H
 
#include <prc/core/basic/sequence.hpp>
#include <prc/core/container/subscripts.hpp>
#include <prc/core/log/log.hpp>
#include <prc/core/tensor/declarations.hpp>
 
namespace pRC::TensorViews
{
    template<class T, class N, class B, class V>
    class Block;
 
    template<class T, Size... Ns, Size... Bs, class V>
        requires IsTensorView<V>
    class Block<T, Sizes<Ns...>, Sizes<Bs...>, V>
        : public Conditional<IsAssignable<V>,
              Assignable<T, Sizes<Ns...>,
                  Block<T, Sizes<Ns...>, Sizes<Bs...>, V>>,
              View<T, Sizes<Ns...>, Block<T, Sizes<Ns...>, Sizes<Bs...>, V>>>
    {
    private:
        using Base = Conditional<IsAssignable<V>,
            Assignable<T, Sizes<Ns...>, Block>, View<T, Sizes<Ns...>, Block>>;
 
    public:
        template<class X>
            requires IsSame<V, RemoveReference<X>>
        Block(X &&a, Subscripts<Bs...> const &offsets)
            : mA(forward<X>(a))
            , mOffsets(offsets)
        {
            if constexpr(cDebugLevel >= DebugLevel::Low)
            {
                if(mOffsets.isOutOfRange())
                {
                    Logging::error("Block offset indices out of range.");
                }
            }
        }
 
        using Base::operator=;
 
        template<IsConvertible<Index>... Is>
            requires(sizeof...(Is) == Base::Dimension)
        constexpr decltype(auto) operator()(Is const... indices)
        {
            return this->call(indices...);
        }
 
        template<IsConvertible<Index>... Is>
            requires(sizeof...(Is) == Base::Dimension)
        constexpr decltype(auto) operator()(Is const... indices) const
        {
            return this->call(indices...);
        }
 
        constexpr decltype(auto) operator()(
            typename Base::Subscripts const &subscripts)
        {
            return expand(makeSeriesFor<Index, Ns...>(),
                [this, &subscripts](auto const... seq) -> decltype(auto)
                {
                    return mA((mOffsets[seq] * Ns + subscripts[seq])...);
                });
        }
 
        constexpr decltype(auto) operator()(
            typename Base::Subscripts const &subscripts) const
        {
            return expand(makeSeriesFor<Index, Ns...>(),
                [this, &subscripts](auto const... seq) -> decltype(auto)
                {
                    return mA((mOffsets[seq] * Ns + subscripts[seq])...);
                });
        }
 
        constexpr decltype(auto) operator[](Index const index) = delete;
        constexpr decltype(auto) operator[](Index const index) const = delete;
 
    private:
        V mA;
        BEGIN_IGNORE_DIAGNOSTIC_GCC("-Wpedantic")
        Subscripts<Bs...> const mOffsets;
        END_IGNORE_DIAGNOSTIC_GCC
    };
}
#endif // pRC_CORE_TENSOR_VIEWS_BLOCK_H