dyn_pool_2d_expr.hpp

//=======================================================================
// Copyright (c) 2014-2023 Baptiste Wicht
// Distributed under the terms of the MIT License.
// (See accompanying file LICENSE or copy at
//  http://opensource.org/licenses/MIT)
//=======================================================================

#pragma once

#include "etl/expr/base_temporary_expr.hpp"

namespace etl {

template <etl_expr A, typename Impl>
struct dyn_pool_2d_expr : base_temporary_expr_un<dyn_pool_2d_expr<A, Impl>, A, false> {
    using value_type = value_t<A>;                                  
    using this_type  = dyn_pool_2d_expr<A, Impl>;                   
    using base_type  = base_temporary_expr_un<this_type, A, false>; 
    using sub_traits = decay_traits<A>;                             

    static constexpr auto storage_order = sub_traits::storage_order; 

    static constexpr bool gpu_computable = Impl::template gpu_computable<A>;

    const size_t c1; 
    const size_t c2; 
    const size_t s1; 
    const size_t s2; 
    const size_t p1; 
    const size_t p2; 

    explicit dyn_pool_2d_expr(A a, size_t c1, size_t c2, size_t s1, size_t s2, size_t p1, size_t p2)
            : base_type(a), c1(c1), c2(c2), s1(s1), s2(s2), p1(p1), p2(p2) {
        //Nothing else to init
    }

    // Assignment functions

    template <same_dimensions<A> C>
    void assign_to(C&& c) const {
        inc_counter("temp:assign");

        auto& a = this->a();

        Impl::template apply<>(a, c, c1, c2, s1, s2, p1, p2);
    }

    template <same_dimensions<A> L>
    void assign_add_to(L&& lhs) const {
        std_add_evaluate(*this, lhs);
    }

    template <same_dimensions<A> L>
    void assign_sub_to(L&& lhs) const {
        std_sub_evaluate(*this, lhs);
    }

    template <same_dimensions<A> L>
    void assign_mul_to(L&& lhs) const {
        std_mul_evaluate(*this, lhs);
    }

    template <same_dimensions<A> L>
    void assign_div_to(L&& lhs) const {
        std_div_evaluate(*this, lhs);
    }

    template <same_dimensions<A> L>
    void assign_mod_to(L&& lhs) const {
        std_mod_evaluate(*this, lhs);
    }

    friend std::ostream& operator<<(std::ostream& os, const dyn_pool_2d_expr& expr) {
        return os << "pool2(" << expr._a << ")";
    }
};

template <typename A, typename Impl>
struct etl_traits<etl::dyn_pool_2d_expr<A, Impl>> {
    using expr_t     = etl::dyn_pool_2d_expr<A, Impl>; 
    using sub_expr_t = std::decay_t<A>;                
    using sub_traits = etl_traits<sub_expr_t>;         
    using value_type = value_t<A>;                     

    static constexpr size_t D = sub_traits::dimensions(); 

    static constexpr bool is_etl         = true;                                 
    static constexpr bool is_transformer = false;                                
    static constexpr bool is_view        = false;                                
    static constexpr bool is_magic_view  = false;                                
    static constexpr bool is_fast        = false;                                
    static constexpr bool is_linear      = false;                                
    static constexpr bool is_thread_safe = true;                                 
    static constexpr bool is_value       = false;                                
    static constexpr bool is_direct      = true;                                 
    static constexpr bool is_generator   = false;                                
    static constexpr bool is_padded      = false;                                
    static constexpr bool is_aligned     = true;                                 
    static constexpr bool is_temporary   = true;                                 
    static constexpr order storage_order = sub_traits::storage_order;            
    static constexpr bool gpu_computable = is_gpu_t<value_type> && cuda_enabled; 

    template <vector_mode_t V>
    static constexpr bool vectorizable = true;

    static size_t dim(const expr_t& e, size_t d) {
        if (d == D - 2) {
            return (etl::dim(e._a, d) - e.c1 + 2 * e.p1) / e.s1 + 1;
        } else if (d == D - 1) {
            return (etl::dim(e._a, d) - e.c2 + 2 * e.p2) / e.s2 + 1;
        } else {
            return etl::dim(e._a, d);
        }
    }

    static size_t size(const expr_t& e) {
        size_t acc = 1;
        for (size_t i = 0; i < D; ++i) {
            acc *= dim(e, i);
        }
        return acc;
    }

    static constexpr size_t dimensions() {
        return D;
    }

    static constexpr int complexity() noexcept {
        return -1;
    }
};

template <etl_expr E>
dyn_pool_2d_expr<detail::build_type<E>, impl::max_pool_2d> max_pool_2d(E&& value, size_t c1, size_t c2) {
    return dyn_pool_2d_expr<detail::build_type<E>, impl::max_pool_2d>{value, c1, c2, c1, c2, 0, 0};
}

template <etl_expr E>
dyn_pool_2d_expr<detail::build_type<E>, impl::max_pool_2d> max_pool_2d(E&& value, size_t c1, size_t c2, size_t s1, size_t s2, size_t p1 = 0, size_t p2 = 0) {
    return dyn_pool_2d_expr<detail::build_type<E>, impl::max_pool_2d>{value, c1, c2, s1, s2, p1, p2};
}

template <etl_expr E>
dyn_pool_2d_expr<detail::build_type<E>, impl::avg_pool_2d> avg_pool_2d(E&& value, size_t c1, size_t c2) {
    return dyn_pool_2d_expr<detail::build_type<E>, impl::avg_pool_2d>{value, c1, c2, c1, c2, 0, 0};
}

template <etl_expr E>
dyn_pool_2d_expr<detail::build_type<E>, impl::avg_pool_2d> avg_pool_2d(E&& value, size_t c1, size_t c2, size_t s1, size_t s2, size_t p1 = 0, size_t p2 = 0) {
    return dyn_pool_2d_expr<detail::build_type<E>, impl::avg_pool_2d>{value, c1, c2, s1, s2, p1, p2};
}

template <etl_expr E>
dyn_pool_2d_expr<detail::build_type<E>, impl::standard::dyn_pmp_p_impl> p_max_pool_p(E&& value, size_t c1, size_t c2) {
    validate_pmax_pooling(value, c1, c2);
    return dyn_pool_2d_expr<detail::build_type<E>, impl::standard::dyn_pmp_p_impl>{value, c1, c2, c1, c2, 0, 0};
}

} //end of namespace etl