Function test

Tests/ArrayND/main.cpp:82–122

template <typename T, int N, bool C>
void test (ArrayND<T,N,C>& a)
{
    auto sz = a.size();
    auto* p = a.dataPtr();
    amrex::ParallelForRNG(sz, [=] AMREX_GPU_DEVICE (std::size_t i,
                                                    RandomEngine const& eng)
    {
        p[i] = amrex::Random_int(1000000, eng);
    });
    T tot = Reduce::Sum<T>(sz, p, 0);

    if constexpr (std::remove_reference_t<decltype(a)>::IsArray4_v) {
        test_array4(a, tot);
    }

    if constexpr (std::remove_reference_t<decltype(a)>::IsLastDimComponent_v) {
        test_comp(a, tot);
    }

    {
        ReduceOps<ReduceOpSum> reduce_op;
        ReduceData<T> reduce_data(reduce_op);
        using ReduceTuple = typename decltype(reduce_data)::Type;
        BoxND<N> box(a.begin, a.end-1);
        reduce_op.eval(box, reduce_data,
        [=] AMREX_GPU_DEVICE (IntVectND<N> const& iv) -> ReduceTuple
        {
            auto v0 = a(iv);
            auto* p0 = a.ptr(iv);
            auto v1 = amrex::Apply([&] (auto&&... i) {
                return a(i...);
            }, iv);
            auto* p1 = amrex::Apply([&] (auto&&... i) {
                return a.ptr(i...);
            }, iv);
            return (v0+v1+*p0+*p1)/4;
        });
        ReduceTuple hv = reduce_data.value(reduce_op);
        AMREX_ALWAYS_ASSERT(tot == amrex::get<0>(hv));
    }
}

int main (int argc, char* argv[])
{